Data assimilation(DA)and uncertainty quantification(UQ)are extensively used in analysing and reducing error propagation in high-dimensional spatial-temporal dynamics.Typical applications span from computational fluid ...Data assimilation(DA)and uncertainty quantification(UQ)are extensively used in analysing and reducing error propagation in high-dimensional spatial-temporal dynamics.Typical applications span from computational fluid dynamics(CFD)to geoscience and climate systems.Recently,much effort has been given in combining DA,UQ and machine learning(ML)techniques.These research efforts seek to address some critical challenges in high-dimensional dynamical systems,including but not limited to dynamical system identification,reduced order surrogate modelling,error covariance specification and model error correction.A large number of developed techniques and methodologies exhibit a broad applicability across numerous domains,resulting in the necessity for a comprehensive guide.This paper provides the first overview of state-of-the-art researches in this interdisciplinary field,covering a wide range of applications.This review is aimed at ML scientists who attempt to apply DA and UQ techniques to improve the accuracy and the interpretability of their models,but also at DA and UQ experts who intend to integrate cutting-edge ML approaches to their systems.Therefore,this article has a special focus on how ML methods can overcome the existing limits of DA and UQ,and vice versa.Some exciting perspectives of this rapidly developing research field are also discussed.Index Terms-Data assimilation(DA),deep learning,machine learning(ML),reduced-order-modelling,uncertainty quantification(UQ).展开更多
Failure of irregular rock samples may provide implications in the rapid estimation of rock strength,which is imperative in rock engineering practice.In this work,analytical,experimental and numerical investigations we...Failure of irregular rock samples may provide implications in the rapid estimation of rock strength,which is imperative in rock engineering practice.In this work,analytical,experimental and numerical investigations were carried out to study the mechanical properties and failure characteristics of rock spheres under paired point loads.Analytical solutions indicted that with the increase in sample size(contact angle)and decrease in Poisson’s ratio,the uneven tensile stress in theta direction decreased.Then laboratory experiments were carried out to investigate the load characteristics and failure mode of spherical marble samples with different sizes subjected to a pair of diametral point loads.The discrete element method(DEM)was adopted to study the failure process of rock spheres.The effect of the sphere diameter on the point load contact angle was examined in terms of peak load,crushed zone distribution and energy dissipation.Experimental and numerical results showed that the samples primarily fail in tension,with crushed zones formed at both loading points.With increase in the sample size,the contact angle,crushed area and total work increase.As the specimen diameter increases from 30 mm to 50 mm,the peak load on the specimen increases from 3.6 kN to 8.8 kN,and the percentage of crushed zone(ratio of crushing zone to sample radius,d/r)increased from 0.191 to 0.262.The results of the study have implications for understanding the failure of irregular rock specimens under point loading conditions and their size effects.展开更多
In order to explore the pressure relief and structure stability mechanism of lateral cantilever structure in the stope under the direct coverage of thick hard roof and its impact on the gob-side entry retaining, a lat...In order to explore the pressure relief and structure stability mechanism of lateral cantilever structure in the stope under the direct coverage of thick hard roof and its impact on the gob-side entry retaining, a lateral cantilever fractured structural mechanical model was established on the basis of clarification for the stress environment of gob-side entry retaining, and the equation of roof given deformation and the balance judgment for fracture block were obtained. The optimal cantilever length was proposed based on the comparison of roof structural characteristics and the stress, deformation law of surrounding rocks under six different cantilever lengths by numerical simulation method. Double stress peaks exist on the sides of gob-side entry retaining and the entry located in the low stress area. The pressure of gob-side entry retaining can be relieved by reducing the cantilever length. However, due to the impact of arch structure of rock beam, unduly short cantilever would result in insufficient pressure relief and unduly long cantilever would bring larger roof stress which results in intense deformation. Therefore, there is optimal cantilever length, which was 7-8 m in this project that enables to achieve the minimum deformation and the most stabilized rock structure of entry retaining. An engineering case of gob-side entry retaining with the direct coverage of 10 m thick hard limestone roof was put forward, and the measured data verified the reasonability of conclusion.展开更多
The flow properties of geomechanically generated discrete fracture networks are examined in the context of channelling.Fracture networks are generated by growing fractures in tension,modelling the low permeability roc...The flow properties of geomechanically generated discrete fracture networks are examined in the context of channelling.Fracture networks are generated by growing fractures in tension,modelling the low permeability rock as a linear elastic material.Fractures are modelled as discrete surfaces which grow quasi-statically within a three-dimensional(3D)volume.Fractures may have their locations specified as a simulation input,or be generated as a function of damage,quantified using the local variation in equivalent strain.The properties of the grown networks are shown to be a product of in situ stress,relative orientation of initial flaws,and competitive process of fracture interaction and growth.Fractures grow preferentially in the direction perpendicular to the direction of maximum tension and may deviate from this path due to mechanical fracture interaction.Flow is significantly channelled through a subset of the fractures in the full domain,consistent with observations of other real and simulated fractures.As the fracture networks grow,small changes in the geometry of the fractures lead to large changes in the locations and scale of primary flow channels.The flow variability and formation of channels are examined for two growing networks,one with a fixed amount of fractures,and another with nucleating fractures.The interaction between fractures is shown to modify the local stress field,and in turn the aperture of the fractures.Pathways for single-phase flow are the results of hydro-mechanical effects in fracture networks during growth.These are the results of changes to the topology of the network as well as the result of mechanical self-organisation which occurs during interaction leading to growth and intersection.展开更多
The stable isotope systems of the transition metals potentially provide constraints on the current and past operation of the biological pump,and on the state of ocean redox in Earth history.Here we focus on two exempl...The stable isotope systems of the transition metals potentially provide constraints on the current and past operation of the biological pump,and on the state of ocean redox in Earth history.Here we focus on two exemplar metals,nickel(Ni)and zinc(Zn).The oceanic dissolved pool of both elements is isotopically heavier than the known inputs,implying an output with light isotope compositions.The modern oceanic cycle of both these elements is dominated by biological uptake into photosynthesised organic matter and output to sediment.It is increasingly clear,however,that such uptake is associated with only very minor isotope fractionation.We suggest that the isotopic balance is instead closed by the sequestration of light isotopes to sulphide in anoxic and organic-rich sediments,so that it is ocean chemistry that controls these isotope systems,and suggesting a different but equally interesting array of questions in Earth history that can be addressed with these systems.展开更多
Deep coal seams are one of the world’s most widespread deposits for carbon dioxide(C02)disposal and are generally located near large point sources of CO_(2)emissions.The injection of CO_(2)into coal seams has great p...Deep coal seams are one of the world’s most widespread deposits for carbon dioxide(C02)disposal and are generally located near large point sources of CO_(2)emissions.The injection of CO_(2)into coal seams has great potential to sequester CO_(2)while simultaneously enhancing coalbed methane(CO_(2)-ECBM)recovery.Pilot tests of CO_(2)-ECBM have been conducted in coal seams worldwide with favorable early results.However,one of the main technical barriers in coal seams needs to be resolved:Injecting CO_(2)reduces coal permeability and well injectivity.Here,using in situ synchrotron X-ray microtomography,we provide the first observational evidence that injecting nitrogen(N_(2))can reverse much of this lost permeability by reopening fractures that have closed due to coal swelling induced by CO_(2)adsorption.Our findings support the notion that injecting minimally treated flue gas-a mixture of mainly N_(2) and CO_(2)-is an attractive alternative for ECBM recovery instead of pure CO_(2)injection in deep coal seams.Firstly,flue gas produced by power plants could be directly injected after particulate removal,thus avoiding high CO_(2)-separation costs.Secondly,the presence of N_(2)makes it possible to maintain a sufficiently high level of coal permeability.These results suggest that flue-gas ECBM for deep coal seams may provide a promising path toward net-zero emissions from coal mines.展开更多
Mixed ionic electronic conductors(MIECs)have attracted increasing attention as anode materials for solid oxide fuel cells(SOFCs)and they hold great promise for lowering the operation temperature of SOFCs.However,there...Mixed ionic electronic conductors(MIECs)have attracted increasing attention as anode materials for solid oxide fuel cells(SOFCs)and they hold great promise for lowering the operation temperature of SOFCs.However,there has been a lack of understanding of the performance-limiting factors and guidelines for rational design of composite metal-MIEC electrodes.Using a newly-developed approach based on 3 D-tomography and electrochemical impedance spectroscopy,here for the first time we quantify the contribution of the dual-phase boundary(DPB)relative to the three-phase boundary(TPB)reaction pathway on real MIEC electrodes.A new design strategy is developed for Ni/gadolinium doped ceria(CGO)electrodes(a typical MIEC electrode)based on the quantitative analyses and a novel Ni/CGO fiber-matrix structure is proposed and fabricated by combining electrospinning and tape-casting methods using commercial powders.With only 11.5 vol%nickel,the designer Ni/CGO fiber-matrix electrode shows 32%and 67%lower polarization resistance than a nano-Ni impregnated CGO scaffold electrode and conventional cermet electrode respectively.The results in this paper demonstrate quantitatively using real electrode structures that enhancing DPB and hydrogen kinetics are more efficient strategies to enhance electrode performance than simply increasing TPB.展开更多
We measured Electron Spin Resonance(ESR) signal intensity and Crystallinity Index(CI) of fine-(<16 m) and coarse-grained(>16 m) quartz in surface samples collected from the Taklimakan desert in western China, th...We measured Electron Spin Resonance(ESR) signal intensity and Crystallinity Index(CI) of fine-(<16 m) and coarse-grained(>16 m) quartz in surface samples collected from the Taklimakan desert in western China, the Badain Juran, Tengger and Mu Us deserts in northern China, and the Gobi desert in southern Mongolia. The aim was to test if these two parameters can be used as the accurate provenance tracers of Asian dust. Coarse-grained quartz displays spatial variability of the ESR signal intensity within each desert and of the CI within the Taklimakan and the Mongolian Gobi deserts. Coarse-grained quartz derived from the Mongolian Gobi and northern China deserts can be differentiated from the Taklimakan desert using the ESR signal intensity. By contrast, variability of ESR signal intensity and CI of fine-grained quartz particles is distinct within the Taklimakan and the Mongolian Gobi deserts, but less significant in three deserts of northern China. Spatially, fine-grained dust originating from three major Asian dust sources(i.e., the Gobi-sandy deserts in western China, northern China and southern Mongolia) can be differentiated using the combination of the ESR and CI signals. Our results suggest that combination of the ESR and CI results can discriminate the sources of fine-grained dust better than coarse-grained dust, providing an effective approach to trace the provenance of fine-grained dust depositions on the land and in the ocean.展开更多
In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation ...In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation in inversion. A feasible way to avoid the excessive storage demand is to reconstruct the source wavefield backward in time by storing the entire history of the wavefield in perfectly matched layers. In this paper, we make full use of the elementwise global property of the Laplace operator of the spectral element method (SEM) and propose an efficient source wavefield reconstruction method at the cost of storing the wavefield history only at single boundary layer nodes. Numerical experiments indicate that the accuracy of the proposed method is identical to that of the conventional method and is independent of the order of the Lagrange polynomials, the element type, and the temporal discretization method. In contrast, the memory-saving ratios of the conventional method versus our method is at least N when using either quadrilateral or hexahedron elements, respectively, where N is the order of the Lagrange polynomials used in the SEM. A higher memorysaving ratio is achieved with triangular elements versus quadrilaterals. The new method is applied to reverse time migration by considering the Marmousi model as a benchmark. Numerical results demonstrate that the method is able to provide the same result as the conventional method but with about 1/25 times lower storage demand. With the proposed wavefield reconstruction method, the storage demand is dramatically reduced;therefore, in-core memory storage is feasible even for large-scale three-dimensional adjoint inversion problems.展开更多
1.Historical context Flow in porous media is encountered in a wide variety of natural and engineered settings,including gas exchange in lungs and blood vessels,water extraction from underground aquifers,hydrocarbon pr...1.Historical context Flow in porous media is encountered in a wide variety of natural and engineered settings,including gas exchange in lungs and blood vessels,water extraction from underground aquifers,hydrocarbon production,and microfluidics in drug delivery and food manufacturing.Underground rock is a porous medium that holds water and hydrocarbons,and is a potential storage site for carbon dioxide(CO_(2))and hydrogen.Many manufactured devices also contain porous media to allow the transport of gases and water,including fuel cells,electrolyzers,and electrocatalysts for CO_(2) reduction.Their performance is limited by how rapidly two fluid phases can flow through a porous layer.展开更多
The creep phenomenon of rocks is quite complex and the creep mechanisms are far from being well understood.Although laboratory creep tests have been carried out to determine the creep deformation of various rocks,thes...The creep phenomenon of rocks is quite complex and the creep mechanisms are far from being well understood.Although laboratory creep tests have been carried out to determine the creep deformation of various rocks,these tests are expensive and time-consuming.Nanoindentation creep tests,as an alternative method,can be performed to investigate the mechanical and viscoelastic properties of granite samples.In this study,the reduced Young’s modulus,hardness,fracture toughness,creep strain rate,stress exponent,activation volume and maximum creep displacement of common rock-forming minerals of granite were calculated from nanoindentation results.It was found that the hardness decreases with the increase of holding time and the initial decrease in hardness was swift,and then it decreased slowly.The stress exponent values obtained were in the range from 4.5 to 22.9,which indicates that dislocation climb is the creep deformation mechanism.In addition,fracture toughness of granite’s rock-forming minerals was calculated using energy-based method and homogenization method was adopted to upscale the micro-scale mechanical properties to macro-scale mechanical properties.Last but not least,both three-element Voigt model and Burgers model fit the nanoindentation creep curves well.This study is beneficial to the understanding of the long-term mechanical properties of rock samples from a microscale perspective,which is of great significance to the understanding of localized deformation processes of rocks.展开更多
Decarbonizing power systems is crucial to mitigating climate change impacts and achieving carbon neutrality.Increasing renewable energy supply can reduce greenhouse gas emissions and accelerate the decarbonization pro...Decarbonizing power systems is crucial to mitigating climate change impacts and achieving carbon neutrality.Increasing renewable energy supply can reduce greenhouse gas emissions and accelerate the decarbonization process.However,renewable energy sources(RESs)such as wind and solar power are characterized by intermittency and often non-dispatchability,significantly challenging their high-level integration into power systems.Energy storage is acknowledged as a vital indispensable solution for mitigating the intermittency of renewables such as wind and solar power and boosting the penetrations of renewables.In the CSEE JPES Forum,five well-known experts were invited to give keynote speeches,and the participating experts and scholars had comprehensive exchanges and discussions on energy storage technologies.Specifically,the views on the design,control,performance,and applications of new energy storage technologies,such as the fuel cell vehicle,water electrolysis,and flow battery,in the coordination and operation of power and energy systems were analyzed.The experts also provided experience that could be used to develop energy storage for constructing and decarbonizing new power systems.展开更多
Wave breaking plays an important role in wave-structure interaction. A novel control volume finite element method with adaptive unstructured meshes is employed here to study 3-D breaking waves. The numerical framework...Wave breaking plays an important role in wave-structure interaction. A novel control volume finite element method with adaptive unstructured meshes is employed here to study 3-D breaking waves. The numerical framework consists of a "volume of fluid" type method for the interface capturing and adaptive unstructured meshes to improve computational efficiency. The numerical model is validated against experimental measurements of breaking wave over a sloping beach and is then used to study the breaking wave impact on a vertical circular cylinder on a slope. Detailed complex interfacial structures during wave impact, such as plunging jet formation and splash-up are captured in the simulation, demonstrating the capability of the present method.展开更多
Shale gas reservoirs have fine-grained textures and high organic contents,leading to complex pore structures.Therefore,accurate well-log derived pore size distributions are difficult to acquire for this unconventional...Shale gas reservoirs have fine-grained textures and high organic contents,leading to complex pore structures.Therefore,accurate well-log derived pore size distributions are difficult to acquire for this unconventional reservoir type,despite their importance.However,nuclear magnetic resonance(NMR)logging can in principle provide such information via hydrogen relaxation time measurements.Thus,in this paper,NMR response curves(of shale samples)were rigorously mathematically analyzed(with an Expectation Maximization algorithm)and categorized based on the NMR data and their geology,respectively.Thus the number of the NMR peaks,their relaxation times and amplitudes were analyzed to characterize pore size distributions and lithofacies.Seven pore size distribution classes were distinguished;these were verified independently with Pulsed-Neutron Spectrometry(PNS)well-log data.This study thus improves the interpretation of well log data in terms of pore structure and mineralogy of shale reservoirs,and consequently aids in the optimization of shale gas extraction from the subsurface.展开更多
The Kuqa Depression in Tarim Basin develops fan deltaic,braided river deltaic,and lacustrine sedi-mentation,espeially extensive braided river deltaic sedimentation,in the Cretaceous-Paleogene Sys-tems.However,fine sed...The Kuqa Depression in Tarim Basin develops fan deltaic,braided river deltaic,and lacustrine sedi-mentation,espeially extensive braided river deltaic sedimentation,in the Cretaceous-Paleogene Sys-tems.However,fine sedimentary pattern of the Kuqa Depression and North Tarim Uplift belt since the Middle Cenozoic still needs to be studied.As a combination zone of the Kuqa Depression and North Tarim Uplift,the sedimentary characteristics of the DB area in the southern margin of the east of the Kuqa depression have been paid more attention.To better understand the sedimentary framework in the Kuqa Depression and North Tarim Uplift in the Mesozoic and Cenozoic,through joint geologic-geophysical study,sedime ntation of sandstone at the bottom of the Paleogene Suweiyi Formation in the DB area and its peripheral area is investigated.Sedimentary facies and sandstone distribution in the area are identified through core observation,component analysis,logging interpretation and seismic inversion.Based on seismic facies analysis,sedimentary facies distribution in the area is delineated.Results show that sandstone at the bottom of the Paleogene Suweiyi Fommation in the DB area and its peripheral area is developed in shore shallow lacustrine beach-bar facies.The beach-bar arenaceous sediments are mainly distnibuted in the southern DB area and the shallow lacustrine mudstone is developed in northern DB area,showing obvious north south di ferentiation charactenistics.展开更多
Because of the combination of optimization algorithms and full wave equations, full-waveform inversion(FWI) has become the frontier of the study of seismic exploration and is gradually becoming one of the essential to...Because of the combination of optimization algorithms and full wave equations, full-waveform inversion(FWI) has become the frontier of the study of seismic exploration and is gradually becoming one of the essential tools for obtaining the Earth interior information. However, the application of conventional FWI to pure reflection data in the absence of a highly accurate starting velocity model is difficult. Compared to other types of seismic waves, reflections carry the information of the deep part of the subsurface. Reflection FWI, therefore, is able to improve the accuracy of imaging the Earth interior further. Here, we demonstrate a means of achieving this successfully by interleaving least-squares RTM with a version of reflection FWI in which the tomographic gradient that is required to update the background macro-model is separated from the reflectivity gradient using the Born approximation during forward modeling. This provides a good update to the macro-model. This approach is then followed by conventional FWI to obtain a final high-fidelity high-resolution result from a poor starting model using only reflection data.Further analysis reveals the high-resolution result is achieved due to a deconvolution imaging condition implicitly used by FWI.展开更多
A least-squares reverse-time migration scheme is presented for reflectivity imaging. Based on an accurate reflection modeling formula, this scheme produces amplitude-preserved stacked reflectivity images with zero pha...A least-squares reverse-time migration scheme is presented for reflectivity imaging. Based on an accurate reflection modeling formula, this scheme produces amplitude-preserved stacked reflectivity images with zero phase. Spatial preconditioning, weighting and the Barzilai-Borwein method are applied to speed up the convergence of the least-squares inversion. In addition, this scheme compensates the effect of ghost waves to broaden the bandwidth of the reflectivity images. Furthermore, roughness penalty constraint is used to regularize the inversion, which in turn stabilizes inversion and removes high-wavenumber artifacts and mitigates spatial aliasing. The examples of synthetic and field datasets demonstrate the scheme can generate zerophase reflectivity images with broader bandwidth, higher resolution, fewer artifacts and more reliable amplitudes than conventional reverse-time migration.展开更多
基金the support of the Leverhulme Centre for Wildfires,Environment and Society through the Leverhulme Trust(RC-2018-023)Sibo Cheng,César Quilodran-Casas,and Rossella Arcucci acknowledge the support of the PREMIERE project(EP/T000414/1)+5 种基金the support of EPSRC grant:PURIFY(EP/V000756/1)the Fundamental Research Funds for the Central Universitiesthe support of the SASIP project(353)funded by Schmidt Futures–a philanthropic initiative that seeks to improve societal outcomes through the development of emerging science and technologiesDFG for the Heisenberg Programm Award(JA 1077/4-1)the National Natural Science Foundation of China(61976120)the Natural Science Key Foundat ion of Jiangsu Education Department(21KJA510004)。
文摘Data assimilation(DA)and uncertainty quantification(UQ)are extensively used in analysing and reducing error propagation in high-dimensional spatial-temporal dynamics.Typical applications span from computational fluid dynamics(CFD)to geoscience and climate systems.Recently,much effort has been given in combining DA,UQ and machine learning(ML)techniques.These research efforts seek to address some critical challenges in high-dimensional dynamical systems,including but not limited to dynamical system identification,reduced order surrogate modelling,error covariance specification and model error correction.A large number of developed techniques and methodologies exhibit a broad applicability across numerous domains,resulting in the necessity for a comprehensive guide.This paper provides the first overview of state-of-the-art researches in this interdisciplinary field,covering a wide range of applications.This review is aimed at ML scientists who attempt to apply DA and UQ techniques to improve the accuracy and the interpretability of their models,but also at DA and UQ experts who intend to integrate cutting-edge ML approaches to their systems.Therefore,this article has a special focus on how ML methods can overcome the existing limits of DA and UQ,and vice versa.Some exciting perspectives of this rapidly developing research field are also discussed.Index Terms-Data assimilation(DA),deep learning,machine learning(ML),reduced-order-modelling,uncertainty quantification(UQ).
文摘Failure of irregular rock samples may provide implications in the rapid estimation of rock strength,which is imperative in rock engineering practice.In this work,analytical,experimental and numerical investigations were carried out to study the mechanical properties and failure characteristics of rock spheres under paired point loads.Analytical solutions indicted that with the increase in sample size(contact angle)and decrease in Poisson’s ratio,the uneven tensile stress in theta direction decreased.Then laboratory experiments were carried out to investigate the load characteristics and failure mode of spherical marble samples with different sizes subjected to a pair of diametral point loads.The discrete element method(DEM)was adopted to study the failure process of rock spheres.The effect of the sphere diameter on the point load contact angle was examined in terms of peak load,crushed zone distribution and energy dissipation.Experimental and numerical results showed that the samples primarily fail in tension,with crushed zones formed at both loading points.With increase in the sample size,the contact angle,crushed area and total work increase.As the specimen diameter increases from 30 mm to 50 mm,the peak load on the specimen increases from 3.6 kN to 8.8 kN,and the percentage of crushed zone(ratio of crushing zone to sample radius,d/r)increased from 0.191 to 0.262.The results of the study have implications for understanding the failure of irregular rock specimens under point loading conditions and their size effects.
基金Project(51404251)supported by the National Natural Science Foundation of ChinaProject(BK20140198)supported by the Natural Science Foundation of Jiangsu Province of China+1 种基金Project(2014XT01)supported by the Fundamental Research Funds for the Central UniversitiesChina
文摘In order to explore the pressure relief and structure stability mechanism of lateral cantilever structure in the stope under the direct coverage of thick hard roof and its impact on the gob-side entry retaining, a lateral cantilever fractured structural mechanical model was established on the basis of clarification for the stress environment of gob-side entry retaining, and the equation of roof given deformation and the balance judgment for fracture block were obtained. The optimal cantilever length was proposed based on the comparison of roof structural characteristics and the stress, deformation law of surrounding rocks under six different cantilever lengths by numerical simulation method. Double stress peaks exist on the sides of gob-side entry retaining and the entry located in the low stress area. The pressure of gob-side entry retaining can be relieved by reducing the cantilever length. However, due to the impact of arch structure of rock beam, unduly short cantilever would result in insufficient pressure relief and unduly long cantilever would bring larger roof stress which results in intense deformation. Therefore, there is optimal cantilever length, which was 7-8 m in this project that enables to achieve the minimum deformation and the most stabilized rock structure of entry retaining. An engineering case of gob-side entry retaining with the direct coverage of 10 m thick hard limestone roof was put forward, and the measured data verified the reasonability of conclusion.
基金the Natural Environment Research Council for the funding received for project NE/R018065/1the Royal Society for funding this research through fellowship UF160443the Swedish Nuclear Fuel and Waste Management Company(SKB)for partially funding this research。
文摘The flow properties of geomechanically generated discrete fracture networks are examined in the context of channelling.Fracture networks are generated by growing fractures in tension,modelling the low permeability rock as a linear elastic material.Fractures are modelled as discrete surfaces which grow quasi-statically within a three-dimensional(3D)volume.Fractures may have their locations specified as a simulation input,or be generated as a function of damage,quantified using the local variation in equivalent strain.The properties of the grown networks are shown to be a product of in situ stress,relative orientation of initial flaws,and competitive process of fracture interaction and growth.Fractures grow preferentially in the direction perpendicular to the direction of maximum tension and may deviate from this path due to mechanical fracture interaction.Flow is significantly channelled through a subset of the fractures in the full domain,consistent with observations of other real and simulated fractures.As the fracture networks grow,small changes in the geometry of the fractures lead to large changes in the locations and scale of primary flow channels.The flow variability and formation of channels are examined for two growing networks,one with a fixed amount of fractures,and another with nucleating fractures.The interaction between fractures is shown to modify the local stress field,and in turn the aperture of the fractures.Pathways for single-phase flow are the results of hydro-mechanical effects in fracture networks during growth.These are the results of changes to the topology of the network as well as the result of mechanical self-organisation which occurs during interaction leading to growth and intersection.
文摘The stable isotope systems of the transition metals potentially provide constraints on the current and past operation of the biological pump,and on the state of ocean redox in Earth history.Here we focus on two exemplar metals,nickel(Ni)and zinc(Zn).The oceanic dissolved pool of both elements is isotopically heavier than the known inputs,implying an output with light isotope compositions.The modern oceanic cycle of both these elements is dominated by biological uptake into photosynthesised organic matter and output to sediment.It is increasingly clear,however,that such uptake is associated with only very minor isotope fractionation.We suggest that the isotopic balance is instead closed by the sequestration of light isotopes to sulphide in anoxic and organic-rich sediments,so that it is ocean chemistry that controls these isotope systems,and suggesting a different but equally interesting array of questions in Earth history that can be addressed with these systems.
文摘Deep coal seams are one of the world’s most widespread deposits for carbon dioxide(C02)disposal and are generally located near large point sources of CO_(2)emissions.The injection of CO_(2)into coal seams has great potential to sequester CO_(2)while simultaneously enhancing coalbed methane(CO_(2)-ECBM)recovery.Pilot tests of CO_(2)-ECBM have been conducted in coal seams worldwide with favorable early results.However,one of the main technical barriers in coal seams needs to be resolved:Injecting CO_(2)reduces coal permeability and well injectivity.Here,using in situ synchrotron X-ray microtomography,we provide the first observational evidence that injecting nitrogen(N_(2))can reverse much of this lost permeability by reopening fractures that have closed due to coal swelling induced by CO_(2)adsorption.Our findings support the notion that injecting minimally treated flue gas-a mixture of mainly N_(2) and CO_(2)-is an attractive alternative for ECBM recovery instead of pure CO_(2)injection in deep coal seams.Firstly,flue gas produced by power plants could be directly injected after particulate removal,thus avoiding high CO_(2)-separation costs.Secondly,the presence of N_(2)makes it possible to maintain a sufficiently high level of coal permeability.These results suggest that flue-gas ECBM for deep coal seams may provide a promising path toward net-zero emissions from coal mines.
基金the financial support from EPSRC(EP/P024807/1,EP/M014045/1,EP/S000933/1 and EP/N009924/1)by the EPSRC energy storage for low carbon grids project(EP/K002252/1)+3 种基金the EPSRC Joint UK-India Clean Energy center(JUICE)(EP/P003605/1)the Integrated Development of Low-Carbon Energy Systems(IDLES)project(EP/R045518/1)the Innovate UK BAFTA project,the Innovate UK for Advanced Battery Lifetime Extension(ABLE)project for funding underthe China Scholarship Council。
文摘Mixed ionic electronic conductors(MIECs)have attracted increasing attention as anode materials for solid oxide fuel cells(SOFCs)and they hold great promise for lowering the operation temperature of SOFCs.However,there has been a lack of understanding of the performance-limiting factors and guidelines for rational design of composite metal-MIEC electrodes.Using a newly-developed approach based on 3 D-tomography and electrochemical impedance spectroscopy,here for the first time we quantify the contribution of the dual-phase boundary(DPB)relative to the three-phase boundary(TPB)reaction pathway on real MIEC electrodes.A new design strategy is developed for Ni/gadolinium doped ceria(CGO)electrodes(a typical MIEC electrode)based on the quantitative analyses and a novel Ni/CGO fiber-matrix structure is proposed and fabricated by combining electrospinning and tape-casting methods using commercial powders.With only 11.5 vol%nickel,the designer Ni/CGO fiber-matrix electrode shows 32%and 67%lower polarization resistance than a nano-Ni impregnated CGO scaffold electrode and conventional cermet electrode respectively.The results in this paper demonstrate quantitatively using real electrode structures that enhancing DPB and hydrogen kinetics are more efficient strategies to enhance electrode performance than simply increasing TPB.
文摘We measured Electron Spin Resonance(ESR) signal intensity and Crystallinity Index(CI) of fine-(<16 m) and coarse-grained(>16 m) quartz in surface samples collected from the Taklimakan desert in western China, the Badain Juran, Tengger and Mu Us deserts in northern China, and the Gobi desert in southern Mongolia. The aim was to test if these two parameters can be used as the accurate provenance tracers of Asian dust. Coarse-grained quartz displays spatial variability of the ESR signal intensity within each desert and of the CI within the Taklimakan and the Mongolian Gobi deserts. Coarse-grained quartz derived from the Mongolian Gobi and northern China deserts can be differentiated from the Taklimakan desert using the ESR signal intensity. By contrast, variability of ESR signal intensity and CI of fine-grained quartz particles is distinct within the Taklimakan and the Mongolian Gobi deserts, but less significant in three deserts of northern China. Spatially, fine-grained dust originating from three major Asian dust sources(i.e., the Gobi-sandy deserts in western China, northern China and southern Mongolia) can be differentiated using the combination of the ESR and CI signals. Our results suggest that combination of the ESR and CI results can discriminate the sources of fine-grained dust better than coarse-grained dust, providing an effective approach to trace the provenance of fine-grained dust depositions on the land and in the ocean.
基金国家自然科学基金项目(41176052、41576035、41276050)南方海洋科学与工程广东省实验室(广州)人才团队引进重大专项(GML2019ZD0104)+2 种基金广东省自然科学基金面上项目(S2013010014523)中国科学院南海生态环境工程创新研究院创新发展基金项目(ISEE2018PY02)the U.S.Department of Energy(DE-AC02-05CH11231)~~
基金financial support for this work contributed by the National Key Research and Development Program of China (grant numbers 2016YFC0600101 and 2016YFC 0600201)the National Natural Science Foundation of China (grant numbers 41874065, 41604076, 41674102, 41674095, 41522401, 41574082, and 41774097)
文摘In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation in inversion. A feasible way to avoid the excessive storage demand is to reconstruct the source wavefield backward in time by storing the entire history of the wavefield in perfectly matched layers. In this paper, we make full use of the elementwise global property of the Laplace operator of the spectral element method (SEM) and propose an efficient source wavefield reconstruction method at the cost of storing the wavefield history only at single boundary layer nodes. Numerical experiments indicate that the accuracy of the proposed method is identical to that of the conventional method and is independent of the order of the Lagrange polynomials, the element type, and the temporal discretization method. In contrast, the memory-saving ratios of the conventional method versus our method is at least N when using either quadrilateral or hexahedron elements, respectively, where N is the order of the Lagrange polynomials used in the SEM. A higher memorysaving ratio is achieved with triangular elements versus quadrilaterals. The new method is applied to reverse time migration by considering the Marmousi model as a benchmark. Numerical results demonstrate that the method is able to provide the same result as the conventional method but with about 1/25 times lower storage demand. With the proposed wavefield reconstruction method, the storage demand is dramatically reduced;therefore, in-core memory storage is feasible even for large-scale three-dimensional adjoint inversion problems.
文摘1.Historical context Flow in porous media is encountered in a wide variety of natural and engineered settings,including gas exchange in lungs and blood vessels,water extraction from underground aquifers,hydrocarbon production,and microfluidics in drug delivery and food manufacturing.Underground rock is a porous medium that holds water and hydrocarbons,and is a potential storage site for carbon dioxide(CO_(2))and hydrogen.Many manufactured devices also contain porous media to allow the transport of gases and water,including fuel cells,electrolyzers,and electrocatalysts for CO_(2) reduction.Their performance is limited by how rapidly two fluid phases can flow through a porous layer.
基金financial support provided by the China Scholarship Council(No.201706420055)funded by various projects,namely,China Postdoctoral Science Foundation(No.2020M673142)+3 种基金Chongqing Science and Technology Bureau Foundation(No.cstc2020jcyj-bsh0059)Chongqing Human Resources and Social Security Bureau Foundation(No.cx2019100)Chongqing Science and Technology Bureau Foundation(No.cstc2020jcyj-zdxm0122)National Natural Science Foundation of China(No.U19B2009)。
文摘The creep phenomenon of rocks is quite complex and the creep mechanisms are far from being well understood.Although laboratory creep tests have been carried out to determine the creep deformation of various rocks,these tests are expensive and time-consuming.Nanoindentation creep tests,as an alternative method,can be performed to investigate the mechanical and viscoelastic properties of granite samples.In this study,the reduced Young’s modulus,hardness,fracture toughness,creep strain rate,stress exponent,activation volume and maximum creep displacement of common rock-forming minerals of granite were calculated from nanoindentation results.It was found that the hardness decreases with the increase of holding time and the initial decrease in hardness was swift,and then it decreased slowly.The stress exponent values obtained were in the range from 4.5 to 22.9,which indicates that dislocation climb is the creep deformation mechanism.In addition,fracture toughness of granite’s rock-forming minerals was calculated using energy-based method and homogenization method was adopted to upscale the micro-scale mechanical properties to macro-scale mechanical properties.Last but not least,both three-element Voigt model and Burgers model fit the nanoindentation creep curves well.This study is beneficial to the understanding of the long-term mechanical properties of rock samples from a microscale perspective,which is of great significance to the understanding of localized deformation processes of rocks.
文摘Decarbonizing power systems is crucial to mitigating climate change impacts and achieving carbon neutrality.Increasing renewable energy supply can reduce greenhouse gas emissions and accelerate the decarbonization process.However,renewable energy sources(RESs)such as wind and solar power are characterized by intermittency and often non-dispatchability,significantly challenging their high-level integration into power systems.Energy storage is acknowledged as a vital indispensable solution for mitigating the intermittency of renewables such as wind and solar power and boosting the penetrations of renewables.In the CSEE JPES Forum,five well-known experts were invited to give keynote speeches,and the participating experts and scholars had comprehensive exchanges and discussions on energy storage technologies.Specifically,the views on the design,control,performance,and applications of new energy storage technologies,such as the fuel cell vehicle,water electrolysis,and flow battery,in the coordination and operation of power and energy systems were analyzed.The experts also provided experience that could be used to develop energy storage for constructing and decarbonizing new power systems.
基金the financial support by the National Natural Science Foundation of China (Grant No. 51490673)the Open Awards of the State Key Laboratory of Coastal and Offshore Engineering+1 种基金funded by the EPSRC MEMPHIS multiphase Programme (Grant No. EP/K003976/1)funding from the European Union Seventh Framework Programme (FP7/20072013) under grant agreement No. 603663 for the research project PEARL (Preparing for Extreme and Rare events in coasta L regions)
文摘Wave breaking plays an important role in wave-structure interaction. A novel control volume finite element method with adaptive unstructured meshes is employed here to study 3-D breaking waves. The numerical framework consists of a "volume of fluid" type method for the interface capturing and adaptive unstructured meshes to improve computational efficiency. The numerical model is validated against experimental measurements of breaking wave over a sloping beach and is then used to study the breaking wave impact on a vertical circular cylinder on a slope. Detailed complex interfacial structures during wave impact, such as plunging jet formation and splash-up are captured in the simulation, demonstrating the capability of the present method.
基金National Natural Science Foundation of China(41902145)Natural Science Basic Research Plan in Shaanxi Province of China(2020JQ-594)+2 种基金Young Talent fund of University Association for Science and Technology in Shaanxi,China(20180701)National and Local Joint Engineering Research Center for Carbon Capture Utilization and Sequestration at Northwest University in China.The measurements were performed using the mCT system of the National Geosequestration Laboratory(NGL)of Australia.Funding for the facilities was provided by the Australian Federal Governmentsupported by the Pawsey Supercomputing Centre,who provided the Avizo 9.5 image processing software and workstation,with funding from the Australian Government and the Government of Western Australia.
文摘Shale gas reservoirs have fine-grained textures and high organic contents,leading to complex pore structures.Therefore,accurate well-log derived pore size distributions are difficult to acquire for this unconventional reservoir type,despite their importance.However,nuclear magnetic resonance(NMR)logging can in principle provide such information via hydrogen relaxation time measurements.Thus,in this paper,NMR response curves(of shale samples)were rigorously mathematically analyzed(with an Expectation Maximization algorithm)and categorized based on the NMR data and their geology,respectively.Thus the number of the NMR peaks,their relaxation times and amplitudes were analyzed to characterize pore size distributions and lithofacies.Seven pore size distribution classes were distinguished;these were verified independently with Pulsed-Neutron Spectrometry(PNS)well-log data.This study thus improves the interpretation of well log data in terms of pore structure and mineralogy of shale reservoirs,and consequently aids in the optimization of shale gas extraction from the subsurface.
文摘The Kuqa Depression in Tarim Basin develops fan deltaic,braided river deltaic,and lacustrine sedi-mentation,espeially extensive braided river deltaic sedimentation,in the Cretaceous-Paleogene Sys-tems.However,fine sedimentary pattern of the Kuqa Depression and North Tarim Uplift belt since the Middle Cenozoic still needs to be studied.As a combination zone of the Kuqa Depression and North Tarim Uplift,the sedimentary characteristics of the DB area in the southern margin of the east of the Kuqa depression have been paid more attention.To better understand the sedimentary framework in the Kuqa Depression and North Tarim Uplift in the Mesozoic and Cenozoic,through joint geologic-geophysical study,sedime ntation of sandstone at the bottom of the Paleogene Suweiyi Formation in the DB area and its peripheral area is investigated.Sedimentary facies and sandstone distribution in the area are identified through core observation,component analysis,logging interpretation and seismic inversion.Based on seismic facies analysis,sedimentary facies distribution in the area is delineated.Results show that sandstone at the bottom of the Paleogene Suweiyi Fommation in the DB area and its peripheral area is developed in shore shallow lacustrine beach-bar facies.The beach-bar arenaceous sediments are mainly distnibuted in the southern DB area and the shallow lacustrine mudstone is developed in northern DB area,showing obvious north south di ferentiation charactenistics.
基金supported by the National Natural Science Foundation of China(Grant No.41504106&41274099)the Science Foundation of China University of Petroleum(Beijing)(Grant No.2462015YJRC012)State Laboratory of Petroleum Resource and Prospecting(Grant No.PRP/indep-3-1508)
文摘Because of the combination of optimization algorithms and full wave equations, full-waveform inversion(FWI) has become the frontier of the study of seismic exploration and is gradually becoming one of the essential tools for obtaining the Earth interior information. However, the application of conventional FWI to pure reflection data in the absence of a highly accurate starting velocity model is difficult. Compared to other types of seismic waves, reflections carry the information of the deep part of the subsurface. Reflection FWI, therefore, is able to improve the accuracy of imaging the Earth interior further. Here, we demonstrate a means of achieving this successfully by interleaving least-squares RTM with a version of reflection FWI in which the tomographic gradient that is required to update the background macro-model is separated from the reflectivity gradient using the Born approximation during forward modeling. This provides a good update to the macro-model. This approach is then followed by conventional FWI to obtain a final high-fidelity high-resolution result from a poor starting model using only reflection data.Further analysis reveals the high-resolution result is achieved due to a deconvolution imaging condition implicitly used by FWI.
基金partly supported by the National Naural Science Foundation of China(Grant No.41272099)the Science Foundation of China University of Petroleum,Beijing(Grant No.2462015YJRC012)
文摘A least-squares reverse-time migration scheme is presented for reflectivity imaging. Based on an accurate reflection modeling formula, this scheme produces amplitude-preserved stacked reflectivity images with zero phase. Spatial preconditioning, weighting and the Barzilai-Borwein method are applied to speed up the convergence of the least-squares inversion. In addition, this scheme compensates the effect of ghost waves to broaden the bandwidth of the reflectivity images. Furthermore, roughness penalty constraint is used to regularize the inversion, which in turn stabilizes inversion and removes high-wavenumber artifacts and mitigates spatial aliasing. The examples of synthetic and field datasets demonstrate the scheme can generate zerophase reflectivity images with broader bandwidth, higher resolution, fewer artifacts and more reliable amplitudes than conventional reverse-time migration.