Three-dimensional gravity inversion based on sparse recovery iteration using approximate zero norm

This research proposes a novel three-dimensional gravity inversion based on sparse recovery in compress sensing. Zero norm is selected as the objective function, which is then iteratively solved by the approximate zero norm solution. The inversion approach mainly employs forward modeling; a depth weight function is introduced into the objective function of the zero norms. Sparse inversion results are obtained by the corresponding optimal mathematical method. To achieve the practical geophysical and geological significance of the results, penalty function is applied to constrain the density values. Results obtained by proposed provide clear boundary depth and density contrast distribution information. The method's accuracy, validity, and reliability are verified by comparing its results with those of synthetic models. To further explain its reliability, a practical gravity data is obtained for a region in Texas, USA is applied. Inversion results for this region are compared with those of previous studies, including a research of logging data in the same area. The depth of salt dome obtained by the inversion method is 4.2 km, which is in good agreement with the 4.4 km value from the logging data. From this, the practicality of the inversion method is also validated.

Three-dimensional frequency-domain full waveform inversion based on the nearly-analytic discrete method

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.

Electrical Structure of Wulingshan and Middle Jiangnan Orogen by Three-Dimensional Magnetotelluric Data Inversion

Funded by The National Key Research and Development Program of China,China Deep Exploration(Sinoprobe)and The China Geological Suvery Project on 2009–2019,a large scale magnetotelluric sounding(MT)survey grid(Fig.1)has covered whole south China.

Three-dimensional gravity inversion based on 3D U-Net++

The gravity inversion is to restore genetic density distribution of the underground target to be explored for explaining the internal structure and distribution of the Earth.In this paper,we propose a new 3D gravity inversion method based on 3D U-Net++.Compared with two-dimensional gravity inversion,three-dimensional(3D)gravity inversion can more precisely describe the density distribution of underground space.However,conventional 3D gravity inversion method input is two-dimensional,the input and output of the network proposed in our method are three-dimensional.In the training stage,we design a large number of diversifi ed simulation model-data pairs by using the random walk method to improve the generalization ability of the network.In the test phase,we verify the network performance by using the model-data pairs generated by the simulation.To further illustrate the eff ectiveness of the algorithm,we apply the method to the inversion of the San Nicolas mining area,and the inversion results are basically consistent with the borehole measurement results.Moreover,the results of the 3D U-Net++inversion and the 3D U-Net inversion are compared.The density models of the 3D U-Net++inversion have higher resolution,more concentrated inversion results,and a clearer boundary of the density model.

Inversion of source mechanism of 1989 Loma Prieta earthquake by three-dimensional FEM Green's function

Based on three-dimensional joint finite element, this paper discusses the theory and methodology about inversionof geodetic data. The FEM and inversion formula is given in detail; also a related code is developed. By use of theGreen's function about 3-D FEM, we invert geodetic measurementS of coseismic deformation of the 1989 Ms=7. 1Loma Prieta earthquake to datermine itS source mechanism. The result indicates that the slip on the fault plane isvery heterogeneous. The maximum slip and shear stress are located about 10 kin to northwest of the eathquakesource, the stress drop is about more than 1 MPa.

Concealed porphyry delineation based on nonlinear three-dimensional densitydifference inversion: An example in the Beiya mine area, Western Yunnan, China

Intermediate acid-complex rock masses with low-density characteristics are the most important prospecting sign in the Beiya area, of western Yunnan province, and provide a physical basis for good gravity exploration. It is usually difficult to obtaining solutions in connection with actual geological situations due to the ambiguity of the conventional gravity-processing results and lack of deep constraints. Thus, the three-dimensional (3D) inversion technology is considered as the main channel for reducing the number of solutions and improving the vertical resolution at the current stage. The current study starts from a model test and performs nonlinear 3D density-difference inversion called “model likelihood exploration”, which performs 3D inversion imaging and inversion of the known model while considering the topographic effects. The inversion results are highly consistent with those of the known models. Simultaneously, we consider the Beiya gold mine in Yunnan as an example. The nonlinear 3D densitydifference inversion technology, which is restricted by geological information, is explored to obtain the 3D density body structure below 5 km in the mine area, and the 3D structure of the deep and concealed rock masses are obtained using the density constraints of the intermediate-acid-complex rock masses. The results are well consistent with the surface geological masses and drilling-controlled deep geological masses. The model test and examples both show that the 3D density-difference nonlinear inversion technology can reduce inversion ambiguity, improve resolution, optimize the inversion results, and realize “transparency” in deeply concealed rock masses in ore-concentrated areas,which is useful in guiding the deep ore prospecting.

