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.

The Forward and Inverse Problem Based on Magneto-Acoustic Tomography with Current Injection

The Magneto-acoustic Tomography with Current Injection (MAT-CI) is a new biological electrical impedance imaging technique that combines Electrical Impedance Tomography (EIT) with Ultrasonic Imaging (UI), which possesses the non-invasive and high-contrast of the EIT and the high-resolution of the UI. The MAT-CI is expected to acquire high quality image and embraces a wide application. Its principle is to put the conductive sample in the Static Magnetic Field(SMF) and inject a time-varying current, during which the SMF and the current interact and generate the Lorentz Force that inspire ultrasonic signal received by the ultrasonic transducers positioned around the sample. And then according to related reconstruction algorithm and ultrasonic signal, electrical conductivity image is obtained. In this paper, a forward problem mathematical model of the MAT-CI has been set up to deduce the theoretical equation of the electromagnetic field and solve the sound source distribution by Green’s function. Secondly, a sound field restoration by Wiener filtering and reconstruction of current density by time-rotating method have deduced the Laplace’s equation that caters to the current density to further acquire the electrical conductivity distribution image of the sample through iteration method. In the end, double-loop coils experiments have been conducted to verify its feasibility.

2D numerical study on the effect of conductor between the transmitter and survey area in CSEM exploration

In CSEM exploration, the receivers are generally located about three to five times the skin depth from the transmitter. In this paper, we study the effect of a conductor between the transmitter and the survey area on the target conductor response using forward modeling and inversion. The 2D forward finite element calculations show that the conductor mainly affects the response at middle and low frequencies. The lower the resistivity and the larger the conductor, the larger the effect and the effect increases with decreasing frequency. The inversion results indicate that the calculated position of the target body can move towards the source, leading to an incorrect interpretation without considering the conductor. In order to reduce the effect of a conductor between the source and the survey area, CSEM acquisition should be conducted in three dimensions using multiple sources and 3D inversion should be used during interpretation.

Three-dimensional forward modeling and inversion of borehole-to-surface electrical imaging with different power sources

Borehole-to-surface electrical imaging （BSEI） uses a line source and a point source to generate a stable electric field in the ground. In order to study the surface potential of anomalies, three-dimensional forward modeling of point and line sources was conducted by using the finite-difference method and the incomplete Cholesky conjugate gradient （ICCG） method. Then, the damping least square method was used in the 3D inversion of the formation resistivity data. Several geological models were considered in the forward modeling and inversion. The forward modeling results suggest that the potentials generated by the two sources have different surface signatures. The inversion data suggest that the low- resistivity anomaly is outlined better than the high-resistivity anomaly. Moreover, when the point source is under the anomaly, the resistivity anomaly boundaries are better outlined than when using a line source.

AN INVERSE DESIGN OF ELECTROSTATIC FOCUSING FIELD FOR ELECTROSTATIC AND MAGNETIC IMAGING

An inverse design of electrostatic focusing field for electrostatic and magneticimaging is investigated.Using the potential superimposition theorem of electrostatic field inmulti-electrode system,a mathematical model has been developed and an optimization methodhas been introduced into computation for designing the electrostatic focusing field of the imagingsystem.

Real-time forward modeling and inversion of logging-while-drilling electromagnetic measurements in horizontal wells

Based on the pseudo-analytical equation of electromagnetic log for layered formation,an optimal boundary match method is proposed to adaptively truncate the encountered formation structures.An efficient integral method is put forward to significantly accelerate the convergence of Sommerfeld integral.By asymptotically approximating and subtracting the first reflection/transmission waves from the scattered field,the new Sommerfeld integral method has addressed difficulties encountered by the traditional digital filtering method,such as low computational precision and limited operating range,and realized the acceleration of the computation speed of logging-while-drilling electromagnetic measurements(LWD EM).By making use of the priori information from the offset/pilot wells and interactively adjusting the formation model,the optimum initial guesses of the inversion model is determined in order to predict the nearby formation boundaries.The gradient optimization algorithm is developed and an interactive inversion system for the LWD EM data from the horizontal wells is established.The inverted results of field data demonstrated that the real-time interactive inversion method is capable of providing the accurate boundaries of layers around the wellbore from the LWD EM,and it will benefit the wellbore trajectory optimization and reservoir interpretation.

