Inversion of Reservoir Parameters by Ultra-High Temperature and High Pressure CO2 Gas Reservoirs with Nuclear Logging

2021 ◽  
Author(s):  
Hengrong Zhang ◽  
Lizhi Xiao ◽  
Wensheng Wu ◽  
Xinyue Fu ◽  
Shenglin He
Keyword(s):  
Joint Inversion ◽  
Simulation Method ◽  
Inversion Method ◽  
Burial Depth ◽  
Gas Reservoirs ◽  
Co2 Gas ◽  
Carbonate Cement ◽  
Response Equation ◽  
Temperature And Pressure ◽  
The Difference

Abstract The Yinggehai basin is located in the western part of the South China Sea, the burial depth of the Huangliu and Meishan formations in the target layer is close to 4000 meters, the formation temperature is close to 200 degrees Celsius, and the formation pressure is up to 100 MPa. The reservoir is characterized by low porosity-ultra-low permeability, heavy carbonate cement, complex CO2 content, this leads to complex neutron and density logging effects. The solubility of CO2 Above CH4, the solubility change with temperature and pressure is different from CH4, which makes it difficult to identify the CO2 gas layer. In this paper, based on the difference in the physical characteristics of CO2 and CH4, the Boltzmann equation combined with MCNP software was used to simulate the neutron and density logging responses under different CO2 saturations. Environmental factors such as temperature and pressure, carbonate cement, mud content and pores were studied To measure the effect of logging response, the LM inversion method is used to jointly invert CO2 saturation of density and neutron logs. The purpose of the inversion is to reduce the non-uniqueness of the evaluation of porosity and CO2 saturation. By introducing the Levenberg-Marquardt (LM) method, the neutron logging response equation of the porosity, argillaceous content, CO2, CH4 in the rock and the corresponding temperature and pressure is solved, and also the response equation of above parameters to density logging, where porosity and CO2 content are the key parameters, and the calculation results prove the effectiveness of the method by comparing the sampling data. The results show that the accuracy of the estimated CO2 saturation is increased by 10% compared with the conventional interpretation method, and the new simulation method improves the calculation speed several times compared to the MCNP software. The joint inversion method has been successfully applied to field data, which has greatly improved the saturation evaluation results of traditional logging interpretation methods, can be extended to other fields of nuclear logging simulation and inversion.

scholarly journals Imaging the subsurface using induced seismicity and ambient noise: 3-D tomographic Monte Carlo joint inversion of earthquake body wave traveltimes and surface wave dispersion

2020 ◽  
Vol 222 (3) ◽  
pp. 1639-1655
Author(s):  
Xin Zhang ◽  
Corinna Roy ◽  
Andrew Curtis ◽  
Andy Nowacki ◽  
Brian Baptie
Keyword(s):  
Surface Wave ◽  
Ambient Noise ◽  
Velocity Structure ◽  
Joint Inversion ◽  
Body Wave ◽  
Source Parameters ◽  
Wave Dispersion ◽  
Inversion Method ◽  
Surface Wave Dispersion ◽  
Wave Data

SUMMARY Seismic body wave traveltime tomography and surface wave dispersion tomography have been used widely to characterize earthquakes and to study the subsurface structure of the Earth. Since these types of problem are often significantly non-linear and have non-unique solutions, Markov chain Monte Carlo methods have been used to find probabilistic solutions. Body and surface wave data are usually inverted separately to produce independent velocity models. However, body wave tomography is generally sensitive to structure around the subvolume in which earthquakes occur and produces limited resolution in the shallower Earth, whereas surface wave tomography is often sensitive to shallower structure. To better estimate subsurface properties, we therefore jointly invert for the seismic velocity structure and earthquake locations using body and surface wave data simultaneously. We apply the new joint inversion method to a mining site in the United Kingdom at which induced seismicity occurred and was recorded on a small local network of stations, and where ambient noise recordings are available from the same stations. The ambient noise is processed to obtain inter-receiver surface wave dispersion measurements which are inverted jointly with body wave arrival times from local earthquakes. The results show that by using both types of data, the earthquake source parameters and the velocity structure can be better constrained than in independent inversions. To further understand and interpret the results, we conduct synthetic tests to compare the results from body wave inversion and joint inversion. The results show that trade-offs between source parameters and velocities appear to bias results if only body wave data are used, but this issue is largely resolved by using the joint inversion method. Thus the use of ambient seismic noise and our fully non-linear inversion provides a valuable, improved method to image the subsurface velocity and seismicity.