Creation of targeted inversion mutations in plants using an RNA-guided endonuclease

Inversions are DNA rearrangements that are essential for plant gene evolution and adaptation to environmental changes. We demonstrate the creation of targeted inversions and previously reported targeted deletion mutations via delivery of a pair of RNA-guided endonucleases(RGENs) of CRISPR/Cas9. The efficiencies of the targeted inversions were2.6% and 2.2% in the Arabidopsis FLOWERING TIME(At FT) and TERMINAL FLOWER 1(At TFL1)loci, respectively. Thus, we successfully established an approach that can potentially be used to introduce targeted DNA inversions of interest for functional studies and crop improvement.

Three-dimensional suboptimal guidance law for fly-over and shoot-down smart ammunition via virtual target

An optimal burst height is required for the fly-over and shoot-down smart ammunition with an EFP warhead at the instant of explosion which brings a special requirement to the miss distance of the terminal guidance law. In this paper, a guidance law based on the virtual target scheme is proposed. First, the practical pursuit-evasion issue between the ammunition and the target with specific miss distance is transformed into a virtuai pursuit-evasion problem with zero miss distance. Secondly, a complete three-dimensional pursuit-evasion kinematics model is established without any simplifications. And then, a suboptimal guidance law is designed based on the θ-D method which has constraints of the elevation and azimuth angular velocity of the virtual line of sight （LOS）. Finally, in order to verify the performance of the proposed guidance law, three test cases are conducted. Numericai results show that under the proposed terminal guidance law, the smart ammunition not only can fly above the target with an optimal burst height but also have a smaller normal acceleration on the terminal trajectory.

The distribution of deep source rocks in the GS Sag:joint MT-gravity modeling and constrained inversion

The coal-bearing strata of the deep Upper Paleozoic in the GS Sag have high hydrocarbon potential. Because of the absence of seismic data, we use electromagnetic （MT） and gravity data jointly to delineate the distribution of deep targets based on well logging and geological data. First, a preliminary geological model is established by using three-dimensional （3D） MT inversion results. Second, using the formation density and gravity anomalies, the preliminary geological model is modified by interactive inversion of the gravity data. Then, we conduct MT-constrained inversion based on the modified model to obtain an optimal geological model until the deviations at all stations are minimized. Finally, the geological model and a seismic profile in the middle of the sag is analysed. We determine that the deep reflections of the seismic profile correspond to the Upper Paleozoic that reaches thickness up to 800 m. The processing of field data suggests that the joint MT-gravity modeling and constrained inversion can reduce the multiple solutions for single geophysical data and thus improve the recognition of deep formations. The MT-constrained inversion is consistent with the geological features in the seismic section. This suggests that the joint MT and gravity modeling and constrained inversion can be used to delineate deep targets in similar basins.

3D parallel inversion of time-domain airborne EM data

To improve the inversion accuracy of time-domain airborne electromagnetic data, we propose a parallel 3D inversion algorithm for airborne EM data based on the direct Gauss-Newton optimization. Forward modeling is performed in the frequency domain based on the scattered secondary electrical field. Then, the inverse Fourier transform and convolution of the transmitting waveform are used to calculate the EM responses and the sensitivity matrix in the time domain for arbitrary transmitting waves. To optimize the computational time and memory requirements, we use the EM ＂footprint＂ concept to reduce the model size and obtain the sparse sensitivity matrix. To improve the 3D inversion, we use the OpenMP library and parallel computing. We test the proposed 3D parallel inversion code using two synthetic datasets and a field dataset. The time-domain airborne EM inversion results suggest that the proposed algorithm is effective, efficient, and practical.

Micro-motion effect in inverse synthetic aperture radar imaging of ballistic mid-course targets

For ballistic mid-course targets,in addition to constant orbital motion,the target or any structure on the target undergoes micro-motion dynamics,such as spin,precession and tumbling.The micro-motion characteristics of the ballistic mid-course targets were discussed.The target motion model and inverse synthetic aperture radar(ISAR) imaging model for this kind of targets were built.Then,the influence of micro-motion on ISAR imaging based on the established imaging model was presented.The computer simulation to get mid-course target echoes from static darkroom electromagnetic scattering data based on the established target motion model was realized.The imaging results of computer simulation show the validity of ISAR imaging analysis for micro-motion targets.