The “Inverse Hall-Petch” effect on the impact response of single crystal copper

Based on the available experimental and compu- tational capabilities, a phenomenological approach has been proposed to formulate a hypersurface in both spatial and temporal domains to predict combined specimen size and load- ing rate effects on the material properties [ 1-2]. A systematic investigation is being performed to understand the combined size, rate and thermal effects on the properties and deformation patterns of representative materials with different nanos- tructures and under various types of loading conditions [3- 16]. The recent study on the single crystal copper response to impact loading has revealed the size-dependence of the Hugoniot curve. In this paper, the ＂inverse Hall-Petch＂ behavior as observed in the impact response of single crystal copper, which has not been reported in the open literature, is investigated by performing molecular dynamics simulations of the response of copper nanobeam targets subjected to impacts by copper nanobeam flyers with different impact velocities. It appears from the preliminary results that the ＂inverse Hall-Petch＂ behavior in single crystal copper is mainly due to the formation and evolution of disordered atoms and the interaction between ordered and disordered atoms, as compared with the physics behind the ＂inverse Hall-Petch＂behavior as observed in nanocrystalline materials

Effi cient modeling of the gravity anomaly caused by a sedimentary basin with lateral variable density contrast and its application in basement relief estimation

The forward calculation of gravity anomalies is a non-negligible aspect contributing to the time consumption of the entire process of basement relief estimation.In this study,we develop a fast hybrid computing scheme to compute the gravity anomaly of a basement.We use the vertical prism source equation in a given region R centered at a certain gravity observation point and the vertical line source equation outside R to derive the gravity anomaly.We observe that the computation with the vertical line source equation is much faster than that of the vertical prism source equation,but the former is slightly inaccurate.Therefore,our method is highly effi cient and able to avoid the errors caused by the low accuracy of the vertical line source equation near the observation point.We then derive the general principle of choosing the size of R via a series of prism model tests.Our tests on the gravity anomaly over the Los Angeles Basin confirm the correctness of our proposed forward strategy.We modify Bott’s method with an accelerating factor to expedite the inversion procedure and presume that the density contrast between the sediments and the basement in a sedimentary basin varies laterally and can be obtained using the equivalent equation.Synthetic data and real data applications in the Weihe Basin illustrate that our proposed method can accurately and effi ciently estimate the basement relief of sedimentary basins.

Performance of Gas-Steam Combined Cycle Cogeneration Units Influenced by Heating Network Terminal Steam Parameters

The determination of source-side extracted heating parameters is of great significance to the economic operation of cogeneration systems.This paper investigated the coupling performance of a cogeneration heating and power system multidimensionally based on the operating characteristics of the cogeneration units,the hydraulic and thermodynamic characteristics of the heating network,and the energy loads.Taking a steam network supported by a gas-steam combined cycle cogeneration system as the research case,the interaction effect among the source-side prime movers,the heating networks,and the terminal demand thermal parameters were investigated based on the designed values,the plant testing data,and the validated simulation.The operating maps of the gas-steam combined cycle cogeneration units were obtained using THERMOFLEX,and the minimum source-side steam parameters of the steam network were solved using an inverse solution procedure based on the hydro-thermodynamic coupling model.The cogeneration operating maps indicate that the available operating domain considerably narrows with the rise of the extraction steam pressure and flow rate.The heating network inverse solution demonstrates that the source-side steam pressure and temperature can be optimized from the originally designed 1.11 MPa and 238.8°C to 1.074 MPa and 191.15°C,respectively.Under the operating strategy with the minimum source-side heating parameters,the power peak regulation depth remarkably increases to 18.30%whereas the comprehensive thermal efficiency decreases.The operation under the minimum source-side heating steam parameters can be superior to the originally designed one in the economy at a higher price of the heating steam.At a fuel price of$0.38/kg and the power to fuel price of 0.18 kg/(kW·h),the critical price ratio of heating steam to fuel is 119.1 kg/t.The influence of the power-fuel price ratio on the economic deviation appears relatively weak.

Characteristics of Zhangsanying-Tongshanzi aeromagnetic anomaly zone and prospecting potential of iron deposits in northern Hebei,China

