首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
This paper presents an application of the normalized surface magnetic charge (NSMC) model to discriminate objects of interest, such as unexploded ordnance (UXO), from innocuous items in cases when UXO electromagnetic induction (EMI) responses are contaminated by signals from other objects. Over the entire EMI spectrum considered here (tens of Hertz up to several hundreds of kHz), the scattered magnetic field outside the object can be produced mathematically by equivalent magnetic charges. The amplitudes of these charges are determined from measurement data and normalized by the excitation field. The model takes into account the scatterer's heterogeneity and near- and far-field effects. For classification algorithms, the frequency spectrum of the total NSMC is proposed and investigated as a discriminant. The NSMC is combined with the differential evolution (DE) algorithm in a two-step inversion procedure. To illustrate the applicability of the DE–NSMC algorithm, blind test data are processed and analyzed for cases in which signals from nearby objects frequently overlap. The method was highly successful in distinguishing UXO from accompanying clutter.  相似文献   

2.
Statistical signal processing techniques have shown progress in discriminating UXO from clutter when the objects occur in isolation. Under this condition, only a single object contributes to the measured sensor data. For multiple closely spaced subsurface objects, however, the unprocessed sensor measurement is a mixture of the responses from several objects. Consequently, the unprocessed measurements cannot be used directly to discriminate UXO from clutter. In this paper, we implement independent component analysis (ICA), a well-established blind source separation (BSS) technique, to recover the unobserved object signatures from the mixed measurement data obtained by simulating electromagnetic induction (EMI) sensor data, and then use the recovered signatures for UXO/clutter discrimination. Discrimination performance depends on multiple factors, including the number of clutter objects in proximity to the UXO, the separation distance between the UXO and clutter, and the number of mixed measurements available. Simulation results are presented illustrating the impact of these factors on discrimination performance.  相似文献   

3.
The time-domain EM induction response of non-magnetic and magnetic targets can be approximated using a conductive permeable prism composed of six faces of conductive plates, each face being composed of a set of conductive ribbons. The effect of magnetic permeability is included by the use of two “apparent flux gathering” coefficients, and two “effective magnetic permeability” coefficients, in the axial and transverse directions. These four magnetic property coefficients are a function of physical properties and geometry of the target, but are independent of prism orientation relative to a transmitter. The approximation algorithm is computationally fast, allowing inversions for target parameters to be achieved in seconds. The model is tested on profiles acquired with a Geonics EM63 time-domain EM metal detector over a non-magnetic copper pipe target, and a steel artillery shell in horizontal and vertical orientations. Results show that this approximation to a permeable prism has a capability of fitting geometric, conductivity and magnetic parameters at both early and late sample times. The magnetic parameters show strong change from early to late times on the EMI decay curve, indicating that the magnetic properties of the target have non-linear characteristics. It is proposed that these magnetic parameters and the nature of their non-linearity may carry additional discrimination information for distinguishing between intact munitions and scrap in UXO studies.  相似文献   

4.
Fluxgate vector magnetometers can be powerful instruments in magnetic mapping and in the detection of unexploded ordnance (UXO). Being lightweight and having low power requirements, they allow the development of multi-sensor use. The main problems in using fluxgate magnetometers arise from calibration errors and drift but these can be overcome using a quick and simple method of calibration in the field. This method also has the advantage of compensating permanent and induced magnetic fields generated by magnetized objects carried with the sensors. Measurement accuracy is similar to that obtained with scalar magnetometers. Multi-magnetometer profiling allows direct inversion of raw magnetic data along profiles in order to locate and characterize dipoles typically generated by UXO. The method used is a non-linear inverse procedure that estimates the three coordinates and magnetization vector of the dipole. Application of the technique to both a synthetic case and a field example illustrates the advantages of the method compared with the use of the analytic signal.  相似文献   

