Modeling of multiphase flow in low permeability porous media:Effect of wettability and pore structure properties

Multiphase flow in low permeability porous media is involved in numerous energy and environmental applications.However,a complete description of this process is challenging due to the limited modeling scale and the effects of complex pore structures and wettability.To address this issue,based on the digital rock of low permeability sandstone,a direct numerical simulation is performed considering the interphase drag and boundary slip to clarify the microscopic water-oil displacement process.In addition,a dual-porosity pore network model(PNM)is constructed to obtain the water-oil relative permeability of the sample.The displacement efficiency as a recovery process is assessed under different wetting and pore structure properties.Results show that microscopic displacement mechanisms explain the corresponding macroscopic relative permeability.The injected water breaks through the outlet earlier with a large mass flow,while thick oil films exist in rough hydrophobic surfaces and poorly connected pores.The variation of water-oil relative permeability is significant,and residual oil saturation is high in the oil-wet system.The flooding is extensive,and the residual oil is trapped in complex pore networks for hydrophilic pore surfaces;thus,water relative permeability is lower in the water-wet system.While the displacement efficiency is the worst in mixed-wetting systems for poor water connectivity.Microporosity negatively correlates with invading oil volume fraction due to strong capillary resistance,and a large microporosity corresponds to low residual oil saturation.This work provides insights into the water-oil flow from different modeling perspectives and helps to optimize the development plan for enhanced recovery.

Pore-pressure and stress-coupled creep behavior in deep coal:Insights from real-time NMR analysis

Understanding the variations in microscopic pore-fracture structures(MPFS) during coal creep under pore pressure and stress coupling is crucial for coal mining and effective gas treatment. In this manuscript, a triaxial creep test on deep coal at various pore pressures using a test system that combines in-situ mechanical loading with real-time nuclear magnetic resonance(NMR) detection was conducted.Full-scale quantitative characterization, online real-time detection, and visualization of MPFS during coal creep influenced by pore pressure and stress coupling were performed using NMR and NMR imaging(NMRI) techniques. The results revealed that seepage pores and microfractures(SPM) undergo the most significant changes during coal creep, with creep failure gradually expanding from dense primary pore fractures. Pore pressure presence promotes MPFS development primarily by inhibiting SPM compression and encouraging adsorption pores(AP) to evolve into SPM. Coal enters the accelerated creep stage earlier at lower stress levels, resulting in more pronounced creep deformation. The connection between the micro and macro values was established, demonstrating that increased porosity at different pore pressures leads to a negative exponential decay of the viscosity coefficient. The Newton dashpot in the ideal viscoplastic body and the Burgers model was improved using NMR experimental results, and a creep model that considers pore pressure and stress coupling using variable-order fractional operators was developed. The model’s reasonableness was confirmed using creep experimental data. The damagestate adjustment factors ω and β were identified through a parameter sensitivity analysis to characterize the effect of pore pressure and stress coupling on the creep damage characteristics(size and degree of difficulty) of coal.

A hybrid machine learning optimization algorithm for multivariable pore pressure prediction

Pore pressure is essential data in drilling design,and its accurate prediction is necessary to ensure drilling safety and improve drilling efficiency.Traditional methods for predicting pore pressure are limited when forming particular structures and lithology.In this paper,a machine learning algorithm and effective stress theorem are used to establish the transformation model between rock physical parameters and pore pressure.This study collects data from three wells.Well 1 had 881 data sets for model training,and Wells 2 and 3 had 538 and 464 data sets for model testing.In this paper,support vector machine(SVM),random forest(RF),extreme gradient boosting(XGB),and multilayer perceptron(MLP)are selected as the machine learning algorithms for pore pressure modeling.In addition,this paper uses the grey wolf optimization(GWO)algorithm,particle swarm optimization(PSO)algorithm,sparrow search algorithm(SSA),and bat algorithm(BA)to establish a hybrid machine learning optimization algorithm,and proposes an improved grey wolf optimization(IGWO)algorithm.The IGWO-MLP model obtained the minimum root mean square error(RMSE)by using the 5-fold cross-validation method for the training data.For the pore pressure data in Well 2 and Well 3,the coefficients of determination(R^(2))of SVM,RF,XGB,and MLP are 0.9930 and 0.9446,0.9943 and 0.9472,0.9945 and 0.9488,0.9949 and 0.9574.MLP achieves optimal performance on both training and test data,and the MLP model shows a high degree of generalization.It indicates that the IGWO-MLP is an excellent predictor of pore pressure and can be used to predict pore pressure.

