Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userFocusing inversion technique applied to radar tomographic data
168
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Traveltime tomography is a very effective tool to reconstruct acoustic, seismic or electromagnetic wave speed distribution. To infer the velocity image of the medium from the measurements of first arrivals is a typical example of ill-posed problem. In the framework of Tikhonov regularization theory, in order to replace an ill-posed problem by a well-posed one and to get a unique and stable solution, a stabilizing functional (stabilizer) has to be introduced. The stabilizer selects the desired solution from a class of solutions with a specific physical and/or geometrical property; e.g., the existence of sharp boundaries separating media with different petrophysical parameters. Usually stabilizers based on maximum smoothness criteria are used during the inversion process; in these cases the solutions provide smooth images which, in many situations, do not describe the examined objects properly. Recently a new algorithm of direct minimization of the Tikhonov parametric functional with minimum support stabilizer has been introduced; it produces clear and focused images of targets with sharp boundaries. In this research we apply this new technique to real radar tomographic data and we compare the obtained result with the solution generated by the more traditional minimum norm stabilizer.
rate research
Read More
We present an algorithm for focusing inversion of electrical resistivity tomography (ERT) data. ERT is a typical example of ill-posed problem. Regularization is the most common way to face this kind of problems; it basically consists in using a priori information about targets to reduce the ambiguity and the instability of the solution. By using the minimum gradient support (MGS) stabilizing functional, we introduce the following geometrical prior information in the reconstruction process: anomalies have sharp boundaries. The presented work is embedded in a project (L.A.R.A.) which aims at the estimation of hydrogeological properties from geophysical investigations. L.A.R.A. facilities include a simulation tank (4 m x 8 m x 1.35 m); 160 electrodes are located all around the tank and used for 3-D ERT. Because of the large number of electrodes and their dimensions, it is important to model their effect in order to correctly evaluate the electrical system response. The forward modelling in the presented algorithm is based on the so-called complete electrode model that takes into account the presence of the electrodes and their contact impedances. In this paper, we compare the results obtained with different regularizing functionals applied on a synthetic model.
Seismic full-waveform inversion (FWI), which uses iterative methods to estimate high-resolution subsurface models from seismograms, is a powerful imaging technique in exploration geophysics. In recent years, the computational cost of FWI has grown exponentially due to the increasing size and resolution of seismic data. Moreover, it is a non-convex problem and can encounter local minima due to the limited accuracy of the initial velocity models or the absence of low frequencies in the measurements. To overcome these computational issues, we develop a multiscale data-driven FWI method based on fully convolutional networks (FCN). In preparing the training data, we first develop a real-time style transform method to create a large set of synthetic subsurface velocity models from natural images. We then develop two convolutional neural networks with encoder-decoder structure to reconstruct the low- and high-frequency components of the subsurface velocity models, separately. To validate the performance of our data-driven inversion method and the effectiveness of the synthesized training set, we compare it with conventional physics-based waveform inversion approaches using both synthetic and field data. These numerical results demonstrate that, once our model is fully trained, it can significantly reduce the computation time, and yield more accurate subsurface velocity models in comparison with conventional FWI.
Multistatic ground-penetrating radar (GPR) signals can be imaged tomographically to produce three-dimensional distributions of image intensities. In the absence of objects of interest, these intensities can be considered to be estimates of clutter. These clutter intensities spatially vary over several orders of magnitude, and vary across different arrays, which makes direct comparison of these raw intensities difficult. However, by gathering statistics on these intensities and their spatial variation, a variety of metrics can be determined. In this study, the clutter distribution is found to fit better to a two-parameter Weibull distribution than Gaussian or lognormal distributions. Based upon the spatial variation of the two Weibull parameters, scale and shape, more information may be gleaned from these data. How well the GPR array is illuminating various parts of the ground, in depth and cross-track, may be determined from the spatial variation of the Weibull scale parameter, which may in turn be used to estimate an effective attenuation coefficient in the soil. The transition in depth from clutter-limited to noise-limited conditions (which is one possible definition of GPR penetration depth) can be estimated from the spatial variation of the Weibull shape parameter. Finally, the underlying clutter distributions also provide an opportunity to standardize image intensities to determine when a statistically significant deviation from background (clutter) has occurred, which is convenient for buried threat detection algorithm development which needs to be robust across multiple different arrays.
The three electromagnetic properties appearing in Maxwells equations are dielectric permittivity, electrical conductivity and magnetic permeability. The study of point diffractors in a homogeneous, isotropic, linear medium suggests the use of logarithms to describe the variations of electromagnetic properties in the earth. A small anomaly in electrical properties (permittivity and conductivity) responds to an incident electromagnetic field as an electric dipole, whereas a small anomaly in the magnetic property responds as a magnetic dipole. Neither property variation can be neglected without justification. Considering radiation patterns of the different diffracting points, diagnostic interpretation of electric and magnetic variations is theoretically feasible but is not an easy task using Ground Penetrating Radar. However, using an effective electromagnetic impedance and an effective electromagnetic velocity to describe a medium, the radiation patterns of a small anomaly behave completely differently with source-receiver offset. Zero-offset reflection data give a direct image of impedance variations while large-offset reflection data contain information on velocity variations.
We introduce a generalization of time-domain wavefield reconstruction inversion to anisotropic acoustic modeling. Wavefield reconstruction inversion has been extensively researched in recent years for its ability to mitigate cycle skipping. The original method was formulated in the frequency domain with acoustic isotropic physics. However, frequency-domain modeling requires sophisticated iterative solvers that are difficult to scale to industrial-size problems and more realistic physical assumptions, such as tilted transverse isotropy, object of this study. The work presented here is based on a recently proposed dual formulation of wavefield reconstruction inversion, which allows time-domain propagator that are suitable to both large scales and more accurate physics.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