Origin and distribution of dolomite in Permian Rotliegend siliciclastic sandstones (Dutch Southern Permian Basin)

2019 ◽  
Vol 89 (10) ◽  
pp. 1055-1073 ◽  
Author(s):  
Nicolaas Molenaar ◽  
Marita Felder
Keyword(s):  
Large Scale ◽  
Permian Basin ◽  
Advective Flow ◽  
Gas Reservoirs ◽  
Carbonate Cement ◽  
Thin Section Petrography ◽  
The Uk ◽  
Pressure Dissolution ◽  
Dolomite Cement ◽  
Detrital Carbonate

ABSTRACT Dolomite is a common and volumetrically important mineral in many siliciclastic sandstones, including Permian Rotliegend sandstones (the Slochteren Formation). These sandstones form extensive gas reservoirs in the Southern Permian Basin in the Netherlands, Germany, Poland, and the UK. The reservoir quality of these sandstones is negatively influenced by the content and distribution of dolomite. The origin and the stratigraphic distribution of the dolomite is not yet fully understood. The aim of this study is to identify the origin of carbonate. The main methods used to achieve those aims are a combination of thin-section petrography, scanning electron microscopy (SEM and EDX), and XRD analyses. The present study shows that the typical dispersed occurrence of the dolomite is a consequence of dispersed detrital carbonate grains that served both as nuclei and source for authigenic dolomite cement. The dolomite cement formed syntaxial outgrowths and overgrowths around detrital carbonate grains. The study also shows that dolomite cement, often in combination with ankerite and siderite, precipitated during burial after mechanical compaction. Most of the carbonate grains consisted of dolomite before deposition. The carbonate grains were affected by compaction and pressure dissolution, and commonly have no well-defined outlines anymore. The distribution of dolomite cement in the Rotliegend sandstones was controlled by the presence of stable carbonate grains. Due to the restricted and variable content of carbonate grains and their dispersed occurrence, the cement is also dispersed and the degree of cementation heterogeneous. Our findings have important implications on diagenesis modeling. The presence of detrital carbonate excludes the need for external supply by any large-scale advective flow of diagenetic fluids. By knowing that the carbonate source is local and related to detrital grains instead of being externally derived from an unknown source, the presence of carbonate cement can be linked to a paleogeographic and sedimentological model.

Self-diffusion of lignite/water under different temperatures and pressure: A molecular dynamics study

2016 ◽  
Vol 30 (01) ◽  
pp. 1550253 ◽  
Author(s):  
Xinjian Liu ◽  
Yu Jin ◽  
Congliang Huang ◽  
Jingfeng He ◽  
Zhonghao Rao ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Water System ◽  
Simulation Method ◽  
The Self ◽  
Dynamics Simulation ◽  
Low Rank ◽  
Self Diffusion ◽  
Change Mechanism ◽  
Different Temperatures ◽  
Temperature And Pressure

Temperature and pressure have direct and remarkable implications for drying and dewatering effect of low rank coals such as lignite. To understand the microenergy change mechanism of lignite, the molecular dynamics simulation method was performed to study the self-diffusion of lignite/water under different temperatures and pressure. The results showed that high temperature and high pressure can promote the diffusion of lignite/water system, which facilitates the drying and dewatering of lignite. The volume and density of lignite/water system will increase and decrease with temperature increasing, respectively. Though the pressure within simulation range can make lignite density increase, the increasing pressure showed a weak impact on variation of density.

scholarly journals Simulation method for determining traction power of ATN–PRT vehicle

Transport ◽  
2016 ◽  
Vol 33 (2) ◽  
pp. 335-343 ◽  
Author(s):  
Maciej Kozłowski
Keyword(s):  
Aerodynamic Force ◽  
Dynamic Properties ◽  
Simulation Method ◽  
Drive System ◽  
Constant Speed ◽  
Technical Parameters ◽  
Personal Rapid Transit ◽  
Motion Resistance ◽  
The Difference ◽  
Traction Power