Experimental investigation on coal pore-fracture variation and fractal characteristics synergistically affected by solvents for improving clean gas extraction

Chemical solvents instead of pure water being as hydraulic fracturing fluid could effectively increase permeability and improve clean methane extraction efficiency.However,pore-fracture variation features of lean coal synergistically affected by solvents have not been fully understood.Ultrasonic testing,nuclear magnetic resonance analysis,liquid phase mass spectrometry was adopted to comprehensively analyze pore-fracture change characteristics of lean coal treated by combined solvent(NMP and CS_(2)).Meanwhile,quantitative characterization of above changing properties was conducted using geometric fractal theory.Relationship model between permeability,fractal dimension and porosity were established.Results indicate that the end face fractures of coal are well developed after CS2and combined solvent treatments,of which,end face box-counting fractal dimensions range from 1.1227 to 1.4767.Maximum decreases in ultrasonic longitudinal wave velocity of coal affected by NMP,CS_(2)and combined solvent are 2.700%,20.521%,22.454%,respectively.Solvent treatments could lead to increasing amount of both mesopores and macropores.Decrease ratio of fractal dimension Dsis 0.259%–2.159%,while permeability increases ratio of NMR ranges from 0.1904 to 6.4486.Meanwhile,combined solvent could dissolve coal polar and non-polar small molecules and expand flow space.Results could provide reference for solvent selection and parameter optimization of permeability-enhancement technology.

The Negative Thermal Expansion Property of NdMnO_(3) Based on Pores Effect and Phase Transition

A novel negative thermal expansion(NTE) material NdMnO_(3) was synthesized by solid-state method at 1 523 K. The crystal structure, phase transition, pores effect and negative expansion properties of NdMnO_(3) were investigated by variable temperature X-ray diffraction(XRD), scanning electron microscope(SEM) and variable temperature Raman spectra. The compound exhibits NTE properties in the orderly O' phase crystal structure. When the temperature is from 293 to 759 K, the ceramic NdMnO_(3) shows negative thermal expansion of-4.7×10^(-6)/K. As temperature increases, the ceramic NdMnO_(3) presents NTE property range from 759 to 1 007 K. The average linear expansion coefficient is-18.88×10^(-6)/K. The physical mechanism of NTE is discussed and clarified through experiments.

A novel box-counting method for quantitative fractal analysis of threedimensional pore characteristics in sandstone

Fractal theory offers a powerful tool for the precise description and quantification of the complex pore structures in reservoir rocks,crucial for understanding the storage and migration characteristics of media within these rocks.Faced with the challenge of calculating the three-dimensional fractal dimensions of rock porosity,this study proposes an innovative computational process that directly calculates the three-dimensional fractal dimensions from a geometric perspective.By employing a composite denoising approach that integrates Fourier transform(FT)and wavelet transform(WT),coupled with multimodal pore extraction techniques such as threshold segmentation,top-hat transformation,and membrane enhancement,we successfully crafted accurate digital rock models.The improved box-counting method was then applied to analyze the voxel data of these digital rocks,accurately calculating the fractal dimensions of the rock pore distribution.Further numerical simulations of permeability experiments were conducted to explore the physical correlations between the rock pore fractal dimensions,porosity,and absolute permeability.The results reveal that rocks with higher fractal dimensions exhibit more complex pore connectivity pathways and a wider,more uneven pore distribution,suggesting that the ideal rock samples should possess lower fractal dimensions and higher effective porosity rates to achieve optimal fluid transmission properties.The methodology and conclusions of this study provide new tools and insights for the quantitative analysis of complex pores in rocks and contribute to the exploration of the fractal transport properties of media within rocks.

Effects of confining pressure and pore pressure on multipole borehole acoustic field in fluid-saturated porous media

In-situ stress is a common stress in the exploration and development of oil reservoirs. Therefore, it is of great significance to study the propagation characteristics of borehole acoustic waves in fluid-saturated porous media under stress.Based on the acoustoelastic theory of fluid-saturated porous media, the field equation of fluid-saturated porous media under the conditions of confining pressure and pore pressure and the acoustic field formula of multipole source excitation in open hole are given. The influences of pore pressure and confining pressure on guided waves of multipole borehole acoustic field in fluid-saturated porous media are investigated. The numerical results show that the phase velocity and excitation intensity of guided wave increase significantly under the confining pressure. For a given confining pressure, the phase velocity of the guided wave decreases with pore pressure increasing. The excitation intensity of guided wave increases at low frequency and then decreases at high frequency with pore pressure increasing, except for that of Stoneley wave which decreases in the whole frequency range. These results will help us get an insight into the influences of confining pressure and pore pressure on the acoustic field of multipole source in borehole around fluid-saturated porous media.

