An operating simulation mode for multireach canal system in series under gate regulating is established. The discharge feedforward plus water lever feedback PID controller was adopted. And the performance indicators a...An operating simulation mode for multireach canal system in series under gate regulating is established. The discharge feedforward plus water lever feedback PID controller was adopted. And the performance indicators are introduced to evaluate the control effects of canal operation. Influence on dynamic response process and dynamic performance as well as the control actions of canal system are studied according to the variation of the gate deadband and water level deadband. The results showed that the larger the gate deadband, the worse, the control effect of canal system over water level, and so done the stable process of gate discharge, yet the stable process of gate opening was getting better, while the upstream gate regulation of each canal pool was less at the time when canal stabilizes; the closer to downstream end of canal, the lager steady-state error of this canal pool was, and the lager influence of gate deadband on the steady-state error as well as response time of this canal pool was; the canal system can not acquire superiority in both the control performance and control actions when only the water level deadband was established.展开更多
The measured data of vibrations induced by excavation of deep-buried cavern and open pit with method of bench blasting were analyzed by Fourier Transform and Wavelet Transform,and the characteristics of vibrations ind...The measured data of vibrations induced by excavation of deep-buried cavern and open pit with method of bench blasting were analyzed by Fourier Transform and Wavelet Transform,and the characteristics of vibrations induced under these two circumstances were studied.It is concluded that with the similar rock condition and drilling-blasting parameters,vibration induced by bench blasting in deep-buried cavern has a higher main frequency and more scattered energy distribution than that in open pit.The vibration induced by bench blasting in open pit is mainly originated from the blast load,while the vibration induced by bench blasting in deep-buried cavern is the superposition of vibrations induced by blast load and transient release of in-situ stress.The vibration induced by transient release of in-situ stress increases with the stress level.展开更多
This study investigates the factors affecting the rock-breaking efficiency of the TBM disc cutter in deep rock excavation,including confining pressure,penetration,cutter spacing,and revolution speed.The finite element...This study investigates the factors affecting the rock-breaking efficiency of the TBM disc cutter in deep rock excavation,including confining pressure,penetration,cutter spacing,and revolution speed.The finite element method is employed to formulate a rock-breaking model of the rotary disc cutters and a numerical simulation is also implemented.The rock breaking effect,rock breaking volume,and rock breaking specific energy consumption under different combinations of the factors are investigated.An orthogonal test of four factors at four levels was constructed.Based on the test results and range analysis in the process of deep rock mass breaking,the order of sensitivity of each influencing factor with respect to the rock breaking specific energy for the disc cutter is cutter spacing>revolution speed>penetration>confining pressure.By constructing a numerical simulation comparison scheme,the orthogonal test results are analyzed and corroborated,and the rock breaking law and rock breaking efficiency under different influencing factors are derived.Finally,the sensitivity of different influencing factors on the rock-breaking efficiency is verified.展开更多
The hydration and thermal properties of cement-based materials containing various proportions of limestone powder as a partial replacement for ordinary Portland cement, were investigated and reported. Both compressive...The hydration and thermal properties of cement-based materials containing various proportions of limestone powder as a partial replacement for ordinary Portland cement, were investigated and reported. Both compressive and flexural strengths of cement mortar with various contents of limestone powder were tested to study the influence of limestone powder on the strength development of resulting mixtures. The hydration heat and its rate of evolution were also tested, which clearly showed that the replacement percentage of limestone powder had significant effects on the total hydration heat but only a modest influence on the rate of heat evolution of cement-limestone binder. Importantly, the reduction coefficient of limestone powder on the hydration heat, needed for estimation of adiabatic temperature rise of cement-limestone binder, was found to be approximately 0.51. Fundamental thermal properties of these concrete mixtures containing limestone powder were also studied. Increasing the percentage of limestone powder resulted in a significant reduction in the adiabatic temperature rise but only a slight increase in other thermal properties such as thermal conductivity, thermal diffusivity and specific heat. In addition, thermal analysis using finite-element modelling indicated that inclusion of limestone powder did not significantly affect the rate of temperature rise nor the occurrence time of the highest temperature at early ages.展开更多
The controlled volume method of operation is especially suitable for large-scale water delivery canal system with complex operation requirements. An operating simulation model based on the storage volume control metho...The controlled volume method of operation is especially suitable for large-scale water delivery canal system with complex operation requirements. An operating simulation model based on the storage volume control method for multi-reach canal system in series was established. In allusion to the deficiency of existing controlled volume algorithm, the improved controlled volume algorithm of the whole canal pools was proposed, and the simulation results indicated that the storage volume and water level of each canal pool could be accurately controlled after the improved algorithm had been adopted. However, for some typical discharge demand operating conditions, if the previously mentioned algorithm was adopted, then it certainly would cause some unnecessary gate adjustments, and consequently the disturbed canal pools would be increased. Therefore, the idea of controlled volume operation method of continuous canal pools was proposed, and corresponding algorithm was designed. Through simulating practical project, the results indicated that the new controlled volume algorithm proposed for typical operating conditions could comparatively and obviously reduce the number of regulated check gates and disturb canal pools for some typical discharge demand operating conditions, thus the control efficiency of canal system could be improved.展开更多
Geotechnical stability analyses based on classical continuum may lead to remarkable underestimations on geotechnical safety.To attain better estimations on geotechnical stability,the micro-polar continuum is employed ...Geotechnical stability analyses based on classical continuum may lead to remarkable underestimations on geotechnical safety.To attain better estimations on geotechnical stability,the micro-polar continuum is employed so that its internal characteristic length(lc)can be utilized to model the shear band width.Based on two soil slope examples,the role of internal characteristic length in modeling the shear band width of geomaterial is investigated by the second-order cone programming optimized micro-polar continuum finite element method.It is recognized that the underestimation on factor of safety(FOS)calculated from the classical continuum tends to be more pronounced with the increase of lc.When the micro-polar continuum is applied,the shear band dominated by lc is almost kept unaffected as long as the adopted meshes are fine enough,but it does not generally present a slip surface like in the cases from the classical continuum,indicating that the micro-polar continuum is capable of capturing the non-local geotechnical failure characteristic.Due to the coupling effects of lc and strain softening,softening behavior of geomaterial tends to be postponed.Additionally,the bearing capacity of a geotechnical system may be significantly underestimated,if the effects of lc are not modeled or considered in numerical analyses.展开更多
Uniaxial or triaxial compression test of cylindrical rock specimens using rock mechanics testing machine is a basic experimental method to study the strength and deformation characteristics of rock and the development...Uniaxial or triaxial compression test of cylindrical rock specimens using rock mechanics testing machine is a basic experimental method to study the strength and deformation characteristics of rock and the development process of rock fracture. Extensive literature review has been conducted on this issue;meanwhile, experimental and numerical studies have been conducted on the stress-drop effect on the brittleness of rock materials. A plastic flow factor of λ is proposed to describe the stress-drop effect. Evaluation methods of the factor λ corresponding to the four yield criteria of rock mass are proposed. Those four yield criteria are Tresca criterion, von-Mises criterion, Mohr-Coulomb criterion and Drucker-Prager criterion. For simplicity purposes, an engineering approximation approach has been proposed to evaluate the stress-drop with a non-zero strain increment. Numerical simulation results validated the effectiveness of the plastic flow factors λ as well as the engineering approximation approach. Based on the results in this study, finite element code can be programmed for brittle materials with stress-drop, which has the potential to be readily incorporated in finite element codes.展开更多
The hydration mechanism of low quality fly ash in cement-based materials was investigated. The hydration heat of the composite cementitious materials was determined by isothermal calorimetry, and the hydration product...The hydration mechanism of low quality fly ash in cement-based materials was investigated. The hydration heat of the composite cementitious materials was determined by isothermal calorimetry, and the hydration products, quantity, pore structure and morphology were measured by X-ray diffraction(XRD), thermalgravity-differential thermal analysis(TG-DTA), mercury intrusion porosimetry(MIP) and scanning electron microscopy(SEM), respectively. The results indicate that grinding could not only improve the physical properties of the low quality fly ash on particle effect, but also improve hydration properties of the cementitious system from various aspects compared with raw low quality fly ash(RLFA). At the early stage of hydration, the low quanlity fly ash acts as almost inert material; but then at the later stage, high chemical activity, especially for ground low quality fly ash(GLFA), could be observed. It can accelerate the formation of hydration products containing more chemical bonded water, resulting in higher degree of cement hydration, thus denser microstructure and more reasonable pore size distribution, but the hydration heat in total is reduced. It can also delay the induction period, but the accelerating period is shortened and there is little influence on the second exothermic peak.展开更多
A finite volume,multiphase solver in the framework of OpenFOAM was used to calculate the flow field of the cavitating flow over the Clark-Y hydrofoil. This solver used Transport Based Equation Model(TEM) to solve the ...A finite volume,multiphase solver in the framework of OpenFOAM was used to calculate the flow field of the cavitating flow over the Clark-Y hydrofoil. This solver used Transport Based Equation Model(TEM) to solve the liquid volume fraction,and utilized volume of fluid(VOF) technique to predict the interface between liquid and vapor phases. The simulation was designed to study the cavitation shedding and different fluid characteristics in the cloud cavitation regime when adopting two different Large Eddy Simulation(LES) models,namely,one equation eddy viscosity(one EqEddy) model and Smagorinsky model. It is shown that these two models can be used to study the cavitation shedding dynamics and predict the velocity profiles.展开更多
Regarding excavation-damaged zone （EDZ） around underground opening as non-homogeneous rockmass with spatial deterioration effect on stuffiness and strength, a parametric model of EDZ using radius-displacement-depend...Regarding excavation-damaged zone （EDZ） around underground opening as non-homogeneous rockmass with spatial deterioration effect on stuffiness and strength, a parametric model of EDZ using radius-displacement-dependent deformation modulus （RDDM） was proposed. Considering the nonlinearity characteristic of deformation and locality otherness of surrounding rock, deterioration parameter field of deformation modulus of rockmass around opening was quantitatively calculated through a given function. Applicability for multi-cavern condition and parameter sensibility of the model was analyzed by numerical experiments using synthetic data. Furthermore, the model was applied to identify EDZ of underground caverns of Pubugou hydropower station by calculating deterioration parameter field. Based on the parametric analysis of spatial effect and geological investigation, it is recognized that large radial deformation of deep fractured rock at the spandrel position and insufficient supporting bolts mainly result in great deformation pressure to act on the shotcrete and cause partial crack and spalling. It is shown that deterioration parameter field along the longitudinal axis of main powerhouse is evidently non-homogeneous in space and distributes exponentially along the radius from the opening. The model provides a simple and convenient way to identify the EDZ in the working state for rapid construction feedback analysis and support optimization of underground cavem from quantitative point of view and also aids in interpreting monitoring displacements and estimating support requirements.展开更多
The seepage under a dam foundation is mainly controlled by the performance of the curtain.Its anti-seepage behavior may be weakened by the long term physic-chemical actions from groundwater.According to seepage hydrau...The seepage under a dam foundation is mainly controlled by the performance of the curtain.Its anti-seepage behavior may be weakened by the long term physic-chemical actions from groundwater.According to seepage hydraulics and geochemistry theory,a coupling model for assessing the behavior of the curtain under a dam foundation is set up,which consists of seepage module,solute transport module,geochemistry module and curtain erosion module,solved by FEM.A case study was carried out.The result shows that the curtain efficiency is weakened all the time.Aqueous calcium from the curtain is always in dissolution during the stress period for simulation,which leads to the increasing amount in groundwater reaching 846.35-865.312 g/m3.Within the domain,reaction extent differs in different parts of the curtain.The dissolution of Ca(OH)2 accounts to 877.884 g/m3 near the bottom and is much higher than that of the other parts.The erosion is much more serious near the bottom of the curtain than the other parts,which is the same to the upstream and downstream.Calcium dissolution is mainly controlled by hydraulic condition and dispersion,and it varies in a non-linear way within the domain.展开更多
The bubble dynamic near a rigid wall with a wall jet was investigated by codynamics(CFD)method with the volume of fluid(VOF)model,which had been validated by vious experimental data.The effects of different velocities...The bubble dynamic near a rigid wall with a wall jet was investigated by codynamics(CFD)method with the volume of fluid(VOF)model,which had been validated by vious experimental data.The effects of different velocities of the wall jet and ditances on the bubble dynamics were studied.The results show that the bubble is squjet due to more force added on the bubble.When the velocity of the wall jet increa,the wall anthe pressure along the wall at collapse time increase because of the extra push indAs the stand-off distance increases,the pressure along the wall first increases then decrethe distance from the bubble to the wall increases.展开更多
基金Supported by the National Natural Science Foundation of China (59879016)The Key Projects in the National Science & Technology Pillar Program during the 11th Five-Year Plan Period of China (2006BAB04A12)
文摘An operating simulation mode for multireach canal system in series under gate regulating is established. The discharge feedforward plus water lever feedback PID controller was adopted. And the performance indicators are introduced to evaluate the control effects of canal operation. Influence on dynamic response process and dynamic performance as well as the control actions of canal system are studied according to the variation of the gate deadband and water level deadband. The results showed that the larger the gate deadband, the worse, the control effect of canal system over water level, and so done the stable process of gate discharge, yet the stable process of gate opening was getting better, while the upstream gate regulation of each canal pool was less at the time when canal stabilizes; the closer to downstream end of canal, the lager steady-state error of this canal pool was, and the lager influence of gate deadband on the steady-state error as well as response time of this canal pool was; the canal system can not acquire superiority in both the control performance and control actions when only the water level deadband was established.
基金Project(2010CB732003) supported by the National Basic Research Program of ChinaProjects(50725931,50779050 and 50909077) supported by the National Natural Science Foundation of China
文摘The measured data of vibrations induced by excavation of deep-buried cavern and open pit with method of bench blasting were analyzed by Fourier Transform and Wavelet Transform,and the characteristics of vibrations induced under these two circumstances were studied.It is concluded that with the similar rock condition and drilling-blasting parameters,vibration induced by bench blasting in deep-buried cavern has a higher main frequency and more scattered energy distribution than that in open pit.The vibration induced by bench blasting in open pit is mainly originated from the blast load,while the vibration induced by bench blasting in deep-buried cavern is the superposition of vibrations induced by blast load and transient release of in-situ stress.The vibration induced by transient release of in-situ stress increases with the stress level.
基金Project(51979156)supported by the National Natural Science Foundation of ChinaProject(tsqn202103087)supported by the Young Taishan Scholars,ChinaProject(2019KJG015)supported by the Youth Innovation Technology Project of Higher School in Shandong Province,China。
文摘This study investigates the factors affecting the rock-breaking efficiency of the TBM disc cutter in deep rock excavation,including confining pressure,penetration,cutter spacing,and revolution speed.The finite element method is employed to formulate a rock-breaking model of the rotary disc cutters and a numerical simulation is also implemented.The rock breaking effect,rock breaking volume,and rock breaking specific energy consumption under different combinations of the factors are investigated.An orthogonal test of four factors at four levels was constructed.Based on the test results and range analysis in the process of deep rock mass breaking,the order of sensitivity of each influencing factor with respect to the rock breaking specific energy for the disc cutter is cutter spacing>revolution speed>penetration>confining pressure.By constructing a numerical simulation comparison scheme,the orthogonal test results are analyzed and corroborated,and the rock breaking law and rock breaking efficiency under different influencing factors are derived.Finally,the sensitivity of different influencing factors on the rock-breaking efficiency is verified.
基金Project(51579192) supported by the National Natural Science Foundation of ChinaProject(2013BC0359001) supported by the National Basic Research Program of ChinaProject(201506270058) supported by China Scholarship Council
文摘The hydration and thermal properties of cement-based materials containing various proportions of limestone powder as a partial replacement for ordinary Portland cement, were investigated and reported. Both compressive and flexural strengths of cement mortar with various contents of limestone powder were tested to study the influence of limestone powder on the strength development of resulting mixtures. The hydration heat and its rate of evolution were also tested, which clearly showed that the replacement percentage of limestone powder had significant effects on the total hydration heat but only a modest influence on the rate of heat evolution of cement-limestone binder. Importantly, the reduction coefficient of limestone powder on the hydration heat, needed for estimation of adiabatic temperature rise of cement-limestone binder, was found to be approximately 0.51. Fundamental thermal properties of these concrete mixtures containing limestone powder were also studied. Increasing the percentage of limestone powder resulted in a significant reduction in the adiabatic temperature rise but only a slight increase in other thermal properties such as thermal conductivity, thermal diffusivity and specific heat. In addition, thermal analysis using finite-element modelling indicated that inclusion of limestone powder did not significantly affect the rate of temperature rise nor the occurrence time of the highest temperature at early ages.
基金Supported by the Governmental Public Industry Research Special Funds for Projects of MWR (200901002,200901003,200901006)Key Projects in the National Science & Technology Pillar Program During the 11th Five-year Plan Period of China (2006BAB04A12)
文摘The controlled volume method of operation is especially suitable for large-scale water delivery canal system with complex operation requirements. An operating simulation model based on the storage volume control method for multi-reach canal system in series was established. In allusion to the deficiency of existing controlled volume algorithm, the improved controlled volume algorithm of the whole canal pools was proposed, and the simulation results indicated that the storage volume and water level of each canal pool could be accurately controlled after the improved algorithm had been adopted. However, for some typical discharge demand operating conditions, if the previously mentioned algorithm was adopted, then it certainly would cause some unnecessary gate adjustments, and consequently the disturbed canal pools would be increased. Therefore, the idea of controlled volume operation method of continuous canal pools was proposed, and corresponding algorithm was designed. Through simulating practical project, the results indicated that the new controlled volume algorithm proposed for typical operating conditions could comparatively and obviously reduce the number of regulated check gates and disturb canal pools for some typical discharge demand operating conditions, thus the control efficiency of canal system could be improved.
基金Projects(2017YFC0804602,2017YFC0404802)supported by the National Key R&D Program of ChinaProject(2019JBM092)supported by the Fundamental Research Funds for the Central Universities,China。
文摘Geotechnical stability analyses based on classical continuum may lead to remarkable underestimations on geotechnical safety.To attain better estimations on geotechnical stability,the micro-polar continuum is employed so that its internal characteristic length(lc)can be utilized to model the shear band width.Based on two soil slope examples,the role of internal characteristic length in modeling the shear band width of geomaterial is investigated by the second-order cone programming optimized micro-polar continuum finite element method.It is recognized that the underestimation on factor of safety(FOS)calculated from the classical continuum tends to be more pronounced with the increase of lc.When the micro-polar continuum is applied,the shear band dominated by lc is almost kept unaffected as long as the adopted meshes are fine enough,but it does not generally present a slip surface like in the cases from the classical continuum,indicating that the micro-polar continuum is capable of capturing the non-local geotechnical failure characteristic.Due to the coupling effects of lc and strain softening,softening behavior of geomaterial tends to be postponed.Additionally,the bearing capacity of a geotechnical system may be significantly underestimated,if the effects of lc are not modeled or considered in numerical analyses.
基金Projects(51678083,41302226)supported by the National Natural Science Foundation of China
文摘Uniaxial or triaxial compression test of cylindrical rock specimens using rock mechanics testing machine is a basic experimental method to study the strength and deformation characteristics of rock and the development process of rock fracture. Extensive literature review has been conducted on this issue;meanwhile, experimental and numerical studies have been conducted on the stress-drop effect on the brittleness of rock materials. A plastic flow factor of λ is proposed to describe the stress-drop effect. Evaluation methods of the factor λ corresponding to the four yield criteria of rock mass are proposed. Those four yield criteria are Tresca criterion, von-Mises criterion, Mohr-Coulomb criterion and Drucker-Prager criterion. For simplicity purposes, an engineering approximation approach has been proposed to evaluate the stress-drop with a non-zero strain increment. Numerical simulation results validated the effectiveness of the plastic flow factors λ as well as the engineering approximation approach. Based on the results in this study, finite element code can be programmed for brittle materials with stress-drop, which has the potential to be readily incorporated in finite element codes.
基金Project(51208391) supported by the National Natural Science Foundation of China
文摘The hydration mechanism of low quality fly ash in cement-based materials was investigated. The hydration heat of the composite cementitious materials was determined by isothermal calorimetry, and the hydration products, quantity, pore structure and morphology were measured by X-ray diffraction(XRD), thermalgravity-differential thermal analysis(TG-DTA), mercury intrusion porosimetry(MIP) and scanning electron microscopy(SEM), respectively. The results indicate that grinding could not only improve the physical properties of the low quality fly ash on particle effect, but also improve hydration properties of the cementitious system from various aspects compared with raw low quality fly ash(RLFA). At the early stage of hydration, the low quanlity fly ash acts as almost inert material; but then at the later stage, high chemical activity, especially for ground low quality fly ash(GLFA), could be observed. It can accelerate the formation of hydration products containing more chemical bonded water, resulting in higher degree of cement hydration, thus denser microstructure and more reasonable pore size distribution, but the hydration heat in total is reduced. It can also delay the induction period, but the accelerating period is shortened and there is little influence on the second exothermic peak.
基金The research was supported by the National Natural Science Foundation of China(No.51422906).References
文摘A finite volume,multiphase solver in the framework of OpenFOAM was used to calculate the flow field of the cavitating flow over the Clark-Y hydrofoil. This solver used Transport Based Equation Model(TEM) to solve the liquid volume fraction,and utilized volume of fluid(VOF) technique to predict the interface between liquid and vapor phases. The simulation was designed to study the cavitation shedding and different fluid characteristics in the cloud cavitation regime when adopting two different Large Eddy Simulation(LES) models,namely,one equation eddy viscosity(one EqEddy) model and Smagorinsky model. It is shown that these two models can be used to study the cavitation shedding dynamics and predict the velocity profiles.
基金Project(2010CB732005) supported by the National Basic Research Program of ChinaProjects(51279136, 51209164) supported by the National Natural Science Foundation of China
文摘Regarding excavation-damaged zone （EDZ） around underground opening as non-homogeneous rockmass with spatial deterioration effect on stuffiness and strength, a parametric model of EDZ using radius-displacement-dependent deformation modulus （RDDM） was proposed. Considering the nonlinearity characteristic of deformation and locality otherness of surrounding rock, deterioration parameter field of deformation modulus of rockmass around opening was quantitatively calculated through a given function. Applicability for multi-cavern condition and parameter sensibility of the model was analyzed by numerical experiments using synthetic data. Furthermore, the model was applied to identify EDZ of underground caverns of Pubugou hydropower station by calculating deterioration parameter field. Based on the parametric analysis of spatial effect and geological investigation, it is recognized that large radial deformation of deep fractured rock at the spandrel position and insufficient supporting bolts mainly result in great deformation pressure to act on the shotcrete and cause partial crack and spalling. It is shown that deterioration parameter field along the longitudinal axis of main powerhouse is evidently non-homogeneous in space and distributes exponentially along the radius from the opening. The model provides a simple and convenient way to identify the EDZ in the working state for rapid construction feedback analysis and support optimization of underground cavem from quantitative point of view and also aids in interpreting monitoring displacements and estimating support requirements.
基金Project(50139030) supported by the National Natural Science Foundation of ChinaProject(501072) supported by the Scientific Research Foundation for the Returned Overseas Scholars of the Ministry of Education of China
文摘The seepage under a dam foundation is mainly controlled by the performance of the curtain.Its anti-seepage behavior may be weakened by the long term physic-chemical actions from groundwater.According to seepage hydraulics and geochemistry theory,a coupling model for assessing the behavior of the curtain under a dam foundation is set up,which consists of seepage module,solute transport module,geochemistry module and curtain erosion module,solved by FEM.A case study was carried out.The result shows that the curtain efficiency is weakened all the time.Aqueous calcium from the curtain is always in dissolution during the stress period for simulation,which leads to the increasing amount in groundwater reaching 846.35-865.312 g/m3.Within the domain,reaction extent differs in different parts of the curtain.The dissolution of Ca(OH)2 accounts to 877.884 g/m3 near the bottom and is much higher than that of the other parts.The erosion is much more serious near the bottom of the curtain than the other parts,which is the same to the upstream and downstream.Calcium dissolution is mainly controlled by hydraulic condition and dispersion,and it varies in a non-linear way within the domain.
基金National Natural Science Foundation of China(51422906,51609177)Key Projects in the National Science & Technology Pillar Program During the Twelfth Five-Year Plan Period(2012BAD08B03)
文摘The bubble dynamic near a rigid wall with a wall jet was investigated by codynamics(CFD)method with the volume of fluid(VOF)model,which had been validated by vious experimental data.The effects of different velocities of the wall jet and ditances on the bubble dynamics were studied.The results show that the bubble is squjet due to more force added on the bubble.When the velocity of the wall jet increa,the wall anthe pressure along the wall at collapse time increase because of the extra push indAs the stand-off distance increases,the pressure along the wall first increases then decrethe distance from the bubble to the wall increases.
基金国家自然科学基金项目(41176052、41576035、41276050)南方海洋科学与工程广东省实验室(广州)人才团队引进重大专项(GML2019ZD0104)+2 种基金广东省自然科学基金面上项目(S2013010014523)中国科学院南海生态环境工程创新研究院创新发展基金项目(ISEE2018PY02)the U.S.Department of Energy(DE-AC02-05CH11231)~~