5.
Advanced discrimination methods and careful optimization of operational procedures are critical for efficient remediation of unexploded ordnance (UXO) contaminated sites. In this paper, we report on our experiences with a 200 acre magnetic survey that was collected and processed under production survey conditions at Chevallier Ranch, Montana. All anomalies with fitted moments above 0.05 Am2 were excavated. During the survey the magnetic remanence metric was predicted but not used to guide the discrimination. The retrospective analysis presented here reveals that discrimination using remanence would have significantly reduced the total number of anomalies (with good dipolar fits) that needed to be excavated, from 524 to 290 while still recovering all 69 UXO. The false alarm rate (FAR = number of non-UXOs excavated divided / number of UXO found) was reduced from 6.3 to 2.9. At a cut-off of 75% remanence, 77% of anomalies due to shrapnel and metallic debris and 64% of geological anomalies were rejected.Geological anomalies due to variations in magnetite concentration introduced a significant human-element into the interpretation process. Three different interpreters added a total of 305 additional anomalies that were not fit with a dipole model and which were later found to be non-UXO. Between 40 and 50% of anomalies picked by the two relatively inexperienced interpreters who analyzed the data turned out to be geology, as compared to 14% for an experienced interpreter. Critical analysis of results, operator training and feedback from the UXO technicians validating the anomaly are essential components towards improving the quality and consistency of the anomaly interpretations. This is consistent with the tenants of Total Quality Management (TQM). We compare the actual FAR that resulted during the survey when there was little feedback between UXO technician validation results, to a hypothetical result that could have been achieved had there been a constant feedback system in place at the onset of operations. Feedback would have significantly reduced the number of geological anomalies and decreased the FAR from 10.7 to 4.0.The hypothetical results presented here demonstrate the value of using TQM principles to guide the UXO remediation process. They further show that improvements in the efficiency and costs of UXO remediation require both technological advances and operational optimization of the technology when implemented in a production setting. Furthermore, by treating geophysical modeling and UXO validation as separate entities, both with respect to contracting and operational reporting, there is little incentive for the geophysicist to leave an anomaly off the dig-sheet. Only potential negative consequences will result if that anomaly is later found to be a UXO. An incentive based mechanism that rewards the geophysicist for reductions in follow-on costs would have a strong potential to reduce the number of unnecessary excavations, and hence reduce the total cost of the UXO remediation effort.  相似文献   

6.
Hydrothermal systems are often studied by collecting thermal gradient data and temperature/depth curves. These data contain important information about the flow field, the evolution of the hydrothermal system, and the location and nature of the ultimate heat sources. Thermal data are interpreted by the “forward” method; the thermal field is calculated based on selected initial conditions and boundary conditions such as temperature and permeability distributions. If the calculated thermal field matches the data, the chosen conditions are inferred to be possibly correct. Because many sets of initial conditions may produce similar thermal fields, users of the “forward” method may inadvertently miss the correct set of initial conditions. Analytical methods for “inverting” data also allow the determination of all the possible solutions consistent with the definition of the problem. In this paper we suggest an approach for inverting thermal data from a hydrothermal system, and compare it to the more conventional approach. We illustrate the difference in the methods by comparing their application to the Salton Sea Geothermal Field by Lau (1980a) and Kasameyer, et al. (1984). In this particular example, the inverse method was used to draw conclusions about the age and total rate of fluid flow into the hydrothermal system.  相似文献   

7.
This work adopts a continuation approach, based on path tracking in model space, to solve the non-linear least-squares problem for discrimination of unexploded ordnance (UXO) using multi-receiver electromagnetic induction (EMI) data. The forward model corresponds to a stretched-exponential decay of eddy currents induced in a magnetic spheroid. We formulate an over-determined, or under-parameterized, inverse problem. An example using synthetic multi-receiver EMI responses illustrates the efficiency of the method. The fast inversion of actual field multi-receiver EMI responses of inert, buried ordnances is also shown. Software based on the continuation method could be installed within a multi-receiver EMI sensor and used for near-real-time UXO decision-making purposes without the need for a highly-trained operator.  相似文献   

