UAV-based transient electromagnetic 3D forward modeling and inversion and analysis of exploration capability

Unmanned aerial vehicle transient electromagnetic(UAV-TEM)is a novel airborne exploration method that offers advantages such as low cost,simple operation,high exploration efficiency and suitability for near-surface exploration in complex terrain areas.To improve the accuracy of data interpretation in this method,the authors conducted a systematic three-dimensional(3D)forward modeling and inversion of the UAV-TEM.This study utilized the finite element method based on unstructured tetrahedral elements and employed the second-order backward Euler method for time discretization.This allowed for accurate 3D modeling and accounted for the effects of complex terrain.Based on these,the influence characteristics of flight altitudes and the sizes,burial depths,and resistivities of anomalies are compared and analyzed to explore the UAV-TEM systems’exploration capability.Lastly,four typical geoelectrical models of landslides are designed,and the inversion method based on the Gauss-Newton optimization method is used to image the landslide models and analyze the imaging effect of the UAV-TEM method on landslide geohazards.Numerical results showed that UAV-TEM could have better exploration resolution and fine imaging of nearsurface structures,providing important technical support for monitoring,early warning,and preventing landslides and other geological hazards.

3D time-domain forward modeling of airborne transient electromagnetism considering superparamagnetic effect

The superparamagnetic effect arises from the superparamagnetism exhibited by a multitude of nano-sized magnetic mineral particles under an external electric field.This phenomenon manifests in transient electromagnetic data primarily as a deceleration in the attenuation rate of late-stage signals,a characteristic difficult to discern directly from airborne transient electromagnetic signals,consequently leading to significant misinterpretations of subterranean electrical structures.This study embarks on 3D forward modeling of airborne electromagnetic responses in the frequency domain,accounting for the superparamagnetic effect,utilizing an unstructured finite element method.Superparamagnetic responses in the time domain were obtained through frequency-time conversion.This investigation explores the influence of various parameters-such as magnetic susceptibility,time constants,and flight altitude-on the superparamagnetic effect by examining the response characteristics of typical targets.Findings indicate that in its late stages,the superparamagnetic effect can induce a relative anomaly of up to 300%.There is a positive correlation between magnetic susceptibility and the strength of the superparamagnetic effect.The influence of the time constant's upper and lower limits on the superparamagnetic effect is minimal;however,the range between these limits significantly affects the effect,showing a negative correlation with its intensity.Higher flight altitudes weaken the superparamagnetic signal.The impact is most pronounced when superparamagnetic minerals are shallowly buried,effectively shielding the underlying geology with the characteristics of a good conductivity anomaly,but this effect diminishes with greater depth.The insights from this study provide a theoretical framework for a deeper understanding of the superparamagnetic effect in transient electromagnetic signals and for more accurate interpretations of subterranean geological and electrical structures.

Meshing effects of the 3-D FEM numerical modeling in seismo-electromagnetics:An application in selectivity of seismic electric signal (SES)

We investigated how density and quality of mesh around interest domain affect electromagnetic （EM） responses of 3D Earth layered media using finite element method （FEM）. Effect of different mesh shapes was also investigated using a method of mixing structured and unstructured mesh. As a case study, we estimated the effects of meshing on selectivity phenomenon of seismic electric signal （SES）. Our results suggest that the relative errors resulting from mesh effects may not be negligible, which may lead to some unconvincing explanation of the SES selectivity based on the numerical modeling results.

2.5-D modeling of cross-hole electromagnetic measurement by nite element method

A finite element method is developed for simulating frequency domain electromagnetic responses due to a dipole source in the 2-D conductive structures. Computing costs are considerably minimized by reducing the full three-dimensional problem to a series of two-dimensional problems. This is accomplished by transforming the problem into y-wave number （Ky） domain using Fourier transform and the y-axis is parallel to the structural strike. In the Ky domain, two coupled partial differential equations for magnetic field Hy and electric field Ey are derived. For a specific value of Ky, the coupled equations are solved by the finite element method with isoparametric elements in the x-z plane. Application of the inverse Fourier transform to the Ky, domain provides the electric and magnetic fields in real space. The equations derived can be applied to general complex two-dimensional structures containing either electric or magnetic dipole source in any direction. In the modeling of the electromagnetic measurement, we adopted a pseudo-delta function to distribute the dipole source current and circumvent the problem of singularity at the source point. Moreover, the suggested method used isoparametric finite elements to accommodate the complex subsurface formation. For the large scale linear system derived from the discretization of the Maxwell＇s equations, several iterative solvers were used and compared to select the optimal one. A quantitative test of accuracy was presented which compared the finite element results with analytic solutions for a dipole source in homogeneous space for different ranges and different wave numbers Ky. to validate the addressed the effects of the distribution range τ of the homogeneous medium. code and check its effectiveness. In addition, we pseudo-delta function on the numerical results in

