To comprehensively utilize the valuable geological map,exploration profile,borehole,and geochemical logging data and the knowledge on the formation of the Jinshan Ag-Au deposit for forecasting the exploration targets ...To comprehensively utilize the valuable geological map,exploration profile,borehole,and geochemical logging data and the knowledge on the formation of the Jinshan Ag-Au deposit for forecasting the exploration targets of concealed ore bodies,three-dimensional Mineral Prospectivity Modeling(MPM)of the deposit has been conducted using the weights-of-evidence(WofE)method.Conditional independence between evidence layers was tested,and the outline results using the prediction-volume(P-V)and Student's t-statistic methods for delineating favorable mineralization areas from continuous posterior probability map were critically compared.Four exploration targets delineated ultimately by the Student's t-statistic method for the discovery of minable ore bodies in each of the target areas were discussed in detail.The main conclusions include:(1)three-dimensional modeling of a deposit using multi-source reconnaissance data is useful for MPM in interpreting their relationships with known ore bodies;(2)WofE modeling can be used as a straightforward tool for integrating deposit model and reconnaissance data in MPM;(3)the Student's t-statistic method is more applicable in binarizing the continuous prospectivity map for exploration targeting than the PV approach;and(4)two target areas within high potential to find undiscovered ore bodies were diagnosed to guide future near-mine exploration activities of the Jinshan deposit.展开更多
Satellite Component Layout Optimization(SCLO) is crucial in satellite system design.This paper proposes a novel Satellite Three-Dimensional Component Assignment and Layout Optimization(3D-SCALO) problem tailored to en...Satellite Component Layout Optimization(SCLO) is crucial in satellite system design.This paper proposes a novel Satellite Three-Dimensional Component Assignment and Layout Optimization(3D-SCALO) problem tailored to engineering requirements, aiming to optimize satellite heat dissipation while considering constraints on static stability, 3D geometric relationships between components, and special component positions. The 3D-SCALO problem is a challenging bilevel combinatorial optimization task, involving the optimization of discrete component assignment variables in the outer layer and continuous component position variables in the inner layer,with both influencing each other. To address this issue, first, a Mixed Integer Programming(MIP) model is proposed, which reformulates the original bilevel problem into a single-level optimization problem, enabling the exploration of a more comprehensive optimization space while avoiding iterative nested optimization. Then, to model the 3D geometric relationships between components within the MIP framework, a linearized 3D Phi-function method is proposed, which handles non-overlapping and safety distance constraints between cuboid components in an explicit and effective way. Subsequently, the Finite-Rectangle Method(FRM) is proposed to manage 3D geometric constraints for complex-shaped components by approximating them with a finite set of cuboids, extending the applicability of the geometric modeling approach. Finally, the feasibility and effectiveness of the proposed MIP model are demonstrated through two numerical examples"and a real-world engineering case, which confirms its suitability for complex-shaped components and real engineering applications.展开更多
3-D geological modeling plays an increasingly important role in Petroleum Geology, Mining Geology and Engineering Geology. The complexity of geological conditions requires different modeling methods in different situa...3-D geological modeling plays an increasingly important role in Petroleum Geology, Mining Geology and Engineering Geology. The complexity of geological conditions requires different modeling methods in different situations. This paper summarizes the general concept of geological modeling; compares the characteristics of borehole-based modeling, cross-section based modeling and multi- source interactive modeling; analyses key techniques in 3-D geological modeling; and highlights the main difficulties and directions of future studies.展开更多
The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete elem...The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete element method.A randomly generating algorithm was proposed to capture the three-dimensional irregular shape of coarse aggregate.And then,modeling algorithm and method for graded aggregates were built.Based on the combination of modeling of coarse aggregates,asphalt mastic and air voids,three-dimensional virtual sample of asphalt mixture was modeled by using PFC^(3D).Virtual tests for penetration test of aggregate and uniaxial creep test of asphalt mixture were built and conducted by using PFC^(3D).By comparison of the testing results between virtual tests and actual laboratory tests,the validity of the microstructure modeling and virtual test built in this study was verified.Additionally,compared with laboratory test,the virtual test is easier to conduct and has less variability.It is proved that microstructure modeling and virtual test based on three-dimensional discrete element method is a promising way to conduct research of asphalt mixture.展开更多
A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been deve...A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.展开更多
A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forwar...A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.展开更多
The Horne deposit with rich Cu and Au in Noranda region of Black River Group in Quebec has high economic significance.Current researches on Horne deposit are mostly based on two-dimensional maps and statistical data.I...The Horne deposit with rich Cu and Au in Noranda region of Black River Group in Quebec has high economic significance.Current researches on Horne deposit are mostly based on two-dimensional maps and statistical data.It is hard to reflect the spatial structure and characteristics of Horne orebody directly.In this paper,GIS was used to digitize the mining plan-view maps at different depths,stope maps,the boundary of the massive sulfide in drilling trajectories as well as the grade data of Au and Cu of Horne deposit.Meanwhile,the authors established the grade attribute database.Subsequently the three-dimensional(3D)geological model and grade attribute model of Horne orebody were established by Geological Object Computer Aided Design(GOCAD).Positions of two vents and directions of hydrothermal alteration in Horne deposit were inferred based on the property of the major fault,characteristics of hydrothermal alteration,the enrichment morphology and spatial distribution of high-grade Cu in the Cu attribute model.展开更多
With drilling and seismic data of Transtensional(strike-slip)Fault System in the Ziyang area of the central Sichuan Basin,SW China plane-section integrated structural interpretation,3-D fault framework model building,...With drilling and seismic data of Transtensional(strike-slip)Fault System in the Ziyang area of the central Sichuan Basin,SW China plane-section integrated structural interpretation,3-D fault framework model building,fault throw analyzing,and balanced profile restoration,it is pointed out that the transtensional fault system in the Ziyang 3-D seismic survey consists of the northeast-trending F_(I)19 and F_(I)20 fault zones dominated by extensional deformation,as well as 3 sets of northwest-trending en echelon normal faults experienced dextral shear deformation.Among them,the F_(I)19 and F_(I)20 fault zones cut through the Neoproterozoic to Lower Triassic Jialingjiang Formation,presenting a 3-D structure of an“S”-shaped ribbon.And before Permian and during the Early Triassic,the F_(I)19 and F_(I)20 fault zones underwent at least two periods of structural superimposition.Besides,the 3 sets of northwest-trending en echelon normal faults are composed of small normal faults arranged in pairs,with opposite dip directions and partially left-stepped arrangement.And before Permian,they had formed almost,restricting the eastward growth and propagation of the F_(I)19 fault zone.The F_(I)19 and F_(I)20 fault zones communicate multiple sets of source rocks and reservoirs from deep to shallow,and the timing of fault activity matches well with oil and gas generation peaks.If there were favorable Cambrian-Triassic sedimentary facies and reservoirs developing on the local anticlinal belts of both sides of the F_(I)19 and F_(I)20 fault zones,the major reservoirs in this area are expected to achieve breakthroughs in oil and gas exploration.展开更多
A new back-analysis method of ground stress is proposed with comprehensive consideration of influence of topography, geology and nonlinear physical mechanical properties of rock on ground stress. This method based on ...A new back-analysis method of ground stress is proposed with comprehensive consideration of influence of topography, geology and nonlinear physical mechanical properties of rock on ground stress. This method based on non-uniform rational B-spline (NURBS) technology provides the means to build a refined three-dimensional finite element model with more accurate meshing under complex terrain and geological conditions. Meanwhile, this method is a back-analysis of ground stress with combination of multivariable linear regression model and neural network (ANN) model. Firstly, the regression model is used to fit approximately boundary loads. Regarding the regressed loads as mean value, some sets of boundary loads with the same interval are constructed according to the principle of orthogonal design, to calculate the corresponding ground stress at the observation positions using finite element method. The results (boundary loads and the corresponding ground stress) are added to the samples for ANN training. And on this basis, an ANN model is established to implement higher precise back-analysis of initial ground stress. A practical application case shows that the relative error between the inversed ground stress and observed value is mostly less than 10 %, which can meet the need of engineering design and construction requirements.展开更多
To dynamically update the shape of orebody according to the knowledge of a structural geologist’s insight,an approach of orebody implicit modeling from raw drillhole data using the generalized radial basis function i...To dynamically update the shape of orebody according to the knowledge of a structural geologist’s insight,an approach of orebody implicit modeling from raw drillhole data using the generalized radial basis function interpolant was presented.A variety of constraint rules,including geology trend line,geology constraint line,geology trend surface,geology constraint surface and anisotropy,which can be converted into interpolation constraints,were developed to dynamically control the geology trends.Combined with the interactive tools of constraint rules,this method can avoid the shortcomings of the explicit modeling method based on the contour stitching,such as poor model quality,and is difficult to update dynamically,and simplify the modeling process of orebody.The results of numerical experiments show that the 3D ore body model can be reconstructed quickly,accurately and dynamically by the implicit modeling method.展开更多
Three-dimensional geological modeling （3DGM） assists geologists to quantitatively study in three-dimensional （3D） space structures that define temporal and spatial relationships between geological objects. The 3D ...Three-dimensional geological modeling （3DGM） assists geologists to quantitatively study in three-dimensional （3D） space structures that define temporal and spatial relationships between geological objects. The 3D property model can also be used to infer or deduce causes of geological objects. 3DGM technology provides technical support for extraction of diverse geoscience information, 3D modeling, and quantitative calculation of mineral resources. Based on metallogenic concepts and an ore deposit model, 3DGM technology is applied to analyze geological characteristics of the Tongshan Cu deposit in order to define a metallogenic model and develop a virtual borehole technology; a BP neural network and a 3D interpolation technique were combined to integrate multiple geoscience information in a 3D environment. The results indicate： （1） on basis of the concept of magmatic-hydrothermal Cu polymetallic mineraliza- tion and a porphyry Cu deposit model, a spatial relational database of multiple geoscience information for mineralization in the study area （geology, geophysics, geochemistry, borehole, and cross-section data） was established, and 3D metallogenic geological objects including mineralization stratum, granodiorite, alteration rock, and magnetic anomaly were constructed; （2） on basis of the 3D ore deposit model, 23,800 effective surveys from 94 boreholes and 21 sections were applied to establish 3D orebody models with a kriging interpolation method; （3） combined 23,800 surveys involving 21 sections, using VC＋＋ and OpenGL platform, virtual borehole and virtual section with BP network, and an improved inverse distance interpolation （IDW） method were used to predict and delineate mineralization potential targets （Cu-grade of cell not less than 0.1%）; （4） comparison of 3D ore bodies, metallogenic geological objects of mineralization, and potential targets of mineralization models in the study area, delineated the 3D spatial and temporal relationship and causal processes among the ore bodies, alteration rock, metallo- genic stratum, intrusive rock, and the Tongshan Fault. This study provides important technical support and a scientific basis for assessment of the Tongshan Cu deposit and surrounding exploration and mineral resources.展开更多
The development of digital twins for geotechnical structures necessitates the real-time updates of threedimensional(3D)virtual models(e.g.numerical finite element method(FEM)model)to accurately predict time-varying ge...The development of digital twins for geotechnical structures necessitates the real-time updates of threedimensional(3D)virtual models(e.g.numerical finite element method(FEM)model)to accurately predict time-varying geotechnical responses(e.g.consolidation settlement)in a 3D spatial domain.However,traditional 3D numerical model updating approaches are computationally prohibitive and therefore difficult to update the 3D responses in real time.To address these challenges,this study proposes a novel machine learning framework called sparse dictionary learning(T-3D-SDL)for real-time updating of time-varying 3D geotechnical responses.In T-3D-SDL,a concerned dataset(e.g.time-varying 3D settlement)is approximated as a linear superposition of dictionary atoms generated from 3D random FEM analyses.Field monitoring data are then used to identify non-trivial atoms and estimate their weights within a Bayesian framework for model updating and prediction.The proposed approach enables the real-time update of temporally varying settlements with a high 3D spatial resolution and quantified uncertainty as field monitoring data evolve.The proposed approach is illustrated using an embankment construction project.The results show that the proposed approach effectively improves settlement predictions along temporal and 3D spatial dimensions,with minimal latency(e.g.within minutes),as monitoring data appear.In addition,the proposed approach requires only a reasonably small number of 3D FEM model evaluations,avoids the use of widely adopted yet often criticized surrogate models,and effectively addresses the limitations(e.g.computational inefficiency)of existing 3D model updating approaches.展开更多
Taking hundreds of pieces of hazardous geological maps （1 ： 10 000） of Three Gorges res-ervoir area （3GR） as background, we establish regional three-dimensional （3D） geo-hazard modelusing DEM （digital elevatio...Taking hundreds of pieces of hazardous geological maps （1 ： 10 000） of Three Gorges res-ervoir area （3GR） as background, we establish regional three-dimensional （3D） geo-hazard modelusing DEM （digital elevation model） superposed surface images and geo-hazards elements. Based on landslides and other geo-hazard survey data,using improved B-REP（boundary representa-tion）entity data structure （two-body 3D data structure）, we set up 3D solid models for each hazardous bodies in each hazardous geological maps. Then we integrate the two types of 3D models with different scales from area to point, which are the regional geo-hazard 3D model and the solid models of each disaster body, in order to provide a visual processing and analysis plat-form for danger partition, stability evaluation, disaster prevention and control, early warning and command.展开更多
This paper focuses on the progress in geomechanical modeling associated with carbon dioxide(CO2)geological storage.The detailed review of some geomechanical aspects,including numerical methods,stress analysis,ground d...This paper focuses on the progress in geomechanical modeling associated with carbon dioxide(CO2)geological storage.The detailed review of some geomechanical aspects,including numerical methods,stress analysis,ground deformation,fault reactivation,induced seismicity and crack propagation,is presented.It is indicated that although all the processes involved are not fully understood,integration of all available data,such as ground survey,geological conditions,microseismicity and ground level deformation,has led to many new insights into the rock mechanical response to CO2injection.The review also shows that in geomechanical modeling,continuum modeling methods are predominant compared with discontinuum methods.It is recommended to develop continuum-discontinuum numerical methods since they are more convenient for geomechanical modeling of CO2geological storage,especially for fracture propagation simulation.The Mohr-Coulomb criterion is widely used in prediction of rock mass mechanical behavior.It would be better to use a criterion considering the effect of the intermediate principal stress on rock mechanical behavior,especially for the stability analysis of deeply seated rock engineering.Some challenges related to geomechanical modeling of CO2geological storage are also discussed.展开更多
Mine safety have top-five disasters,which including the water,gas,fire,dust and geological dynamic disaster.The coal mine water disaster is one of the important factors which restricted the development of China’s coa...Mine safety have top-five disasters,which including the water,gas,fire,dust and geological dynamic disaster.The coal mine water disaster is one of the important factors which restricted the development of China’s coal production.It is showed by statistics that 60%of mine accidents are affected by groundwater,which not only result in the production losses,casualties and a variety of展开更多
An orthogonal 2D training image is constructed from the geological analysis results of well logs and sedimentary facies;the 2 D probabilities in three directions are obtained through linear pooling method and then agg...An orthogonal 2D training image is constructed from the geological analysis results of well logs and sedimentary facies;the 2 D probabilities in three directions are obtained through linear pooling method and then aggregated by the logarithmic linear pooling to determine the 3 D multi-point pattern probabilities at the unknown points,to realize the reconstruction of a 3 D model from 2D cross-section.To solve the problems of reducing pattern variability in the 2 D training image and increasing sampling uncertainty,an adaptive spatial sampling method is introduced,and an iterative simulation strategy is adopted,in which sample points from the region with higher reliability of the previous simulation results are extracted to be additional condition points in the following simulation to improve the pattern probability sampling stability.The comparison of lateral accretion layer conceptual models shows that the reconstructing algorithm using self-adaptive spatial sampling can improve the accuracy of pattern sampling and rationality of spatial structure characteristics,and accurately reflect the morphology and distribution pattern of the lateral accretion layer.Application of the method in reconstructing the meandering river reservoir of the Cretaceous McMurray Formation in Canada shows that the new method can accurately reproduce the shape,spatial distribution pattern and development features of complex lateral accretion layers in the meandering river reservoir under tide effect.The test by sparse wells shows that the simulation accuracy is above 85%,and the coincidence rate of interpretation and prediction results of newly drilled horizontal wells is up to 80%.展开更多
The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling ...The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.展开更多
The essential for microseismic monitoring is fast and accurate calculation of seismic wave source location. The precision of most traditional microseismic monitoring processes of mines, using TDOA location method in t...The essential for microseismic monitoring is fast and accurate calculation of seismic wave source location. The precision of most traditional microseismic monitoring processes of mines, using TDOA location method in two-dimensional space to position the microseismic events, as well as the accuracy of positioning microseismic events, may be reduced by the two-dimensional model and simple method, and ill-conditioned equations produced by TDOA location method will increase the positioning error. This article, based on inversion theory, studies the mathematical model of TDOA location method, polariza- tion analysis location method, and comprehensive difference location method of adding angle factor in the traditional TDOA location method. The feasibility of three methods is verified by numerical simulation and analysis of the positioning error of them. The results show that the comprehensive location method of adding angle difference has strong positioning stability and high positioning accuracy, and it may reduce the impact effectively about ill-conditioned equations to positioning results. Comprehensive location method with the data of actual measure may get better positioning results.展开更多
The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the tim...The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.展开更多
In the present study, considering the transport and transformation processes of variables, a threedimensional water quality model for the river system was established, which coupled the volume of fluid(VOF) method wit...In the present study, considering the transport and transformation processes of variables, a threedimensional water quality model for the river system was established, which coupled the volume of fluid(VOF) method with the k-ε turbulence mathematical model. Then, the water hydrodynamic characteristics and transport processes for BOD_5, NH_(3^-)N and TP were analyzed. The results showed that the water surface of convex bank was a little lower than that of concave bank due to the centrifugal force near the bend, and most concentrations were inferior to the type Ⅴ standard indexes of surface water environmental quality. The model validation indicated that the errors between the simulated and monitored values were comparatively small, satisfying the application demands and providing scientific basis and decision support for the restoration and protection of water quality.展开更多
基金financially supported by the Ministry of Science and Technology of China(Nos.2022YFF0801201,2021YFC2900300)the National Natural Science Foundation of China(Nos.41872245,U1911202)the Guangdong Basic and Applied Basic Research Foundation(No.2020A1515010666)。
文摘To comprehensively utilize the valuable geological map,exploration profile,borehole,and geochemical logging data and the knowledge on the formation of the Jinshan Ag-Au deposit for forecasting the exploration targets of concealed ore bodies,three-dimensional Mineral Prospectivity Modeling(MPM)of the deposit has been conducted using the weights-of-evidence(WofE)method.Conditional independence between evidence layers was tested,and the outline results using the prediction-volume(P-V)and Student's t-statistic methods for delineating favorable mineralization areas from continuous posterior probability map were critically compared.Four exploration targets delineated ultimately by the Student's t-statistic method for the discovery of minable ore bodies in each of the target areas were discussed in detail.The main conclusions include:(1)three-dimensional modeling of a deposit using multi-source reconnaissance data is useful for MPM in interpreting their relationships with known ore bodies;(2)WofE modeling can be used as a straightforward tool for integrating deposit model and reconnaissance data in MPM;(3)the Student's t-statistic method is more applicable in binarizing the continuous prospectivity map for exploration targeting than the PV approach;and(4)two target areas within high potential to find undiscovered ore bodies were diagnosed to guide future near-mine exploration activities of the Jinshan deposit.
基金supported by the National Natural Science Foundation of China(No.92371206)the Postgraduate Scientific Research Innovation Project of Hunan Province,China(No.CX2023063).
文摘Satellite Component Layout Optimization(SCLO) is crucial in satellite system design.This paper proposes a novel Satellite Three-Dimensional Component Assignment and Layout Optimization(3D-SCALO) problem tailored to engineering requirements, aiming to optimize satellite heat dissipation while considering constraints on static stability, 3D geometric relationships between components, and special component positions. The 3D-SCALO problem is a challenging bilevel combinatorial optimization task, involving the optimization of discrete component assignment variables in the outer layer and continuous component position variables in the inner layer,with both influencing each other. To address this issue, first, a Mixed Integer Programming(MIP) model is proposed, which reformulates the original bilevel problem into a single-level optimization problem, enabling the exploration of a more comprehensive optimization space while avoiding iterative nested optimization. Then, to model the 3D geometric relationships between components within the MIP framework, a linearized 3D Phi-function method is proposed, which handles non-overlapping and safety distance constraints between cuboid components in an explicit and effective way. Subsequently, the Finite-Rectangle Method(FRM) is proposed to manage 3D geometric constraints for complex-shaped components by approximating them with a finite set of cuboids, extending the applicability of the geometric modeling approach. Finally, the feasibility and effectiveness of the proposed MIP model are demonstrated through two numerical examples"and a real-world engineering case, which confirms its suitability for complex-shaped components and real engineering applications.
文摘3-D geological modeling plays an increasingly important role in Petroleum Geology, Mining Geology and Engineering Geology. The complexity of geological conditions requires different modeling methods in different situations. This paper summarizes the general concept of geological modeling; compares the characteristics of borehole-based modeling, cross-section based modeling and multi- source interactive modeling; analyses key techniques in 3-D geological modeling; and highlights the main difficulties and directions of future studies.
基金Project(51378006) supported by National Natural Science Foundation of ChinaProject(141076) supported by Huoyingdong Foundation of the Ministry of Education of China+1 种基金Project(2242015R30027) supported by Excellent Young Teacher Program of Southeast University,ChinaProject(BK20140109) supported by the Natural Science Foundation of Jiangsu Province,China
文摘The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete element method.A randomly generating algorithm was proposed to capture the three-dimensional irregular shape of coarse aggregate.And then,modeling algorithm and method for graded aggregates were built.Based on the combination of modeling of coarse aggregates,asphalt mastic and air voids,three-dimensional virtual sample of asphalt mixture was modeled by using PFC^(3D).Virtual tests for penetration test of aggregate and uniaxial creep test of asphalt mixture were built and conducted by using PFC^(3D).By comparison of the testing results between virtual tests and actual laboratory tests,the validity of the microstructure modeling and virtual test built in this study was verified.Additionally,compared with laboratory test,the virtual test is easier to conduct and has less variability.It is proved that microstructure modeling and virtual test based on three-dimensional discrete element method is a promising way to conduct research of asphalt mixture.
基金Projects(41674080,41674079)supported by the National Natural Science Foundation of China
文摘A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.
基金Project(60672042) supported by the National Natural Science Foundation of China
文摘A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.
文摘The Horne deposit with rich Cu and Au in Noranda region of Black River Group in Quebec has high economic significance.Current researches on Horne deposit are mostly based on two-dimensional maps and statistical data.It is hard to reflect the spatial structure and characteristics of Horne orebody directly.In this paper,GIS was used to digitize the mining plan-view maps at different depths,stope maps,the boundary of the massive sulfide in drilling trajectories as well as the grade data of Au and Cu of Horne deposit.Meanwhile,the authors established the grade attribute database.Subsequently the three-dimensional(3D)geological model and grade attribute model of Horne orebody were established by Geological Object Computer Aided Design(GOCAD).Positions of two vents and directions of hydrothermal alteration in Horne deposit were inferred based on the property of the major fault,characteristics of hydrothermal alteration,the enrichment morphology and spatial distribution of high-grade Cu in the Cu attribute model.
基金Supported by the Key Project of National Natural Science Foundation of China(42330810).
文摘With drilling and seismic data of Transtensional(strike-slip)Fault System in the Ziyang area of the central Sichuan Basin,SW China plane-section integrated structural interpretation,3-D fault framework model building,fault throw analyzing,and balanced profile restoration,it is pointed out that the transtensional fault system in the Ziyang 3-D seismic survey consists of the northeast-trending F_(I)19 and F_(I)20 fault zones dominated by extensional deformation,as well as 3 sets of northwest-trending en echelon normal faults experienced dextral shear deformation.Among them,the F_(I)19 and F_(I)20 fault zones cut through the Neoproterozoic to Lower Triassic Jialingjiang Formation,presenting a 3-D structure of an“S”-shaped ribbon.And before Permian and during the Early Triassic,the F_(I)19 and F_(I)20 fault zones underwent at least two periods of structural superimposition.Besides,the 3 sets of northwest-trending en echelon normal faults are composed of small normal faults arranged in pairs,with opposite dip directions and partially left-stepped arrangement.And before Permian,they had formed almost,restricting the eastward growth and propagation of the F_(I)19 fault zone.The F_(I)19 and F_(I)20 fault zones communicate multiple sets of source rocks and reservoirs from deep to shallow,and the timing of fault activity matches well with oil and gas generation peaks.If there were favorable Cambrian-Triassic sedimentary facies and reservoirs developing on the local anticlinal belts of both sides of the F_(I)19 and F_(I)20 fault zones,the major reservoirs in this area are expected to achieve breakthroughs in oil and gas exploration.
基金Innovative Research Groups of the National Natural Science Foundation of China (No.51021004)National Science Foundation of China (No. 51079096)Program for New Century Excellent Talents in University (No. NCET-08-0391)
文摘A new back-analysis method of ground stress is proposed with comprehensive consideration of influence of topography, geology and nonlinear physical mechanical properties of rock on ground stress. This method based on non-uniform rational B-spline (NURBS) technology provides the means to build a refined three-dimensional finite element model with more accurate meshing under complex terrain and geological conditions. Meanwhile, this method is a back-analysis of ground stress with combination of multivariable linear regression model and neural network (ANN) model. Firstly, the regression model is used to fit approximately boundary loads. Regarding the regressed loads as mean value, some sets of boundary loads with the same interval are constructed according to the principle of orthogonal design, to calculate the corresponding ground stress at the observation positions using finite element method. The results (boundary loads and the corresponding ground stress) are added to the samples for ANN training. And on this basis, an ANN model is established to implement higher precise back-analysis of initial ground stress. A practical application case shows that the relative error between the inversed ground stress and observed value is mostly less than 10 %, which can meet the need of engineering design and construction requirements.
文摘To dynamically update the shape of orebody according to the knowledge of a structural geologist’s insight,an approach of orebody implicit modeling from raw drillhole data using the generalized radial basis function interpolant was presented.A variety of constraint rules,including geology trend line,geology constraint line,geology trend surface,geology constraint surface and anisotropy,which can be converted into interpolation constraints,were developed to dynamically control the geology trends.Combined with the interactive tools of constraint rules,this method can avoid the shortcomings of the explicit modeling method based on the contour stitching,such as poor model quality,and is difficult to update dynamically,and simplify the modeling process of orebody.The results of numerical experiments show that the 3D ore body model can be reconstructed quickly,accurately and dynamically by the implicit modeling method.
基金supported by the National Basic Research Program of China(Grant No.1212010881001 )the National Scicnce of the 12th "Five-Year Technology Support Program"(Grant No.2010BAE00281-6)+1 种基金the National Natural Science Foundation of China(Grant Nos.40772157,40972232, 41072070)the State Key Laboratory of Geological Processes and Mineral Resources(Grant Nos.GPMR0941,200624)
文摘Three-dimensional geological modeling （3DGM） assists geologists to quantitatively study in three-dimensional （3D） space structures that define temporal and spatial relationships between geological objects. The 3D property model can also be used to infer or deduce causes of geological objects. 3DGM technology provides technical support for extraction of diverse geoscience information, 3D modeling, and quantitative calculation of mineral resources. Based on metallogenic concepts and an ore deposit model, 3DGM technology is applied to analyze geological characteristics of the Tongshan Cu deposit in order to define a metallogenic model and develop a virtual borehole technology; a BP neural network and a 3D interpolation technique were combined to integrate multiple geoscience information in a 3D environment. The results indicate： （1） on basis of the concept of magmatic-hydrothermal Cu polymetallic mineraliza- tion and a porphyry Cu deposit model, a spatial relational database of multiple geoscience information for mineralization in the study area （geology, geophysics, geochemistry, borehole, and cross-section data） was established, and 3D metallogenic geological objects including mineralization stratum, granodiorite, alteration rock, and magnetic anomaly were constructed; （2） on basis of the 3D ore deposit model, 23,800 effective surveys from 94 boreholes and 21 sections were applied to establish 3D orebody models with a kriging interpolation method; （3） combined 23,800 surveys involving 21 sections, using VC＋＋ and OpenGL platform, virtual borehole and virtual section with BP network, and an improved inverse distance interpolation （IDW） method were used to predict and delineate mineralization potential targets （Cu-grade of cell not less than 0.1%）; （4） comparison of 3D ore bodies, metallogenic geological objects of mineralization, and potential targets of mineralization models in the study area, delineated the 3D spatial and temporal relationship and causal processes among the ore bodies, alteration rock, metallo- genic stratum, intrusive rock, and the Tongshan Fault. This study provides important technical support and a scientific basis for assessment of the Tongshan Cu deposit and surrounding exploration and mineral resources.
基金supported by a grant from the Research Grant Council of Hong Kong Special Administrative Region(Project No.11207724).
文摘The development of digital twins for geotechnical structures necessitates the real-time updates of threedimensional(3D)virtual models(e.g.numerical finite element method(FEM)model)to accurately predict time-varying geotechnical responses(e.g.consolidation settlement)in a 3D spatial domain.However,traditional 3D numerical model updating approaches are computationally prohibitive and therefore difficult to update the 3D responses in real time.To address these challenges,this study proposes a novel machine learning framework called sparse dictionary learning(T-3D-SDL)for real-time updating of time-varying 3D geotechnical responses.In T-3D-SDL,a concerned dataset(e.g.time-varying 3D settlement)is approximated as a linear superposition of dictionary atoms generated from 3D random FEM analyses.Field monitoring data are then used to identify non-trivial atoms and estimate their weights within a Bayesian framework for model updating and prediction.The proposed approach enables the real-time update of temporally varying settlements with a high 3D spatial resolution and quantified uncertainty as field monitoring data evolve.The proposed approach is illustrated using an embankment construction project.The results show that the proposed approach effectively improves settlement predictions along temporal and 3D spatial dimensions,with minimal latency(e.g.within minutes),as monitoring data appear.In addition,the proposed approach requires only a reasonably small number of 3D FEM model evaluations,avoids the use of widely adopted yet often criticized surrogate models,and effectively addresses the limitations(e.g.computational inefficiency)of existing 3D model updating approaches.
基金supported by the 3D Model Library of Geo-hazards in the 3GR (No. SXJC-3ZH1A7)the software development of 3D area disaster geology map in the 3GR (No. SXJC-3ZH1A6)+1 种基金survey data acquisition and geologic map CAD system in the 3GR (No. SXKY4-02)985 Platform Projects,3D modeling and space analysis system of geo-hazards and the National Natural Science Foundation of China (No. 41172300)
文摘Taking hundreds of pieces of hazardous geological maps （1 ： 10 000） of Three Gorges res-ervoir area （3GR） as background, we establish regional three-dimensional （3D） geo-hazard modelusing DEM （digital elevation model） superposed surface images and geo-hazards elements. Based on landslides and other geo-hazard survey data,using improved B-REP（boundary representa-tion）entity data structure （two-body 3D data structure）, we set up 3D solid models for each hazardous bodies in each hazardous geological maps. Then we integrate the two types of 3D models with different scales from area to point, which are the regional geo-hazard 3D model and the solid models of each disaster body, in order to provide a visual processing and analysis plat-form for danger partition, stability evaluation, disaster prevention and control, early warning and command.
基金finically supported by the National Natural Science Foundation of China(Grant Nos.41272349 and 51322906)Key Research Program of Frontier Sciences,Chinese Academy of Sciences(Grant No.QYZDB-SSW-DQC029)
文摘This paper focuses on the progress in geomechanical modeling associated with carbon dioxide(CO2)geological storage.The detailed review of some geomechanical aspects,including numerical methods,stress analysis,ground deformation,fault reactivation,induced seismicity and crack propagation,is presented.It is indicated that although all the processes involved are not fully understood,integration of all available data,such as ground survey,geological conditions,microseismicity and ground level deformation,has led to many new insights into the rock mechanical response to CO2injection.The review also shows that in geomechanical modeling,continuum modeling methods are predominant compared with discontinuum methods.It is recommended to develop continuum-discontinuum numerical methods since they are more convenient for geomechanical modeling of CO2geological storage,especially for fracture propagation simulation.The Mohr-Coulomb criterion is widely used in prediction of rock mass mechanical behavior.It would be better to use a criterion considering the effect of the intermediate principal stress on rock mechanical behavior,especially for the stability analysis of deeply seated rock engineering.Some challenges related to geomechanical modeling of CO2geological storage are also discussed.
文摘Mine safety have top-five disasters,which including the water,gas,fire,dust and geological dynamic disaster.The coal mine water disaster is one of the important factors which restricted the development of China’s coal production.It is showed by statistics that 60%of mine accidents are affected by groundwater,which not only result in the production losses,casualties and a variety of
基金Supported by the China National Science and Technology Major Project(2017ZX05005-004-002,2016ZX05031-002-001)National Natural Science Foundation of China(41872138)Open Foundation of Top Disciplines in Yangtze University(2019KFJJ0818029)。
文摘An orthogonal 2D training image is constructed from the geological analysis results of well logs and sedimentary facies;the 2 D probabilities in three directions are obtained through linear pooling method and then aggregated by the logarithmic linear pooling to determine the 3 D multi-point pattern probabilities at the unknown points,to realize the reconstruction of a 3 D model from 2D cross-section.To solve the problems of reducing pattern variability in the 2 D training image and increasing sampling uncertainty,an adaptive spatial sampling method is introduced,and an iterative simulation strategy is adopted,in which sample points from the region with higher reliability of the previous simulation results are extracted to be additional condition points in the following simulation to improve the pattern probability sampling stability.The comparison of lateral accretion layer conceptual models shows that the reconstructing algorithm using self-adaptive spatial sampling can improve the accuracy of pattern sampling and rationality of spatial structure characteristics,and accurately reflect the morphology and distribution pattern of the lateral accretion layer.Application of the method in reconstructing the meandering river reservoir of the Cretaceous McMurray Formation in Canada shows that the new method can accurately reproduce the shape,spatial distribution pattern and development features of complex lateral accretion layers in the meandering river reservoir under tide effect.The test by sparse wells shows that the simulation accuracy is above 85%,and the coincidence rate of interpretation and prediction results of newly drilled horizontal wells is up to 80%.
基金The Construction S&T Project of the Department of Transportation of Sichuan Province(Grant No.2023A02)the National Natural Science Foundation of China(No.52109135).
文摘The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.
文摘The essential for microseismic monitoring is fast and accurate calculation of seismic wave source location. The precision of most traditional microseismic monitoring processes of mines, using TDOA location method in two-dimensional space to position the microseismic events, as well as the accuracy of positioning microseismic events, may be reduced by the two-dimensional model and simple method, and ill-conditioned equations produced by TDOA location method will increase the positioning error. This article, based on inversion theory, studies the mathematical model of TDOA location method, polariza- tion analysis location method, and comprehensive difference location method of adding angle factor in the traditional TDOA location method. The feasibility of three methods is verified by numerical simulation and analysis of the positioning error of them. The results show that the comprehensive location method of adding angle difference has strong positioning stability and high positioning accuracy, and it may reduce the impact effectively about ill-conditioned equations to positioning results. Comprehensive location method with the data of actual measure may get better positioning results.
基金supported by the Joint Fund of Seismological Science(Grant No.U1839206)the National R&D Program on Monitoring,Early Warning and Prevention of Major Natural Disaster(Grant No.2017YFC1500301)+2 种基金supported by IGGCAS Research Start-up Funds(Grant No.E0515402)National Natural Science Foundation of China(Grant No.E1115401)supported by National Natural Science Foundation of China(Grant No.11971258).
文摘The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.
基金Supported by the Innovative Research Groups of National Natural Science Foundation of China(No.51321065)the Major Science and Technology Program for Water Pollution Control and Treatment(2012ZX07101-008)the National Natural Science Foundation of China(No.51439005)
文摘In the present study, considering the transport and transformation processes of variables, a threedimensional water quality model for the river system was established, which coupled the volume of fluid(VOF) method with the k-ε turbulence mathematical model. Then, the water hydrodynamic characteristics and transport processes for BOD_5, NH_(3^-)N and TP were analyzed. The results showed that the water surface of convex bank was a little lower than that of concave bank due to the centrifugal force near the bend, and most concentrations were inferior to the type Ⅴ standard indexes of surface water environmental quality. The model validation indicated that the errors between the simulated and monitored values were comparatively small, satisfying the application demands and providing scientific basis and decision support for the restoration and protection of water quality.