Based on the latest high-precision aeromagnetic data,an aeromagnetic anomaly zone is identified at Zhangsanying--Tongshanzi in northern Hebei Province.By the potential field conversion processing,including the reduction to the pole,vertical derivative,upward continuation and residual anomaly,the authors analyzed the characteristics of three typical aeromagnetic anomalies in Zhangsanying--Tongshanzi aeromagnetic anomaly zone and their geological origin.The methods include the forward and inversion methods,such as 2.5D optimization fitting and Euler deconvolution.Moreover,combined with the geological outcrop,known iron deposits,ground magnetic survey and verification,the authors studied the relationship between the aeromagnetic anomalies and iron deposits.The result shows that the Zhangsanying--Tongshanzi aeromagnetic anomaly zone is composed of 10 large magnetic anomalies with high amplitude and clear boundary.The aeromagnetic anomalies are comparable and intrinsically related to the ground magnetic anomalies and IP anomalies,indicating that the anomalies are caused by magnetite deposits.It has good magnetite prospecting potential in the Zhangsanying--Tongshanzi aeromagnetic anomaly zone.

Modified Cooperative Subchannel Allocation Algorithms and PSO Based Power Allocation for an Alamouti Decode and Forward Relaying Protocol in Multiuser OFDMA Systems

The present work is a discussion on the performance analysis of Modified Cooperative Subchannel Allocation (CSA) Algorithms which is used in Alamouti Decoded and Forward (Alamouti DF) Relaying Protocol for wireless multi-user Orthogonal Frequency Division Multiplexing Access (OFDMA) systems. In addition, the performance of approximate Symbol Error Rate (SER) for the Alamouti DF Relaying Protocol with the Cooperative Maximum Ratio Combining Technique (C-MRC) is analyzed and compared with SER upper bound. The approximate SER is asymptotically tight bound at higher Signal-to-Noise Ratio (SNR). From the asymptotic tight bound approximate SER, Particle Swarm Optimization (PSO) based Power Allocation (PA) is determined for the Alamouti DF Relaying Protocol. The simulation results suggested that the Modified Throughput based Subchannel Allocation Algorithm achieved an improved throughput of 6% to 33% compared to that of existing cooperative diversity protocol. Further, the Modified Fairness based Subchannel Allocation Algorithm rendered fairness of 7.2% to 17% among the multiuser against the existing cooperative diversity protocol.

ON THE INVERTIBILITY OF HILBERT SPACE IDEMPOTENTS

This note is to present some results on the group invertibility of linear combina- tions of idempotents when the difference of two idempotents is group invertible.

Application of Reconstruction and Optimization Algorithms in Optical Tomography

Optical tomography is a non-invasive technique that uses visible or near infrared radiation to analyze biological tissues. Researchers take immense attention towards advancement in optical tomography because of its low cost and an advantage of providing anatomical information. Based on the information of optical characteristics, forward and inverse problem of tomography are solved. In this research, finite element method is employed for forward problem and gradient-based optimization algorithm is developed for inverse problem of optical tomography. It is found from simulations that information about imaging is processed more distinctly and in less computational time. Normal and abnormal conditions in imaging are readily distinguished.?Simulations are carried out in Matlab. Different scenarios are developed and are simulated to validate the performance of reconstruction and optimization algorithms in optical tomography.

Broadband diffuse optical spectroscopy of two-layered scattering media containing oxyhemoglobin,deoxyhemoglobin,water,and lipids

We investigated the relationship between chromophore concentrations in two-layered scattering media and the apparent chromophore concentrations measured with broadband optical spectroscopy in conjunction with commonly used homogeneous medium inverse models.We used diffusion theory to generate optical data from a two-layered distribution of relevant tissue absorbers,namely,oxyhemoglobin,deoxyhemoglobin,water,and lipids,with a top-layer thickness in the range 1–15 mm.The generated data consisted of broadband continuous-wave(CW)diffuse reflectance in the wavelength range 650–1024 nm,and frequency-domain(FD)diffuse reflectance at 690 and 830 nm;two source-detector distances of 25 and 35mm were used to simulate a dual-slope technique.The data were inverted using diffusion theory for a semi-infinite homogeneous medium to generate reduced scattering coeffcients at 690 and 830nm(from FD data)and effective absorption spectra in the range 650–1024nm(from CW data).The absorption spectra were then converted into effective total concentration and oxygen saturation of hemoglobin,as well as water and lipid concentrations.For absolute values,it was found that the effective hemoglobin parameters are typically representative of the bottom layer,whereas water and lipid represent some average of the respective concentrations in the two layers.For concentration changes,lipid showed a significant cross-talk with other absorber concentrations,thus indicating that lipid dynamics obtained in these conditions may not be reliable.These systematic simulations of broadband spectroscopy of two-layered media provide guidance on how to interpret effective optical properties measured with similar instrumental setups under the assumption of medium homogeneity.

