Geological conditions and reservoir characteristics of various shales in major shalehosted regions of China

China is home to shales of three facies:Marine shale,continental shale,and marine-continental transitional shale.Different types of shale gas are associated with significantly different formation conditions and major controlling factors.This study compared the geological characteristics of various shales and analyzed the influences of different parameters on the formation and accumulation of shale gas.In general,shales in China’s several regions exhibit high total organic carbon(TOC)contents,which lays a sound material basis for shale gas generation.Marine strata generally show high degrees of thermal evolution.In contrast,continental shales manifest low degrees of thermal evolution,necessitating focusing on areas with relatively high degrees of thermal evolution in the process of shale gas surveys for these shales.The shales of the Wufeng and Silurian formations constitute the most favorable shale gas reservoirs since they exhibit the highest porosity among the three types of shales.These shales are followed by those in the Niutitang and Longtan formations.In contrast,the shales of the Doushantuo,Yanchang,and Qingshankou formations manifest low porosities.Furthermore,the shales of the Wufeng and Longmaxi formations exhibit high brittle mineral contents.Despite a low siliceous mineral content,the shales of the Doushantuo Formation feature a high carbonate mineral content,which can increase the shales’brittleness to some extent.For marine-continental transitional shales,where thin interbeds of tight sandstone with unequal thicknesses are generally found,it is recommended that fracturing combined with drainage of multiple sets of lithologic strata should be employed to enhance their shale gas production.

A Mathematical Modeling of 3D Cubical Geometry Hypothetical Reservoir under the Effect of Nanoparticles Flow Rate,Porosity,and Relative Permeability

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.

Modularized and Parametric Modeling Technology for Finite Element Simulations of Underground Engineering under Complicated Geological Conditions

The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.

Reactive transport modeling constraints on the complex genesis of a lacustrine dolomite reservoir:A case from the Eocene Qaidam Basin,China

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.

Changes in filtration and capacitance properties of highly porous reservoir in underground gas storage:CT-based and geomechanical modeling

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.

Reservoir Stochastic Modeling Constrained by Quantitative Geological Conceptual Patterns

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.

Thermo-hydro-poro-mechanical responses of a reservoir-induced landslide tracked by high-resolution fiber optic sensing nerves

Thermo-poro-mechanical responses along sliding zone/surface have been extensively studied.However,it has not been recognized that the potential contribution of other crucial engineering geological interfaces beyond the slip surface to progressive failure.Here,we aim to investigate the subsurface multiphysics of reservoir landslides under two extreme hydrologic conditions(i.e.wet and dry),particularly within sliding masses.Based on ultra-weak fiber Bragg grating(UWFBG)technology,we employ specialpurpose fiber optic sensing cables that can be implanted into boreholes as“nerves of the Earth”to collect data on soil temperature,water content,pore water pressure,and strain.The Xinpu landslide in the middle reach of the Three Gorges Reservoir Area in China was selected as a case study to establish a paradigm for in situ thermo-hydro-poro-mechanical monitoring.These UWFBG-based sensing cables were vertically buried in a 31 m-deep borehole at the foot of the landslide,with a resolution of 1 m except for the pressure sensor.We reported field measurements covering the period 2021 and 2022 and produced the spatiotemporal profiles throughout the borehole.Results show that wet years are more likely to motivate landslide motions than dry years.The annual thermally active layer of the landslide has a critical depth of roughly 9 m and might move downward in warmer years.The dynamic groundwater table is located at depths of 9e15 m,where the peaked strain undergoes a periodical response of leap and withdrawal to annual hydrometeorological cycles.These interface behaviors may support the interpretation of the contribution of reservoir regulation to slope stability,allowing us to correlate them to local damage events and potential global destabilization.This paper also offers a natural framework for interpreting thermo-hydro-poro-mechanical signatures from creeping reservoir bank slopes,which may form the basis for a landslide monitoring and early warning system.

Projecting Spring Consecutive Rainfall Events in the Three Gorges Reservoir Based on Triple-Nested Dynamical Downscaling

Spring consecutive rainfall events(CREs) are key triggers of geological hazards in the Three Gorges Reservoir area(TGR), China. However, previous projections of CREs based on the direct outputs of global climate models(GCMs) are subject to considerable uncertainties, largely caused by their coarse resolution. This study applies a triple-nested WRF(Weather Research and Forecasting) model dynamical downscaling, driven by a GCM, MIROC6(Model for Interdisciplinary Research on Climate, version 6), to improve the historical simulation and reduce the uncertainties in the future projection of CREs in the TGR. Results indicate that WRF has better performances in reproducing the observed rainfall in terms of the daily probability distribution, monthly evolution and duration of rainfall events, demonstrating the ability of WRF in simulating CREs. Thus, the triple-nested WRF is applied to project the future changes of CREs under the middle-of-the-road and fossil-fueled development scenarios. It is indicated that light and moderate rainfall and the duration of continuous rainfall spells will decrease in the TGR, leading to a decrease in the frequency of CREs. Meanwhile, the duration, rainfall amount, and intensity of CREs is projected to regional increase in the central-west TGR. These results are inconsistent with the raw projection of MIROC6. Observational diagnosis implies that CREs are mainly contributed by the vertical moisture advection. Such a synoptic contribution is captured well by WRF, which is not the case in MIROC6,indicating larger uncertainties in the CREs projected by MIROC6.

