利用不同深度场源磁异常分离的方法解译隐伏岩体

为了查明安徽省殷家汇地区隐伏岩体的赋存状况,初步圈定侵入性中酸性岩体的范围并推测有利赋矿位置,对地面实测磁异常进行化极、延拓、不同深度场源磁异常分离、垂向导数等数据处理,提取利于识别的剩余磁异常及垂向一阶导数磁异常,结合地质资料解译磁异常特征。结果表明:区内NE向剩余磁异常带为沿断裂构造侵入的脉岩引起,NE向构造为区内主要的控岩和控矿构造;推测引起局部剩余磁异常的磁性地质体为隐伏或半隐伏中酸性岩体。该研究为长江中下游铁铜金成矿带中的安庆—盱眙地质矿产调查区成矿远景区划分及找矿靶区初选提供参考。

Wireless location algorithm using digital broadcasting signals based on neural network

In order to enhance the accuracy and reliability of wireless location under non-line-of-sight （NLOS） environments,a novel neural network （NN） location approach using the digital broadcasting signals is presented. By the learning ability of the NN and the closely approximate unknown function to any degree of desired accuracy,the input-output mapping relationship between coordinates and the measurement data of time of arrival （TOA） and time difference of arrival （TDOA） is established. A real-time learning algorithm based on the extended Kalman filter （EKF） is used to train the multilayer perceptron （MLP） network by treating the linkweights of a network as the states of the nonlinear dynamic system. Since the EKF-based learning algorithm approximately gives the minimum variance estimate of the linkweights,the convergence is improved in comparison with the backwards error propagation （BP） algorithm. Numerical results illustrate thatthe proposedalgorithmcanachieve enhanced accuracy,and the performance ofthe algorithmis betterthanthat of the BP-based NN algorithm and the least squares （LS） algorithm in the NLOS environments. Moreover,this location method does not depend on a particular distribution of the NLOS error and does not need line-of-sight （ LOS ） or NLOS identification.

A Multi-level Approach for Complex Fault Isolation Based on Structured Residuals

In industrial processes,there exist faults that have complex effect on process variables.Complex and simple faults are defined according to their effect dimensions.The conventional approaches based on structured residuals cannot isolate complex faults.This paper presents a multi-level strategy for complex fault isolation.An extraction procedure is employed to reduce the complex faults to simple ones and assign them to several levels.On each level,faults are isolated by their different responses in the structured residuals.Each residual is obtained insensitive to one fault but more sensitive to others.The faults on different levels are verified to have different residual responses and will not be confused.An entire incidence matrix containing residual response characteristics of all faults is obtained,based on which faults can be isolated.The proposed method is applied in the Tennessee Eastman process example,and the effectiveness and advantage are demonstrated.

TWO-LAYERED IMAGE COMPRESSION USING 2G BANDELETS AND WAVELETS

An image can be decomposed into the structural component and the geometric texturalcomponent.Based on this idea,an efficient two-layered compressing algorithm is proposed,which uses2nd generation bandelets and wavelets.First,an original image is decomposed into the structuralcomponent and the textural component,and then these two components are compressed using waveletsand 2nd generation bandelets respectively.Numerical tests show that the proposed method worksbetter than the bandelets and JPEG2000 in some specific SAR scene.

Numerical modeling and dynamic analysis of the 2017 Xinmo landslide in Maoxian County, China

A catastrophic landslide occurred at Xinmo village in Maoxian County, Sichuan Province,China, on June 24, 2017. A 2.87×106 m3 rock mass collapsed and entrained the surface soil layer along the landslide path. Eighty-three people were killed or went missing and more than 103 houses were destroyed. In this paper, the geological conditions of the landslide are analyzed via field investigation and high-resolution imagery. The dynamic process and runout characteristics of the landslide are numerically analyzed using a depth-integrated continuum method and Mac Cormack-TVD finite difference algorithm.Computational results show that the evaluated area of the danger zone matchs well with the results of field investigation. It is worth noting that soil sprayed by the high-speed blast needs to be taken into account for such kind of large high-locality landslide. The maximum velocity is about 55 m/s, which is consistent with most cases. In addition, the potential danger zone of an unstable block is evaluated. The potential risk area evaluated by the efficient depthintegrated continuum method could play a significant role in disaster prevention and secondary hazard avoidance during rescue operations.

Borehole-GPR numerical simulation of full wave field based on convolutional perfect matched layer boundary