3D forward modeling and response characteristics of low-resistivity overburden of the CFS-PML absorbing boundary for ground-well transient electromagnetic method

This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of complex frequency-shifted perfectly matched layer(CFS-PML)was used for truncation so that the low-frequency electromagnetic wave can be better absorbed at the model boundary.A typical three-dimensional(3D)homogeneous half-space model was established and a low-resistivity cube model was analyzed under the half-space condition.The response patterns and drivers of the low-resistivity cube model were discussed under the influence of a low-resistivity overburden.The absorption boundary conditions of CFS-PML significantly affected the low-frequency electromagnetic waves.For a low-resistivity cube around the borehole,its response curve exhibited a single-peak,and the extreme point of the curve corresponded to the center of the low-resistivity body.When the low-resistivity cube was directly below the borehole,the response curve showed three extreme values(two high and one low),with the low corresponding to the center of the low-resistivity body.The total field response of the low-resistivity overburden was stronger than that of the uniform half-space model due to the low-resistivity shielding effect of electromagnetic waves.When the receiving-transmitting distance gradually increased,the effect of the low-resistivity overburden was gradually weakened,and the response of the low-resistivity cube was strengthened.It was affected by the ratio of the overburden resistivity to the resistivity of the low-resistivity body.

3-D Geological Modeling-Concept,Methods and Key Techniques

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.

Application of 3-D Geoscience Modeling Technology for the Estimation of Solid Mineral Reserves

Applying new approaches, methods, and technologies for the estimation of reserves can effectively improve the efficiency and accuracy of assessments of solid mineral resources. After analyzing the development of 3-D geoscience modeling technology （3-D GMT）, this paper discusses the application of 3-D GMT for the estimation of solid mineral reserves, emphatically introducing its workflow and two key technologies, 3-D orebody surface modeling, and property modeling. Moreover, the paper analyzes the limitations of traditional methods, such as the section method and geological block method, and points out the advantages of 3-D GMT： building more accurate 3-D orebody models, expressing the internal inhomogeneous attributes of an orebody, reducing the potential for errors in the estimation of reserves, and implementing dynamic estimations of reserves.

3-D modeling of rock burst in pillar No. 19 of Fetr6 chromite mine

Fetr6 is an underground mine in which chromite is extracted using stope and pillar mining method. Despite of all improving works such as roof supporting and replacing of ore pillars with concrete pillars, pillar No. 19 failed and other pillars failed progressively as a domino effect and 4000 m2 of mine collapsed within a few minutes, consequently. For detail investigation, two 3-D numerical models were developed by 3Dec. The first, a base model, was used for estimation of stress on pillars just before failure and the other for investigation of rock burst in pillar No. 19. The results show that discontinuity parameters such as friction angle and shear stiffness is critical parameters in this pillar failure. In addition, it indicates that W/H ratio equal 0.3, the lack of ore extraction strategy and inadequate roof support are the major reasons for this failure. In this paper, the procedure of study was described.

MODELING OF 3-D IMAGES AND SPACE MARKOV CUBIC MESH MODELS

Three-dimensional(3-D)Markov cubic random mesh models are presented andproved in the form of two theorems in details.Its applications to the modeling and description of3-D images are described.The model presented here is a appropriate mathematical tool for thesegmentation,modeling,classification and other processing.Finally,an example is given.

Numerical study on residual current and its impact on mass transport in the Hangzhou Bay and the Changjiang Estuary I. A3-D joint model of the Hangzhou Bay and the Changjiang Estuary

A 3-D numerical model is set up in a large domain covering the Hangzhou Bay and the Changjiang Estuary based on the ECOM model in orthogonal curvilinear coordinates.The numerical schemes for baroclinic pressure gradient (BPG)terms and convective terms are improved in the paper according to the characteristics of velocity field and mass transport in the area.The model is validated by the simulations of residual current and salinity transport in the Hangzhou Bay and the Changjiang Estuary.