Pore Connectivity of Deep Lacustrine Shale and its Effect on Gas-bearing Characteristics in the Songliao Basin:Implications from Continental Scientific Drilling

The lacustrine shale of deep Shahezi Formation in the Songliao basin has great gas potential,but its pore evolution,heterogeneity,and connectivity characteristics remain unclear.In this work,total organic carbon analysis,rock pyrolysis,X-ray diffraction field emission scanning electron microscopy,the particle and crack analysis system software,low-temperature nitrogen adsorption experiment,fractal theory,high-pressure mercury injection experiment and nuclear magnetic resonance experiment were used to study the Shahezi shale from Well SK-2.The result indicated that the organic pores in Shahezi shale are not developed,and the intergranular and intragranular pores are mainly formed by illitedominated clay.As the burial depth increases,the pore size and slit-shaped pores formed by clay decrease,and dissolved pores in the feldspar and carbonate minerals and dissolved fractures in the quartz increase.The pore evolution is affected by clay,compaction,and high-temperature corrosion.Based on the pore structure characteristics reflected by the pore size distribution and pore structure parameters obtained by multiple experimental methods,the pore development and evolution are divided into three stages.During stageⅠandⅡ,the pore heterogeneity of the shale reservoirs increases with the depth,the physical properties and pore connectivity deteriorate,but the gas-bearing property is good.In stageⅢ,the pore heterogeneity is the highest,its gas generation and storage capacity are low,but the increase of micro-fractures makes pore connectivity and gas-bearing better.

Capacitive energy storage from single pore to porous electrode identified by frequency response analysis

Rate capability,peak power,and energy density are of vital importance for the capacitive energy storage(CES)of electrochemical energy devices.The frequency response analysis(FRA)is regarded as an efficient tool in studying the CES.In the present work,a bi-scale impedance transmission line model(TLM)is firstly developed for a single pore to a porous electrode.Not only the TLM of the single pore is reparameterized but also the particle packing compactness is defined in the bi-scale.Subsequently,the CES properties are identified by FRA,focused on rate capability vs.characteristic frequency,peak power vs.equivalent series resistance,and energy density vs.low frequency limiting capacitance for a single pore to a porous electrode.Based on these relationships,the CES properties are numerically simulated and theoretically predicted for a single pore to a porous electrode in terms of intra-particle pore length,intra-particle pore diameter,inter-particle pore diameter,electrolyte conductivity,interfacial capacitance&exponent factor,electrode thickness,electrode apparent surface area,and particle packing compactness.Finally,the experimental diagnosis of four supercapacitors(SCs)with different electrode thicknesses is conducted for validating the bi-scale TLM and gaining an insight into the CES properties for a porous electrode to a single pore.The calculating results suggest,to some extent,the inter-particle pore plays a more critical role than the intra-particle pore in the CES properties such as the rate capability and the peak power density for a single pore to a porous electrode.Hence,in order to design a better porous electrode,more attention should be given to the inter-particle pore.

Preparation of Chiral Silica Nanostructures with Radial Pores through Single-templating Approach

A chiral low-molecular-weight gelator(LMWG) L-16Ala5PyPF6 was synthesized from L-alanine, which can cause physical gel in n-propanol, ethyl acetate, butylene oxide, water, benzene, 1,4-dioxane and chloroform. The sol-gel reactions were carried out in a mixture of stronger ammonia water and n-propanol at the volume ratio of 2:8. Single-handed twisted silica nanostructures with pore channels vertical to the wall surfaces were first prepared through a single-templating approach comparing with the reported double template method. The formation mechanism of radial pore structure was studied by transmission electron microscopy at different reaction time intervals, which indicated that the radial pore structure was formed via a structural transition in the sol-gel transcription process.

Influence of Acid Treatment on Pore Structure and Fractal Characterization of a Tight Sandstone:A Case Study from Wudun Sag,Dunhuang Basin