8.
We present two case studies in the night and evening sides of the auroral oval, based on plasma and field measurements made at low altitudes by the AUREOL-3 satellite, during a long period of stationary magnetospheric convection (SMC) on November 24, 1981. The basic feature of both oval crossings was an evident double oval pattern, including (1) a weak arc-type structure at the equatorial edge of the oval/polar edge of the diffuse auroral band, collocated with an upward field-aligned current (FAC) sheet of ≈1.0 μA m−2, (2) an intermediate region of weaker precipitation within the oval, (3) a more intense auroral band at the polar oval boundary, and (4) polar diffuse auroral zone near the polar cap boundary. These measurements are compared with the published magnetospheric data during this SMC period, accumulated by Yahnin et al. and Sergeev et al., including a semi-empirical radial magnetic field profile BZ in the near-Earth neutral sheet, with a minimum at about 10–14 RE. Such a radial BZ profile appears to be very similar to that assumed in the “minimum B/cross-tail line current” model by Galperin et al. (GVZ92) as the “root of the arc”, or the arc generic region. This model considers a FAC generator mechanism by Grad-Vasyliunas-Boström-Tverskoy operating in the region of a narrow magnetic field minimum in the near-Earth neutral sheet, together with the concept of ion non-adiabatic scattering in the “wall region”. The generated upward FAC branch of the double sheet current structure feeds the steady auroral arc/inverted-V at the equatorial border of the oval. When the semi-empirical BZ profile is introduced in the GVZ92 model, a good agreement is found between the modelled current and the measured characteristics of the FACs associated with the equatorial arc. Thus the main predictions of the GVZ92 model concerning the “minimum-B” region are consistent with these data, while some small-scale features are not reproduced. Implications of the GVZ92 model are discussed, particularly concerning the necessary conditions for a substorm onset that were not fulfilled during the SMC period.  相似文献   

9.
Supported by the recent advancement of experimental test methods, numerical simulation, and high‐speed communication networks, it is possible to distribute geographically the testing of structural systems using hybrid experimental–computational simulation. One of the barriers for this advanced testing is the lack of flexible software for hybrid simulation using heterogeneous experimental equipment. To address this need, an object‐oriented software framework is designed, developed, implemented, and demonstrated for distributed experimental–computational simulation of structural systems. The software computes the imposed displacements for a range of test methods and co‐ordinates the control of local and distributed configurations of experimental equipment. The object‐oriented design of the software promotes the sharing of modules for experimental equipment, test set‐ups, simulation models, and test methods. The communication model for distributed hybrid testing is similar to that used for parallel computing to solve structural simulation problems. As a demonstration, a distributed pseudodynamic test was conducted using a client–server approach, in which the server program controlled the test equipment in Japan and the client program performed the computational simulation in the United States. The distributed hybrid simulation showed that the software framework is flexible and reliable. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

10.
A case study of the dayside cusp/cleft region during an interval of stationary magnetospheric convection (SMC) on November, 24, 1981 is presented, based on detailed measurements made by the AUREOL-3 satellite. Layered small-scale field-aligned current sheets, or loops, superimposed to a narrow V-shaped ion dispersion structure, were observed just equatorward from the region of the “cusp proper”. The equatorward sheet was accompanied by a very intense and short (less than 1 s) ion intensity spike at 100 eV. No major differences were noted of the characteristics of the LLBL, or “boundary cusp”, and plasma mantle precipitation during this SMC period from those typical of the cusp/cleft region for similar IMF conditions. Simultaneous NOAA-6 and NOAA-7 measurements described in Despirak et al. were used to estimate the average extent of the “cusp proper” (defined by dispersed precipitating ions with the energy flux exceeding 10−3 erg cm−2 s−1) during the SMC period, as ≈0.73∼ ILAT width, 2.6–3.4 h in MLT, and thus the recently merged magnetic flux, 0.54–0.70 × 107 Wb. This, together with the average drift velocity across the cusp at the convection throat, ≈0.5 km s−1, allowed to evaluate the cusp merging contribution to the total cross-polar cap potential difference, ≈33.8–43.8 kV. It amounts to a quite significant part of the total cross-polar cap potential difference evaluated from other data. A “shutter” scenario is suggested for the ion beam injection/penetration through the stagnant plasma region in the outer cusp to explain the pulsating nature of the particle injections in the low- and medium-altitude cusp region.  相似文献   