Multi-component joint inversion of gravity gradient based on fast forward calculation

With the development of gravity gradient full tensor measurement technique,three-dimensional( 3D) inversion based on gravity gradient tensor can provide more accurate information. But the forward calculation of 3D full tensor sensitivity matrix is very time-consuming,which restricts its development and application.According to the symmetry of the kernel function,the authors reconstruct the underground source of geological body to avoid repeat computation of the same value,and work out the corresponding relationship between the response of geological body to the observation point and the response of reconstructed geological body to the observation point. According to the relationship,rapid calculation of full tensor gravity sensitivity matrix can be achieved. The model calculation shows that this method can increase the speed of 30-45 times compared with the traditional calculation method. The sensitivity matrix is applied to the multi-component inversion of gravity gradient. The application of this method on the measured data provides the basis for the promotion of the method.

Simulation research and application on response characteristics of detecting water-filled goaf by transient electromagnetic method

Water inrush disasters poses a great threat to the safe exploitation of coal resources.To solve this problem,the transient electromagnetic method(TEM)was proposed to accurately detect the water accumulation in the goaf.The electromagnetic response characteristics of diferent water-flled goaves were studied by electromagnetic feld theory,numerical simulation and feld verifcation.Through the models of 100%water accumulation,50%water accumulation,0%water accumulation,100%water accumulation with collapsed rock,50%water accumulation with collapsed rock and 0%water accumulation with collapsed rock goaf,the characteristics of induced voltage attenuation curves were studied.Meanwhile,the relationship between the attenuation voltage value and area of the transmitting coil,the depth of the goaf,the background resistivity,and the delay time were also simulated.The results illustrate that the attenuation curve of induced voltage presented a regular exponential decay form in the 0%water accumulation model but existed abnormal exaltation for voltage in water-flled model.Through the linear ftting curve,it can be seen that the abnormal intensity of the induced voltage becomes stronger as the distance between the measuring point and the center of the target decrement.Moreover,the abnormal amplitude of the induced voltage increases with the rise of the water accumulation and collapsed rock will weakly reduce the low-resistivity anomalous efect on the water-accumulated goaf.In addition,the response value of the attenuation voltage increased as the area of the transmitting coil increases,but decreased with increasing delay time and increasing background resistivity and depth of the target body.The feld detection results of the Majiliang coal mine also confrmed the theoretical analysis and the numerical simulation.

Breakthrough and significance of technology on internal multiple recognition and suppression: A case study of Ordovician Dengying Formation in Central Sichuan Basin, SW China

Aimed at the serious mismatch between the synthetic seismogram and the real data of the Sinian Dengying Formation and the Lower Cambrian Qiongzhusi Formation in the GS1 well area, Sichuan Basin, four aspects of internal multiples identification and suppression were studied. Firstly, a forward modeling method of internal multiple based on reflectivity method was developed. Through eight means such as post-stack and pre-stack forward modeling of internal multiple, and combined with VSP data, it was demonstrated that well-seismic mismatching in this area is mainly caused by the internal multiples. Secondly, the simulation results of internal multiple forward modeling using the stripping method combined with the internal multiple periodicity analysis showed that four groups of overburden layers with velocity inversion were the main sources of the internal multiples. Thirdly, by identifying internal multiples accurately and using suppression technology based on pattern recognition, an effective and replicable suppression scheme suitable for these formations was established, overcoming the difficulty of the small speed difference between internal multiple and primary reflection wave which makes the current methods ineffective. Fourthly, an evaluation index of internal multiple intensity was proposed, and the internal multiple intensity distribution diagram of the fourth member of Dengying Formation(Deng-4 Member) in Gaoshiti-Moxi area was compiled. This scheme greatly improved the well-seismic match, and the strata sedimentary features are clearer on the new seismic profiles with higher lateral resolution, with which smaller faults and geological anomalies can be identified and a series of the bead reflections in the Dengying Formation are first discovered. The coincidence rate of reservoir prediction based on seismic waveform classification was increased from 60% to 90%, and that of hydrocarbon detection based on dual phase medium theory was increased from 70% to 100%.

