To speed up three-dimensional （3D） DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method （AGMG）. We first discretize the differential equation of the secondar...To speed up three-dimensional （3D） DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method （AGMG）. We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner （AGMG-CG） to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling （3DDCXH）. In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity.展开更多
A three-dimensional path-planning approach has been developed to coordinate multiple fixed-wing unmanned aerial vehicles(UAVs)while avoiding collisions.The hierarchical path-planning architecture that divides the path...A three-dimensional path-planning approach has been developed to coordinate multiple fixed-wing unmanned aerial vehicles(UAVs)while avoiding collisions.The hierarchical path-planning architecture that divides the path-planning process into two layers is proposed by designing the velocityobstacle strategy for satisfying timeliness and effectiveness.The upper-level layer focuses on creating an efficient Dubins initial path considering the dynamic constraints of the fixed wing.Subsequently,the lower-level layer detects potential collisions and adjusts its flight paths to avoid collisions by using the threedimensional velocity obstacle method,which describes the maneuvering space of collision avoidance as the intersection space of half space.To further handle the dynamic and collisionavoidance constraints,a priority mechanism is designed to ensure that the adjusted path is still feasible for fixed-wing UAVs.Simulation experiments demonstrate the effectiveness of the proposed method.展开更多
A parallel finite element scheme for 3D resistivity method forward modeling is introduced in this article.The domain decomposition algorithm,along with a message passing interface,is used to implement parallelism.The ...A parallel finite element scheme for 3D resistivity method forward modeling is introduced in this article.The domain decomposition algorithm,along with a message passing interface,is used to implement parallelism.The computational domain is divided into subdomains,and mesh partitioning is combined with load balancing.Unstructured meshes and local mesh refinement strategies are used to realize high precision for complex topography models.Furthermore,an improved linear solver for multi-electrode resistivity method modeling is adopted.Recycling preconditioned conjugate gradient,which is a linear solver,is based on the similarity of linear systems between point sources.The multiple right-hand-side linear systems corresponding to different point source positions are constructed,and the accelerated convergence is obtained through recycling subspace using the linear solver.The computational accuracy and efficiency of the forward scheme for complex topography models are verified using the numerical test results.展开更多
In the leakage detection of reservoir dam bodies,traditional geophysical methods can only achieve one-time detection.Meanwhile,due to the non-uniqueness of geophysicalin version,how to improve the fidelity of geophysi...In the leakage detection of reservoir dam bodies,traditional geophysical methods can only achieve one-time detection.Meanwhile,due to the non-uniqueness of geophysicalin version,how to improve the fidelity of geophysical detection inversion profiles has become a key challenge in the industry.This study aims to construct a long-term real-time monitoring system using the time-lapse resistivity method,reveal the spatiotemporal evolution law of dam leakage,and provide technical support for accurate treatment.By integrating the Internet of Things(IoT),5G technology,and AI technology,real-time data acquisition,real-time transmission,and automatic inversionare realized.Through dynamic imaging analysis of the electrical anomaly characteristics of the leakage area and comparison between corresponding rainfall events,the leakage range and resistivity changes,the reliability and efficiency of dam leakagedete ction are significantly improved.This achieves long-termdynamic monitoring of dam leakage and provides a new perspective for the safe operation and maintenance of reservoirs.展开更多
In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow e...In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow equation.The velocity and pressure are computed simultaneously.The accuracy of velocity is improved one order.The concentration equation is solved by using mixed finite element,multi-step difference and upwind approximation.A multi-step method is used to approximate time derivative for improving the accuracy.The upwind approximation and an expanded mixed finite element are adopted to solve the convection and diffusion,respectively.The composite method could compute the diffusion flux and its gradient.It possibly becomes an eficient tool for solving convection-dominated diffusion problems.Firstly,the conservation of mass holds.Secondly,the multi-step method has high accuracy.Thirdly,the upwind approximation could avoid numerical dispersion.Using numerical analysis of a priori estimates and special techniques of differential equations,we give an error estimates for a positive definite problem.Numerical experiments illustrate its computational efficiency and feasibility of application.展开更多
Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell＇s equations for an anisotropic mo...Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell＇s equations for an anisotropic model, we introduce scalar potentials based on the divergence-free characteristic of the electric and magnetic （EM） fields. We then continue the EM fields down into the deep earth and upward into the seawater and couple them at the ocean bottom to the transmitting source. By studying both the DC apparent resistivity curves and their polar plots, we can resolve the anisotropy of the ocean bottom. Forward modeling of a high-resistivity thin layer in an anisotropic half-space demonstrates that the marine DC resistivity method in shallow water is very sensitive to the resistive reservoir but is not influenced by airwaves. As such, it is very suitable for oil and gas exploration in shallowwater areas but, to date, most modeling algorithms for studying marine DC resistivity are based on isotropic models. In this paper, we investigate one-dimensional anisotropic forward modeling for marine DC resistivity method, prove the algorithm to have high accuracy, and thus provide a theoretical basis for 2D and 3D forward modeling.展开更多
The hydrodynamic characteristics of a rigid, single, circular cylinder in a three dimensional, incompressible, uniform cross flow were calculated using the large-eddy simulation method of CFX5. Solutions to the three ...The hydrodynamic characteristics of a rigid, single, circular cylinder in a three dimensional, incompressible, uniform cross flow were calculated using the large-eddy simulation method of CFX5. Solutions to the three dimensional N-S equations were obtained by the finite volume method. The focus of this numerical simulation was to research the characteristics of pressure distribution （drag and litt forces） and vortex tubes at high Reynolds numbers. The results of the calculations showed that the forces at every section in the spanwise direction of the cylinder were symmetrical about the middle section and smaller than the forces calculated in two dimensional cases. Moreover, the flow around the cylinder obviously presents three dimensional characteristics.展开更多
The research about subsurface characteristics by using transient electromagnetic method（TEM） and high density resistivity method（HDRM） were already conducted in Ordos. The objective of this research is to detect c...The research about subsurface characteristics by using transient electromagnetic method（TEM） and high density resistivity method（HDRM） were already conducted in Ordos. The objective of this research is to detect coalmine goaf areas based on rock resistivity. The data processing using wavelet transform, three point smoothing, RES2 DINV and Maxwell processing software to obtain 2D resistivity structure. The results showed that the layers with maximum resistivity values（30e33 U m on Line 1, 30e31 U m on Line 2, 32e40 U m on Line3） are founded at station 1e7, and 14e20 on Line 1,13e18 on Line 2, and 8e13 and 16e20 on Line 3 which is predicted as goaf layer, and the minimum resistivity values（20e26 U m of TEM, 45e75 U m of HDRM） at the other layers. This resistivity difference was caused by the geology and characteristics of the study area which is located close by the cleugh with rich coal, so the goaf area distinguishable with aquifer layer and coal seam. The results were also significant accidents and serious destruction of ecological environment.展开更多
Electrical measurement was employed to investigate the early hydration characteristics of cement pastes with different dosages of superplasticizer in the same W/C ratio. The hyperbolic method was applied to analyze th...Electrical measurement was employed to investigate the early hydration characteristics of cement pastes with different dosages of superplasticizer in the same W/C ratio. The hyperbolic method was applied to analyze the electrical resistivity development. The peak point （Ph） on the hyperbolic curve could be easily read. The time （th） to reach the point Ph had strong relations with the setting time. th was delayed with the increment of the dosage of superplasticizer. The time th was used to plot the relationship between the initial setting time and final setting time. The hyperbolic equation was established to predict the ultimate resistivity. The retardation effect of the superplasticizer was confirmed in the same W/C ratio by setting time and isothermal heat evolution.展开更多
Microseismic/acoustic emission(MS/AE)source localization method is crucial for predicting and controlling of potentially dangerous sources of complex structures.However,the locating errors induced by both the irregula...Microseismic/acoustic emission(MS/AE)source localization method is crucial for predicting and controlling of potentially dangerous sources of complex structures.However,the locating errors induced by both the irregular structure and pre-measured velocity are poorly understood in existing methods.To meet the high-accuracy locating requirements in complex three-dimensional hole-containing structures,a velocity-free MS/AE source location method is developed in this paper.It avoids manual repetitive training by using equidistant grid points to search the path,which introduces A*search algorithm and uses grid points to accommodate complex structures with irregular holes.It also takes advantage of the velocity-free source location method.To verify the validity of the proposed method,lead-breaking tests were performed on a cubic concrete test specimen with a size of 10 cm10 cm10 cm.It was cut out into a cylindrical empty space with a size of/6cm10 cm.Based on the arrivals,the classical Geiger method and the proposed method are used to locate lead-breaking sources.Results show that the locating error of the proposed method is 1.20 cm,which is less than 2.02 cm of the Geiger method.Hence,the proposed method can effectively locate sources in the complex three-dimensional structure with holes and achieve higher precision requirements.展开更多
The estimation of shear strength of rock mass discontinuity is always a focal, but difficult, problem in the field of geotechnical engineering. Considering the disadvantages and limitation of exist- ing estimation met...The estimation of shear strength of rock mass discontinuity is always a focal, but difficult, problem in the field of geotechnical engineering. Considering the disadvantages and limitation of exist- ing estimation methods, a new approach based on the shadow area percentage （SAP） that can be used to quantify surface roughness is proposed in this article. Firstly, by the help of laser scanning technique, the three-dimensional model of the surface of rock discontinuity was established. Secondly, a light source was simulated, and there would be some shadows produced on the model surface. Thirdly, to obtain the value of SAP of each specimen, the shadow detection technique was introduced for use. Fourthly, compared with the result from direct shear testing and based on statistics, an empirical for- mula was found among SAP, normal stress, and shear strength. Data of Yujian （~ River were used as an example, and the following conclusions have been made. （1） In the case of equal normal stress, the peak shear stress is positively proportional to the SAP. （2） The formula for estimating was derived, and the predictions of peak-shear strength made with this equation well agreed with the experimental re- suits obtained in laboratory tests.展开更多
Based on the particle-in-cell technology and the secondary electron emission theory, a three-dimensional simulation method for multipactor is presented in this paper. By combining the finite difference time domain met...Based on the particle-in-cell technology and the secondary electron emission theory, a three-dimensional simulation method for multipactor is presented in this paper. By combining the finite difference time domain method and the panicle tracing method, such an algorithm is self-consistent and accurate since the interaction between electromagnetic fields and particles is properly modeled. In the time domain aspect, the generation of multipactor can be easily visualized, which makes it possible to gain a deeper insight into the physical mechanism of this effect. In addition to the classic secondary electron emission model, the measured practical secondary electron yield is used, which increases the accuracy of the algorithm. In order to validate the method, the impedance transformer and ridge waveguide filter are studied. By analyzing the evolution of the secondaries obtained by our method, multipactor thresholds of these components are estimated, which show good agreement with the experimental results. Furthermore, the most sensitive positions where multipactor occurs are determined from the phase focusing phenomenon, which is very meaningful for multipactor analysis and design.展开更多
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.展开更多
The objective of this work is to develop a novel methodology for determining real resistivity of conductive asphalt concrete based on two-electrode method.Due to an influence of contact resistance,the measured resisti...The objective of this work is to develop a novel methodology for determining real resistivity of conductive asphalt concrete based on two-electrode method.Due to an influence of contact resistance,the measured resistivity is always not equal to the real resistivity.To determine the real resistivity,a linear relationship of the measured resistivity,contact resistance and the real resistivity was established.Then experiments for six specimens with varying graphite contents were designed and performed to validate the formulation.Results of experiments demonstrate that the slope of the line represents contact resistance,and the intercept indicates the real resistivity.The effects of graphite content on contact resistance and real resistivity are also revealed.Finally,results show that the influence of contact resistance on accuracy of resisitvity measurement becomes more serious if graphite content is beyond 3%.Hence,it is the time to choose this novel methodology to determine the real resistivity of asphalt concrete by taking account of contact resistance.展开更多
The detecting examples using the high density resistivity method, about the evaluation of containing water characteristic from the floor rock and the height of overburden failure, were given. It expresses that the hig...The detecting examples using the high density resistivity method, about the evaluation of containing water characteristic from the floor rock and the height of overburden failure, were given. It expresses that the high density resistivity method has good effect with strong maneuverability and continuous observing section during the prevention and cure for mine water disaster. At the same time, the article pointed out that the study of space data inversion and dynamic watching technology about the high density resistivity method must be enhanced in the future because of special condition of data collecting in mine.展开更多
The compressive strength development of Portland cement pastes was investigated by the electrical resistivity method and the maturity method.The experiments were carried out on the cement pastes with different water-c...The compressive strength development of Portland cement pastes was investigated by the electrical resistivity method and the maturity method.The experiments were carried out on the cement pastes with different water-cement ratios at different curing temperatures.The results show that the application of the maturity method has limitation to obtain the strength.It is found that both of the compressive strength and the electrical resistivity follow hyperbolic trend for all the mixes.The hyperbolic equation of each mix is obtained to estimate the ultimate resistivity value which can probably be reached.The relationship between electrical resistivity and compressive strength of the cement pastes is established based on the test results and interpreted by the empirical Archie equation and a strength-porosity equation.The relationship between the electrical resistivity after temperature correction and the compressive strength was linear and independent of curing temperature and water-cement ratio.展开更多
A new method,the stream surface strip element method,for simulating the three-dimensional deformation of plate and strip rolling process was proposed.The rolling deformation zone was divided into a number of stream su...A new method,the stream surface strip element method,for simulating the three-dimensional deformation of plate and strip rolling process was proposed.The rolling deformation zone was divided into a number of stream surface(curved surface)strip elements along metal flow traces,and the stream surface strip elements were mapped into the corresponding plane strip elements for analysis and computation.The longitudinal distributions of the lateral displacement and the altitudinal displacement of metal were respectively constructed to be a quartic curve and a quadratic curve,of which the lateral distributions were expressed as the third-power spline function,and the altitudinal distributions were fitted in the quadratic curve.From the flow theory of plastic mechanics,the mathematical models of the three-dimensional deformations and stresses of the deformation zone were constructed.Compared with the streamline strip element method proposed by the first author of this paper,the stream surface strip element method takes into account the uneven distributions of stresses and deformations along altitudinal direction,and realizes the precise three-dimensional analysis and computation.The simulation example of continuous hot rolled strip indicates that the method and the model accord with facts and provide a new reliable engineering-computation method for the three-dimensional mechanics simulation of plate and strip rolling process.展开更多
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.展开更多
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.展开更多
A meshless approach, called the rigid-plastic reproducing kernel particle method （RKPM）, is presented for three-dimensional （3D） bulk metal forming simulation. The approach is a combination of RKPM with the flow t...A meshless approach, called the rigid-plastic reproducing kernel particle method （RKPM）, is presented for three-dimensional （3D） bulk metal forming simulation. The approach is a combination of RKPM with the flow theory of 3D rigid-plastic mechanics. For the treatments of essential boundary conditions and incompressibility constraint, the boundary singular kernel method and the modified penalty method are utilized, respectively. The arc-tangential friction model is employed to treat the contact conditions. The compression of rectangular blocks, a typical 3D upsetting operation, is analyzed for different friction conditions and the numerical results are compared with those obtained using commercial rigid-plastic FEM （finite element method） software Deform^3D. As results show, when handling 3D plastic deformations, the proposed approach eliminates the need of expensive meshing and remeshing procedures which are unavoidable in conventional FEM and can provide results that are in good agreement with finite element predictions.展开更多
基金supported by the Natural Science Foundation of China(Nos.41404057,41674077 and 411640034)the Nuclear Energy Development Project of China,and the‘555’Project of Gan Po Excellent People
文摘To speed up three-dimensional （3D） DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method （AGMG）. We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner （AGMG-CG） to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling （3DDCXH）. In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity.
基金supported by the National Science Fund for Distinguished Young Scholars(52425211)BIT Research Fund Program for Young Scholars(XSQD-202201005).
文摘A three-dimensional path-planning approach has been developed to coordinate multiple fixed-wing unmanned aerial vehicles(UAVs)while avoiding collisions.The hierarchical path-planning architecture that divides the path-planning process into two layers is proposed by designing the velocityobstacle strategy for satisfying timeliness and effectiveness.The upper-level layer focuses on creating an efficient Dubins initial path considering the dynamic constraints of the fixed wing.Subsequently,the lower-level layer detects potential collisions and adjusts its flight paths to avoid collisions by using the threedimensional velocity obstacle method,which describes the maneuvering space of collision avoidance as the intersection space of half space.To further handle the dynamic and collisionavoidance constraints,a priority mechanism is designed to ensure that the adjusted path is still feasible for fixed-wing UAVs.Simulation experiments demonstrate the effectiveness of the proposed method.
基金supported by the National Natural Science Foundation of China(No:42274182)Guangxi Natural Science Foundation(No:2020 GXNSFAA297079)。
文摘A parallel finite element scheme for 3D resistivity method forward modeling is introduced in this article.The domain decomposition algorithm,along with a message passing interface,is used to implement parallelism.The computational domain is divided into subdomains,and mesh partitioning is combined with load balancing.Unstructured meshes and local mesh refinement strategies are used to realize high precision for complex topography models.Furthermore,an improved linear solver for multi-electrode resistivity method modeling is adopted.Recycling preconditioned conjugate gradient,which is a linear solver,is based on the similarity of linear systems between point sources.The multiple right-hand-side linear systems corresponding to different point source positions are constructed,and the accelerated convergence is obtained through recycling subspace using the linear solver.The computational accuracy and efficiency of the forward scheme for complex topography models are verified using the numerical test results.
文摘In the leakage detection of reservoir dam bodies,traditional geophysical methods can only achieve one-time detection.Meanwhile,due to the non-uniqueness of geophysicalin version,how to improve the fidelity of geophysical detection inversion profiles has become a key challenge in the industry.This study aims to construct a long-term real-time monitoring system using the time-lapse resistivity method,reveal the spatiotemporal evolution law of dam leakage,and provide technical support for accurate treatment.By integrating the Internet of Things(IoT),5G technology,and AI technology,real-time data acquisition,real-time transmission,and automatic inversionare realized.Through dynamic imaging analysis of the electrical anomaly characteristics of the leakage area and comparison between corresponding rainfall events,the leakage range and resistivity changes,the reliability and efficiency of dam leakagedete ction are significantly improved.This achieves long-termdynamic monitoring of dam leakage and provides a new perspective for the safe operation and maintenance of reservoirs.
基金supported by the Natural Science Foundation of Shandong Province(ZR2021MA019)the National Natural Science Foundation of China(11871312)。
文摘In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow equation.The velocity and pressure are computed simultaneously.The accuracy of velocity is improved one order.The concentration equation is solved by using mixed finite element,multi-step difference and upwind approximation.A multi-step method is used to approximate time derivative for improving the accuracy.The upwind approximation and an expanded mixed finite element are adopted to solve the convection and diffusion,respectively.The composite method could compute the diffusion flux and its gradient.It possibly becomes an eficient tool for solving convection-dominated diffusion problems.Firstly,the conservation of mass holds.Secondly,the multi-step method has high accuracy.Thirdly,the upwind approximation could avoid numerical dispersion.Using numerical analysis of a priori estimates and special techniques of differential equations,we give an error estimates for a positive definite problem.Numerical experiments illustrate its computational efficiency and feasibility of application.
基金financially supported by the National Hi-tech Research and Development Program of China(863 Program)(No.2012AA09A20103)
文摘Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell＇s equations for an anisotropic model, we introduce scalar potentials based on the divergence-free characteristic of the electric and magnetic （EM） fields. We then continue the EM fields down into the deep earth and upward into the seawater and couple them at the ocean bottom to the transmitting source. By studying both the DC apparent resistivity curves and their polar plots, we can resolve the anisotropy of the ocean bottom. Forward modeling of a high-resistivity thin layer in an anisotropic half-space demonstrates that the marine DC resistivity method in shallow water is very sensitive to the resistive reservoir but is not influenced by airwaves. As such, it is very suitable for oil and gas exploration in shallowwater areas but, to date, most modeling algorithms for studying marine DC resistivity are based on isotropic models. In this paper, we investigate one-dimensional anisotropic forward modeling for marine DC resistivity method, prove the algorithm to have high accuracy, and thus provide a theoretical basis for 2D and 3D forward modeling.
文摘The hydrodynamic characteristics of a rigid, single, circular cylinder in a three dimensional, incompressible, uniform cross flow were calculated using the large-eddy simulation method of CFX5. Solutions to the three dimensional N-S equations were obtained by the finite volume method. The focus of this numerical simulation was to research the characteristics of pressure distribution （drag and litt forces） and vortex tubes at high Reynolds numbers. The results of the calculations showed that the forces at every section in the spanwise direction of the cylinder were symmetrical about the middle section and smaller than the forces calculated in two dimensional cases. Moreover, the flow around the cylinder obviously presents three dimensional characteristics.
基金supported by the Institute of Seismology Foundation, China Earthquake Administration (201326126)
文摘The research about subsurface characteristics by using transient electromagnetic method（TEM） and high density resistivity method（HDRM） were already conducted in Ordos. The objective of this research is to detect coalmine goaf areas based on rock resistivity. The data processing using wavelet transform, three point smoothing, RES2 DINV and Maxwell processing software to obtain 2D resistivity structure. The results showed that the layers with maximum resistivity values（30e33 U m on Line 1, 30e31 U m on Line 2, 32e40 U m on Line3） are founded at station 1e7, and 14e20 on Line 1,13e18 on Line 2, and 8e13 and 16e20 on Line 3 which is predicted as goaf layer, and the minimum resistivity values（20e26 U m of TEM, 45e75 U m of HDRM） at the other layers. This resistivity difference was caused by the geology and characteristics of the study area which is located close by the cleugh with rich coal, so the goaf area distinguishable with aquifer layer and coal seam. The results were also significant accidents and serious destruction of ecological environment.
基金the National Natural Science Foundation of China(No.50778078)
文摘Electrical measurement was employed to investigate the early hydration characteristics of cement pastes with different dosages of superplasticizer in the same W/C ratio. The hyperbolic method was applied to analyze the electrical resistivity development. The peak point （Ph） on the hyperbolic curve could be easily read. The time （th） to reach the point Ph had strong relations with the setting time. th was delayed with the increment of the dosage of superplasticizer. The time th was used to plot the relationship between the initial setting time and final setting time. The hyperbolic equation was established to predict the ultimate resistivity. The retardation effect of the superplasticizer was confirmed in the same W/C ratio by setting time and isothermal heat evolution.
基金The authors wish to acknowledge financial support from the National Natural Science Foundation of China(51822407 and 51774327)Natural Science Foundation of Hunan Province in China(2018JJ1037)Innovation Driven project of Central South University(2020CX014).
文摘Microseismic/acoustic emission(MS/AE)source localization method is crucial for predicting and controlling of potentially dangerous sources of complex structures.However,the locating errors induced by both the irregular structure and pre-measured velocity are poorly understood in existing methods.To meet the high-accuracy locating requirements in complex three-dimensional hole-containing structures,a velocity-free MS/AE source location method is developed in this paper.It avoids manual repetitive training by using equidistant grid points to search the path,which introduces A*search algorithm and uses grid points to accommodate complex structures with irregular holes.It also takes advantage of the velocity-free source location method.To verify the validity of the proposed method,lead-breaking tests were performed on a cubic concrete test specimen with a size of 10 cm10 cm10 cm.It was cut out into a cylindrical empty space with a size of/6cm10 cm.Based on the arrivals,the classical Geiger method and the proposed method are used to locate lead-breaking sources.Results show that the locating error of the proposed method is 1.20 cm,which is less than 2.02 cm of the Geiger method.Hence,the proposed method can effectively locate sources in the complex three-dimensional structure with holes and achieve higher precision requirements.
基金supported by the China Geological Survey (No.1212011014030)the Major State Basic Research Development Program of China (973 Program) (No.2011CB710600)
文摘The estimation of shear strength of rock mass discontinuity is always a focal, but difficult, problem in the field of geotechnical engineering. Considering the disadvantages and limitation of exist- ing estimation methods, a new approach based on the shadow area percentage （SAP） that can be used to quantify surface roughness is proposed in this article. Firstly, by the help of laser scanning technique, the three-dimensional model of the surface of rock discontinuity was established. Secondly, a light source was simulated, and there would be some shadows produced on the model surface. Thirdly, to obtain the value of SAP of each specimen, the shadow detection technique was introduced for use. Fourthly, compared with the result from direct shear testing and based on statistics, an empirical for- mula was found among SAP, normal stress, and shear strength. Data of Yujian （~ River were used as an example, and the following conclusions have been made. （1） In the case of equal normal stress, the peak shear stress is positively proportional to the SAP. （2） The formula for estimating was derived, and the predictions of peak-shear strength made with this equation well agreed with the experimental re- suits obtained in laboratory tests.
基金Project supported by the National Key Laboratory Foundation,China(Grant No.9140C530103110C5301)
文摘Based on the particle-in-cell technology and the secondary electron emission theory, a three-dimensional simulation method for multipactor is presented in this paper. By combining the finite difference time domain method and the panicle tracing method, such an algorithm is self-consistent and accurate since the interaction between electromagnetic fields and particles is properly modeled. In the time domain aspect, the generation of multipactor can be easily visualized, which makes it possible to gain a deeper insight into the physical mechanism of this effect. In addition to the classic secondary electron emission model, the measured practical secondary electron yield is used, which increases the accuracy of the algorithm. In order to validate the method, the impedance transformer and ridge waveguide filter are studied. By analyzing the evolution of the secondaries obtained by our method, multipactor thresholds of these components are estimated, which show good agreement with the experimental results. Furthermore, the most sensitive positions where multipactor occurs are determined from the phase focusing phenomenon, which is very meaningful for multipactor analysis and design.
基金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.
基金Project(51178348)supported by the National Natural Science Foundation of China
文摘The objective of this work is to develop a novel methodology for determining real resistivity of conductive asphalt concrete based on two-electrode method.Due to an influence of contact resistance,the measured resistivity is always not equal to the real resistivity.To determine the real resistivity,a linear relationship of the measured resistivity,contact resistance and the real resistivity was established.Then experiments for six specimens with varying graphite contents were designed and performed to validate the formulation.Results of experiments demonstrate that the slope of the line represents contact resistance,and the intercept indicates the real resistivity.The effects of graphite content on contact resistance and real resistivity are also revealed.Finally,results show that the influence of contact resistance on accuracy of resisitvity measurement becomes more serious if graphite content is beyond 3%.Hence,it is the time to choose this novel methodology to determine the real resistivity of asphalt concrete by taking account of contact resistance.
文摘The detecting examples using the high density resistivity method, about the evaluation of containing water characteristic from the floor rock and the height of overburden failure, were given. It expresses that the high density resistivity method has good effect with strong maneuverability and continuous observing section during the prevention and cure for mine water disaster. At the same time, the article pointed out that the study of space data inversion and dynamic watching technology about the high density resistivity method must be enhanced in the future because of special condition of data collecting in mine.
基金Funding by the National Natural Science Foundation of China (Nos.50778078 and 51178202)the Doctoral Research Fund from Wuhan Institute of Technology
文摘The compressive strength development of Portland cement pastes was investigated by the electrical resistivity method and the maturity method.The experiments were carried out on the cement pastes with different water-cement ratios at different curing temperatures.The results show that the application of the maturity method has limitation to obtain the strength.It is found that both of the compressive strength and the electrical resistivity follow hyperbolic trend for all the mixes.The hyperbolic equation of each mix is obtained to estimate the ultimate resistivity value which can probably be reached.The relationship between electrical resistivity and compressive strength of the cement pastes is established based on the test results and interpreted by the empirical Archie equation and a strength-porosity equation.The relationship between the electrical resistivity after temperature correction and the compressive strength was linear and independent of curing temperature and water-cement ratio.
基金Sponsored by National Natural Science Foundation of China(50175095)Provincial Natural Science Foundation of Hebei of China(502173)
文摘A new method,the stream surface strip element method,for simulating the three-dimensional deformation of plate and strip rolling process was proposed.The rolling deformation zone was divided into a number of stream surface(curved surface)strip elements along metal flow traces,and the stream surface strip elements were mapped into the corresponding plane strip elements for analysis and computation.The longitudinal distributions of the lateral displacement and the altitudinal displacement of metal were respectively constructed to be a quartic curve and a quadratic curve,of which the lateral distributions were expressed as the third-power spline function,and the altitudinal distributions were fitted in the quadratic curve.From the flow theory of plastic mechanics,the mathematical models of the three-dimensional deformations and stresses of the deformation zone were constructed.Compared with the streamline strip element method proposed by the first author of this paper,the stream surface strip element method takes into account the uneven distributions of stresses and deformations along altitudinal direction,and realizes the precise three-dimensional analysis and computation.The simulation example of continuous hot rolled strip indicates that the method and the model accord with facts and provide a new reliable engineering-computation method for the three-dimensional mechanics simulation of plate and strip rolling process.
基金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.
基金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.
基金This work was supported by the National Natural Science Foundation of China （No. 50275094）.
文摘A meshless approach, called the rigid-plastic reproducing kernel particle method （RKPM）, is presented for three-dimensional （3D） bulk metal forming simulation. The approach is a combination of RKPM with the flow theory of 3D rigid-plastic mechanics. For the treatments of essential boundary conditions and incompressibility constraint, the boundary singular kernel method and the modified penalty method are utilized, respectively. The arc-tangential friction model is employed to treat the contact conditions. The compression of rectangular blocks, a typical 3D upsetting operation, is analyzed for different friction conditions and the numerical results are compared with those obtained using commercial rigid-plastic FEM （finite element method） software Deform^3D. As results show, when handling 3D plastic deformations, the proposed approach eliminates the need of expensive meshing and remeshing procedures which are unavoidable in conventional FEM and can provide results that are in good agreement with finite element predictions.