In this study,X-ray diffraction,N_(2)adsorption(N_(2)A),and mercury intrusion(MI)experiments were used to investigate the influence of acid treatment on pore structure and fractal characterization of tight sandstones.The results showed that acid treatment generated a certain number of ink-bottle pores in fine sandstone,aggravated the ink-bottle effect in the sandy mudstone,and transformed some smaller pores into larger ones.After the acid treatment,both the pore volume in the range of 2–11 nm and 0.271–8μm for the fine sandstone and the entire pore size range for the sandy mudstone significantly increased.The dissolution of sandstone cement causes the fine sandstone particles to fall off and fill the pores;the porosity increased at first but then decreased with acid treatment time.The fractal dimension obtained using the Frenkel-Halsey-Hill model was positively correlated with acid treatment time.However,the total fractal dimensions obtained by MI tests showed different changes with acid treatment time in fine sandstone and sandy mudstone.These results provide good guiding significance for reservoir acidification stimulation.

Rock physics model for velocity–pressure relations and its application to shale pore pressure estimation

Acoustic wave velocity has been commonly utilized to predict subsurface geopressure using empirical relations.Acoustic wave velocity is, however, affected by many factors. To estimate pore pressure accurately, we here propose to use elastic rock physics models to understand and analyze quantitatively the various contributions from these different factors affecting wave velocity. We report a closed-form relationship between the frame flexibility factor(γ) in a rock physics model and differential pressure, which presents the major control of pressure on elastic properties such as bulk modulus and compressional wave velocity. For a gas-bearing shale with abundant micro-cracks and fractures, its bulk modulus is much lower at abnormally high pore pressure(high γ values) where thin cracks and flat pores are open than that at normal hydrostatic pressure(low γ values) where pores are more rounded on average. The developed relations between bulk modulus and differential pressure have been successfully applied to the Upper Ordovician Wufeng and Lower Silurian Longmaxi formations in the Dingshan area of the Sichuan Basin to map the three-dimensional spatial distribution of pore pressure in the shale, integrating core, log and seismic data. The estimated results agree well with field measurements. Pressure coefficient is positively correlated to gas content. The relations and methods reported here could be useful for hydrocarbon exploration, production, and drilling safety in both unconventional and conventional fields.

Pore throat structure heterogeneity and its effect on gas-phase seepage capacity in tight sandstone reservoirs:A case study from the Triassic Yanchang Formation,Ordos Basin

The microscopic heterogeneity of pore-throat structures in tight sandstone is a crucial parameter for understanding the transport mechanism of fluid flow.In this work,we firstly developed the new procedure to characterize the pore size distribution(PSD)and throat size distribution(TSD)by combining the nuclear magnetic resonance(NMR),cast thin section(CTS),and constant-rate mercury injection(CRMI)tests,and used the permeability estimated model to verify the full-scale PSD and TSD.Then,we respectively analyzed the fractal feature of the pore and throat,and characterized the heterogeneity of pores and throats.Finally,we elaborated the effect of the pore and throat heterogeneity on the gas-phase seepage capacity base on the analysis of the simple capillary tube model and gas-flooding experiment.The results showed that(1)The PSD and TSD of the tight sandstone sample ranged from 0.01 to 10 mm and from 0.1 to 57 mm,respectively,mainly contributed by the micropores and mesopores.Meanwhile,the permeability estimated by the PSD and TSD was consistent with the experimental permeability,and relative error was lower than 8%.(2)The PSD and TSD exhibited multifractal characteristics,and singularity strength range,Δα,could be used as the indicator for characterizing the heterogeneity of pore and throat.Furthermore,the throat of the sample showed stronger heterogeneity than that the pore.(3)The throats played an important role for the fluid transport in the tight sandstone,and the effect of the throat heterogeneity on the gas-phase seepage capacity was different under the lower and higher injection pressure.The macropores and micropores maybe respectively become the preferential migration pathways at the lower and higher injection pressure.In the end,the identification plate was established in our paper,and could be described the relationship among the throat heterogeneity,injection pressure,permeability and flow path of the gas phase in the tight sandstone.

Effect of seepage-induced erosion on soil macropore structure

Internal erosion is one of the important factors causing geological disasters.The microstructure of soil can change with seepage erosion,resulting in changes in the hydraulic and mechanical properties of the soil.The evolution of seepage erosion is investigated with X-ray computed tomography(CT)in this study.The change in macropore structure characteristics during the seepage erosion test is quantified and the influence of seepage erosion on soil deformation is analyzed.Moreover,a pore network model(PNM)is established for the specimens and the evolution of the connected pore size characteristics is assessed.The results show that the macropore structure is significantly affected by seepage erosion,especially in terms of the porosity and pore geometry characteristics.The changes in macropore structure characteristics are most obvious in the lower part of the specimen.The influence of seepage erosion on the pore size distribution(PSD)and soil deformation is heterogeneous and closely dependent on the spatial location of the soil.Moreover,seepage erosion enhances macropore connectivity and has a directional impact on macropore orientation.These findings can provide a reference for the theoretical modeling and numerical simulation of the seepage erosion and improve the understanding of the seepage erosion evolution in engineering practice.