11.
The “fluid-flow tomography”, an advanced technique for geoelectrical survey based on the conventional mise-à-la-masse measurement, has been developed by Exploration Geophysics Laboratory at the Kyushu University. This technique is proposed to monitor fluid-flow behavior during water injection and production in a geothermal field. However data processing of this technique is very costly. In this light, this paper will discuss the solution to cost reduction by applying a neural network in the data processing. A case study in the Takigami geothermal field in Japan will be used to illustrate this. The achieved neural network in this case study is three-layered and feed-forward. The most successful learning algorithm in this network is the Resilient Propagation (RPROP). Consequently, the study advances the pragmatism of the “fluid-flow tomography” technique which can be widely used for geothermal fields. Accuracy of the solution is then verified by using root mean square (RMS) misfit error as an indicator.  相似文献   

12.
Analysis of Pc3 observational data along the 210° magnetic meridian showed a complicated frequency-latitude structure at middle latitudes. The observed period-latitude distributions vary between events with a “noisy source”: the D component has a colored-noise spectrum, while the spectrum of H component exhibits regular peaks that vary with latitude, and events with a “band-limited source”: the spectral power density of the D component is enhanced at certain frequencies throughout the network. For most ULF events a local gap of the H component amplitude has been exhibited at both conjugate stations at L ≃ 2.1. A quantitative interpretation has been given assuming that band-limited MHD emission from an extra-magnetospheric source is distorted by local field line resonances. Resonant frequencies had been singled out with the use of the asymmetry between spectra of H and D components. Additionally, a local resonant frequency at L ≃ 1.6 was determined by the quasi-gradient method using the data from nearly conjugate stations. The experimentally determined local resonance frequencies agree satisfactorily with those obtained from a numerical model of the Alfven resonator with the equatorial plasma density taken by extrapolation of Carpenter-Anderson model. We demonstrate how simple methods of hydromagnetic spectroscopy enable us to monitor simultaneously both the magnitude of the IMF and the magnetospheric plasma density from ULF data.  相似文献   

13.
电磁场数值模拟的背景场/异常场算法是三维正演的有效策略之一,优点为采用解析法计算电磁场背景场代替场源项、克服了场源奇异性,缺点为不适用于发射源布置于起伏地表或背景模型复杂的情形.总场算法是直接对电磁场总场开展数值模拟,其难点是有效加载场源、保证近区与过渡区数值解精度.本文以水平电偶源形式分段加载接地长导线源,并以电场总场Helmholtz方程为矢量有限元法控制方程,实现了基于非结构化四面体网格剖分的接地长导线源频率域电磁法三维正演.通过与均匀全空间中水平电偶源产生的电场解析解对比,验证了本文算法的正确性,并分析了四面体外接圆半径与其最短棱边的最大比值和四面体二面角最小值对数值解精度的影响规律.通过与块状高导体地电模型的积分方程法、有限体积法和基于磁矢量势Helmholtz方程的有限元法数值解对比,进一步验证了本文算法正确性,同时说明了非结构化四面体网格能够更加精细地剖分电性异常体,利于获得精确数值解.  相似文献   