ANALYSIS OF SURFACE CRACK ON FORWARD EXTRUDING BAR DURING AXISYMMETRIC CUP-BAR COMBINED EXTRUSION PROCESS

In this paper, the kinematically admissible velocity field with surface crack on forward extruding bar is put forward during the axisymmetric cup-bar combined extrusion process, in accordance with the results of model experiments.On the basis of velocity field, the necessary condition for surface crack formation on the forward extruding bar is derived, with the help of upper bound theorem and the minimum energy principle. Meanwhile, the relationships between surface crack formation and combination of reduction in area for the part of forward and backward extursions relative residual thickness of billet (T/R0),frictional factor (m) or relative land length of ram and chamber are calculated during the extrusion process. Therefore, whether the surface crack on forward exturding bar occurs can be predicted before extruding the lower-plasticity metals for axisymmetric cup-bar combined extrusion process.The analytical results agree very well with experimental results of aluminium alloy LY12 (ASTM 2024) and LC4 (ASTM 7075).

基于2.5D有限元的高密度电法不同装置勘探效果研究

高密度电法是一种适用性广泛的电阻率勘探方法,拥有多种电极排列装置选择。不同装置由于电极分布方式的差异,在针对不同地质目标和测区环境时往往表现出明显不同的探测能力。为探究各排列装置的特点及其适用情况,推导了2.5D电位的变分问题,采用Delaunay三角化算法实现了非结构化网格剖分,进而实现了有限元正演模拟。结合实践,对常见地质情况进行建模,分别使用温纳α、温纳β、施伦贝谢尔、偶极-偶极4种装置开展正反演对比研究。研究发现:(1)在探测未知区域单个异常体时,温纳β、施伦贝谢尔2种排列的效果更好;(2)在探测相邻较近的多个异常体时,应优先采用具有更高水平分辨率的施伦贝谢尔和偶极-偶极排列;(3)对于低阻破碎带的勘探,温纳β和偶极-偶极排列更为适用;(4)当地层具有较好分层界限时,4种排列均可体现良好的分层效果。因此,高密度电法勘探实践中应综合考虑不同排列装置的特点,根据具体情况选用合适的多种排列方式进行测量,并对采集数据进行综合对比分析,以获得更为准确的解释结果。

钻孔瞬变电磁法扫描探测RCQPSO-LMO组合算法2.5D反演

利用钻孔进行超前探测地质构造及含水体是地下开挖工程中的常规手段,如何利用这些钻孔进行钻孔瞬变电磁法扫描探测,从而实现钻孔孔壁外围地质异常体的精细探测,对实现地下工程地质透明化具有重要的指导意义.本文提出钻孔瞬变电磁法扫描探测2.5D反演的数据解译方法,首先针对随机性反演算法时效性低,易陷入局部最优解,而确定性反演算法依赖初始模型的问题,提出了组合策略的量子粒子群优化算法用来随机搜索最优初始模型.在此基础上,利用Levenberg-Marquarat方法求解Occam反演的目标函数,形成了RCQPSO-LMO组合算法进行2.5D反演,通过对比组合算法和单一算法,验证了组合算法具有更精确的反演结果.其次结合屏蔽条件下扫描探测,对比分析了有无屏蔽的2.5D反演结果,通过设定屏蔽系数对非探测方向信号进行部分压制,可以较好地解决钻孔径向扫描探测中对非探测方向信号部分屏蔽下的反演及成像.最后建立三组理论模型进行组合算法2.5D反演,结果表明:组合算法反演结果与理论模型的一致性较好,对低阻异常体的反演精度较高,验证了组合算法对钻孔孔壁外围低阻异常体具有较高的反演精度和分辨能力.