Pore structure characterization based on NMR experiment: A case from the Shanxi Formation tight sandstones in the Daning-Jixian area, eastern Ordos Basin

The study of pore structure requires consideration of important factors including pore throat size,pore radius composition,and pore-throat configuration.As the nuclear magnetic resonance(NMR)experimental results contain rich information about pore structures and fluid occurrence states,this study investigated the pore structures of the tight sandstone reservoirs of the Shanxi Formation in the Daning-Jixian area,eastern Ordos Basin.Firstly,by making the inverse cumulative curve of the NMR T2 spectrum coincide with the capillary pressure curves which were obtained by the mercury injection capillary pressure(MICP)technique,this study derived a conversion coefficient that can be used to convert the NMR T2 spectrum into the pore throat radius distribution curves based on the NMR experimental results.Subsequently,we determined the pore radius intervals corresponding to irreducible water distribution using the NMR-derived pore radius distribution curves.Finally,the NMR T2 distribution curves based on the fractal theory were analyzed and the relationships between fractal dimensions and parameters,including permeability,porosity,reservoir quality index(RQI),flow zone indicator(FZI),irreducible water saturation,RT35,and RT50,were also discussed.The NMR-derived pore throat radius distribution curves of the study area are mainly unimodal,with some curves showing slightly bimodal distributions.The irreducible water mainly occurs in small pores with a pore radius less than 100 nm.As the permeability decreases,the contribution rate of small pores to the irreducible water gradually increases.The NMR-based fractal dimensions of pores show a two-segment distribution.Small pores have small fractal dimensions and are evenly distributed,while large pores have large fractal dimensions and complex pore structures.The fractal dimension of large pores(Dmax)is poorly correlated with porosity but strongly correlated with FZI,RQI,RT35,and RT50.These results indicate that large pores are the main pore zones that determine the seepage capacity of the reservoirs.Additionally,there is a certain correlation between Dmax and the irreducible water saturation.

Multi-scale pore fractal characteristics of differently ranked coal and its impact on gas adsorption

Well-developed pores and cracks in coal reservoirs are the main venues for gas storage and migration.To investigate the multi-scale pore fractal characteristics,six coal samples of different rankings were studied using high-pressure mercury injection(HPMI),low-pressure nitrogen adsorption(LPGA-N2),and scanning electron microscopy(SEM)test methods.Based on the Frankel,Halsey and Hill(FHH)fractal theory,the Menger sponge model,Pores and Cracks Analysis System(PCAS),pore volume complexity(D_(v)),coal surface irregularity(Ds)and pore distribution heterogeneity(D_(p))were studied and evaluated,respectively.The effect of three fractal dimensions on the gas adsorption ability was also analyzed with high-pressure isothermal gas adsorption experiments.Results show that pore structures within these coal samples have obvious fractal characteristics.A noticeable segmentation effect appears in the Dv1and Dv2fitting process,with the boundary size ranging from 36.00 to 182.95 nm,which helps differentiate diffusion pores and seepage fractures.The D values show an asymmetric U-shaped trend as the coal metamorphism increases,demonstrating that coalification greatly affects the pore fractal dimensions.The three fractal dimensions can characterize the difference in coal microstructure and reflect their influence on gas adsorption ability.Langmuir volume(V_(L))has an evident and positive correlation with Dsvalues,whereas Langmuir pressure(P_(L))is mainly affected by the combined action of Dvand Dp.This study will provide valuable knowledge for the appraisal of coal seam gas reservoirs of differently ranked coals.

Probing deactivation by coking in catalyst pellets for dry reforming of methane using a pore network model