14.
Tholeiitic basalts and associated intrusives are the major component of the Karoo igneous province. They are of Mesozoic age and constitute one of the world's classic continental flood basalt (CFB) provinces. It has been argued that most Karoo basalts have not undergone significant contamination with continental crust and that their lithospheric mantle source areas were enriched in incompatible minor and trace elements during the Proterozoic. The only exceptions to this are late-stage MORB-like dolerites near the present-day continental margins which are considered to be of asthenospheric origin.When data for the “southern” Karoo basalts are plotted on many of the geochemical discriminant diagrams which have been used to infer tectonic setting, essentially all of them would be classified as calc-alkali basalts (CAB's) or low-K tholeiites. Virtually none of them plot in the compositional fields designated as characteristic of “within-plate” basalts. There is little likelihood that the compositions of the Karoo basalts can be controlled by active subduction at the time of their eruption and no convincing evidence that a “subduction component” has been added to the subcontinental lithospheric mantle under the entire area in which the basalts crop out. It must be concluded that the mantle source areas for CAB's and the southern Karoo basalts have marked similarities.In contrast, the data for “northern” Karoo basalts largely plot in the “within-plate” field on geochemical discriminant diagrams. Available data suggest that the source composition and/or the restite mineralogy and degree of partial melting are different for southern and northern Karoo basalts. There is no evidence for any difference in tectonic setting between the southern and northern Karoo basalts at the time they were erupted. This appears to be clear evidence that specific mantle source characteristics and/or magmatic processes can vary within a single CFB province to an extent that renders at least some geochemical discriminant diagrams most unreliable for classifying tectonic environment with respect to continental volcanic rocks.  相似文献   

15.
  相似文献   

16.
We study an interval of 56 h on January 16 to 18, 1995, during which the GEOTAIL spacecraft traversed the duskside magnetosheath from X ≅ −15 to −40 RE and the EPIC/ICS and EPIC/STICS sensors sporadically detected tens of energetic particle bursts. This interval coincides with the expansion and growth of a great geomagnetic storm. The flux bursts are strongly dependent on the magnetic field orientation. They switch on whenever the Bz component approaches zero (Bz ≅ 0 nT). We strongly suggest a magnetospheric origin for the energetic ions and electrons streaming along these “exodus channels”. The time profiles for energetic protons and “tracer” O+ ions are nearly identical, which suggests a common source. We suggest that the particles leak out of the magnetosphere all the time and that when the magnetosheath magnetic field connects the spacecraft to the magnetotail, they stream away to be observed by the GEOTAIL sensors. The energetic electron fluxes are not observed as commonly as the ions, indicating that their source is more limited in extent. In one case study the magnetosheath magnetic field lines are draped around the magnetopause within the YZ plane and a dispersed structure for peak fluxes of different species is detected and interpreted as evidence for energetic electrons leaking out from the dawn LLBL and then being channelled along the draped magnetic field lines over the magnetopause. Protons leak from the equatorial dusk LLBL and this spatial differentiation between electron and proton sources results in the observed dispersion. A gradient of energetic proton intensities toward the ZGSM= 0 plane is inferred. There is a permanent layer of energetic particles adjacent to the magnetosheath during this interval in which the dominant component of the magnetic field was Bz.  相似文献   

17.
A simple library based algorithm for the identification of unexploded ordnance is tested on time domain electromagnetic data. A library of polarization tensors is generated from data acquired on a test stand over a collection of different UXO. The objective of the algorithm is to determine which target of the library is most likely to have produced an observed data anomaly. For each target in our library a non-linear inverse problem is solved for the position and orientation that minimizes the least-squares difference between the observed data anomaly and the data predicted from each target. This technique avoids direct inversion for polarization tensor, making it feasible for cases where sensor data quality may not be sufficient to support confident estimation of model parameters. For cases where the background soil response is significant, we also invert for the t−1 decay characteristic of viscous remnant magnetic soil. We present results of preliminary tests of the library technique to Geonics EM63 time domain electromagnetic data collected on a test plot seeded with UXO. These tests demonstrate an excellent ability to accurately identify isolated targets. Misidentification of single targets occur when data anomalies have low signal to noise ratios or when targets within the library have similar polarization tensors.  相似文献   