Geostatistical seismic inversion and 3D modelling of metric flow units,porosity and permeability in Brazilian presalt reservoir

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.

Construction of carbonate reservoir knowledge base and its application in fracture-cavity reservoir geological modeling

To improve the efficiency and accuracy of carbonate reservoir research,a unified reservoir knowledge base linking geological knowledge management with reservoir research is proposed.The reservoir knowledge base serves high-quality analysis,evaluation,description and geological modeling of reservoirs.The knowledge framework is divided into three categories:technical service standard,technical research method and professional knowledge and cases related to geological objects.In order to build a knowledge base,first of all,it is necessary to form a knowledge classification system and knowledge description standards;secondly,to sort out theoretical understandings and various technical methods for different geologic objects and work out a technical service standard package according to the technical standard;thirdly,to collect typical outcrop and reservoir cases,constantly expand the content of the knowledge base through systematic extraction,sorting and saving,and construct professional knowledge about geological objects.Through the use of encyclopedia based collaborative editing architecture,knowledge construction and sharing can be realized.Geological objects and related attribute parameters can be automatically extracted by using natural language processing(NLP)technology,and outcrop data can be collected by using modern fine measurement technology,to enhance the efficiency of knowledge acquisition,extraction and sorting.In this paper,the geological modeling of fracture-cavity reservoir in the Tarim Basin is taken as an example to illustrate the construction of knowledge base of carbonate reservoir and its application in geological modeling of fracture-cavity carbonate reservoir.

Reservoir characteristics and formation model of Upper Carboniferous bauxite series in eastern Ordos Basin,NW China

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.

Reservoir geological modeling and significance of Cambrian Xiaoerblak Formation in Keping outcrop area, Tarim Basin, NW China

Take the Cambrian Xiaoerblak Formation in the Keping(Kalpin) outcrop area as an example, a 28 km reservoir scale geological model was built based on description of 7 profiles, observation of more than 1000 thin sections, petrophysical analysis of 556 samples and many geochemical tests. The Xiaoerblak Formation, 158–178 m thick, is divided into three members and 5 submembers, and is composed of laminated microbialite dolomite(LMD), thrombolite dolomite(TD), foamy-stromatolite dolomite(FSD), oncolite dolomite(OD), grain dolomite(GD)/crystalline dolomite with grain ghost and micritic dolomite(MD)/argillaceous dolomite. The petrology features show that its sediment sequence is micro-organism layer – microbial mound/shoal – tidal flat in carbonate ramp background from bottom up. The reservoir has 5 types of pores, namely, framework pore, dissolved vug, intergranular and intragranular dissolved pore and intercrystalline dissolved pore, as main reservoir space. It is found that the development of pore has high lithofacies selectivity, FSD has the highest average porosity, TD, OD and GD come second. The reservoir is pore-vug reservoir with medium-high porosity and medium-low permeability. The dolomite of Xiaoerblak Formation was formed in para-syngenetic to early diagenetic stage through dolomitization caused by seawater. The reservoir development is jointly controlled by sedimentary facies, micro-organism type, high frequency sequence interface and early dolomitization. The classⅠand Ⅱ reservoirs, with an average thickness of 41.2 m and average reservoir-stratum ratio of about 25.6%, have significant potential. It is predicted that the microbial mounds and shoals in the middle ramp around the ancient uplift are the favorable zones for reservoir development.

GeoNER:Geological Named Entity Recognition with Enriched Domain Pre-Training Model and Adversarial Training

As important geological data,a geological report contains rich expert and geological knowledge,but the challenge facing current research into geological knowledge extraction and mining is how to render accurate understanding of geological reports guided by domain knowledge.While generic named entity recognition models/tools can be utilized for the processing of geoscience reports/documents,their effectiveness is hampered by a dearth of domain-specific knowledge,which in turn leads to a pronounced decline in recognition accuracy.This study summarizes six types of typical geological entities,with reference to the ontological system of geological domains and builds a high quality corpus for the task of geological named entity recognition(GNER).In addition,Geo Wo BERT-adv BGP(Geological Word-base BERTadversarial training Bi-directional Long Short-Term Memory Global Pointer)is proposed to address the issues of ambiguity,diversity and nested entities for the geological entities.The model first uses the fine-tuned word granularitybased pre-training model Geo Wo BERT(Geological Word-base BERT)and combines the text features that are extracted using the Bi LSTM(Bi-directional Long Short-Term Memory),followed by an adversarial training algorithm to improve the robustness of the model and enhance its resistance to interference,the decoding finally being performed using a global association pointer algorithm.The experimental results show that the proposed model for the constructed dataset achieves high performance and is capable of mining the rich geological information.