The construction of Personal Rapid Transit (PRT) vehicle made within the framework of Eco-Mobility project has been described in the present paper. Key features of the vehicles were identified – e.g. drive with three-phase linear motor with winding on the vehicle and fixed rotor in the road surface, contactless dynamic vehicle powering. Attention was paid to the difference in dynamic properties compared to rail vehicles, related to the lack of the so-called ‘centering mechanism’. A development of a nominal model for the analysis of vehicle drive properties was presented. Results of simulation studies were presented for a vehicle with running-drive system construction, planned for implementation in the city of Rzeszów (Poland). While discussing the problems of building a PRT system, there was a focus on the issue of determining power and traction of the vehicle. A methodology for determining the power and traction energy consumption of the vehicle was presented for assumed conditions of travel on road segments. Input values for the calculation of power are variables describing the curvature (or bends radii) of paths of movement between stops and the course of the current speed. Output values are total traction power or traction energy (where ‘traction’ refers to the power or mechanical work of drive forces). Three basic elements of traction power were isolated: the power of kinetic energy (for acceleration/delay of vehicle movement) basic (to offset the aerodynamic force of motion resistance at constant speed) and additional losses (to offset additional motion resistance forces operating in turns at constant speed). Due to the lack of vehicle prototypes with assumed structure, it was proposed that these components are determined via simulation. The presented results relate to the calculation of demand for power and energy for the planned test section. The scope of further work was indicated: determining the required traction characteristics of electric drive, selecting the best values for supercapacitor’s capacity in the drive system, determining the technical parameters of substation.

Interaction of elements in carbon-graphite—Al—Mg—Zn—Cu melt system under combined action of temperature and pressure

2021 ◽  
pp. 514-518
Author(s):  
N.Yu. Miroshkin ◽  
V.A. Gulevsky ◽  
S.N. Tsurikhin ◽  
A.I. Bogdanov ◽  
L.M. Gurevich ◽  
...  
Keyword(s):  
Matrix Alloy ◽  
Impregnation Process ◽  
Active Elements ◽  
Cu Alloy ◽  
Inner Surface ◽  
The Matrix ◽  
Combined Action ◽  
Temperature And Pressure ◽  
The Difference ◽  
Titanium Compounds

Redistribution of chemically active elements is established on its inner surface of pores and at the interface with the alloy when impregnating carbon-graphite framework with Al—Mg—Zn—Cu alloy at temperature of 800 °C under pressure of up to 3 MPa. In this case, change in the solubility of melt elements in aluminum is possible as result of the combined action of temperature and pressure in the impregnation process, created due to the difference in the coefficients of thermal and thermal expansion of the matrix alloy, and the material of the impregnation device during impregnation. Titanium compounds are found in the pores filled with metal that are not added to the matrix alloy, but are formed as result of the contact of the matrix alloy melt with the walls of the impregnation device.

A new crustal shear-velocity model in Southwest China from joint seismological inversion and its implications for regional crustal dynamics

2019 ◽  
Author(s):  
Yan Yang ◽  
Huajian Yao ◽  
Hanxiao Wu ◽  
Ping Zhang ◽  
Maomao Wang
Keyword(s):  
Tibetan Plateau ◽  
Rayleigh Wave ◽  
Joint Inversion ◽  
Shear Velocity ◽  
P Wave ◽  
Tectonic Deformation ◽  
Inversion Method ◽  
Low Velocity ◽  
Sw China ◽  
S Period

SUMMARY Southwest (SW) China is located in a transition site from the active Tibetan Plateau to the stable Yangtze craton, which has complicated tectonic deformation and severe seismic hazards. We combine data from ambient noise, teleseismic body and surface waves, and petroleum wells to better constrain the crustal shear-velocity structure in SW China. We jointly invert the Rayleigh wave dispersion (5–40 s period), Rayleigh wave ZH ratio (20–60 s period), and P-wave receiver function for 114 permanent stations with a stepwise linearized joint inversion method. Compared to previous tomography results, we observe higher shear velocity in the sedimentary rocks within the Sichuan Basin, which is consistent with sonic logging measurements. Our model reveals widespread low-velocity zones in the mid-lower crust, and their boundaries correlate well with major fault systems. Between two main mid-crustal low-velocity channels, a prominent high-velocity region surrounded by earthquakes is observed in the inner zone of the Emeishan large igneous province (ELIP) and around the Anninghe-Zemuhe fault zone. These observations are comparable to regional tomography results using very dense arrays. Based on the results, we suggest that mid-lower crustal ductile flow and upper-crustal rigid fault movement play equally important roles in controlling the regional deformation styles and earthquake distribution in SW China. Our results also resolve thick crust-mantle transition zones beneath the eastern Tibetan Plateau and the inner zone of the ELIP due to ‘top-down’ and ‘bottom-up’ crust-mantle interactions, respectively. Our new model can serve as a reference crustal model of future high resolution model construction in SW China.