18.
Data collected along transects are becoming more common in environmental studies as indirect measurement devices, such as geophysical sensors, that can be attached to mobile platforms become more prevalent. Because exhaustive sampling is not always possible under constraints of time and costs, geostatistical interpolation techniques are used to estimate unknown values at unsampled locations from transect data. It is known that outlying observations can receive significantly greater ordinary kriging weights than centrally located observations when the data are contiguously aligned along a transect within a finite search window. Deutsch (1994) proposed a kriging algorithm, finite domain kriging, that uses a redundancy measure in place of the covariance function in the data-to-data kriging matrix to address the problem of overweighting the outlying observations. This paper compares the performances of two kriging techniques, ordinary kriging (OK) and finite domain kriging (FDK), on examining unexploded ordnance (UXO) densities by comparing prediction errors at unsampled locations. The impact of sampling design on object count prediction is also investigated using data collected from transects and at random locations. The Poisson process is used to model the spatial distribution of UXO for three 5000 × 5000 m fields; one of which does not have any ordnance target (homogeneous field), while the other two sites have an ordnance target in the center of the site (isotropic and anisotropic fields). In general, for a given sampling transects width, the differences between OK and FDK in terms of the mean error and the mean square error are not significant regardless of the sampled area and the choice of the field. When 20% or more of the site is sampled, the estimation of object counts is unbiased on average for all three fields regardless of the choice of the transect width and the choice of the kriging algorithm. However, for non-homogeneous fields (isotropic and anisotropic fields), the mean error fluctuates considerably when a small number of transects are sampled. The difference between the transect sampling and the random sampling in terms of prediction errors becomes almost negligible if more than 20% of the site is sampled. Overall, FDK is no better than OK in terms of the prediction performances when the transect sampling procedure is used.  相似文献   

19.
A processing strategy and the corresponding software architecture for the processing of GOCE (Gravity field and steady-state Ocean Circulation Explorer) observables is presented and described, with the major objective to compute a high-accuracy, high-resolution spherical harmonic model of the Earth's gravity field. The combination of two numerical solution strategies, i.e. the rigorous solution of the corresponding large normal equation systems applying parallel processing (on a PC cluster) as the core solver, and the fast semianalytic approach as a quick-look gravity field analysis (QL-GFA) tool, is proposed. Such a method fusion benefits from the advantages of the individual components: the rigorous inversion of the system providing also the full variance-covariance information, and the quickness enabling the consecutive production of intermediate gravity field solutions, for the purpose to analyse partial and incomplete data sets and to derive a diagnosis of the performance of the GOCE measurement system. The functionality and operability of the individual components are demonstrated in the framework of a closed loop simulation, which is based on a realistic mission scenario both in terms of the orbit configuration and the coloured measuring noise. Special concern is given to the accuracy of the recovered coefficients, the numerical behaviour, the required computing time, and the particular role of the individual modules within the processing chain. In the case of the core solver, it is demonstrated that the assembling and rigorous solution of large normal equation systems can be handled by using Beowulf clusters within a reasonable computing time. The application of the quick-look tool to partial data sets with short-term data gaps is demonstrated on the basis of several case studies. Additionally, the spectral analysis of the residuals of the adjustment is presented as a valuable tool for the verification of the noise characteristics of the GOCE gradiometer.  相似文献   

20.
In this paper we deal with an indirect measure of the dielectric permittivity of the soil starting from GPR surface data collected on a buried “cooperative” target, meant as an object buried on purpose and whose extent is known a-priori. This target is exploited in order to achieve, from its image obtained from a suitable GPR data processing, an indirect measure of the dielectric permittivity of the embedding soil. GPR data processing is based on a linear microwave tomographic approach funded on the Born Approximation. Using this Born approach on two-dimensional inversion tests, we investigate the effect of the soil's electrical conductivity and permittivity on this indirect measure and demonstrate that the electrical field scattered by a spot-like buried object permits an accurate estimation of the soil permittivity even when no information of the soil conductivity is available.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号