The absorbing boundary is the key in numerical simulation of borehole radar.Perfect match layer(PML) was chosen as the absorbing boundary in numerical simulation of GPR.But CPML(convolutional perfect match layer) approach that we have chosen has the advantage of being media independent.Beginning with the Maxwell equations in a two-dimensional structure,numerical formulas of finite-difference time-domain(FDTD) method with CPML boundary condition for transverse electric(TE) or transverse magnetic(TM) wave are presented in details.Also,there are three models for borehole-GPR simulation.By analyzing the simulation results,the features of targets in GPR are obtained,which can provide a better interpretation of real radar data.The results show that CPML is well suited for the simulation of borehole-GPR.

The study of perfectly matched layer absorbing boundaries for SH wave fields

When modeling wave propagation in infinite space, it is necessary to have stable absorbing boundaries to effectively eliminate spurious reflections from the truncation boundaries. The SH wave equations for Perfectly Matched Layers （PML） are deduced and their Crank-Nicolson scheme are presented in this paper. We use the second-, sixth-, and tenth-order finite difference and pseudo-spectral algorithms to compute the spatial derivatives. Two numerical models, a homogeneous isotropic medium and a multi-layer model with a cave, are designed to investigate how the absorbing boundary width and the algorithms determine PML effects. Numerical results show that, for PML, the low-order finite difference algorithms have fairly good absorbing effects when the absorbing boundary is thin, whereas, high-order algorithms always have good absorption when the boundary is thick. Finally, we discuss the reflection coefficient and point out its shortcomings, which is why we use the SNR to quantitatively scale the PML effects,

WLS filter for reducing atmospheric effects in spaceborne SAR tomography

A new approach was presented to eliminate the atmosphere-induced phase error utilizing only the single look complex(SLC) synthetic aperture radar(SAR) image set. This method exploited the space-invariance characteristic of phase error components contained in image pixels and estimates the phase error using the weighted least-squares(WLS) filter. Actually, this sort of method can be classified as autofocus algorithm which was generally applied in airborne SAR 2-D imaging to compensate the phase error introduced by airplane's nonideal motion. Real data processing, which is relevant to Honda center and Angel stadium of Anaheim test-sites and acquired by Envisat-ASAR during the period from June 2004 to October 2007, was carried out to evaluate this WLS estimation algorithm. Experimental results show that the phase error estimated from WLS filter is very accurate and the focusing quality along NSR dimension is improved prominently via phase correction, which verifies the practicability of this new method.

Surface-borehole electric potential survey in layered media and its applications

The method in which a source is set on the surface and electric potential is received in the borehole is called surface-borehole electric potential technique. Technique of surface-borehole electric potential was employed to study electric response of layered formation. The electric potential was obtained by solving Poisson equation with finite difference method. In the course of calculation, forward modeling wilth finite difference method was realized by adopting bandwidth non-zero storage technique and the incomplete Cholesky conjugate gradient method. The results show that method of surface-borehole can acquire anomalous electric potential corresponding tc geo-electric layers. In addition, application of appropriate mathematical operator can improve the resolution. Moreover, overburden low resistivity layers have severe influence on measuring results of surface-borehole electric potential. However, bottom low resistivity layers play a positive role in the measurement.

RESIDUAL A POSTERIORI ERROR ESTIMATE OF A NEW TWO-LEVEL METHOD FOR STEADY NAVIER-STOKES EQUATIONS

Residual-based a posteriori error estimate for conforming finite element solutions of incompressible Navier-Stokes equations, which is computed with a new two-level method that is different from Volker John, is derived. A posteriori error estimate contains additional terms in comparison to the estimate for the solution obtained by the standard finite element method. The importance of the additional terms in the error estimates is investigated by studying their asymptotic behavior. For optimal scaled meshes, these bounds are not of higher order than the convergence of discrete solution. The two-level method aims to solve the nonlinear problem on a coarse grid with less computational work, then to solve the linear problem on a fine grid, which is superior to the usual finite element method solving a similar nonlinear problem on the fine grid.

An efficient hybrid model for fractured reservoirs

In this study, an efficient hybrid model is proposed to simulate the fluid flow in the reservoirs with multi-scale fractures, which cannot be easily modeled by neither the continuum models nor the discrete fracture models. In the proposed method, the small fractures are modeled by using an improved Multiple Sub-Region method, which can capture the strongly anisotropy of fracture elements and the effects of border region on the transmissibility and provide more accurate results, on the other hand, the large fractures are modeled explicitly as major fluid conduits by the Embedded Discrete Fracture Model. Then, an efficient numerical algorithm based on the Mimetic Finite Difference method is developed to solve the hybrid method. At the end, several numerical examples are carried out to verify the accuracy and applicability of the proposed numerical model.