Study on three dimensional electromagnetic model for permanent magnet linear synchronous motor

In this paper, with the non-salient pole permanent magnet linear synchronous motor (PMLSM)being cited, by using Fourier transform method and "slot-by-slot", "pole -by-pole" current approach, a 3D electromagnetic field model of PMLSM is established. Special attention is paid to its structure and the influence of longitudinal and transverse end effect. The distribution of electromagnetic field of PMLSM can be obtained directly and promptly by using FFT algorithm. It can also be used for the analysis of other LSM.

Spatiotemporal Evolution Analysis of the 3-D Deformation Field of Land Subsidence in Beijing after the South-to-North Water Diversion Project

The implementation of the South-to-North Water Diversion Project (SNWDP) has alleviated groundwater resource pressure in North China to some extent, resulting in a gradual deceleration of land subsidence and even rebound in some areas. To investigate the spatiotemporal evolution characteristics of land subsidence in the eastern plain of Beijing following the SNWDP, this study employs Ascending (ASC) and Descending (DES) InSAR data combined with a Strain Model (SM) to obtain a Three-Dimensional (3-D) deformation field from 2016 to 2018. Through analysis of the 3-D deformation characteristics and spatiotemporal evolution of land subsidence in this region from 2016 to 2018, the results reveal a shift in the distribution of subsiding areas after the South-to-North Water Diversion, with a marked decrease in subsidence rates in certain areas. The maximum subsidence rate in the Beijing area has decreased to 110 mm/yr, accompanied by horizontal deformation at a rate of 12 mm/yr. Additionally, by examining the spatial relationship between major active faults and subsidence deformation in this region, the study further elucidates the influence of fault activity on the spatial distribution of subsidence deformation.

Simulation of subcritical flow pattern in 180° uniform and convergent open-channel bends using SSIIM 3-D model

In meandering rivers, the flow pattern is highly complex, with specific characteristics at bends that are not observed along straight paths. A numerical model can be effectively used to predict such flow fields. Since river bends are not uniform-some are divergent and others convergent-in this study, after the SSIIM 3-D model was calibrated using the result of measurements along a uniform 180° bend with a width of 0.6 m, a similar but convergent 180v bend, 0.6 m to 0.45 m wide, was simulated using the SSI1M 3-D numerical model. Flow characteristics of the convergent 180° bend, including lengthwise and vertical velocity profiles, primary and secondary flows, lengthwise and widtbwise slopes of the water surface, and the helical flow strength, were compared with those of the uniform 180° bend. The verification results of the model show that the numerical model can effectively simulate the flow field in the uniform bend. In addition, this research indicates that, in a convergent channel, the maximum velocity path at a plane near the water surface crosses the channel＇s centerline at about a 30° to 40° cross-section, while in the uniform bend, this occurs at about the 50° cross-section. The varying range of the water surface elevation is wider in the convergent channel than in the uniform one, and the strength of the helical flow is generally greater in the uniform channel than in the convergent one. Also, unlike the uniform bend, the convergent bend exhibits no rotational cell against the main direction of secondary flow rotation at the 135° cross-section.

3-D Magnetic Sensor Module for Locating and Tracking MEMS Swallowable Capsule Based on Scalar Form of Magnetic Dipole Model

MEMS swallowable capsule is a novel technology in the non-invasive surgery. This technology provides a way to diagnose directly into the deep intestinal where the traditional invasive technology implemented, such as X-Ray, endoscopy. It is a key for us to locate and track the position of a MEMS capsule in clinical applications. To solve this problem, we implemented a magnetic sensor module based on the scalar form of the magnetic dipole model,which was designed with very small size （5.2 ＊ 2. 1 ＊ 1.2 em） and easy to assemble to satisfy the system requirement. Here we discuss in detail the principle of magnetic dipole model, rules of selecting sensor and functions of the module. Some trials are established to test the characteristic of the module. The results of the Cm experiment demonstrates that the module follows the rules of the new magnetic dipole model form.

3-D crustal-scale gravity model of the San Rafael Block and Payenia volcanic province in Mendoza,Argentina

