This study aims to formulate a steady-state mathematical model for a three-dimensional permeable enclosure(cavity)to determine the oil extraction rate using three distinct nanoparticles,SiO_(2),Al_(2)O_(3),and Fe_(2)O...This study aims to formulate a steady-state mathematical model for a three-dimensional permeable enclosure(cavity)to determine the oil extraction rate using three distinct nanoparticles,SiO_(2),Al_(2)O_(3),and Fe_(2)O_(3),in unconventional oil reservoirs.The simulation is conducted for different parameters of volume fractions,porosities,and mass flow rates to determine the optimal oil recovery.The impact of nanoparticles on relative permeability(kr)and water is also investigated.The simulation process utilizes the finite volume ANSYS Fluent.The study results showed that when the mass flow rate at the inlet is low,oil recovery goes up.In addition,they indicated that silicon nanoparticles are better at getting oil out of the ground(i.e.,oil reservoir)than Al_(2)O_(3)and Fe_(2)O_(3).Most oil can be extracted from SiO_(2),Al_(2)O_(3),and Fe_(2)O_(3)at a rate of 97.8%,96.5%,and 88%,respectively.展开更多
Reactive transport modeling(RTM)is an emerging method used to address geological issues in diagenesis research.However,the extrapolation of RTM results to practical reservoir prediction is not sufficiently understood....Reactive transport modeling(RTM)is an emerging method used to address geological issues in diagenesis research.However,the extrapolation of RTM results to practical reservoir prediction is not sufficiently understood.This paper presents a case study of the Eocene Qaidam Basin that combines RTM results with petrological and mineralogical evidence.The results show that the Eocene Xiaganchaigou Formation is characterized by mixed siliciclastic-carbonate-evaporite sedimentation in a semiclosed saline lacustrine environment.Periodic evaporation and salinization processes during the syngeneticpenecontemporaneous stage gave rise to the replacive genesis of dolomites and the cyclic enrichment of dolomite in the middle-upper parts of the meter-scale depositional sequences.The successive change in mineral paragenesis from terrigenous clastics to carbonates to evaporites was reconstructed using RTM simulations.Parametric uncertainty analyses further suggest that the evaporation intensity(brine salinity)and particle size of sediments(reactive surface area)were important rate-determining factors in the dolomitization,as shown by the relatively higher reaction rates under conditions of higher brine salinity and fine-grained sediments.Combining the simulation results with measured mineralogical and reservoir physical property data indicates that the preservation of original intergranular pores and the generation of porosity via replacive dolomitization were the major formation mechanisms of the distinctive lacustrine dolomite reservoirs(widespread submicron intercrystalline micropores)in the Eocene Qaidam Basin.The results confirm that RTM can be effectively used in geological studies,can provide a better general understanding of the dolomitizing fluid-rock interactions,and can shed light on the spatiotemporal evolution of mineralogy and porosity during dolomitization and the formation of lacustrine dolomite reservoirs.展开更多
The paper presents the results of comprehensive studies of filtration and capacitance properties of highly porous reservoir rocks of the aquifer of an underground gas storage facility.The geomechanical part of the res...The paper presents the results of comprehensive studies of filtration and capacitance properties of highly porous reservoir rocks of the aquifer of an underground gas storage facility.The geomechanical part of the research included studying the dependence of rock permeability on the stress-strain state in the vicinity of the wells,and physical modeling of the implementation of the method of increasing the permeability of the wellbore zone-the method of directional unloading of the reservoir.The digital part of the research included computed tomography(CT)-based computer analysis of the internal structure,pore space characteristics,and filtration properties before and after the tests.According to the results of physical modeling of deformation and filtration processes,it is found that the permeability of rocks before fracture depends on the stress-strain state insignificantly,and this influence is reversible.However,when downhole pressure reaches 7-8 MPa,macrocracks in the rock begin to grow,accompanied by irreversible permeability increase.Porosity,geodesic tortuosity and permeability values were obtained based on digital studies and numerical modeling.A weak degree of transversal anisotropy of the filtration properties of rocks was detected.Based on the analysis of pore size distribution,pressure field and flow velocities,high homogeneity and connectivity of the rock pore space is shown.The absence of pronounced changes in pore space characteristics and pore permeability after non-uniform triaxial loading rocks was shown.On the basis of geometrical analysis of pore space,the reasons for weak permeability anisotropy were identified.The filtration-capacitance properties obtained from the digital analysis showed very good agreement with the results of field and laboratory measurements.The physical modeling has confirmed the efficiency of application of the directional unloading method for the reservoir under study.The necessary parameters of its application were calculated:bottomhole geometry,stage of operation,stresses and pressure drawdown value.展开更多
Flow units(FU)rock typing is a common technique for characterizing reservoir flow behavior,producing reliable porosity and permeability estimation even in complex geological settings.However,the lateral extrapolation ...Flow units(FU)rock typing is a common technique for characterizing reservoir flow behavior,producing reliable porosity and permeability estimation even in complex geological settings.However,the lateral extrapolation of FU away from the well into the whole reservoir grid is commonly a difficult task and using the seismic data as constraints is rarely a subject of study.This paper proposes a workflow to generate numerous possible 3D volumes of flow units,porosity and permeability below the seismic resolution limit,respecting the available seismic data at larger scales.The methodology is used in the Mero Field,a Brazilian presalt carbonate reservoir located in the Santos Basin,who presents a complex and heterogenic geological setting with different sedimentological processes and diagenetic history.We generated metric flow units using the conventional core analysis and transposed to the well log data.Then,given a Markov chain Monte Carlo algorithm,the seismic data and the well log statistics,we simulated acoustic impedance,decametric flow units(DFU),metric flow units(MFU),porosity and permeability volumes in the metric scale.The aim is to estimate a minimum amount of MFU able to calculate realistic scenarios porosity and permeability scenarios,without losing the seismic lateral control.In other words,every porosity and permeability volume simulated produces a synthetic seismic that match the real seismic of the area,even in the metric scale.The achieved 3D results represent a high-resolution fluid flow reservoir modelling considering the lateral control of the seismic during the process and can be directly incorporated in the dynamic characterization workflow.展开更多
Through core observation,thin section identification,X-ray diffraction analysis,scanning electron microscopy,and low-temperature nitrogen adsorption and isothermal adsorption experiments,the lithology and pore charact...Through core observation,thin section identification,X-ray diffraction analysis,scanning electron microscopy,and low-temperature nitrogen adsorption and isothermal adsorption experiments,the lithology and pore characteristics of the Upper Carboniferous bauxite series in eastern Ordos Basin were analyzed to reveal the formation and evolution process of the bauxite reservoirs.A petrological nomenclature and classification scheme for bauxitic rocks based on three units(aluminum hydroxides,iron minerals and clay minerals)is proposed.It is found that bauxitic mudstone is in the form of dense massive and clastic structures,while the(clayey)bauxite is of dense massive,pisolite,oolite,porous soil and clastic structures.Both bauxitic mudstone and bauxite reservoirs develop dissolution pores,intercrystalline pores,and microfractures as the dominant gas storage space,with the porosity less than 10% and mesopores in dominance.The bauxite series in the North China Craton can be divided into five sections,i.e.,ferrilite(Shanxi-style iron ore,section A),bauxitic mudstone(section B),bauxite(section C),bauxite mudstone(debris-containing,section D)and dark mudstone-coal section(section E).The burrow/funnel filling,lenticular,layered/massive bauxite deposits occur separately in the karst platforms,gentle slopes and low-lying areas.The karst platforms and gentle slopes are conducive to surface water leaching,with strong karstification,well-developed pores,large reservoir thickness and good physical properties,but poor strata continuity.The low-lying areas have poor physical properties but relatively continuous and stable reservoirs.The gas enrichment in bauxites is jointly controlled by source rock,reservoir rock and fractures.This recognition provides geological basis for the exploration and development of natural gas in the Upper Carboniferous in the study area and similar bauxite systems.展开更多
This research uses both three-dimensional(3D)modeling and geologic well control to piece back together the architectural parts of the Late Cretaceous formations.The goal is to figure out the sizes,directions,locations...This research uses both three-dimensional(3D)modeling and geologic well control to piece back together the architectural parts of the Late Cretaceous formations.The goal is to figure out the sizes,directions,locations,and controls of the layers of the fluvial sandstone reservoirs.Sequence stratigraphy is essential for 3D reservoir modeling and petroleum geology understanding in the Bahga oilfield.The purpose of this work is to create a static model that shows the layers and facies distribution in the reservoir interval.We will use data from nine well logs and 22 seismic lines calibrated by the Abu Roash G Member reservoir core intervals to accomplish this.The petrophysical study discovered three parts in the Abu Roash G Member reservoir rock:channel fill that is affected by tides,channel fill that is dominated by tides(intertidal sands),and channel top with lenticular bedded sandstone.The model's findings point to the existence of an NNW-oriented sand body,which could be a prime location to produce hydrocarbons.The original oil in place(OOIP)is about 3,438,279 Stock Tank Barrels(STB),and the oil reserve reaches up to 1,031,484(STB).Sequence stratigraphic analysis using seismic and well log information(SB)reveals that the Upper Cretaceous AR/G reservoir of the Bahga field is characterized by third-and fourth-order stratigraphic sequences,which are constrained by three Maximum Flooding Surfaces(MFS)and two Sequence Boundaries.The integration of the derived geological model and sequence stratigraphic results can lower future extraction risk by identifying the locations and trends of the geologic facies with the necessary petrophysical properties for the hydrocarbon accumulations.展开更多
Using the typical characteristics of multi-layered marine and continental transitional gas reservoirs as a basis,a model is developed to predict the related well production rate.This model relies on the fractal theory...Using the typical characteristics of multi-layered marine and continental transitional gas reservoirs as a basis,a model is developed to predict the related well production rate.This model relies on the fractal theory of tortuous capillary bundles and can take into account multiple gas flow mechanisms at the micrometer and nanometer scales,as well as the flow characteristics in different types of thin layers(tight sandstone gas,shale gas,and coalbed gas).Moreover,a source-sink function concept and a pressure drop superposition principle are utilized to introduce a coupled flow model in the reservoir.A semi-analytical solution for the production rate is obtained using a matrix iteration method.A specific well is selected for fitting dynamic production data,and the calculation results show that the tight sandstone has the highest gas production per unit thickness compared with the other types of reservoirs.Moreover,desorption and diffusion of coalbed gas and shale gas can significantly contribute to gas production,and the daily production of these two gases decreases rapidly with decreasing reservoir pressure.Interestingly,the gas production from fractures exhibits an approximately U-shaped distribution,indicating the need to optimize the spacing between clusters during hydraulic fracturing to reduce the area of overlapping fracture control.The coal matrix water saturation significantly affects the coalbed gas production,with higher water saturation leading to lower production.展开更多
We apply stochastic seismic inversion and Bayesian facies classification for porosity modeling and igneous rock identification in the presalt interval of the Santos Basin. This integration of seismic and well-derived ...We apply stochastic seismic inversion and Bayesian facies classification for porosity modeling and igneous rock identification in the presalt interval of the Santos Basin. This integration of seismic and well-derived information enhances reservoir characterization. Stochastic inversion and Bayesian classification are powerful tools because they permit addressing the uncertainties in the model. We used the ES-MDA algorithm to achieve the realizations equivalent to the percentiles P10, P50, and P90 of acoustic impedance, a novel method for acoustic inversion in presalt. The facies were divided into five: reservoir 1,reservoir 2, tight carbonates, clayey rocks, and igneous rocks. To deal with the overlaps in acoustic impedance values of facies, we included geological information using a priori probability, indicating that structural highs are reservoir-dominated. To illustrate our approach, we conducted porosity modeling using facies-related rock-physics models for rock-physics inversion in an area with a well drilled in a coquina bank and evaluated the thickness and extension of an igneous intrusion near the carbonate-salt interface. The modeled porosity and the classified seismic facies are in good agreement with the ones observed in the wells. Notably, the coquinas bank presents an improvement in the porosity towards the top. The a priori probability model was crucial for limiting the clayey rocks to the structural lows. In Well B, the hit rate of the igneous rock in the three scenarios is higher than 60%, showing an excellent thickness-prediction capability.展开更多
Spontaneous imbibition(SI)is an important mechanism for enhancing oil recovery in low-permeability reservoirs.Due to the strong heterogeneity,and the non-Darcy flow,the construction of SI model for lowpermeability res...Spontaneous imbibition(SI)is an important mechanism for enhancing oil recovery in low-permeability reservoirs.Due to the strong heterogeneity,and the non-Darcy flow,the construction of SI model for lowpermeability reservoirs is extremely challenging.Commonly,traditional SI models based on single or averaged capillary tortuosity ignore the influence of heterogeneity of pore seepage channels and the threshold pressure(TP)on imbibition.Therefore,in this work,based on capillary model and fractal theory,a mathematical model of characterizing SI considering heterogeneity of pore seepage channels is established.On this basis,the threshold pressure was introduced to determine the pore radius at which the wetted phase can displace oil.The proposed new SI model was verified by imbibition experimental data.The study shows that for weakly heterogeneous cores with permeability of 0-1 m D,the traditional SI model can characterize the imbibition process relatively accurately,and the new imbibition model can increase the coefficient of determination by 1.05 times.However,traditional model has serious deviations in predicting the imbibition recovery for cores with permeability of 10-50 m D.The new SI model coupling with heterogeneity of pore seepage channels and threshold pressure effectively solves this problem,and the determination coefficient is increased from 0.344 to 0.922,which is increased by2.68 times.For low-permeability reservoirs,the production of the oil in transitional pores(0.01-0.1μm)and mesopores(0.1-1μm)significantly affects the imbibition recovery,as the research shows that when the heterogeneity of pore seepage channels is ignored,the oil recovery in transitional pores and mesopores decreases by 7.54%and 4.26%,respectively.Sensitivity analysis shows that increasing interfacial tension,decreasing contact angle,oil-water viscosity ratio and threshold pressure will increase imbibition recovery.In addition,there are critical values for the influence of these factors on the imbibition recovery,which provides theoretical support for surfactant optimization.展开更多
From June 2008 to August 2013,approximately 67 kt of CO_(2) was injected into a deep saline formation at the Ketzin pilot CO_(2) storage site.During injection,3D seismic surveys have been performed to monitor the migr...From June 2008 to August 2013,approximately 67 kt of CO_(2) was injected into a deep saline formation at the Ketzin pilot CO_(2) storage site.During injection,3D seismic surveys have been performed to monitor the migration of sequestered CO_(2).Seismic monitoring results are limited by the acquisition and signal-to-noise ratio of the acquired data.The multiphysical reservoir simulation provides information regarding the CO_(2) fluid behavior,and the approximated model should be calibrated with the monitoring results.In this work,property models are delivered from the multiphysical model during 3D repeated seismic surveys.The simulated seismic data based on the models are compared with the real data,and the results validate the effectiveness of the multiphysical inversion method.Time-lapse analysis shows the trend of CO_(2) migration during and after injection.展开更多
Under the policy background and advocacy of carbon capture,utilization,and storage(CCUS),CO_(2)-EOR has become a promising direction in the shale oil reservoir industry.The multi-scale pore structure distribution and ...Under the policy background and advocacy of carbon capture,utilization,and storage(CCUS),CO_(2)-EOR has become a promising direction in the shale oil reservoir industry.The multi-scale pore structure distribution and fracture structure lead to complex multiphase flow,comprehensively considering multiple mechanisms is crucial for development and CO_(2) storage in fractured shale reservoirs.In this paper,a multi-mechanism coupled model is developed by MATLAB.Compared to the traditional Eclipse300 and MATLAB Reservoir Simulation Toolbox(MRST),this model considers the impact of pore structure on fluid phase behavior by the modified Peng—Robinson equation of state(PR-EOS),and the effect simultaneously radiate to Maxwell—Stefan(M—S)diffusion,stress sensitivity,the nano-confinement(NC)effect.Moreover,a modified embedded discrete fracture model(EDFM)is used to model the complex fractures,which optimizes connection types and half-transmissibility calculation approaches between non-neighboring connections(NNCs).The full implicit equation adopts the finite volume method(FVM)and Newton—Raphson iteration for discretization and solution.The model verification with the Eclipse300 and MRST is satisfactory.The results show that the interaction between the mechanisms significantly affects the production performance and storage characteristics.The effect of molecular diffusion may be overestimated in oil-dominated(liquid-dominated)shale reservoirs.The well spacing and injection gas rate are the most crucial factors affecting the production by sensitivity analysis.Moreover,the potential gas invasion risk is mentioned.This model provides a reliable theoretical basis for CO_(2)-EOR and sequestration in shale oil reservoirs.展开更多
Based on the geochemical,seismic,logging and drilling data,the Fuyu reservoirs of the Lower Cretaceous Quantou Formation in northern Songliao Basin are systematically studied in terms of the geological characteristics...Based on the geochemical,seismic,logging and drilling data,the Fuyu reservoirs of the Lower Cretaceous Quantou Formation in northern Songliao Basin are systematically studied in terms of the geological characteristics,the tight oil enrichment model and its major controlling factors.First,the Quantou Formation is overlaid by high-quality source rocks of the Upper Cretaceous Qingshankou Formation,with the development of nose structure around sag and the broad and continuous distribution of sand bodies.The reservoirs are tight on the whole.Second,the configuration of multiple elements,such as high-quality source rocks,reservoir rocks,fault,overpressure and structure,controls the tight oil enrichment in the Fuyu reservoirs.The source-reservoir combination controls the tight oil distribution pattern.The pressure difference between source and reservoir drives the charging of tight oil.The fault-sandbody transport system determines the migration and accumulation of oil and gas.The positive structure is the favorable place for tight oil enrichment,and the fault-horst zone is the key part of syncline area for tight oil exploration.Third,based on the source-reservoir relationship,transport mode,accumulation dynamics and other elements,three tight oil enrichment models are recognized in the Fuyu reservoirs:(1)vertical or lateral migration of hydrocarbon from source rocks to adjacent reservoir rocks,that is,driven by overpressure,hydrocarbon generated is migrated vertically or laterally to and accumulates in the adjacent reservoir rocks;(2)transport of hydrocarbon through faults between separated source and reservoirs,that is,driven by overpressure,hydrocarbon migrates downward through faults to the sandbodies that are separated from the source rocks;and(3)migration of hydrocarbon through faults and sandbodies between separated source and reservoirs,that is,driven by overpressure,hydrocarbon migrates downwards through faults to the reservoir rocks that are separated from the source rocks,and then migrates laterally through sandbodies.Fourth,the differences in oil source conditions,charging drive,fault distribution,sandbody and reservoir physical properties cause the differential enrichment of tight oil in the Fuyu reservoirs.Comprehensive analysis suggests that the Fuyu reservoir in the Qijia-Gulong Sag has good conditions for tight oil enrichment and has been less explored,and it is an important new zone for tight oil exploration in the future.展开更多
In heterogeneous natural gas reservoirs, gas is generally present as small patchlike pockets embedded in the water-saturated host matrix. This type of heterogeneity, also called "patchy saturation", causes s...In heterogeneous natural gas reservoirs, gas is generally present as small patchlike pockets embedded in the water-saturated host matrix. This type of heterogeneity, also called "patchy saturation", causes significant seismic velocity dispersion and attenuation. To establish the relation between seismic response and type of fluids, we designed a rock physics model for carbonates. First, we performed CT scanning and analysis of the fluid distribution in the partially saturated rocks. Then, we predicted the quantitative relation between the wave response at different frequency ranges and the basic lithological properties and pore fluids. A rock physics template was constructed based on thin section analysis of pore structures and seismic inversion. This approach was applied to the limestone gas reservoirs of the right bank block of the Amu Darya River. Based on poststack wave impedance and prestack elastic parameter inversions, the seismic data were used to estimate rock porosity and gas saturation. The model results were in good agreement with the production regime of the wells.展开更多
Based on the percolation network model characterizing reservoir rock's pore structure and fluid characteristics, this paper qualitatively studies the effects of pore size, pore shape, pore connectivity, and the amoun...Based on the percolation network model characterizing reservoir rock's pore structure and fluid characteristics, this paper qualitatively studies the effects of pore size, pore shape, pore connectivity, and the amount of micropores on the I - Sw curve using numerical modeling. The effects of formation water salinity on the electrical resistivity of the rock are discussed. Then the relative magnitudes of the different influencing factors are discussed. The effects of the different factors on the I - Sw curve are analyzed by fitting simulation results. The results show that the connectivity of the void spaces and the amount of micropores have a large effect on the I - S, curve, while the other factors have little effect. The formation water salinity has a large effect on the absolute resistivity values. The non-Archie phenomenon is prevalent, which is remarkable in rocks with low permeability.展开更多
Geothermal energy is a kind of renewable,sustainable and clean energy resource.Geothermal energy is abundant in carbonate reservoirs.However,low matrix permeability limits its exploitation.The super-critical carbon di...Geothermal energy is a kind of renewable,sustainable and clean energy resource.Geothermal energy is abundant in carbonate reservoirs.However,low matrix permeability limits its exploitation.The super-critical carbon dioxide(SC-CO_(2))jet fracturing is expected to efficiently stimulate the carbonate geothermal reservoirs and achieve the storage of CO_(2) simultaneously.In this paper,we established a transient seepage and fluid-thermo-mechanical coupled model to analyze the impact performance of sc-CO_(2) jet fracturing.The mesh-based parallel code coupling interface was employed to couple the fluid and solid domains by exchanging the data through the mesh interface.The physical properties change of sC-CO_(2) with temperature were considered in the numerical model.Results showed that SC-CO_(2) jet frac-turing is superior to water-jet fracturing with respect to jetting velocity,particle trajectory and pene-trability.Besides,stress distribution on the carbonate rock showed that the tensile and shear failure would more easily occur by SC-CO_(2) jet than that by water jet.Moreover,pressure and temperature control the jet field and seepage field of sC-CO_(2) simultaneously.Increasing the jet temperature can effectively enhance the impingement effect and seepage process by decreasing the viscosity and density of SC-CO_(2).The key findings are expected to provide a theoretical basis and design reference for applying SC-CO_(2) jet fracturing in carbonate geothermal reservoirs.展开更多
Hydrocarbon production from shale has attracted much attention in the recent years. When applied to this prolific and hydrocarbon rich resource plays, our understanding of the complexities of the flow mechanism(sorpt...Hydrocarbon production from shale has attracted much attention in the recent years. When applied to this prolific and hydrocarbon rich resource plays, our understanding of the complexities of the flow mechanism(sorption process and flow behavior in complex fracture systems- induced or natural) leaves much to be desired. In this paper, we present and discuss a novel approach to modeling, history matching of hydrocarbon production from a Marcellus shale asset in southwestern Pennsylvania using advanced data mining, pattern recognition and machine learning technologies. In this new approach instead of imposing our understanding of the flow mechanism, the impact of multi-stage hydraulic fractures, and the production process on the reservoir model, we allow the production history, well log, completion and hydraulic fracturing data to guide our model and determine its behavior. The uniqueness of this technology is that it incorporates the so-called "hard data" directly into the reservoir model, so that the model can be used to optimize the hydraulic fracture process. The "hard data" refers to field measurements during the hydraulic fracturing process such as fluid and proppant type and amount, injection pressure and rate as well as proppant concentration. This novel approach contrasts with the current industry focus on the use of "soft data"(non-measured, interpretive data such as frac length, width,height and conductivity) in the reservoir models. The study focuses on a Marcellus shale asset that includes 135 wells with multiple pads, different landing targets, well length and reservoir properties. The full field history matching process was successfully completed using this data driven approach thus capturing the production behavior with acceptable accuracy for individual wells and for the entire asset.展开更多
Due to inherent limits of data acquisition and geophysical data resolution, there are large uncertainties in the characterization of subsurface fractures. However, outcrop analogies can provide qualitative and quantit...Due to inherent limits of data acquisition and geophysical data resolution, there are large uncertainties in the characterization of subsurface fractures. However, outcrop analogies can provide qualitative and quantitative information on a large number of fractures, based on which the accuracy of subsurface fracture characterization can be improved. Here we take the tectonic fracture modeling of an ultra-low permeability sandstone reservoir based on an outcrop analogy, a case study of the Chang6t~ Formation of the Upper Triassic Yanchang Group of the Wangyao Oilfield in the Ordos Basin of China. An outcrop at the edge of the basin is a suitable analog for the reservoir, but the prerequisite is that they must have equivalent previous stress fields, similar final structural characteristics, relative timing and an identical depositional environment and diagenesis. The relationship among fracture density, rock type and bed thickness based on the outcrop is one of the most important fracture distribution models, and can be used to interpret fracture density in individual wells quantitatively. Fracture orientation, dip, geometry and scale, also should be described and measured in the outcrop, and can be used together with structure restoration and single well fracture density interpretation to guide fracture intensity prediction on bed surfaces and to constrain the construction of the 3D fracture geometry model of the subsurface reservoir. The application of the above principles shows the outcrop-based tectonic fracture models of the target ultra-low permeability sandstone reservoir are consistent with fractures inferred from microseismic interpretation and tracer tests. This illustrated that the fracture modeling based on the outcrop analogy is reliable and can reduce the uncertainty in stochastic fracture modeling.展开更多
This paper discusses the principles of geologic constraints on reservoir stochastic modeling. By using the system science theory, two kinds of uncertainties, including random uncertainty and fuzzy uncertainty, are rec...This paper discusses the principles of geologic constraints on reservoir stochastic modeling. By using the system science theory, two kinds of uncertainties, including random uncertainty and fuzzy uncertainty, are recognized. In order to improve the precision of stochastic modeling and reduce the uncertainty in realization, the fuzzy uncertainty should be stressed, and the "geological genesis-controlled modeling" is conducted under the guidance of a quantitative geological pattern. An example of the Pingqiao horizontal-well division of the Ansai Oilfield in the Ordos Basin is taken to expound the method of stochastic modeling.展开更多
Fractured reservoirs are an important target for oil and gas exploration in the Tarim Basin and the prediction of this type of reservoir is challenging.Due to the complicated fracture system in the Tarim Basin,the con...Fractured reservoirs are an important target for oil and gas exploration in the Tarim Basin and the prediction of this type of reservoir is challenging.Due to the complicated fracture system in the Tarim Basin,the conventional AVO inversion method based on HTI theory to predict fracture development will result in some errors.Thus,an integrated research concept for fractured reservoir prediction is put forward in this paper.Seismic modeling plays a bridging role in this concept,and the establishment of an anisotropic fracture model by Discrete Fracture Network (DFN) is the key part.Because the fracture system in the Tarim Basin shows complex anisotropic characteristics,it is vital to build an effective anisotropic model.Based on geological,well logging and seismic data,an effective anisotropic model of complex fracture systems can be set up with the DFN method.The effective elastic coefficients,and the input data for seismic modeling can be calculated.Then seismic modeling based on this model is performed,and the seismic response characteristics are analyzed.The modeling results can be used in the following AVO inversion for fracture detection.展开更多
Porous carbonate reservoirs,prevalent in the Middle East,are lithologically dominated by bioclastic limestones,exhibiting high porosity,low permeability,intricate pore structure,and strong heterogeneity.Waterflooding ...Porous carbonate reservoirs,prevalent in the Middle East,are lithologically dominated by bioclastic limestones,exhibiting high porosity,low permeability,intricate pore structure,and strong heterogeneity.Waterflooding through horizontal wells is commonly used for exploiting these reservoirs.However,challenges persist,such as significant uncertainty and complex operational procedures regarding adjustment effects during the exploitation.The USH reservoir of the Cretaceous D oilfield,Oman exemplifies typical porous carbonate reservoirs.It initially underwent depletion drive using vertical wells,followed by horizontal well waterflooding in the late stage.Currently,the oilfield is confronted with substantial developmental challenges,involving the understanding of residual oil distribution,effective water cut control,and sustaining oil production since it has entered the late development stage.Employing a microscopic visualization displacement system equipped with electrodes,this study elucidated the waterflooding mechanisms and residual oil distribution in the late-stage development of the USH reservoir.The results reveal that the reservoir's vertical stacking patterns act as the main factor affecting the horizontal well waterflooding efficacy.Distinct water flow channels emerge under varying reservoir stacking patterns,with post-waterflooding residual oil predominantly distributed at the reservoir's top and bottom.The oil recovery can be enhanced by adjusting the waterflooding's flow line and intensity.The findings of this study will provide theoretical insights of waterflooding mechanisms and injection-production adjustments for exploiting other porous carbonate reservoirs in the Middle East through horizontal wells.展开更多
基金The APC of this article is covered by Research Grant YUTP 015LCO-526。
文摘This study aims to formulate a steady-state mathematical model for a three-dimensional permeable enclosure(cavity)to determine the oil extraction rate using three distinct nanoparticles,SiO_(2),Al_(2)O_(3),and Fe_(2)O_(3),in unconventional oil reservoirs.The simulation is conducted for different parameters of volume fractions,porosities,and mass flow rates to determine the optimal oil recovery.The impact of nanoparticles on relative permeability(kr)and water is also investigated.The simulation process utilizes the finite volume ANSYS Fluent.The study results showed that when the mass flow rate at the inlet is low,oil recovery goes up.In addition,they indicated that silicon nanoparticles are better at getting oil out of the ground(i.e.,oil reservoir)than Al_(2)O_(3)and Fe_(2)O_(3).Most oil can be extracted from SiO_(2),Al_(2)O_(3),and Fe_(2)O_(3)at a rate of 97.8%,96.5%,and 88%,respectively.
文摘Reactive transport modeling(RTM)is an emerging method used to address geological issues in diagenesis research.However,the extrapolation of RTM results to practical reservoir prediction is not sufficiently understood.This paper presents a case study of the Eocene Qaidam Basin that combines RTM results with petrological and mineralogical evidence.The results show that the Eocene Xiaganchaigou Formation is characterized by mixed siliciclastic-carbonate-evaporite sedimentation in a semiclosed saline lacustrine environment.Periodic evaporation and salinization processes during the syngeneticpenecontemporaneous stage gave rise to the replacive genesis of dolomites and the cyclic enrichment of dolomite in the middle-upper parts of the meter-scale depositional sequences.The successive change in mineral paragenesis from terrigenous clastics to carbonates to evaporites was reconstructed using RTM simulations.Parametric uncertainty analyses further suggest that the evaporation intensity(brine salinity)and particle size of sediments(reactive surface area)were important rate-determining factors in the dolomitization,as shown by the relatively higher reaction rates under conditions of higher brine salinity and fine-grained sediments.Combining the simulation results with measured mineralogical and reservoir physical property data indicates that the preservation of original intergranular pores and the generation of porosity via replacive dolomitization were the major formation mechanisms of the distinctive lacustrine dolomite reservoirs(widespread submicron intercrystalline micropores)in the Eocene Qaidam Basin.The results confirm that RTM can be effectively used in geological studies,can provide a better general understanding of the dolomitizing fluid-rock interactions,and can shed light on the spatiotemporal evolution of mineralogy and porosity during dolomitization and the formation of lacustrine dolomite reservoirs.
基金supported by the Russian Science Foundation(Grant No.22-11-00273).
文摘The paper presents the results of comprehensive studies of filtration and capacitance properties of highly porous reservoir rocks of the aquifer of an underground gas storage facility.The geomechanical part of the research included studying the dependence of rock permeability on the stress-strain state in the vicinity of the wells,and physical modeling of the implementation of the method of increasing the permeability of the wellbore zone-the method of directional unloading of the reservoir.The digital part of the research included computed tomography(CT)-based computer analysis of the internal structure,pore space characteristics,and filtration properties before and after the tests.According to the results of physical modeling of deformation and filtration processes,it is found that the permeability of rocks before fracture depends on the stress-strain state insignificantly,and this influence is reversible.However,when downhole pressure reaches 7-8 MPa,macrocracks in the rock begin to grow,accompanied by irreversible permeability increase.Porosity,geodesic tortuosity and permeability values were obtained based on digital studies and numerical modeling.A weak degree of transversal anisotropy of the filtration properties of rocks was detected.Based on the analysis of pore size distribution,pressure field and flow velocities,high homogeneity and connectivity of the rock pore space is shown.The absence of pronounced changes in pore space characteristics and pore permeability after non-uniform triaxial loading rocks was shown.On the basis of geometrical analysis of pore space,the reasons for weak permeability anisotropy were identified.The filtration-capacitance properties obtained from the digital analysis showed very good agreement with the results of field and laboratory measurements.The physical modeling has confirmed the efficiency of application of the directional unloading method for the reservoir under study.The necessary parameters of its application were calculated:bottomhole geometry,stage of operation,stresses and pressure drawdown value.
文摘Flow units(FU)rock typing is a common technique for characterizing reservoir flow behavior,producing reliable porosity and permeability estimation even in complex geological settings.However,the lateral extrapolation of FU away from the well into the whole reservoir grid is commonly a difficult task and using the seismic data as constraints is rarely a subject of study.This paper proposes a workflow to generate numerous possible 3D volumes of flow units,porosity and permeability below the seismic resolution limit,respecting the available seismic data at larger scales.The methodology is used in the Mero Field,a Brazilian presalt carbonate reservoir located in the Santos Basin,who presents a complex and heterogenic geological setting with different sedimentological processes and diagenetic history.We generated metric flow units using the conventional core analysis and transposed to the well log data.Then,given a Markov chain Monte Carlo algorithm,the seismic data and the well log statistics,we simulated acoustic impedance,decametric flow units(DFU),metric flow units(MFU),porosity and permeability volumes in the metric scale.The aim is to estimate a minimum amount of MFU able to calculate realistic scenarios porosity and permeability scenarios,without losing the seismic lateral control.In other words,every porosity and permeability volume simulated produces a synthetic seismic that match the real seismic of the area,even in the metric scale.The achieved 3D results represent a high-resolution fluid flow reservoir modelling considering the lateral control of the seismic during the process and can be directly incorporated in the dynamic characterization workflow.
基金Supported by the PetroChina Science and Technology Innovation Fund Project(2021DQ02-1003)Basic Research Project for Central Universities(2022JCCXDC02).
文摘Through core observation,thin section identification,X-ray diffraction analysis,scanning electron microscopy,and low-temperature nitrogen adsorption and isothermal adsorption experiments,the lithology and pore characteristics of the Upper Carboniferous bauxite series in eastern Ordos Basin were analyzed to reveal the formation and evolution process of the bauxite reservoirs.A petrological nomenclature and classification scheme for bauxitic rocks based on three units(aluminum hydroxides,iron minerals and clay minerals)is proposed.It is found that bauxitic mudstone is in the form of dense massive and clastic structures,while the(clayey)bauxite is of dense massive,pisolite,oolite,porous soil and clastic structures.Both bauxitic mudstone and bauxite reservoirs develop dissolution pores,intercrystalline pores,and microfractures as the dominant gas storage space,with the porosity less than 10% and mesopores in dominance.The bauxite series in the North China Craton can be divided into five sections,i.e.,ferrilite(Shanxi-style iron ore,section A),bauxitic mudstone(section B),bauxite(section C),bauxite mudstone(debris-containing,section D)and dark mudstone-coal section(section E).The burrow/funnel filling,lenticular,layered/massive bauxite deposits occur separately in the karst platforms,gentle slopes and low-lying areas.The karst platforms and gentle slopes are conducive to surface water leaching,with strong karstification,well-developed pores,large reservoir thickness and good physical properties,but poor strata continuity.The low-lying areas have poor physical properties but relatively continuous and stable reservoirs.The gas enrichment in bauxites is jointly controlled by source rock,reservoir rock and fractures.This recognition provides geological basis for the exploration and development of natural gas in the Upper Carboniferous in the study area and similar bauxite systems.
文摘This research uses both three-dimensional(3D)modeling and geologic well control to piece back together the architectural parts of the Late Cretaceous formations.The goal is to figure out the sizes,directions,locations,and controls of the layers of the fluvial sandstone reservoirs.Sequence stratigraphy is essential for 3D reservoir modeling and petroleum geology understanding in the Bahga oilfield.The purpose of this work is to create a static model that shows the layers and facies distribution in the reservoir interval.We will use data from nine well logs and 22 seismic lines calibrated by the Abu Roash G Member reservoir core intervals to accomplish this.The petrophysical study discovered three parts in the Abu Roash G Member reservoir rock:channel fill that is affected by tides,channel fill that is dominated by tides(intertidal sands),and channel top with lenticular bedded sandstone.The model's findings point to the existence of an NNW-oriented sand body,which could be a prime location to produce hydrocarbons.The original oil in place(OOIP)is about 3,438,279 Stock Tank Barrels(STB),and the oil reserve reaches up to 1,031,484(STB).Sequence stratigraphic analysis using seismic and well log information(SB)reveals that the Upper Cretaceous AR/G reservoir of the Bahga field is characterized by third-and fourth-order stratigraphic sequences,which are constrained by three Maximum Flooding Surfaces(MFS)and two Sequence Boundaries.The integration of the derived geological model and sequence stratigraphic results can lower future extraction risk by identifying the locations and trends of the geologic facies with the necessary petrophysical properties for the hydrocarbon accumulations.
文摘Using the typical characteristics of multi-layered marine and continental transitional gas reservoirs as a basis,a model is developed to predict the related well production rate.This model relies on the fractal theory of tortuous capillary bundles and can take into account multiple gas flow mechanisms at the micrometer and nanometer scales,as well as the flow characteristics in different types of thin layers(tight sandstone gas,shale gas,and coalbed gas).Moreover,a source-sink function concept and a pressure drop superposition principle are utilized to introduce a coupled flow model in the reservoir.A semi-analytical solution for the production rate is obtained using a matrix iteration method.A specific well is selected for fitting dynamic production data,and the calculation results show that the tight sandstone has the highest gas production per unit thickness compared with the other types of reservoirs.Moreover,desorption and diffusion of coalbed gas and shale gas can significantly contribute to gas production,and the daily production of these two gases decreases rapidly with decreasing reservoir pressure.Interestingly,the gas production from fractures exhibits an approximately U-shaped distribution,indicating the need to optimize the spacing between clusters during hydraulic fracturing to reduce the area of overlapping fracture control.The coal matrix water saturation significantly affects the coalbed gas production,with higher water saturation leading to lower production.
基金Equinor for financing the R&D projectthe Institute of Science and Technology of Petroleum Geophysics of Brazil for supporting this research。
文摘We apply stochastic seismic inversion and Bayesian facies classification for porosity modeling and igneous rock identification in the presalt interval of the Santos Basin. This integration of seismic and well-derived information enhances reservoir characterization. Stochastic inversion and Bayesian classification are powerful tools because they permit addressing the uncertainties in the model. We used the ES-MDA algorithm to achieve the realizations equivalent to the percentiles P10, P50, and P90 of acoustic impedance, a novel method for acoustic inversion in presalt. The facies were divided into five: reservoir 1,reservoir 2, tight carbonates, clayey rocks, and igneous rocks. To deal with the overlaps in acoustic impedance values of facies, we included geological information using a priori probability, indicating that structural highs are reservoir-dominated. To illustrate our approach, we conducted porosity modeling using facies-related rock-physics models for rock-physics inversion in an area with a well drilled in a coquina bank and evaluated the thickness and extension of an igneous intrusion near the carbonate-salt interface. The modeled porosity and the classified seismic facies are in good agreement with the ones observed in the wells. Notably, the coquinas bank presents an improvement in the porosity towards the top. The a priori probability model was crucial for limiting the clayey rocks to the structural lows. In Well B, the hit rate of the igneous rock in the three scenarios is higher than 60%, showing an excellent thickness-prediction capability.
基金supported by China Natural Science Foundation(Grant No.52274053)Beijing Natural Science Foundation(Grant No.3232028)Open Fund of State Key Laboratory of Offshore Oil Exploitation(Grant No.CCL2021RCPS0515KQN)。
文摘Spontaneous imbibition(SI)is an important mechanism for enhancing oil recovery in low-permeability reservoirs.Due to the strong heterogeneity,and the non-Darcy flow,the construction of SI model for lowpermeability reservoirs is extremely challenging.Commonly,traditional SI models based on single or averaged capillary tortuosity ignore the influence of heterogeneity of pore seepage channels and the threshold pressure(TP)on imbibition.Therefore,in this work,based on capillary model and fractal theory,a mathematical model of characterizing SI considering heterogeneity of pore seepage channels is established.On this basis,the threshold pressure was introduced to determine the pore radius at which the wetted phase can displace oil.The proposed new SI model was verified by imbibition experimental data.The study shows that for weakly heterogeneous cores with permeability of 0-1 m D,the traditional SI model can characterize the imbibition process relatively accurately,and the new imbibition model can increase the coefficient of determination by 1.05 times.However,traditional model has serious deviations in predicting the imbibition recovery for cores with permeability of 10-50 m D.The new SI model coupling with heterogeneity of pore seepage channels and threshold pressure effectively solves this problem,and the determination coefficient is increased from 0.344 to 0.922,which is increased by2.68 times.For low-permeability reservoirs,the production of the oil in transitional pores(0.01-0.1μm)and mesopores(0.1-1μm)significantly affects the imbibition recovery,as the research shows that when the heterogeneity of pore seepage channels is ignored,the oil recovery in transitional pores and mesopores decreases by 7.54%and 4.26%,respectively.Sensitivity analysis shows that increasing interfacial tension,decreasing contact angle,oil-water viscosity ratio and threshold pressure will increase imbibition recovery.In addition,there are critical values for the influence of these factors on the imbibition recovery,which provides theoretical support for surfactant optimization.
基金supported by the National Natural Science Foundation of China(Grant No.42025403)the Youth Innovation Promotion Association,Chinese Academy of Sciences(Grant No.2023074).
文摘From June 2008 to August 2013,approximately 67 kt of CO_(2) was injected into a deep saline formation at the Ketzin pilot CO_(2) storage site.During injection,3D seismic surveys have been performed to monitor the migration of sequestered CO_(2).Seismic monitoring results are limited by the acquisition and signal-to-noise ratio of the acquired data.The multiphysical reservoir simulation provides information regarding the CO_(2) fluid behavior,and the approximated model should be calibrated with the monitoring results.In this work,property models are delivered from the multiphysical model during 3D repeated seismic surveys.The simulated seismic data based on the models are compared with the real data,and the results validate the effectiveness of the multiphysical inversion method.Time-lapse analysis shows the trend of CO_(2) migration during and after injection.
基金supported by the National Natural Science Foundation of China(No.52174038 and No.52004307)China Petroleum Science and Technology Project-Major Project-Research on Tight Oil-Shale Oil Reservoir Engineering Methods and Key Technologies in Ordos Basin(No.ZLZX2020-02-04)Science Foundation of China University of Petroleum,Beijing(No.2462018YJRC015)。
文摘Under the policy background and advocacy of carbon capture,utilization,and storage(CCUS),CO_(2)-EOR has become a promising direction in the shale oil reservoir industry.The multi-scale pore structure distribution and fracture structure lead to complex multiphase flow,comprehensively considering multiple mechanisms is crucial for development and CO_(2) storage in fractured shale reservoirs.In this paper,a multi-mechanism coupled model is developed by MATLAB.Compared to the traditional Eclipse300 and MATLAB Reservoir Simulation Toolbox(MRST),this model considers the impact of pore structure on fluid phase behavior by the modified Peng—Robinson equation of state(PR-EOS),and the effect simultaneously radiate to Maxwell—Stefan(M—S)diffusion,stress sensitivity,the nano-confinement(NC)effect.Moreover,a modified embedded discrete fracture model(EDFM)is used to model the complex fractures,which optimizes connection types and half-transmissibility calculation approaches between non-neighboring connections(NNCs).The full implicit equation adopts the finite volume method(FVM)and Newton—Raphson iteration for discretization and solution.The model verification with the Eclipse300 and MRST is satisfactory.The results show that the interaction between the mechanisms significantly affects the production performance and storage characteristics.The effect of molecular diffusion may be overestimated in oil-dominated(liquid-dominated)shale reservoirs.The well spacing and injection gas rate are the most crucial factors affecting the production by sensitivity analysis.Moreover,the potential gas invasion risk is mentioned.This model provides a reliable theoretical basis for CO_(2)-EOR and sequestration in shale oil reservoirs.
基金Supported by the PetroChina Science and Technology Major Project(2016E0201)。
文摘Based on the geochemical,seismic,logging and drilling data,the Fuyu reservoirs of the Lower Cretaceous Quantou Formation in northern Songliao Basin are systematically studied in terms of the geological characteristics,the tight oil enrichment model and its major controlling factors.First,the Quantou Formation is overlaid by high-quality source rocks of the Upper Cretaceous Qingshankou Formation,with the development of nose structure around sag and the broad and continuous distribution of sand bodies.The reservoirs are tight on the whole.Second,the configuration of multiple elements,such as high-quality source rocks,reservoir rocks,fault,overpressure and structure,controls the tight oil enrichment in the Fuyu reservoirs.The source-reservoir combination controls the tight oil distribution pattern.The pressure difference between source and reservoir drives the charging of tight oil.The fault-sandbody transport system determines the migration and accumulation of oil and gas.The positive structure is the favorable place for tight oil enrichment,and the fault-horst zone is the key part of syncline area for tight oil exploration.Third,based on the source-reservoir relationship,transport mode,accumulation dynamics and other elements,three tight oil enrichment models are recognized in the Fuyu reservoirs:(1)vertical or lateral migration of hydrocarbon from source rocks to adjacent reservoir rocks,that is,driven by overpressure,hydrocarbon generated is migrated vertically or laterally to and accumulates in the adjacent reservoir rocks;(2)transport of hydrocarbon through faults between separated source and reservoirs,that is,driven by overpressure,hydrocarbon migrates downward through faults to the sandbodies that are separated from the source rocks;and(3)migration of hydrocarbon through faults and sandbodies between separated source and reservoirs,that is,driven by overpressure,hydrocarbon migrates downwards through faults to the reservoir rocks that are separated from the source rocks,and then migrates laterally through sandbodies.Fourth,the differences in oil source conditions,charging drive,fault distribution,sandbody and reservoir physical properties cause the differential enrichment of tight oil in the Fuyu reservoirs.Comprehensive analysis suggests that the Fuyu reservoir in the Qijia-Gulong Sag has good conditions for tight oil enrichment and has been less explored,and it is an important new zone for tight oil exploration in the future.
基金sponsored by the NSFC(41104066)973 Program of China(No.2014CB239006)+1 种基金NSTMP of China(Nos.2011ZX05004-003 and 2011ZX05029-003)12th 5-Year Basic Research Program of CNPC(No.2011A-3601)
文摘In heterogeneous natural gas reservoirs, gas is generally present as small patchlike pockets embedded in the water-saturated host matrix. This type of heterogeneity, also called "patchy saturation", causes significant seismic velocity dispersion and attenuation. To establish the relation between seismic response and type of fluids, we designed a rock physics model for carbonates. First, we performed CT scanning and analysis of the fluid distribution in the partially saturated rocks. Then, we predicted the quantitative relation between the wave response at different frequency ranges and the basic lithological properties and pore fluids. A rock physics template was constructed based on thin section analysis of pore structures and seismic inversion. This approach was applied to the limestone gas reservoirs of the right bank block of the Amu Darya River. Based on poststack wave impedance and prestack elastic parameter inversions, the seismic data were used to estimate rock porosity and gas saturation. The model results were in good agreement with the production regime of the wells.
基金This project is sponsored by National Natural Science Foundation of China, No. 40574030.
文摘Based on the percolation network model characterizing reservoir rock's pore structure and fluid characteristics, this paper qualitatively studies the effects of pore size, pore shape, pore connectivity, and the amount of micropores on the I - Sw curve using numerical modeling. The effects of formation water salinity on the electrical resistivity of the rock are discussed. Then the relative magnitudes of the different influencing factors are discussed. The effects of the different factors on the I - Sw curve are analyzed by fitting simulation results. The results show that the connectivity of the void spaces and the amount of micropores have a large effect on the I - S, curve, while the other factors have little effect. The formation water salinity has a large effect on the absolute resistivity values. The non-Archie phenomenon is prevalent, which is remarkable in rocks with low permeability.
基金the National Key R&D Program of China(No.2019YFB1504102).
文摘Geothermal energy is a kind of renewable,sustainable and clean energy resource.Geothermal energy is abundant in carbonate reservoirs.However,low matrix permeability limits its exploitation.The super-critical carbon dioxide(SC-CO_(2))jet fracturing is expected to efficiently stimulate the carbonate geothermal reservoirs and achieve the storage of CO_(2) simultaneously.In this paper,we established a transient seepage and fluid-thermo-mechanical coupled model to analyze the impact performance of sc-CO_(2) jet fracturing.The mesh-based parallel code coupling interface was employed to couple the fluid and solid domains by exchanging the data through the mesh interface.The physical properties change of sC-CO_(2) with temperature were considered in the numerical model.Results showed that SC-CO_(2) jet frac-turing is superior to water-jet fracturing with respect to jetting velocity,particle trajectory and pene-trability.Besides,stress distribution on the carbonate rock showed that the tensile and shear failure would more easily occur by SC-CO_(2) jet than that by water jet.Moreover,pressure and temperature control the jet field and seepage field of sC-CO_(2) simultaneously.Increasing the jet temperature can effectively enhance the impingement effect and seepage process by decreasing the viscosity and density of SC-CO_(2).The key findings are expected to provide a theoretical basis and design reference for applying SC-CO_(2) jet fracturing in carbonate geothermal reservoirs.
基金RPSEA and U.S.Department of Energy for partially funding this study
文摘Hydrocarbon production from shale has attracted much attention in the recent years. When applied to this prolific and hydrocarbon rich resource plays, our understanding of the complexities of the flow mechanism(sorption process and flow behavior in complex fracture systems- induced or natural) leaves much to be desired. In this paper, we present and discuss a novel approach to modeling, history matching of hydrocarbon production from a Marcellus shale asset in southwestern Pennsylvania using advanced data mining, pattern recognition and machine learning technologies. In this new approach instead of imposing our understanding of the flow mechanism, the impact of multi-stage hydraulic fractures, and the production process on the reservoir model, we allow the production history, well log, completion and hydraulic fracturing data to guide our model and determine its behavior. The uniqueness of this technology is that it incorporates the so-called "hard data" directly into the reservoir model, so that the model can be used to optimize the hydraulic fracture process. The "hard data" refers to field measurements during the hydraulic fracturing process such as fluid and proppant type and amount, injection pressure and rate as well as proppant concentration. This novel approach contrasts with the current industry focus on the use of "soft data"(non-measured, interpretive data such as frac length, width,height and conductivity) in the reservoir models. The study focuses on a Marcellus shale asset that includes 135 wells with multiple pads, different landing targets, well length and reservoir properties. The full field history matching process was successfully completed using this data driven approach thus capturing the production behavior with acceptable accuracy for individual wells and for the entire asset.
基金supported by Open Fund (PLC201203) of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation (Chengdu University of Technology)Major Project of Education Department in Sichuan Province (13ZA0177)
文摘Due to inherent limits of data acquisition and geophysical data resolution, there are large uncertainties in the characterization of subsurface fractures. However, outcrop analogies can provide qualitative and quantitative information on a large number of fractures, based on which the accuracy of subsurface fracture characterization can be improved. Here we take the tectonic fracture modeling of an ultra-low permeability sandstone reservoir based on an outcrop analogy, a case study of the Chang6t~ Formation of the Upper Triassic Yanchang Group of the Wangyao Oilfield in the Ordos Basin of China. An outcrop at the edge of the basin is a suitable analog for the reservoir, but the prerequisite is that they must have equivalent previous stress fields, similar final structural characteristics, relative timing and an identical depositional environment and diagenesis. The relationship among fracture density, rock type and bed thickness based on the outcrop is one of the most important fracture distribution models, and can be used to interpret fracture density in individual wells quantitatively. Fracture orientation, dip, geometry and scale, also should be described and measured in the outcrop, and can be used together with structure restoration and single well fracture density interpretation to guide fracture intensity prediction on bed surfaces and to constrain the construction of the 3D fracture geometry model of the subsurface reservoir. The application of the above principles shows the outcrop-based tectonic fracture models of the target ultra-low permeability sandstone reservoir are consistent with fractures inferred from microseismic interpretation and tracer tests. This illustrated that the fracture modeling based on the outcrop analogy is reliable and can reduce the uncertainty in stochastic fracture modeling.
文摘This paper discusses the principles of geologic constraints on reservoir stochastic modeling. By using the system science theory, two kinds of uncertainties, including random uncertainty and fuzzy uncertainty, are recognized. In order to improve the precision of stochastic modeling and reduce the uncertainty in realization, the fuzzy uncertainty should be stressed, and the "geological genesis-controlled modeling" is conducted under the guidance of a quantitative geological pattern. An example of the Pingqiao horizontal-well division of the Ansai Oilfield in the Ordos Basin is taken to expound the method of stochastic modeling.
基金co-supported by the National Basic Research Program of China(Grant No.2011CB201103)the National Science and Technology Major Project(GrantNo.2011ZX05004003)
文摘Fractured reservoirs are an important target for oil and gas exploration in the Tarim Basin and the prediction of this type of reservoir is challenging.Due to the complicated fracture system in the Tarim Basin,the conventional AVO inversion method based on HTI theory to predict fracture development will result in some errors.Thus,an integrated research concept for fractured reservoir prediction is put forward in this paper.Seismic modeling plays a bridging role in this concept,and the establishment of an anisotropic fracture model by Discrete Fracture Network (DFN) is the key part.Because the fracture system in the Tarim Basin shows complex anisotropic characteristics,it is vital to build an effective anisotropic model.Based on geological,well logging and seismic data,an effective anisotropic model of complex fracture systems can be set up with the DFN method.The effective elastic coefficients,and the input data for seismic modeling can be calculated.Then seismic modeling based on this model is performed,and the seismic response characteristics are analyzed.The modeling results can be used in the following AVO inversion for fracture detection.
基金funded by a Major Science and Technology Project of China National Petroleum Corporation(CNPC)entitled Research on Key Technologies for Efficient Production of Overseas Large Carbonate Reservoir”(2023ZZ19-08).
文摘Porous carbonate reservoirs,prevalent in the Middle East,are lithologically dominated by bioclastic limestones,exhibiting high porosity,low permeability,intricate pore structure,and strong heterogeneity.Waterflooding through horizontal wells is commonly used for exploiting these reservoirs.However,challenges persist,such as significant uncertainty and complex operational procedures regarding adjustment effects during the exploitation.The USH reservoir of the Cretaceous D oilfield,Oman exemplifies typical porous carbonate reservoirs.It initially underwent depletion drive using vertical wells,followed by horizontal well waterflooding in the late stage.Currently,the oilfield is confronted with substantial developmental challenges,involving the understanding of residual oil distribution,effective water cut control,and sustaining oil production since it has entered the late development stage.Employing a microscopic visualization displacement system equipped with electrodes,this study elucidated the waterflooding mechanisms and residual oil distribution in the late-stage development of the USH reservoir.The results reveal that the reservoir's vertical stacking patterns act as the main factor affecting the horizontal well waterflooding efficacy.Distinct water flow channels emerge under varying reservoir stacking patterns,with post-waterflooding residual oil predominantly distributed at the reservoir's top and bottom.The oil recovery can be enhanced by adjusting the waterflooding's flow line and intensity.The findings of this study will provide theoretical insights of waterflooding mechanisms and injection-production adjustments for exploiting other porous carbonate reservoirs in the Middle East through horizontal wells.