The excavation of deep tunnels crossing faults is highly prone to triggering rockburst disasters,which has become a significant engineering issue.In this study,taking the fault-slip rockbursts from a deep tunnel in so...The excavation of deep tunnels crossing faults is highly prone to triggering rockburst disasters,which has become a significant engineering issue.In this study,taking the fault-slip rockbursts from a deep tunnel in southwestern China as the engineering prototype,large-scale three-dimensional(3D)physical model tests were conducted on a 3D-printed complex geological model containing two faults.Based on the selfdeveloped 3D loading system and excavation device,the macroscopic failure of fault-slip rockbursts was simulated indoors.The stress,strain,and fracturing characteristics of the surrounding rock near the two faults were systematically evaluated during excavation and multistage loading.The test results effectively revealed the evolution and triggering mechanism of fault-slip rockbursts.After the excavation of a highstress tunnel,stress readjustment occurred.Owing to the presence of these two faults,stress continued to accumulate in the rock mass between them,leading to the accumulation of fractures.When the shear stress on a fault surface exceeded its shear strength,sudden fault slip and dislocation occurred,thus triggering rockbursts.Rockbursts occurred twice in the vault between the two faults,showing obvious intermittent characteristics.The rockburst pit was controlled by two faults.When the faults remained stable,tensile failure predominated in the surrounding rock.However,when the fault slip was triggered,shear failure in the surrounding rock increased.These findings provide valuable insights for enhancing the comprehension of fault-slip rockbursts.展开更多
The workload of the 3D magnetotelluric forward modeling algorithm is so large that the traditional serial algorithm costs an extremely large compute time. However, the 3D forward modeling algorithm can process the dat...The workload of the 3D magnetotelluric forward modeling algorithm is so large that the traditional serial algorithm costs an extremely large compute time. However, the 3D forward modeling algorithm can process the data in the frequency domain, which is very suitable for parallel computation. With the advantage of MPI and based on an analysis of the flow of the 3D magnetotelluric serial forward algorithm, we suggest the idea of parallel computation and apply it. Three theoretical models are tested and the execution efficiency is compared in different situations. The results indicate that the parallel 3D forward modeling computation is correct and the efficiency is greatly improved. This method is suitable for large size geophysical computations.展开更多
This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of com...This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of complex frequency-shifted perfectly matched layer(CFS-PML)was used for truncation so that the low-frequency electromagnetic wave can be better absorbed at the model boundary.A typical three-dimensional(3D)homogeneous half-space model was established and a low-resistivity cube model was analyzed under the half-space condition.The response patterns and drivers of the low-resistivity cube model were discussed under the influence of a low-resistivity overburden.The absorption boundary conditions of CFS-PML significantly affected the low-frequency electromagnetic waves.For a low-resistivity cube around the borehole,its response curve exhibited a single-peak,and the extreme point of the curve corresponded to the center of the low-resistivity body.When the low-resistivity cube was directly below the borehole,the response curve showed three extreme values(two high and one low),with the low corresponding to the center of the low-resistivity body.The total field response of the low-resistivity overburden was stronger than that of the uniform half-space model due to the low-resistivity shielding effect of electromagnetic waves.When the receiving-transmitting distance gradually increased,the effect of the low-resistivity overburden was gradually weakened,and the response of the low-resistivity cube was strengthened.It was affected by the ratio of the overburden resistivity to the resistivity of the low-resistivity body.展开更多
Magnetic field gradient tensor technique provides abundant data for delicate inversion of subsurface magnetic susceptibility distribution. Large scale magnetic data inversion imaging requires high speed and accuracy f...Magnetic field gradient tensor technique provides abundant data for delicate inversion of subsurface magnetic susceptibility distribution. Large scale magnetic data inversion imaging requires high speed and accuracy for forward modeling. For arbitrarily distributed susceptibility data on an undulated surface, we propose a fast 3D forward modeling method in the wavenumber domain based on(1) the wavenumber-domain expression of the prism combination model and the Gauss–FFT algorithm and(2) cubic spline interpolation. We apply the proposed 3D forward modeling method to synthetic data and use weighting coefficients in the wavenumber domain to improve the modeling for multiple observation surfaces, and also demonstrate the accuracy and efficiency of the proposed method.展开更多
Unmanned aerial vehicle transient electromagnetic(UAV-TEM)is a novel airborne exploration method that offers advantages such as low cost,simple operation,high exploration efficiency and suitability for near-surface ex...Unmanned aerial vehicle transient electromagnetic(UAV-TEM)is a novel airborne exploration method that offers advantages such as low cost,simple operation,high exploration efficiency and suitability for near-surface exploration in complex terrain areas.To improve the accuracy of data interpretation in this method,the authors conducted a systematic three-dimensional(3D)forward modeling and inversion of the UAV-TEM.This study utilized the finite element method based on unstructured tetrahedral elements and employed the second-order backward Euler method for time discretization.This allowed for accurate 3D modeling and accounted for the effects of complex terrain.Based on these,the influence characteristics of flight altitudes and the sizes,burial depths,and resistivities of anomalies are compared and analyzed to explore the UAV-TEM systems’exploration capability.Lastly,four typical geoelectrical models of landslides are designed,and the inversion method based on the Gauss-Newton optimization method is used to image the landslide models and analyze the imaging effect of the UAV-TEM method on landslide geohazards.Numerical results showed that UAV-TEM could have better exploration resolution and fine imaging of nearsurface structures,providing important technical support for monitoring,early warning,and preventing landslides and other geological hazards.展开更多
The superparamagnetic effect arises from the superparamagnetism exhibited by a multitude of nano-sized magnetic mineral particles under an external electric field.This phenomenon manifests in transient electromagnetic...The superparamagnetic effect arises from the superparamagnetism exhibited by a multitude of nano-sized magnetic mineral particles under an external electric field.This phenomenon manifests in transient electromagnetic data primarily as a deceleration in the attenuation rate of late-stage signals,a characteristic difficult to discern directly from airborne transient electromagnetic signals,consequently leading to significant misinterpretations of subterranean electrical structures.This study embarks on 3D forward modeling of airborne electromagnetic responses in the frequency domain,accounting for the superparamagnetic effect,utilizing an unstructured finite element method.Superparamagnetic responses in the time domain were obtained through frequency-time conversion.This investigation explores the influence of various parameters-such as magnetic susceptibility,time constants,and flight altitude-on the superparamagnetic effect by examining the response characteristics of typical targets.Findings indicate that in its late stages,the superparamagnetic effect can induce a relative anomaly of up to 300%.There is a positive correlation between magnetic susceptibility and the strength of the superparamagnetic effect.The influence of the time constant's upper and lower limits on the superparamagnetic effect is minimal;however,the range between these limits significantly affects the effect,showing a negative correlation with its intensity.Higher flight altitudes weaken the superparamagnetic signal.The impact is most pronounced when superparamagnetic minerals are shallowly buried,effectively shielding the underlying geology with the characteristics of a good conductivity anomaly,but this effect diminishes with greater depth.The insights from this study provide a theoretical framework for a deeper understanding of the superparamagnetic effect in transient electromagnetic signals and for more accurate interpretations of subterranean geological and electrical structures.展开更多
The traditional ground direct current method is not suitable for leakage detection of underground diaphragm walls in foundation pits because of its low accuracy and poor anti-noise ability.Here,we propose a joint surf...The traditional ground direct current method is not suitable for leakage detection of underground diaphragm walls in foundation pits because of its low accuracy and poor anti-noise ability.Here,we propose a joint surface-borehole observation device for leakage electric fi eld detection to achieve rapid measurement of the electric fi eld distribution characteristics at ground level in the foundation pit,thus enabling rapid localization of leakage points.We first establish the mechanism and basic equation of the leakage electric field response by combining the electric field formed by electrokinetic effect(EK)and the stable electric fi eld formed by conduction current in a combined leakage channel.Then,the fi nite–infi nite element coupling method is used to solve the electric fi eld equation to simulate the responses of a three-dimensional foundation pit leakage model.Furthermore,we conduct numerical simulations of diff erent pit models to investigate the infl uencing factors of the detection device and response characteristics of the change in the properties of the leakage channel.The results demonstrate that the proposed joint surface-borehole observation device can effi ciently reveal anomalous potential caused by leakage,and the amplitude of the electric fi eld generated by EK can eff ectively strengthen the leakage electric fi eld signal at the leakage,thus improving detection accuracy and effi ciency.展开更多
Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a...Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a limited range of wavenumbers, and produces large numerical dispersion beyond this range. The optimal FD scheme based on least squares （LS） can guarantee high precision over a larger range of wavenumbers and obtain the best optimization solution at small computational cost. We extend the LS-based optimal FD scheme from two-dimensional （2D） forward modeling to three-dimensional （3D） and develop a 3D acoustic optimal FD method with high efficiency, wide range of high accuracy and adaptability to parallel computing. Dispersion analysis and forward modeling demonstrate that the developed FD method suppresses numerical dispersion. Finally, we use the developed FD method to source wavefield extrapolation and receiver wavefield extrapolation in 3D RTM. To decrease the computation time and storage requirements, the 3D RTM is implemented by combining the efficient boundary storage with checkpointing strategies on GPU. 3D RTM imaging results suggest that the 3D optimal FD method has higher precision than conventional methods.展开更多
Edge reflections are inevitable in numerical modeling of seismic wavefields, and they are usually attenuated by absorbing boundary conditions. However, the commonly used perfectly matched layer （PML） boundary condit...Edge reflections are inevitable in numerical modeling of seismic wavefields, and they are usually attenuated by absorbing boundary conditions. However, the commonly used perfectly matched layer （PML） boundary condition requires special treatment for the absorbing zone, and in three-dimensional （3D） modeling, it has to split each variable into three corresponding variables, which increases the computing time and memory storage. In contrast, the hybrid absorbing boundary condition （HABC） has the advantages such as ease of implementation, less computation time, and near-perfect absorption; it is thus able to enhance the computational efficiency of 3D elastic wave modeling. In this study, a HABC is developed from two-dimensional （2D） modeling into 3D modeling based on the I st Higdon one way wave equations, and a HABC is proposed that is suitable for a 3D elastic wave numerical simulation. Numerical simulation results for a homogenous model and a complex model indicate that the proposed HABC method is more effective and has better absorption than the traditional PML method.展开更多
Interpretation of geophysical material is the prospecting method. Interpretation of Gravity-megnetic data is based on data processing and inversion. When the grid is divided into several million cells, the computing t...Interpretation of geophysical material is the prospecting method. Interpretation of Gravity-megnetic data is based on data processing and inversion. When the grid is divided into several million cells, the computing task is heavy and time-consuming. In order to increase efficiency of the 3D forward modeling, the paper will adopt MPI parallel algorithm and the several processes will deal with data in the method. Finally, we can gather the result. Through comparing the result of sequence algorithm with the result of MPI parallel algorithm, we can see the result is the same. When the number of processes is 2 to 8, the speed-up ratio is 1.97 to 5. The MPI parallel algorithm is very efficient.展开更多
Acoustic-elastic coupled media is often encountered in most marine explorations, and accurate simulation of acoustic-elastic coupled media is of great significance. At present, the study of acoustic-elastic coupled me...Acoustic-elastic coupled media is often encountered in most marine explorations, and accurate simulation of acoustic-elastic coupled media is of great significance. At present, the study of acoustic-elastic coupled media still assumes that the solid of the acoustic-elastic coupled media is isotropic, but this assumption is not in accordance with the actual situation. In this paper, we derive the solid media of acoustic-elastic coupled media from isotropic media to anisotropic media, and propose an acoustic-elastic coupled medium based ontransverse isotropic media with vertical symmetric axes(VTI) to improve the accuracy of forward modeling. Based on the relationship between the Thomsen parameter and the coefficient matrix of the anisotropic elastic wave equation, we transform the Thomson parameter into a velocity model with anisotropic properties. We use a staggered grid finite difference method to simulate the propagation of a wavefield in a three-dimensional acoustic-elastic coupled media. We obtain the snapshots of the wave field when the solid of the acoustic-elastic coupled media is an isotropic medium and a VTI media. When the solid of the acoustic-elastic coupled media is considered VTI media, we can observe the qP wave and qS wave that cannot be observed in the isotropic medium from the wave field snapshot. We can also find that the seismic records obtained by the method we use are more realistic. The algorithm proposed in this paper is of great significance for high-precision ocean numerical simulation.展开更多
How to identify the nested structure of a three-dimensional(3D)hierarchical groundwater flow system is always a difficult problem puzzling hydrogeologists due to the multiple scales and complexity of the 3D flow field...How to identify the nested structure of a three-dimensional(3D)hierarchical groundwater flow system is always a difficult problem puzzling hydrogeologists due to the multiple scales and complexity of the 3D flow field.The main objective of this study was to develop a quantitative method to partition the nested groundwater flow system into different hierarchies in three dimensions.A 3D numerical model with topography derived from the real geomatic data in Jinan,China was implemented to simulate groundwater flow and residence time at the regional scale while the recharge rate,anisotropic permeability and hydrothermal effect being set as climatic and hydrogeological variables in the simulations.The simulated groundwater residence time distribution showed a favorable consistency with the spatial distribution of flow fields.The probability density function of residence time with discontinuous segments indicated the discrete nature of time domain between different flow hierarchies,and it was used to partition the hierarchical flow system into shallow/intermediate/deep flow compartments.The changes in the groundwater flow system can be quantitatively depicted by the climatic and hydrogeological variables.This study provides new insights and an efficient way to analyze groundwater circulation and evolution in three dimensions from the perspective of time domain.展开更多
A model is presented for the simulation of reactive gas-solids flows in large industrial reactors. Circulating fluidized bed （CFB） combustors with several thousands of cubic meters reaction volume are probably the l...A model is presented for the simulation of reactive gas-solids flows in large industrial reactors. Circulating fluidized bed （CFB） combustors with several thousands of cubic meters reaction volume are probably the largest reactors of this type. A semi-empirical modeling approach has been chosen to model the three-dimensional concentration distributions of gas and solids components and temperatures inside the combustion chamber of such boilers. Two industrial CFB boilers are investigated in detail： the 105 MWe Duisburg combustor in Germany and the 235 MWe Turow combustor in Poland. The semi-empirical model approach is described first. Then the model is used to show how the three-dimensional concentration and temperature fields are formed by the interaction of several local phenomena. Good agreement between simulation and measurements has been achieved.展开更多
Recent advances in generative artiflcial intelligence(AI)technologies have been signiflcantly driven by models such as generative adversarial networks(GANs),variational autoencoders(VAEs),and denoising diffusion proba...Recent advances in generative artiflcial intelligence(AI)technologies have been signiflcantly driven by models such as generative adversarial networks(GANs),variational autoencoders(VAEs),and denoising diffusion probabilistic models(DDPMs).Although architects recognize the potential of generative AI in design,personal barriers often restrict their access to the latest technological developments,thereby causing the application of generative AI in architectural design to lag behind.Therefore,it is essential to comprehend the principles and advancements of generative AI models and analyze their relevance in architecture applications.This paper flrst provides an overview of generative AI technologies,with a focus on probabilistic diffusion models(DDPMs),3D generative models,and foundation models,highlighting their recent developments and main application scenarios.Then,the paper explains how the abovementioned models could be utilized in architecture.We subdivide the architectural design process into six steps and review related research projects in each step from 2020 to the present.Lastly,this paper discusses potential future directions for applying generative AI in the architectural design steps.This research can help architects quickly understand the development and latest progress of generative AI and contribute to the further development of intelligent architecture.展开更多
基金funding support from the National Natural Science Foundation of China(Grant Nos.42177136 and 52309126).
文摘The excavation of deep tunnels crossing faults is highly prone to triggering rockburst disasters,which has become a significant engineering issue.In this study,taking the fault-slip rockbursts from a deep tunnel in southwestern China as the engineering prototype,large-scale three-dimensional(3D)physical model tests were conducted on a 3D-printed complex geological model containing two faults.Based on the selfdeveloped 3D loading system and excavation device,the macroscopic failure of fault-slip rockbursts was simulated indoors.The stress,strain,and fracturing characteristics of the surrounding rock near the two faults were systematically evaluated during excavation and multistage loading.The test results effectively revealed the evolution and triggering mechanism of fault-slip rockbursts.After the excavation of a highstress tunnel,stress readjustment occurred.Owing to the presence of these two faults,stress continued to accumulate in the rock mass between them,leading to the accumulation of fractures.When the shear stress on a fault surface exceeded its shear strength,sudden fault slip and dislocation occurred,thus triggering rockbursts.Rockbursts occurred twice in the vault between the two faults,showing obvious intermittent characteristics.The rockburst pit was controlled by two faults.When the faults remained stable,tensile failure predominated in the surrounding rock.However,when the fault slip was triggered,shear failure in the surrounding rock increased.These findings provide valuable insights for enhancing the comprehension of fault-slip rockbursts.
基金This research is sponsored by the National Natural Science Foundation of China （No. 40374024）.
文摘The workload of the 3D magnetotelluric forward modeling algorithm is so large that the traditional serial algorithm costs an extremely large compute time. However, the 3D forward modeling algorithm can process the data in the frequency domain, which is very suitable for parallel computation. With the advantage of MPI and based on an analysis of the flow of the 3D magnetotelluric serial forward algorithm, we suggest the idea of parallel computation and apply it. Three theoretical models are tested and the execution efficiency is compared in different situations. The results indicate that the parallel 3D forward modeling computation is correct and the efficiency is greatly improved. This method is suitable for large size geophysical computations.
基金This work was supported by China Postdoctoral Science Foundation(No.2022M723391)the Science and Technology Innovation Project of Higher Education in Shanxi Province(No.2019L0754)+1 种基金the Central Guiding Local Science and Technology Development Fund Project(No.YDZJSX2021B021)Shanxi Province Basic Research Plan General Project(No.202203021221294).
文摘This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of complex frequency-shifted perfectly matched layer(CFS-PML)was used for truncation so that the low-frequency electromagnetic wave can be better absorbed at the model boundary.A typical three-dimensional(3D)homogeneous half-space model was established and a low-resistivity cube model was analyzed under the half-space condition.The response patterns and drivers of the low-resistivity cube model were discussed under the influence of a low-resistivity overburden.The absorption boundary conditions of CFS-PML significantly affected the low-frequency electromagnetic waves.For a low-resistivity cube around the borehole,its response curve exhibited a single-peak,and the extreme point of the curve corresponded to the center of the low-resistivity body.When the low-resistivity cube was directly below the borehole,the response curve showed three extreme values(two high and one low),with the low corresponding to the center of the low-resistivity body.The total field response of the low-resistivity overburden was stronger than that of the uniform half-space model due to the low-resistivity shielding effect of electromagnetic waves.When the receiving-transmitting distance gradually increased,the effect of the low-resistivity overburden was gradually weakened,and the response of the low-resistivity cube was strengthened.It was affected by the ratio of the overburden resistivity to the resistivity of the low-resistivity body.
基金supported by the National Special Plan for the 13th Five-Year Plan of China(No.2017YFC0602204-10)Independent Exploration of the Innovation Project for Graduate Students at Central South University(No.2017zzts176)+3 种基金National Natural Science Foundation of China(Nos.41574127,41404106,and 41674075)Postdoctoral Fund Projects of China(No.2017M622608)National Key R&D Program of China(No.2018YFC0603602)Natural Science Youth Fund Project of the Hunan Province,China(No.2018JJ3642)
文摘Magnetic field gradient tensor technique provides abundant data for delicate inversion of subsurface magnetic susceptibility distribution. Large scale magnetic data inversion imaging requires high speed and accuracy for forward modeling. For arbitrarily distributed susceptibility data on an undulated surface, we propose a fast 3D forward modeling method in the wavenumber domain based on(1) the wavenumber-domain expression of the prism combination model and the Gauss–FFT algorithm and(2) cubic spline interpolation. We apply the proposed 3D forward modeling method to synthetic data and use weighting coefficients in the wavenumber domain to improve the modeling for multiple observation surfaces, and also demonstrate the accuracy and efficiency of the proposed method.
基金Supported by Key Research and Development Project of Guangxi Pr ovince(No.AB21196028).
文摘Unmanned aerial vehicle transient electromagnetic(UAV-TEM)is a novel airborne exploration method that offers advantages such as low cost,simple operation,high exploration efficiency and suitability for near-surface exploration in complex terrain areas.To improve the accuracy of data interpretation in this method,the authors conducted a systematic three-dimensional(3D)forward modeling and inversion of the UAV-TEM.This study utilized the finite element method based on unstructured tetrahedral elements and employed the second-order backward Euler method for time discretization.This allowed for accurate 3D modeling and accounted for the effects of complex terrain.Based on these,the influence characteristics of flight altitudes and the sizes,burial depths,and resistivities of anomalies are compared and analyzed to explore the UAV-TEM systems’exploration capability.Lastly,four typical geoelectrical models of landslides are designed,and the inversion method based on the Gauss-Newton optimization method is used to image the landslide models and analyze the imaging effect of the UAV-TEM method on landslide geohazards.Numerical results showed that UAV-TEM could have better exploration resolution and fine imaging of nearsurface structures,providing important technical support for monitoring,early warning,and preventing landslides and other geological hazards.
文摘The superparamagnetic effect arises from the superparamagnetism exhibited by a multitude of nano-sized magnetic mineral particles under an external electric field.This phenomenon manifests in transient electromagnetic data primarily as a deceleration in the attenuation rate of late-stage signals,a characteristic difficult to discern directly from airborne transient electromagnetic signals,consequently leading to significant misinterpretations of subterranean electrical structures.This study embarks on 3D forward modeling of airborne electromagnetic responses in the frequency domain,accounting for the superparamagnetic effect,utilizing an unstructured finite element method.Superparamagnetic responses in the time domain were obtained through frequency-time conversion.This investigation explores the influence of various parameters-such as magnetic susceptibility,time constants,and flight altitude-on the superparamagnetic effect by examining the response characteristics of typical targets.Findings indicate that in its late stages,the superparamagnetic effect can induce a relative anomaly of up to 300%.There is a positive correlation between magnetic susceptibility and the strength of the superparamagnetic effect.The influence of the time constant's upper and lower limits on the superparamagnetic effect is minimal;however,the range between these limits significantly affects the effect,showing a negative correlation with its intensity.Higher flight altitudes weaken the superparamagnetic signal.The impact is most pronounced when superparamagnetic minerals are shallowly buried,effectively shielding the underlying geology with the characteristics of a good conductivity anomaly,but this effect diminishes with greater depth.The insights from this study provide a theoretical framework for a deeper understanding of the superparamagnetic effect in transient electromagnetic signals and for more accurate interpretations of subterranean geological and electrical structures.
基金partially supported by the National Natural Science Foundation of China (Nos. 41864004 and 41674077)Jiangxi Provincial Academic Leaders (Youth) Training Program (No. 20204BCJL23058)Open Fund from Engineering Research Center for Seismic Disaster Prevention and Engineering Geological Disaster Detection of Jiangxi Province (SDGD202102)。
文摘The traditional ground direct current method is not suitable for leakage detection of underground diaphragm walls in foundation pits because of its low accuracy and poor anti-noise ability.Here,we propose a joint surface-borehole observation device for leakage electric fi eld detection to achieve rapid measurement of the electric fi eld distribution characteristics at ground level in the foundation pit,thus enabling rapid localization of leakage points.We first establish the mechanism and basic equation of the leakage electric field response by combining the electric field formed by electrokinetic effect(EK)and the stable electric fi eld formed by conduction current in a combined leakage channel.Then,the fi nite–infi nite element coupling method is used to solve the electric fi eld equation to simulate the responses of a three-dimensional foundation pit leakage model.Furthermore,we conduct numerical simulations of diff erent pit models to investigate the infl uencing factors of the detection device and response characteristics of the change in the properties of the leakage channel.The results demonstrate that the proposed joint surface-borehole observation device can effi ciently reveal anomalous potential caused by leakage,and the amplitude of the electric fi eld generated by EK can eff ectively strengthen the leakage electric fi eld signal at the leakage,thus improving detection accuracy and effi ciency.
基金supported by the National Natural Science Foundation of China(No.41474110)Shell Ph.D. Scholarship to support excellence in geophysical research
文摘Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a limited range of wavenumbers, and produces large numerical dispersion beyond this range. The optimal FD scheme based on least squares （LS） can guarantee high precision over a larger range of wavenumbers and obtain the best optimization solution at small computational cost. We extend the LS-based optimal FD scheme from two-dimensional （2D） forward modeling to three-dimensional （3D） and develop a 3D acoustic optimal FD method with high efficiency, wide range of high accuracy and adaptability to parallel computing. Dispersion analysis and forward modeling demonstrate that the developed FD method suppresses numerical dispersion. Finally, we use the developed FD method to source wavefield extrapolation and receiver wavefield extrapolation in 3D RTM. To decrease the computation time and storage requirements, the 3D RTM is implemented by combining the efficient boundary storage with checkpointing strategies on GPU. 3D RTM imaging results suggest that the 3D optimal FD method has higher precision than conventional methods.
基金supported by the National Natural Science Foundation of China(No.41474110)
文摘Edge reflections are inevitable in numerical modeling of seismic wavefields, and they are usually attenuated by absorbing boundary conditions. However, the commonly used perfectly matched layer （PML） boundary condition requires special treatment for the absorbing zone, and in three-dimensional （3D） modeling, it has to split each variable into three corresponding variables, which increases the computing time and memory storage. In contrast, the hybrid absorbing boundary condition （HABC） has the advantages such as ease of implementation, less computation time, and near-perfect absorption; it is thus able to enhance the computational efficiency of 3D elastic wave modeling. In this study, a HABC is developed from two-dimensional （2D） modeling into 3D modeling based on the I st Higdon one way wave equations, and a HABC is proposed that is suitable for a 3D elastic wave numerical simulation. Numerical simulation results for a homogenous model and a complex model indicate that the proposed HABC method is more effective and has better absorption than the traditional PML method.
文摘Interpretation of geophysical material is the prospecting method. Interpretation of Gravity-megnetic data is based on data processing and inversion. When the grid is divided into several million cells, the computing task is heavy and time-consuming. In order to increase efficiency of the 3D forward modeling, the paper will adopt MPI parallel algorithm and the several processes will deal with data in the method. Finally, we can gather the result. Through comparing the result of sequence algorithm with the result of MPI parallel algorithm, we can see the result is the same. When the number of processes is 2 to 8, the speed-up ratio is 1.97 to 5. The MPI parallel algorithm is very efficient.
基金Supported by Major Project of National Science and Technology of China(No.2016ZX05026-002-003)National Natural Science Foundation of China(No.41374108)
文摘Acoustic-elastic coupled media is often encountered in most marine explorations, and accurate simulation of acoustic-elastic coupled media is of great significance. At present, the study of acoustic-elastic coupled media still assumes that the solid of the acoustic-elastic coupled media is isotropic, but this assumption is not in accordance with the actual situation. In this paper, we derive the solid media of acoustic-elastic coupled media from isotropic media to anisotropic media, and propose an acoustic-elastic coupled medium based ontransverse isotropic media with vertical symmetric axes(VTI) to improve the accuracy of forward modeling. Based on the relationship between the Thomsen parameter and the coefficient matrix of the anisotropic elastic wave equation, we transform the Thomson parameter into a velocity model with anisotropic properties. We use a staggered grid finite difference method to simulate the propagation of a wavefield in a three-dimensional acoustic-elastic coupled media. We obtain the snapshots of the wave field when the solid of the acoustic-elastic coupled media is an isotropic medium and a VTI media. When the solid of the acoustic-elastic coupled media is considered VTI media, we can observe the qP wave and qS wave that cannot be observed in the isotropic medium from the wave field snapshot. We can also find that the seismic records obtained by the method we use are more realistic. The algorithm proposed in this paper is of great significance for high-precision ocean numerical simulation.
基金supported by the National Natural Science Foundation of China(Nos.41807219,41877192,U1906209,42072331)the National Key R&D Program of China(No.2017YFC0505304)the Fundamental Research Funds for Central Public Welfare Research Institutes(Nos.CKSF 2019170/TB,CKSF 2016029/TB)。
文摘How to identify the nested structure of a three-dimensional(3D)hierarchical groundwater flow system is always a difficult problem puzzling hydrogeologists due to the multiple scales and complexity of the 3D flow field.The main objective of this study was to develop a quantitative method to partition the nested groundwater flow system into different hierarchies in three dimensions.A 3D numerical model with topography derived from the real geomatic data in Jinan,China was implemented to simulate groundwater flow and residence time at the regional scale while the recharge rate,anisotropic permeability and hydrothermal effect being set as climatic and hydrogeological variables in the simulations.The simulated groundwater residence time distribution showed a favorable consistency with the spatial distribution of flow fields.The probability density function of residence time with discontinuous segments indicated the discrete nature of time domain between different flow hierarchies,and it was used to partition the hierarchical flow system into shallow/intermediate/deep flow compartments.The changes in the groundwater flow system can be quantitatively depicted by the climatic and hydrogeological variables.This study provides new insights and an efficient way to analyze groundwater circulation and evolution in three dimensions from the perspective of time domain.
文摘A model is presented for the simulation of reactive gas-solids flows in large industrial reactors. Circulating fluidized bed （CFB） combustors with several thousands of cubic meters reaction volume are probably the largest reactors of this type. A semi-empirical modeling approach has been chosen to model the three-dimensional concentration distributions of gas and solids components and temperatures inside the combustion chamber of such boilers. Two industrial CFB boilers are investigated in detail： the 105 MWe Duisburg combustor in Germany and the 235 MWe Turow combustor in Poland. The semi-empirical model approach is described first. Then the model is used to show how the three-dimensional concentration and temperature fields are formed by the interaction of several local phenomena. Good agreement between simulation and measurements has been achieved.
基金supported by the Innovative Research Group Project of the National Natural Science Foundation of China(Grant No.202401-202712)。
文摘Recent advances in generative artiflcial intelligence(AI)technologies have been signiflcantly driven by models such as generative adversarial networks(GANs),variational autoencoders(VAEs),and denoising diffusion probabilistic models(DDPMs).Although architects recognize the potential of generative AI in design,personal barriers often restrict their access to the latest technological developments,thereby causing the application of generative AI in architectural design to lag behind.Therefore,it is essential to comprehend the principles and advancements of generative AI models and analyze their relevance in architecture applications.This paper flrst provides an overview of generative AI technologies,with a focus on probabilistic diffusion models(DDPMs),3D generative models,and foundation models,highlighting their recent developments and main application scenarios.Then,the paper explains how the abovementioned models could be utilized in architecture.We subdivide the architectural design process into six steps and review related research projects in each step from 2020 to the present.Lastly,this paper discusses potential future directions for applying generative AI in the architectural design steps.This research can help architects quickly understand the development and latest progress of generative AI and contribute to the further development of intelligent architecture.