The San Rafael Block(SRB)is part of one of the main retroarc volcanic provinces in southern Central Andes in Mendoza,Argentina.This block is located in the Andean foothills between the orogenic front and foreland basement uplifts of late Miocene age.In order to analyze the geochronological evolution of the Quaternary volcanism in the region,several geologic and geophysical studies have been conducted.Nevertheless,the crust,where the SRB is located,has not been well characterized yet.Based on gravimetric and magnetic data,together with isostatic and elastic thickness analyses,we modeled the crustal structure of the area.Information obtained has allowed us to understand the crust where the SRB and the Payenia volcanic province are located.Bouguer anomalies indicate that the SRB presents higher densities to the North of Cerro Nevado and Moho calculations suggest depths for this block between 40 and 50 km.Determinations of elastic thickness would indicate that the crust supporting the San Rafael Block presents values of approximately 10 km,being enough to support the block loading.However,in the Payenia region,elastic thickness values are close to zero due to the regional temperature increase.

3-D MODELLING OF COMPUTER AIDED GARMENT DESIGN

This paper describes a method of the computer aided garment design,and discusses 3-D humanbody,wire frame modelling,approaches of expressing and a shading model of the 3-D garment.

3-D physical model in strong ground motion numerical simulation:A case study of Kunming basin

Seismic hazard assessment based on urban active faults can provide scientific bases for city planning and project construction, while numerical simulation of strong ground motion is an important method for seismic hazard prediction and assessment. A 3-D physical model in conformity with real strata configuration of （mainly） the Quaternary is a prerequisite to ensure the reliability of the simulation results. In this paper, we give a detailed account of the technical scheme and process for creating a 3-D physical model in Kunming basin. The data used are synthe- sized from seismogeological data, borehole data, topographic data, digital elevation mode （DEM） data, seismic exploration results and wave velocity measurements. Strafigraphic division is based mainly on shear wave velocity, with strata sequence taken into consideration. The model construction is finally accomplished with ArcGIS and many relevant programming techniques via layer-by-layer stacking （in depth direction） of the adjacent medium interfaces （meshes）. Meanwhile, a database of 3-D physical models is set up, which provides model data and parameters for strong ground motion simulation. Some processing methods and significant issues are also addressed in the paper in accordance with different types of exploration and experimental data.

Research of the ATR system based on the 3-D models and L-M BP neural network

Automatic target recognition （ATR） is an important issue for military applications, the topic of the ATR system belongs to the field of pattern recognition and classification. In the paper, we present an approach for building an ATR system with improved artificial neural network to recog- nize and classify the typical targets in the battle field. The invariant features of Hu invariant moments and roundness were selected to be the inputs of the neural network because they have the invari- ances of rotation, translation and scaling. The pictures of the targets are generated by the 3-D mod- els to improve the recognition rate because it is necessary to provide enough pictures for training the artificial neural network. The simulations prove that the approach can be implement ed in the ATR system and it has a high recognition rate and can be applied in real time.

Statistical Model of the 3-D Braided Composites Strength

Based on the statistical model for the tensile statistical strength of unidirectional composite materials and the stress analysis of 3-D braided composites, a new method is proposed to calculate the tensile statistical strength of the 3-D braided compos- ites. With this method, the strength of 3-D braided composites can be calculated with very large accuracy, and the statistical parameters of 3-D braided composites can be determined. The numerical result shows that the tensile statistical strength of 3-D braided composites can be predicted using this method.

Some issues on 2.5-D transient electromagnetic forward

Numerical calculation for two integral transforms in 2.5-D transient electromagnetic forward is a difficult and key task, namely, the inverse Fourier transform and the inverse Laplace transform. Some effective algorithms for them were described. Based on the known algorithms in DC resistivity on wave-number distribution and selection, we proposed a principle on how to choose the least wave-number concerning the central-loop transient electromagnetic method. First, observe the behavior of transformation function curve with regard to wave-number in Fourier domain. In the light of its asymptote, ascertain the coverage scope of wave-number. Compared with analytic solution, the least wave-number in Fourier domain can be derived. Furthermore, the Laplace numerical inversion algorithm which needs only a few Laplace variables in pure real domain was also introduced here. The procedure was applied to forward modeling on transient electromagnetic field of a vertical magnetic dipole over uniform half-space to demonstrate them effectiveness and general applicability.