Dry reforming of methane(DRM) is an attractive technology for utilizing the greenhouse gases(CO_(2) and CH_(4)) to produce syngas. However, the catalyst pellets for DRM are heavily plagued by deactivation by coking, which prevents this technology from commercialization. In this work, a pore network model is developed to probe the catalyst deactivation by coking in a Ni/Al_(2)O_(3) catalyst pellet for DRM. The reaction conditions can significantly change the coking rate and then affect the catalyst deactivation. The catalyst lifetime is higher under lower temperature, pressure, and CH_(4)/CO_(2) molar ratio, but the maximum coke content in a catalyst pellet is independent of these reaction conditions. The catalyst pellet with larger pore diameter, narrower pore size distribution and higher pore connectivity is more robust against catalyst deactivation by coking, as the pores in this pellet are more difficult to be plugged or inaccessible.The maximum coke content is also higher for narrower pore size distribution and higher pore connectivity, as the number of inaccessible pores is lower. Besides, the catalyst pellet radius only slightly affects the coke content, although the diffusion limitation increases with the pellet radius. These results should serve to guide the rational design of robust DRM catalyst pellets against deactivation by coking.

Fibrous MXene Aerogels with Tunable Pore of Contaminated Seawater

The seawater desalination based on solardriven interfacial evaporation has emerged as a promising technique to alleviate the global crisis on freshwater shortage.However,achieving high desalination performance on actual,oil-contaminated seawater remains a critical challenge,because the transport channels and evaporation interfaces of the current solar evaporators are easily blocked by the oil slicks,resulting in undermined evaporation rate and conversion efficiency.Herein,we propose a facile strategy for fabricating a modularized solar evaporator based on flexible MXene aerogels with arbitrarily tunable,highly ordered cellular/lamellar pore structures for high-efficiency oil interception and desalination.The core design is the creation of 1D fibrous MXenes with sufficiently large aspect ratios,whose superior flexibility and plentiful link forms lay the basis for controllable 3D assembly into more complicated pore structures.The cellular pore structure is responsible for effective contaminants rejection due to the multi-sieving effect achieved by the omnipresent,isotropic wall apertures together with underwater superhydrophobicity,while the lamellar pore structure is favorable for rapid evaporation due to the presence of continuous,large-area evaporation channels.The modularized solar evaporator delivers the best evaporation rate(1.48 kg m-2h-1)and conversion efficiency(92.08%)among all MXene-based desalination materials on oil-contaminated seawater.

Pore-scale modeling of pore structure properties and wettability effect on permeability of low-rank coal

Permeability is a key parameter for coalbed methane development.Although the absolute permeability of coal has been extensively studied,wettability and pore structure properties continue to challenge the microscopic description of water-gas flow in coal.For this purpose,we reconstructed the microstructures of low-rank coal using micro-computed tomography(micro-CT)images.Pore geometry and pore-throat parameters are introduced to establish a relationship with absolute permeability.A dual-porosity pore network model is developed to study water-gas displacement under different wetting and pore structure properties.Results show that absolute permeability is significantly affected by pore geometry and can be described using a binary quadratic function of porosity and fractal dimension.Water-gas relative permeability varies significantly and the residual gas saturation is lower;the crossover saturation first decreased and then increased with increasing porosity under hydrophobic conditions.While the water relative permeability is lower and a certain amount of gas is trapped in complex pore-throat networks;the crossover saturation is higher under hydrophilic conditions.Models with large percolating porosity and well-developed pore networks have high displacement efficiency due to low capillary resistance and avoidance of trapping.This work provides a systematic description of absolute permeability and water-gas relative permeability in coal microstructure for enhanced gas recovery.

The mechanism for tuning the corrosion resistance and pore density of plasma electrolytic oxidation(PEO)coatings on Mg alloy with fluoride addition

Here we prepared PEO coatings on Mg alloys in silicate-NaOH-phosphate electrolyte containing different concentrations of NaF addition.The detailed microstructural characterizations combining with potentiodynamic polarization and electrochemical impedance spectra(EIS)were employed to investigate the roles of fluoride in the growth and corrosion properties of PEO coating on Mg.The result shows the introduction of NaF led to a fluoride-containing nanolayer(FNL)formed at the Mg/coating interface.The FNL consists of MgO nanoparticles and insoluble MgF_(2)nanoparticles(containing rutile phase and cubic phase).The increase in the NaF concentration of the electrolyte increases the thickness and the MgF_(2)content in the FNL.When anodized in the electrolyte containing 2 g/L NaF,the formed FNL has the highest thickness of 100-200 nm along with the highest value of x of∼0.6 in(MgO)_(1-x)(MgF_(2))x resulted in the highest corrosion performance of PEO coating.In addition,when anodized in the electrolyte containing a low NaF concentration(0.4-0.8 g/L),the formed FNL was thin and discontinuous,which would decrease the pore density and increase the coating's uniformness simultaneously.