Three-dimensional coordinates test method with uncertain projectile proximity explosion position based on dynamic seven photoelectric detection screen

To objectively obtain the three-dimensional coordinates of the projectile fuze proximity explosion when projectile intersects the head of missile target, we propose a dynamic seven photoelectric detection screen test method, which is made up of six plane detection screens and a flash photoelectric dynamic detection screen. The three-dimensional coordinates calculation model of the projectile proximity explosion position based on seven plane detection screens with dynamic characteristics is established.According to the relation of the dynamic seven photoelectric detection screen planes and the time values,the analytical function of the projectile proximity explosion position parameters under non-linear motion is derived. The projectile signal filtering method based on discrete wavelet transform is explored in this work. Additionally, the projectile signal recognition algorithm using an improved particle swarm is proposed. Based on the characteristics of the time duration and the signal peak error for the projectile passing through the detection screen, the signals attribution of the same projectile passing through six detection screens are analyzed for obtaining precise time values of the same projectile passing through the detection screens. On the basis of the projectile fuze proximity explosion test, the linear motion model and the proposed non-linear motion model are used to calculate and compare the same group of projectiles proximity explosion position parameters. The comparison of test results verifies that the proposed test method and calculation model in this work accurately obtain the actual projectile proximity explosion position parameters.

Three-dimensional lithospheric density distribution of China and surrounding regions

In this paper, we analyze lithospheric density distribution of China and surrounding regions on the basis of 30＇ × 30＇ gravity data and 1°×1 °P-wave velodty data, Firstly, we used the empirical equation be- tween the density and the P-wave velocity difference as the base of the initial model of the Asian lithospheric density. Secondly, we calculated the gravity anomaly, caused by the Moho discontinuity and the sedimentary layer discontinuity, by the Parker formula. Thirdly, the gravity anomaly of the spherical harmonics with 2 40 order for the anomalous body below the lithosphere is calculated based on the model of EGM96. Finally, by using Algebra Reconstruction Techniques （ART）, the inversion of 30＇ ~ 30＇ residual lithospheric Bouguer gravity anomaly caused by the lithosphere yields a rather detailed struc- tural model. The results show that the lithospheric density distribution of China and surrounding regions has a certain connection with the tectonic structure. The density is relatively high in the Philippine Sea plate, Japan Sea, the Indian plate, the Kazakhstan shield and the Western Siberia plain, whereas the Tibetan Plateau has low-density characteristics. The minimum value of density lies in the north of Philippines, in the Taiwan province and in the Ryukyu island arc.

A modified OMP method for multi-orbit three dimensional ISAR imaging of the space target

The conventional two dimensional(2D)inverse synthetic aperture radar(ISAR)imaging fails to provide the targets'three dimensional(3D)information.In this paper,a 3D ISAR imaging method for the space target is proposed based on mutliorbit observation data and an improved orthogonal matching pursuit(OMP)algorithm.Firstly,the 3D scattered field data is converted into a set of 2D matrix by stacking slices of the 3D data along the elevation direction dimension.Then,an improved OMP algorithm is applied to recover the space target's amplitude information via the 2D matrix data.Finally,scattering centers can be reconstructed with specific three dimensional locations.Numerical simulations are provided to demonstrate the effectiveness and superiority of the proposed 3D imaging method.

扩展目标跟踪Student’s t逆Wishart平滑算法

脉冲干扰和离群量测信息等因素通常会导致异常的厚尾噪声,这使得以高斯假设为前提的扩展目标跟踪(ETT)估计器的性能急剧降低,针对该问题该文提出一种基于扩展目标随机矩阵模型(RMM)的Student’s t逆Wishart平滑(StIWS)算法。首先,将目标的运动状态以及过程噪声和量测噪声建模为Student’s t分布以表征异常噪声对扩展目标概率分布的影响,将目标扩展状态建模为服从逆Wishart分布的随机矩阵。然后,在Student’s t贝叶斯平滑框架下,详细推导了能在扩展目标的多重特征动态演变的过程中有效估计目标状态的StIWS算法。最后,通过扩展目标跟踪的仿真实验结果和真实场景实验结果验证了所提算法的有效性。