Recent advances in reservoir characterization:Introduction

Reservoir characterization refers to the process of creating a comprehensive model that characterizes the reservoir based on its ability to store and produce hydrocarbons.This includes analyzing reservoir fluid behavior under various conditions and identifying optimal production techniques to maximize hydrocarbon recovery.For this,a holistic understanding is required that integrates data from geophysics,geostatistics,petrophysics,geology,and reservoir engineering.

A Well Productivity Model for Multi-Layered Marine and Continental Transitional Reservoirs with Complex Fracture Networks

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.

Stochastic seismic inversion and Bayesian facies classification applied to porosity modeling and igneous rock identification

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.

Fractal model of spontaneous imbibition in low-permeability reservoirs coupled with heterogeneity of pore seepage channels and threshold pressure

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.

Understanding poromechanical response of a biogenic coalbed methane reservoir

Biogenic coalbed methane(BCBM)reservoirs aim to produce methane from in situ coal deposits following microbial conversion of coal.Success of BCBM reservoirs requires economic methane production within an acceptable timeframe.The work reported here quantifies the findings of previously published qualitative work,where it was found that bioconversion induces strains in the pore,matrix and bulk scales.Using imaging and dynamic strain monitoring techniques,the bioconversion induced strain is quantified here.To understand the effect of these strains from a reservoir geomechanics perspective,a corresponding poromechanical model is developed.Furthermore,findings of imaging experiments are validated using core-flooding flow experiments.Finally,expected field-scale behavior of the permeability response of a BCBM operation is modeled and analyzed.The results of the study indicated that,for Illinois coals,bioconversion induced strains result in a decrease in fracture porosity,resulting in a detrimental permeability drop in excess of 60%during bioconversion,which festers itself exponentially throughout its producing life.Results indicate that reservoirs with high initial permeability that will support higher Darcian flowrates,would be better suited for coal bioconversion,thereby providing a site-selection criteria for BCBM operations.

Effects of acid-rock reaction on physical properties during CO_(2)-rich industrial waste gas(CO_(2)-rich IWG)injection in shale reservoirs

"Carbon peaking and carbon neutrality"is an essential national strategy,and the geological storage and utilization of CO_(2)is a hot issue today.However,due to the scarcity of pure CO_(2)gas sources in China and the high cost of CO_(2)capture,CO_(2)-rich industrial waste gas(CO_(2)-rich IWG)is gradually emerging into the public's gaze.CO_(2)has good adsorption properties on shale surfaces,but acidic gases can react with shale,so the mechanism of the CO_(2)-rich IWG-water-shale reaction and the change in reservoir properties will determine the stability of geological storage.Therefore,based on the mineral composition of the Longmaxi Formation shale,this study constructs a thermodynamic equilibrium model of water-rock reactions and simulates the regularity of reactions between CO_(2)-rich IWG and shale minerals.The results indicate that CO_(2)consumed 12%after reaction,and impurity gases in the CO_(2)-rich IWG can be dissolved entirely,thus demonstrating the feasibility of treating IWG through water-rock reactions.Since IWG inhibits the dissolution of CO_(2),the optimal composition of CO_(2)-rich IWG is 95%CO_(2)and 5%IWG when CO_(2)geological storage is the main goal.In contrast,when the main goal is the geological storage of total CO_(2)-rich IWG or impurity gas,the optimal CO_(2)-rich IWG composition is 50%CO_(2)and 50%IWG.In the CO_(2)-rich IWG-water-shale reaction,temperature has less influence on the water-rock reaction,while pressure is the most important parameter.SO2 has the greatest impact on water-rock reaction in gas.For minerals,clay minerals such as illite and montmorillonite had a significant effect on water-rock reaction.The overall reaction is dominated by precipitation and the volume of the rock skeleton has increased by 0.74 cm3,resulting in a decrease in shale porosity,which enhances the stability of CO_(2)geological storage to some extent.During the reaction between CO_(2)-rich IWG-water-shale at simulated temperatures and pressures,precipitation is the main reaction,and shale porosity decreases.However,as the reservoir water content increases,the reaction will first dissolve and then precipitate before dissolving again.When the water content is less than 0.0005 kg or greater than 0.4 kg,it will lead to an increase in reservoir porosity,which ultimately reduces the long-term geological storage stability of CO_(2)-rich IWG.

3D seismic forward modeling from the multiphysical inversion at the Ketzin CO_(2) storage site

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.