scholarly journals Study on the Origin and Fluid Identification of Low-Resistance Gas Reservoirs

Geofluids ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-12
Author(s):  
Zhou Yuhui ◽  
Hu Qingxiong ◽  
Liu Wentao ◽  
Wu Zhiqi ◽  
Yan Yule ◽  
...  
Keyword(s):  
Technical Problem ◽  
Laboratory Data ◽  
Gas Content ◽  
Formation Mechanisms ◽  
Gas Reservoirs ◽  
Gas Reservoir ◽  
Fluid Identification ◽  
Low Resistance ◽  
Neutron Porosity ◽  
The Difference

The Wu 2 section of the Ke017 well block is a low-resistance gas reservoir with ultralow porosity and low permeability. The comprehensive analysis of rock lithology, physical properties, sedimentary characteristics, and gas content demonstrated that the development of micropores in illite/smectite dominated clay minerals together with the resulted additional conductivity capability and complex reservoir pore structures, as well as the enrichment of self-generating conductivity minerals like zeolites and pyrite which were the formation mechanisms of low-resistance gas layers in the Wu 2 section. A low-resistance gas reservoir has poor physical property, and it is difficult to distinguish the oil layer from the dry, gas, or water layers. In this paper, based on well/mud logging data and laboratory data, by taking advantages of the “excavation effect” of neutron gas and the dual-lateral resistivity difference between different depths, we successfully established a set of low-contrast log response methods for the identification and evaluation of oil layer and formation fluids. For a gas layer, the difference between neutron porosity and acoustic (or density) porosity is smaller than 0 and the difference in dual-lateral resistivity is greater than 0. For a water layer, the neutron porosity is similar to the acoustic (or density) porosity and the dual-lateral resistivity difference will be less than 0. While for a dry layer or a layer with both gas and water, the difference in porosity as well as dual-lateral resistivity is very small. The proposed method effectively solves the technical problem of oil layer and formation fluid identification in low-resistance gas reservoirs.

scholarly journals Poroelastic property analysis of seismic low-frequency shadows associated with gas reservoirs

2020 ◽  
Vol 17 (3) ◽  
pp. 463-474
Author(s):  
Shengjie Li ◽  
Ying Rao
Keyword(s):  
Mode Conversion ◽  
Low Frequency ◽  
Simulation Method ◽  
Wave Mode ◽  
Pore Fluid ◽  
Fluid Distribution ◽  
Gas Reservoirs ◽  
Tight Sandstone ◽  
Reservoir Properties ◽  
Gas Bearing

Abstract Seismic low-frequency amplitude shadows have been widely used as a hydrocarbon indicator. This study investigates the effect of reservoir properties and seismic wave mode conversion on the characteristics of the low-frequency amplitude shadows in gas-bearing reservoirs. The target gas reservoirs are typically related to the lithology of tight sandstone with strong heterogeneity. Pore-fluid distribution within the reservoirs presents patchy saturation in the vertical and horizontal directions, and this patchy saturation easily induces low-frequency shadows beneath gas-bearing reservoirs. These low-frequency shadows are validated by using a poroelastic simulation method. The results of our field case-based study indicate that pore-fluid property, plus the thickness and heterogeneity of reservoirs are the key elements in the generation of low-frequency shadows. The results also indicate that the poroelastic simulation method can be used to effectively predict the spatial distribution of gas-bearing reservoirs, by directly verifying the low-frequency shadow phenomenon existing in the seismic data.

Sign in / Sign up

Export Citation Format

Share Document