微波视觉与SAR图像智能解译

高分辨率雷达成像技术和人工智能、大数据技术的快速发展,有力促进了雷达图像智能解译技术的进步。由于雷达传感器本身的特殊性和电磁散射成像物理的复杂性,雷达图像的解译缺乏光学图像的直观性,准确迅速识别分类的需求对雷达图像解译提出了迫切的挑战。在借鉴人脑光视觉感知机理和计算机视觉图像处理相关技术基础上,进一步融合电磁散射物理规律及其雷达成像机理,我们提出发展微波域雷达图像解译的“微波视觉”的新交叉领域研究。该文介绍微波视觉的概念与内涵,提出微波视觉认知模型,阐述其基础理论问题与技术路线,最后介绍了作者团队在相关问题上的初步研究进展。

基于多特征重构的三维目标反演算法

为了解决基于深度学习的三维反演方法中存在的内存占用大、训练耗时久的问题,提出基于多特征重构的三维目标反演算法.通过特征分解提取目标的水平区域、中心深度、垂直厚度和剩余密度4类特征,实现对三维模型的压缩,降低内存占用.设计多特征重构反演网络(MRNet),通过不同的Decoder实现对目标4类特征的预测,利用4类特征重构三维模型,实现对三维目标的反演.在网络输入端引入梯度联合实现对目标边界信息的增强.在跨层连接处引入CA注意力机制,实现对Decoder预测功能的分化,优化反演效果.模拟实验结果显示,MRNet的局部相对准确度相对于3D U-Net提升了30%以上,达到88.91%,每轮训练时间仅为3D U-Net的1/13.将MRNet应用于Vinton盐丘地区,较准确地得到了盖岩的分布情况,验证了MRNet具备一定的泛化性.

一种联合TDCM和PFA的空间目标大转角ISAR成像方法

空间目标态势感知是当前航天技术领域的研究热点,逆合成孔径雷达(Inverse Synthetic Aperture Radar, ISAR)是空间目标感知的一种有效手段。为了得到利于目标鉴别的高分辨二维ISAR图像,需要通过增大积累转角来提高方位向分辨率。但是,随着转角的增大,平动与转动之间的耦合增大,使得传统的先平动补偿、后利用转动成像的两步式成像方法不再适用。本文提出了一种结合时域信号共轭相乘(Time Domain Conjugate Multiplication, TDCM)处理与极坐标格式算法(Polar Format Algorithm, PFA)的成像算法,首先利用雷达回波的二维相干性,通过TDCM处理进行平动粗补偿,然后借助PFA的插值操作校正大转角引起的散射点越分辨单元走动(Motion Through Resolution Cell, MTRC)并估计目标相对转动参数。通过迭代处理,持续对平动和转动参数进行优化,提升运动补偿效果。仿真验证了所提成像算法的有效性,与传统算法相比,成像聚焦性更好。

以全球导航卫星系统为辐射源的前向散射雷达发展综述

前向散射雷达(FSR)可获得高水平雷达截面积(RCS)的特性使其在反隐身中占据重要地位。利用全球导航卫星系统(GNSS)作为辐射源,具有全天时全天候全地域覆盖的优势,通过部署多个接收节点可构建地面/海上/空中目标监视网络。该文针对基于GNSS的FSR发展现状,从目标检测、目标参数估计、阴影逆合成孔径雷达(SISAR)成像及目标分类识别等方面对关键技术和现存问题进行概述,并从组网探测、多目标定位、布站优化和极化信息获取等方面对基于GNSS的FSR发展趋势提出展望。

基于OpenMV的五自由度机械臂控制方法研究

以五自由度机械臂为研究对象,建立D-H数学模型,探究几何法与代数法求关节角的过程,使用五次多项式插值的运动轨迹规划算法。通过建立仿真模型验证了运动学求解方法与轨迹规划算法的有效性。采用一种基于机械臂几何模型的逆运动学分析方法,设计OpenMV摄像头的图像识别与处理算法来精准定位目标物,并将目标物的位姿发送给主控系统,实现五自由度机械臂自主抓取目标物的功能。以乒乓球为抓取目标进行实验验证,结果表明:此设计能精准定位目标物,完成有效抓取,控制性能稳定、可靠。