Resistivity inversion in 2-D anisotropic media: numerical experiments
NASA Astrophysics Data System (ADS)
Wiese, Timothy; Greenhalgh, Stewart; Zhou, Bing; Greenhalgh, Mark; Marescot, Laurent
2015-04-01
Many rocks and layered/fractured sequences have a clearly expressed electrical anisotropy although it is rare in practice to incorporate anisotropy into resistivity inversion. In this contribution, we present a series of 2.5-D synthetic inversion experiments for various electrode configurations and 2-D anisotropic models. We examine and compare the image reconstructions obtained using the correct anisotropic inversion code with those obtained using the false but widely used isotropic assumption. Superior reconstruction in terms of reduced data misfit, true anomaly shape and position, and anisotropic background parameters were obtained when the correct anisotropic assumption was employed for medium to high coefficients of anisotropy. However, for low coefficient values the isotropic assumption produced better-quality results. When an erroneous isotropic inversion is performed on medium to high level anisotropic data, the images are dominated by patterns of banded artefacts and high data misfits. Various pole-pole, pole-dipole and dipole-dipole data sets were investigated and evaluated for the accuracy of the inversion result. The eigenvalue spectra of the pseudo-Hessian matrix and the formal resolution matrix were also computed to determine the information content and goodness of the results. We also present a data selection strategy based on high sensitivity measurements which drastically reduces the number of data to be inverted but still produces comparable results to that of the comprehensive data set. Inversion was carried out using transversely isotropic model parameters described in two different co-ordinate frames for the conductivity tensor, namely Cartesian versus natural or eigenframe. The Cartesian frame provided a more stable inversion product. This can be simply explained from inspection of the eigenspectra of the pseudo-Hessian matrix for the two model descriptions.
ELRIS2D: A MATLAB Package for the 2D Inversion of DC Resistivity/IP Data
NASA Astrophysics Data System (ADS)
Akca, Irfan
2016-04-01
ELRIS2D is an open source code written in MATLAB for the two-dimensional inversion of direct current resistivity (DCR) and time domain induced polarization (IP) data. The user interface of the program is designed for functionality and ease of use. All available settings of the program can be reached from the main window. The subsurface is discretized using a hybrid mesh generated by the combination of structured and unstructured meshes, which reduces the computational cost of the whole inversion procedure. The inversion routine is based on the smoothness constrained least squares method. In order to verify the program, responses of two test models and field data sets were inverted. The models inverted from the synthetic data sets are consistent with the original test models in both DC resistivity and IP cases. A field data set acquired in an archaeological site is also used for the verification of outcomes of the program in comparison with the excavation results.
2D resistivity inversion using conjugate gradients for a finite element discretization
NASA Astrophysics Data System (ADS)
Bortolozo, C. A.; Santos, F. M.; Porsani, J. L.
2014-12-01
In this work we present a DC 2D inversion algorithm using conjugate gradients relaxation to solve the maximum likelihood inverse equations. We apply, according to Zhang (1995), the maximum likelihood inverse theory developed by Tarantola and Valette (1982) to our 2D resistivity inversion. This algorithm was chosen to this research because it doesn't need to calculate the field's derivatives. Since conjugate gradient techniques only need the results of the sensitivity matrix Ã or its transpose ÃT multiplying a vector, the actual computation of the sensitivity matrix are not performed, according to the methodology described in Zhang (1995). In Zhang (1995), the terms Ãx and ÃTy, are dependent of the stiffness matrix K and its partial derivative ∂K⁄∂ρ. The inversion methodology described in Zhang (1995) is for the case of 3D electrical resistivity by finite differences discretization. So it was necessary to make a series of adjustments to obtain a satisfactory result for 2D electrical inversion using finite element method. The difference between the modeling of 3D resistivity with finite difference and the 2D finite element method are in the integration variable, used in the 2D case. In the 2D case the electrical potential are initially calculated in the transformed domain, including the stiffness matrix, and only in the end is transformed in Cartesian domain. In the case of 3D, described by Zhang (1995) this is done differently, the calculation is done directly in the Cartesian domain. In the literature was not found any work describing how to deal with this problem. Because the calculations of Ãx and ÃTy must be done without having the real stiffness matrix, the adaptation consist in calculate the stiffness matrix and its partial derivative using a set of integration variables. We transform those matrix in the same form has in the potential case, but with different sets of variables. The results will be presented and are very promising.
Comparison of 1D, 2D and 2.5D Constrained Inversion of Electrical Resistivity Data
NASA Astrophysics Data System (ADS)
Catt, L. M.; West, J.; Clark, R. A.
2007-05-01
Clay-rich till plains cover much of the UK. Such sites are attractive locations for landfills, since the till cover lowers the risk of landfill leachate entering groundwater. However, such tills often contain discrete sand and gravel bodies that can act as leachate flow routes. Such bodies may not be detected by conventional site investigation techniques such as drilling boreholes and trial pitting. A combined geoelectrical survey was carried out at a study site typical of such till plains and close to cliff exposures, which allowed direct mapping of sand bodies. Electrical resistivity tomography (ERT), resistivity cone penetrometry (RCPT) and frequency-domain electromagnetic data were collected. In a previous study, the electromagnetic and RCPT data were used to construct reference models for 2D inversion of the ERT data. The use of these reference models improved the solution models produced by inversion. We showed that the best solution model produced by inversion with a range of reference models could be determined without a priori knowledge of the true geoelectrical structure. This was done by using the area-weighted L2 norm between the solution models and associated reference models as a proxy for the misfit between the solution models and the true geoelectrical structure of the ground. In order to assess the most suitable method for combining invasive and non-invasive measurements, we compare both constrained and unconstrained 1D, 2D and 2.5D inversions of resistivity data collected at the study site. Preliminary results suggest that for 2.5D inversion the true 3D geoelectrical structure of the ground at the field study site is not sufficiently well known for comparison between the solution models and the true geoelectrical structure of the ground to be made. The results of work in progress evaluating layer-depth- constrained 1D inversion will be presented at the meeting.
NASA Astrophysics Data System (ADS)
Yi, Myeong-Jong; Sasaki, Yutaka
2015-11-01
Frequency-domain loop-loop electromagnetic (EM) methods are sensitive to the magnetic susceptibility of the Earth as well as its resistivity. Thus, inversion techniques have been used to simultaneously reconstruct both resistivity and susceptibility models from EM data. However, to take full advantage of inversion methods, calibration errors must be assessed and removed because ignoring them can result in misleading models. We present a multidimensional inversion method that jointly inverts EM and direct current (DC) resistivity data to derive offset errors as well as resistivity and susceptibility models, assuming that calibration errors can be represented by in-phase and quadrature offsets at each frequency. Addition of independent data such as DC data is effective for more accurately estimating the offsets, resulting in more reliable subsurface models. Synthetic examples involving small-loop EM data show that simultaneous inversion for resistivity and susceptibility is not stable, because of strong correlations between in-phase offset parameters and background susceptibility, but that the offsets are well determined when the data misfit is reduced rapidly in the early iteration step. Improvements achieved by joint inversion are mainly on the resistivity model. For airborne electromagnetic (AEM) data, the inversion process is stable, because AEM data are acquired using more loop-loop geometries and a wider range of frequencies. As a result, both the resistivity and susceptibility models are significantly improved by joint inversion.
Ding, Lin; Song, An; Dai, Meng; Xu, Min; Sun, Wanwan; Xu, Baihui; Sun, Jichao; Wang, Tiange; Xu, Yu; Lu, Jieli; Wang, Weiqing; Bi, Yufang; Ning, Guang
2015-01-01
Lipoprotein (a) [Lp(a)], an LDL-like particle, has been proposed as a causal risk factor for CVD among general populations. Meanwhile, both serum Lp(a) and diabetes increase the risk of CVD. However, the relationship between serum Lp(a) and T2D is poorly characterized, especially in the Asian population. Therefore, we conducted a cross-sectional study in 10,122 participants aged 40 years or older in Jiading District, Shanghai, China. Our study found that the prevalence of T2D was decreased from 20.9% to 15.0% from the lowest quartile to the highest quartile of serum Lp(a) concentrations (P for trend <0.0001). Logistic regression analyses showed that the odds ratios and 95% confidence intervals of prevalent T2D for quartiles 2–4 versus quartile 1 were 0.86 (0.73–1.01), 0.88 (0.75–1.04), and 0.76 (0.64–0.90) (P for trend = 0.0002), after adjustment for traditional confounding factors. Moreover, the risks for prevalent prediabetes, insulin resistance, and hyperinsulinemia were also decreased from the lowest to the top quartile. This inverse association between serum Lp(a) and T2D was not appreciably changed after we adjusted hypoglycemic medications or excluded the subjects with hypoglycemic and/or lipid-lowering agents and/or a history of self-reported CVD. PMID:25649924
Stochastic Inversion of 2D Magnetotelluric Data
Energy Science and Technology Software Center (ESTSC)
2010-07-01
The algorithm is developed to invert 2D magnetotelluric (MT) data based on sharp boundary parametrization using a Bayesian framework. Within the algorithm, we consider the locations and the resistivity of regions formed by the interfaces are as unknowns. We use a parallel, adaptive finite-element algorithm to forward simulate frequency-domain MT responses of 2D conductivity structure. Those unknown parameters are spatially correlated and are described by a geostatistical model. The joint posterior probability distribution function ismore » explored by Markov Chain Monte Carlo (MCMC) sampling methods. The developed stochastic model is effective for estimating the interface locations and resistivity. Most importantly, it provides details uncertainty information on each unknown parameter. Hardware requirements: PC, Supercomputer, Multi-platform, Workstation; Software requirements C and Fortan; Operation Systems/version is Linux/Unix or Windows« less
Stochastic Inversion of 2D Magnetotelluric Data
Chen, Jinsong
2010-07-01
The algorithm is developed to invert 2D magnetotelluric (MT) data based on sharp boundary parametrization using a Bayesian framework. Within the algorithm, we consider the locations and the resistivity of regions formed by the interfaces are as unknowns. We use a parallel, adaptive finite-element algorithm to forward simulate frequency-domain MT responses of 2D conductivity structure. Those unknown parameters are spatially correlated and are described by a geostatistical model. The joint posterior probability distribution function is explored by Markov Chain Monte Carlo (MCMC) sampling methods. The developed stochastic model is effective for estimating the interface locations and resistivity. Most importantly, it provides details uncertainty information on each unknown parameter. Hardware requirements: PC, Supercomputer, Multi-platform, Workstation; Software requirements C and Fortan; Operation Systems/version is Linux/Unix or Windows
A new inversion method for (T2, D) 2D NMR logging and fluid typing
NASA Astrophysics Data System (ADS)
Tan, Maojin; Zou, Youlong; Zhou, Cancan
2013-02-01
One-dimensional nuclear magnetic resonance (1D NMR) logging technology has some significant limitations in fluid typing. However, not only can two-dimensional nuclear magnetic resonance (2D NMR) provide some accurate porosity parameters, but it can also identify fluids more accurately than 1D NMR. In this paper, based on the relaxation mechanism of (T2, D) 2D NMR in a gradient magnetic field, a hybrid inversion method that combines least-squares-based QR decomposition (LSQR) and truncated singular value decomposition (TSVD) is examined in the 2D NMR inversion of various fluid models. The forward modeling and inversion tests are performed in detail with different acquisition parameters, such as magnetic field gradients (G) and echo spacing (TE) groups. The simulated results are discussed and described in detail, the influence of the above-mentioned observation parameters on the inversion accuracy is investigated and analyzed, and the observation parameters in multi-TE activation are optimized. Furthermore, the hybrid inversion can be applied to quantitatively determine the fluid saturation. To study the effects of noise level on the hybrid method and inversion results, the numerical simulation experiments are performed using different signal-to-noise-ratios (SNRs), and the effect of different SNRs on fluid typing using three fluid models are discussed and analyzed in detail.
Efficient 2d full waveform inversion using Fortran coarray
NASA Astrophysics Data System (ADS)
Ryu, Donghyun; Kim, ahreum; Ha, Wansoo
2016-04-01
We developed a time-domain seismic inversion program using the coarray feature of the Fortran 2008 standard to parallelize the algorithm. We converted a 2d acoustic parallel full waveform inversion program with Message Passing Interface (MPI) to a coarray program and examined performance of the two inversion programs. The results show that the speed of the waveform inversion program using the coarray is slightly faster than that of the MPI version. The standard coarray lacks features for collective communication; however, it can be improved in following standards since it is introduced recently. The parallel algorithm can be applied for 3D seismic data processing.
An inverse design method for 2D airfoil
NASA Astrophysics Data System (ADS)
Liang, Zhi-Yong; Cui, Peng; Zhang, Gen-Bao
2010-03-01
The computational method for aerodynamic design of aircraft is applied more universally than before, in which the design of an airfoil is a hot problem. The forward problem is discussed by most relative papers, but inverse method is more useful in practical designs. In this paper, the inverse design of 2D airfoil was investigated. A finite element method based on the variational principle was used for carrying out. Through the simulation, it was shown that the method was fit for the design.
Image Appraisal for 2D and 3D Electromagnetic Inversion
Alumbaugh, D.L.; Newman, G.A.
1999-01-28
Linearized methods are presented for appraising image resolution and parameter accuracy in images generated with two and three dimensional non-linear electromagnetic inversion schemes. When direct matrix inversion is employed, the model resolution and posterior model covariance matrices can be directly calculated. A method to examine how the horizontal and vertical resolution varies spatially within the electromagnetic property image is developed by examining the columns of the model resolution matrix. Plotting the square root of the diagonal of the model covariance matrix yields an estimate of how errors in the inversion process such as data noise and incorrect a priori assumptions about the imaged model map into parameter error. This type of image is shown to be useful in analyzing spatial variations in the image sensitivity to the data. A method is analyzed for statistically estimating the model covariance matrix when the conjugate gradient method is employed rather than a direct inversion technique (for example in 3D inversion). A method for calculating individual columns of the model resolution matrix using the conjugate gradient method is also developed. Examples of the image analysis techniques are provided on 2D and 3D synthetic cross well EM data sets, as well as a field data set collected at the Lost Hills Oil Field in Central California.
Constriant inversion of 2D magnetotelluric data with anisotropic conductivities
NASA Astrophysics Data System (ADS)
Chen, X.; Weckmann, U.
2011-12-01
Within the framework of the German - South African geo-scientific research initiative Inkaba yeAfrica a series of magnetotelluric (MT) field experiments were conducted along the Agulhas-Karoo Transect in South Africa. This transect crosses several continental collision zones between the Cape Fold Belt, the Namaqua Natal Mobile Belt and the Kaapvaal Craton. Along the Cape Fold Belt (CFB) profile we can identify areas (>10 km) where MT sites exhibit phases over 90°. This phenomenon usually occurs in presence of electrical anisotropy. Due to the dense site spacing we are able to observe this behaviour consistently at several sites. The anisotropy of electrical conductivity is essentially a scale effect: Even if the conductivity is isotropic on the micro scale, it will become anisotropic on a larger scale if, in the averaging volume, preferred orientation (e.g., layering or lamination) exist. Therefore, it is necessary to understand the electrical anisotropy in more details and furthermore electrical anisotropy offers new degrees of freedom, which should allow a better interpretation of data. In 2D MT case with considering of electrical anisotropy, computing of impedance tensor requires two independent electric field solutions computed for two different source polarisations. Based on the forward problem formulation and its numerical approximation we derive partial differential equations for the sensitivities of the magnetotelluric fields with respect to the elements of the conductivity tensor within the medium. For illustration a sensitivity study for a simple synthetic model is shown. We present an algorithm for the inversion of 2D magnetotelluric data with anisotropic conductivities which is a extension of the well-known NLCG minimization algorithm to anisotropic model. To constrain the structure complexity, a penalty function consists of datamisfit, standard model roughness and quadratic variation of the conductivity tensor elements is minimized. To demonstrate the
Almost but not quite 2D, Non-linear Bayesian Inversion of CSEM Data
NASA Astrophysics Data System (ADS)
Ray, A.; Key, K.; Bodin, T.
2013-12-01
The geophysical inverse problem can be elegantly stated in a Bayesian framework where a probability distribution can be viewed as a statement of information regarding a random variable. After all, the goal of geophysical inversion is to provide information on the random variables of interest - physical properties of the earth's subsurface. However, though it may be simple to postulate, a practical difficulty of fully non-linear Bayesian inversion is the computer time required to adequately sample the model space and extract the information we seek. As a consequence, in geophysical problems where evaluation of a full 2D/3D forward model is computationally expensive, such as marine controlled source electromagnetic (CSEM) mapping of the resistivity of seafloor oil and gas reservoirs, Bayesian studies have largely been conducted with 1D forward models. While the 1D approximation is indeed appropriate for exploration targets with planar geometry and geological stratification, it only provides a limited, site-specific idea of uncertainty in resistivity with depth. In this work, we extend our fully non-linear 1D Bayesian inversion to a 2D model framework, without requiring the usual regularization of model resistivities in the horizontal or vertical directions used to stabilize quasi-2D inversions. In our approach, we use the reversible jump Markov-chain Monte-Carlo (RJ-MCMC) or trans-dimensional method and parameterize the subsurface in a 2D plane with Voronoi cells. The method is trans-dimensional in that the number of cells required to parameterize the subsurface is variable, and the cells dynamically move around and multiply or combine as demanded by the data being inverted. This approach allows us to expand our uncertainty analysis of resistivity at depth to more than a single site location, allowing for interactions between model resistivities at different horizontal locations along a traverse over an exploration target. While the model is parameterized in 2D, we
Development of a Geocryologic Model of Permafrost From 2D Inversion of IP Profiling
NASA Astrophysics Data System (ADS)
Fortier, R.; Leblanc, A.
2004-05-01
Non-invasive investigation of permafrost along a planned route of pipeline, road or airstrip in cold regions involves the use of effective methods for detecting, characterizing, mapping and monitoring permafrost conditions on various spatial and temporal scales. Among the available near-surface geophysical methods, the electrical resistivity imaging is probably the most suitable method since the resistivity contrast between unfrozen and frozen ground can be one or two orders of magnitude. Induced polarization (IP) profiling was carried out to study the spatial distribution of ground ice in two permafrost mounds near Umiujaq in Nunavik, Canada. A dipole-dipole array was used to perform the IP profiling. Pseudo-sections of electrical resistivity and chargeability giving a misrepresented cross-section of the sub-surface were first draw. The inversion of IP profiling was also performed using DCIP2D developed by UBC-GIF for estimating the spatial distribution of electrical properties in the ground to create realistic models of sub-surface resistivity and chargeability cross-section. The inverse models show clearly the presence of ice-rich core in the permafrost mounds. The ice-rich cores are underlined by high resistivity values while the unfrozen zones show low resistivity values. The localisation of the permafrost table is highlighted by a strong contrast of resistivity while the permafrost base is marked by a transitional change in resistivity. In the hollow between the permafrost mounds, the models show low resistivity values characteristic of unfrozen zone. A synthetic resistivity sounding built from the most acceptable inverse model correlates well with electrical resistivity logging carried out in the permafrost mound during cone penetration tests. The inversion of IP profiling is fundamental for defining realistic models of sub-surface resistivity and chargeability. Electrical resistivity imaging is a appropriate near-surface geophysical method for permafrost
Full-waveform inversion in 2D VTI media
NASA Astrophysics Data System (ADS)
Kamath, Nishant
Full-waveform inversion (FWI) is a technique designed to produce a high-resolution model of the subsurface by using information contained in entire seismic waveforms. This thesis presents a methodology for FWI in elastic VTI (transversely isotropic with a vertical axis of symmetry) media and discusses synthetic results for heterogeneous VTI models. First, I develop FWI for multicomponent data from a horizontally layered VTI model. The reflectivity method, which permits computation of only PP reflections or a combination of PP and PSV events, is employed to model the data. The Gauss-Newton technique is used to invert for the interval Thomsen parameters, while keeping the densities fixed at the correct values. Eigenvalue/eigenvector decompostion of the Hessian matrix helps analyze the sensitivity of the objective function to the model parameters. Whereas PP data alone are generally sufficient to constrain all four Thomsen parameters even for conventional spreads, including PS reflections provides better constraints, especially for the deeper part of the model. Next, I derive the gradients of the FWI objective function with respect to the stiffness coefficients of arbitrarily anisotropic media by employing the adjoint-state method. From these expressions, it is straightforward to compute the gradients for parameters of 2D heterogeneous VTI media. FWI is implemented in the time domain with the steepest-descent method used to iteratively update the model. The algorithm is tested on transmitted multicomponent data generated for Gaussian anomalies in Thomsen parameters embedded in homogeneous VTI media. To test the sensitivity of the objective function to different model parameters, I derive an an- alytic expression for the Frechet kernel of FWI for arbitrary anisotropic symmetry by using the Born approximation and asymptotic Green's functions. The amplitude of the kernel, which represents the radiation pattern of a secondary source (that source describes a perturbation
Analytic Differentiation of Barlat's 2D Criteria for Inverse Modeling
Endelt, Benny; Nielsen, Karl Brian; Danckert, Joachim
2005-08-05
The demand for alternative identification schemes for identification of constitutive parameters is getting more pronounced as the complexity of the constitutive equations increases, i.e. the number of parameters subject to identification. A general framework for inverse identification of constitutive parameters associated with sheet metal forming is proposed in the article. The inverse problem is solved, through minimization of the least square error between an experimental punch force sampled from a deep drawing and a predicted punch force produced from a coherent finite element model.
Image appraisal for 2D and 3D electromagnetic inversion
Alumbaugh, D.L.; Newman, G.A.
1998-04-01
Linearized methods are presented for appraising image resolution and parameter accuracy in images generated with two and three dimensional non-linear electromagnetic inversion schemes. When direct matrix inversion is employed, the model resolution and model covariance matrices can be directly calculated. The columns of the model resolution matrix are shown to yield empirical estimates of the horizontal and vertical resolution throughout the imaging region. Plotting the square root of the diagonal of the model covariance matrix yields an estimate of how the estimated data noise maps into parameter error. When the conjugate gradient method is employed rather than a direct inversion technique (for example in 3D inversion), an iterative method can be applied to statistically estimate the model covariance matrix, as well as a regularization covariance matrix. The latter estimates the error in the inverted results caused by small variations in the regularization parameter. A method for calculating individual columns of the model resolution matrix using the conjugate gradient method is also developed. Examples of the image analysis techniques are provided on a synthetic cross well EM data set.
2-D acoustic VTI full waveform inversion for CCS monitoring
NASA Astrophysics Data System (ADS)
KIM, S.; Kim, W. K.; Min, D. J.; Jeong, W.; OH, J. W.
2014-12-01
These days many geophysicists have been working not only for oil and gas exploration but also for CO2 monitoring for CCS (Carbon Capture and storage). When CO2 is injected and stored to the target layer, it changes the physical properties of subsurface media like p-wave velocity, density and so on. Seismic method is one of the most widely used geophysical methods for CO2 monitoring, because it can delineate physical properties of subsurface media. To prevent CO2 from leaking out of reservoirs, most target areas require caprocks, and shale often acts as a caprock. However, shale has a strong anisotropic property. Without considering the anisotropic property of subsurface media, interpretations of seismic monitoring data can distort the CO2distribution or movement in the subsurface media. For computational efficiency, seismic data interpretation based on acoustic VTI (Vertical Transversely Isotropic) wave equations has been commonly done although it does not consider the shear waves. To investigate the importance of considering anisotropic properties in acoustic FWI (full waveform inversion) for CO2 monitoring, we compare results obtained by the acoustic VTI FWI with those of the conventional acoustic FWI for isotropic case in the frequency domain. Both methods are based on the node-based finite-element method. Numerical examples show that neglecting anisotropic properties of subsurface media can distort distribution of CO2 and degrade reliability of subsurface image obtained by FWI. Acknowledgements This work was supported by the Human Resources Development program (No. 20134010200510) of the Korea Institute of Energy Technology Evaluation and Planning (KETEP) grant funded by the Korean government Ministry of Trade, Industry, and Energy and by the "Development of Technology for CO2 Marine Geological Storage" grant funded by the Ministry of Oceans and Fisheries of Korea.
A 2-D dynamical model of mesospheric temperature inversions in winter
Hauchecorne, A.; Maillard, A. )
1990-11-01
A 2-D stratospheric and mesospheric dynamical model including drag and diffusion due to gravity wave breaking is used to simulate winter mesospheric temperature inversions similar to those observed by Rayleigh lidar. It is shown that adiabatic heating associated to descending velocities in the mesosphere is the main mechanism involved in the formation of such inversions. Sensitivity tests are performed with the model and confirm this assumption. It is also explained why other previous similar studies with 2-D models did not show mesospheric inversion layers.
NASA Astrophysics Data System (ADS)
Ray, Anandaroop; Key, Kerry; Bodin, Thomas; Myer, David; Constable, Steven
2014-12-01
We apply a reversible-jump Markov chain Monte Carlo method to sample the Bayesian posterior model probability density function of 2-D seafloor resistivity as constrained by marine controlled source electromagnetic data. This density function of earth models conveys information on which parts of the model space are illuminated by the data. Whereas conventional gradient-based inversion approaches require subjective regularization choices to stabilize this highly non-linear and non-unique inverse problem and provide only a single solution with no model uncertainty information, the method we use entirely avoids model regularization. The result of our approach is an ensemble of models that can be visualized and queried to provide meaningful information about the sensitivity of the data to the subsurface, and the level of resolution of model parameters. We represent models in 2-D using a Voronoi cell parametrization. To make the 2-D problem practical, we use a source-receiver common midpoint approximation with 1-D forward modelling. Our algorithm is transdimensional and self-parametrizing where the number of resistivity cells within a 2-D depth section is variable, as are their positions and geometries. Two synthetic studies demonstrate the algorithm's use in the appraisal of a thin, segmented, resistive reservoir which makes for a challenging exploration target. As a demonstration example, we apply our method to survey data collected over the Scarborough gas field on the Northwest Australian shelf.
Efficiency of Pareto joint inversion of 2D geophysical data using global optimization methods
NASA Astrophysics Data System (ADS)
Miernik, Katarzyna; Bogacz, Adrian; Kozubal, Adam; Danek, Tomasz; Wojdyła, Marek
2016-04-01
Pareto joint inversion of two or more sets of data is a promising new tool of modern geophysical exploration. In the first stage of our investigation we created software enabling execution of forward solvers of two geophysical methods (2D magnetotelluric and gravity) as well as inversion with possibility of constraining solution with seismic data. In the algorithm solving MT forward solver Helmholtz's equations, finite element method and Dirichlet's boundary conditions were applied. Gravity forward solver was based on Talwani's algorithm. To limit dimensionality of solution space we decided to describe model as sets of polygons, using Sharp Boundary Interface (SBI) approach. The main inversion engine was created using Particle Swarm Optimization (PSO) algorithm adapted to handle two or more target functions and to prevent acceptance of solutions which are non - realistic or incompatible with Pareto scheme. Each inversion run generates single Pareto solution, which can be added to Pareto Front. The PSO inversion engine was parallelized using OpenMP standard, what enabled execution code for practically unlimited amount of threads at once. Thereby computing time of inversion process was significantly decreased. Furthermore, computing efficiency increases with number of PSO iterations. In this contribution we analyze the efficiency of created software solution taking under consideration details of chosen global optimization engine used as a main joint minimization engine. Additionally we study the scale of possible decrease of computational time caused by different methods of parallelization applied for both forward solvers and inversion algorithm. All tests were done for 2D magnetotelluric and gravity data based on real geological media. Obtained results show that even for relatively simple mid end computational infrastructure proposed solution of inversion problem can be applied in practice and used for real life problems of geophysical inversion and interpretation.
Study on 2D random medium inversion algorithm based on Fuzzy C-means Clustering theory
NASA Astrophysics Data System (ADS)
Xu, Z.; Zhu, P.; Gu, Y.; Yang, X.; Jiang, J.
2015-12-01
Abstract: In seismic exploration for metal deposits, the traditional seismic inversion method based on layered homogeneous medium theory seems difficult to inverse small scale inhomogeneity and spatial variation of the actual medium. The reason is that physical properties of actual medium are more likely random distribution rather than layered. Thus, it is necessary to investigate a random medium inversion algorithm. The velocity of 2D random medium can be described as a function of five parameters: the background velocity (V0), the standard deviation of velocity (σ), the horizontal and vertical autocorrelation lengths (A and B), and the autocorrelation angle (θ). In this study, we propose an inversion algorithm for random medium based on the Fuzzy C-means Clustering (FCM) theory, whose basic idea is that FCM is used to control the inversion process to move forward to the direction we desired by clustering the estimated parameters into groups. Our method can be divided into three steps: firstly, the three parameters (A, B, θ) are estimated from 2D post-stack seismic data using the non-stationary random medium parameter estimation method, and then the estimated parameters are clustered to different groups according to FCM; secondly, the initial random medium model is constructed with clustered groups and the rest two parameters (V0 and σ) obtained from the well logging data; at last, inversion of the random medium are conducted to obtain velocity, impedance and random medium parameters using the Conjugate Gradient Method. The inversion experiments of synthetic seismic data show that the velocity models inverted by our algorithm are close to the real velocity distribution and the boundary of different media can be distinguished clearly.Key words: random medium, inversion, FCM, parameter estimation
2D and 3D separate and joint inversion of airborne ZTEM and ground AMT data: Synthetic model studies
NASA Astrophysics Data System (ADS)
Sasaki, Yutaka; Yi, Myeong-Jong; Choi, Jihyang
2014-05-01
The ZTEM (Z-axis Tipper Electromagnetic) method measures naturally occurring audio-frequency magnetic fields and obtains the tipper function that defines the relationship among the three components of the magnetic field. Since the anomalous tipper responses are caused by the presence of lateral resistivity variations, the ZTEM survey is most suited for detecting and delineating conductive bodies extending to considerable depths, such as graphitic dykes encountered in the exploration of unconformity type uranium deposit. Our simulations shows that inversion of ZTEM data can detect reasonably well multiple conductive dykes placed 1 km apart. One important issue regarding ZTEM inversion is the effect of the initial model, because homogeneous half-space and (1D) layered structures produce no responses. For the 2D model with multiple conductive dykes, the inversion results were useful for locating the dykes even when the initial model was not close to the true background resistivity. For general 3D structures, however, the resolution of the conductive bodies can be reduced considerably depending on the initial model. This is because the tipper magnitudes from 3D conductors are smaller due to boundary charges than the 2D responses. To alleviate this disadvantage of ZTEM surveys, we combined ZTEM and audio-frequency magnetotelluric (AMT) data. Inversion of sparse AMT data was shown to be effective in providing a good initial model for ZTEM inversion. Moreover, simultaneously inverting both data sets led to better results than the sequential approach by enabling to identify structural features that were difficult to resolve from the individual data sets.
Simultaneous elastic parameter inversion in 2-D/3-D TTI medium combined later arrival times
NASA Astrophysics Data System (ADS)
Bai, Chao-ying; Wang, Tao; Yang, Shang-bei; Li, Xing-wang; Huang, Guo-jiao
2016-04-01
Traditional traveltime inversion for anisotropic medium is, in general, based on a "weak" assumption in the anisotropic property, which simplifies both the forward part (ray tracing is performed once only) and the inversion part (a linear inversion solver is possible). But for some real applications, a general (both "weak" and "strong") anisotropic medium should be considered. In such cases, one has to develop a ray tracing algorithm to handle with the general (including "strong") anisotropic medium and also to design a non-linear inversion solver for later tomography. Meanwhile, it is constructive to investigate how much the tomographic resolution can be improved by introducing the later arrivals. For this motivation, we incorporated our newly developed ray tracing algorithm (multistage irregular shortest-path method) for general anisotropic media with a non-linear inversion solver (a damped minimum norm, constrained least squares problem with a conjugate gradient approach) to formulate a non-linear inversion solver for anisotropic medium. This anisotropic traveltime inversion procedure is able to combine the later (reflected) arrival times. Both 2-D/3-D synthetic inversion experiments and comparison tests show that (1) the proposed anisotropic traveltime inversion scheme is able to recover the high contrast anomalies and (2) it is possible to improve the tomographic resolution by introducing the later (reflected) arrivals, but not as expected in the isotropic medium, because the different velocity (qP, qSV and qSH) sensitivities (or derivatives) respective to the different elastic parameters are not the same but are also dependent on the inclination angle.
NASA Astrophysics Data System (ADS)
Randeria, Mohit; Banerjee, Sumilan; Rowland, James
2015-09-01
Most theoretical studies of chiral magnetism, and the resulting spin textures, have focused on 3D systems with broken bulk inversion symmetry, where skyrmions are stabilized by easy-axis anisotropy. In this talk I will describe our results on 2D and quasi-2D systems with broken surface inversion, where we find [1] that skyrmion crystals are much more stable than in 3D, especially for the case of easy-plane anisotropy. These results are of particular interest for thin films, surfaces, and oxide interfaces [2], where broken surface-inversion symmetry and Rashba spin-orbit coupling naturally lead to both the chiral Dzyaloshinskii-Moriya (DM) interaction and to easy-plane compass anisotropy. I will then turn to systems that break both bulk and surface inversion, resulting in two distinct DM terms arising from Dresselhaus and Rashba spin-orbit coupling. I will describe [3] the evolution of the skyrmion structure and of the phase diagram as a function of the ratio of Dresselhaus and Rashba terms, which can be tuned by varying film thickness and strain. [1] S. Banerjee, J. Rowland, O. Erten, and M. Randeria, PRX 4, 031045 (2014). [2] S. Banerjee, O. Erten, and M. Randeria, Nature Phys. 9, 626 (2013). [3] J. Rowland, S. Banerjee and M. Randeria, (unpublished).
2D Seismic Imaging of Elastic Parameters by Frequency Domain Full Waveform Inversion
NASA Astrophysics Data System (ADS)
Brossier, R.; Virieux, J.; Operto, S.
2008-12-01
Thanks to recent advances in parallel computing, full waveform inversion is today a tractable seismic imaging method to reconstruct physical parameters of the earth interior at different scales ranging from the near- surface to the deep crust. We present a massively parallel 2D frequency-domain full-waveform algorithm for imaging visco-elastic media from multi-component seismic data. The forward problem (i.e. the resolution of the frequency-domain 2D PSV elastodynamics equations) is based on low-order Discontinuous Galerkin (DG) method (P0 and/or P1 interpolations). Thanks to triangular unstructured meshes, the DG method allows accurate modeling of both body waves and surface waves in case of complex topography for a discretization of 10 to 15 cells per shear wavelength. The frequency-domain DG system is solved efficiently for multiple sources with the parallel direct solver MUMPS. The local inversion procedure (i.e. minimization of residuals between observed and computed data) is based on the adjoint-state method which allows to efficiently compute the gradient of the objective function. Applying the inversion hierarchically from the low frequencies to the higher ones defines a multiresolution imaging strategy which helps convergence towards the global minimum. In place of expensive Newton algorithm, the combined use of the diagonal terms of the approximate Hessian matrix and optimization algorithms based on quasi-Newton methods (Conjugate Gradient, LBFGS, ...) allows to improve the convergence of the iterative inversion. The distribution of forward problem solutions over processors driven by a mesh partitioning performed by METIS allows to apply most of the inversion in parallel. We shall present the main features of the parallel modeling/inversion algorithm, assess its scalability and illustrate its performances with realistic synthetic case studies.
2-D traveltime and waveform inversion for improved seismic imaging: Naga Thrust and Fold Belt, India
NASA Astrophysics Data System (ADS)
Jaiswal, Priyank; Zelt, Colin A.; Bally, Albert W.; Dasgupta, Rahul
2008-05-01
Exploration along the Naga Thrust and Fold Belt in the Assam province of Northeast India encounters geological as well as logistic challenges. Drilling for hydrocarbons, traditionally guided by surface manifestations of the Naga thrust fault, faces additional challenges in the northeast where the thrust fault gradually deepens leaving subtle surface expressions. In such an area, multichannel 2-D seismic data were collected along a line perpendicular to the trend of the thrust belt. The data have a moderate signal-to-noise ratio and suffer from ground roll and other acquisition-related noise. In addition to data quality, the complex geology of the thrust belt limits the ability of conventional seismic processing to yield a reliable velocity model which in turn leads to poor subsurface image. In this paper, we demonstrate the application of traveltime and waveform inversion as supplements to conventional seismic imaging and interpretation processes. Both traveltime and waveform inversion utilize the first arrivals that are typically discarded during conventional seismic processing. As a first step, a smooth velocity model with long wavelength characteristics of the subsurface is estimated through inversion of the first-arrival traveltimes. This velocity model is then used to obtain a Kirchhoff pre-stack depth-migrated image which in turn is used for the interpretation of the fault. Waveform inversion is applied to the central part of the seismic line to a depth of ~1 km where the quality of the migrated image is poor. Waveform inversion is performed in the frequency domain over a series of iterations, proceeding from low to high frequency (11-19 Hz) using the velocity model from traveltime inversion as the starting model. In the end, the pre-stack depth-migrated image and the waveform inversion model are jointly interpreted. This study demonstrates that a combination of traveltime and waveform inversion with Kirchhoff pre-stack depth migration is a promising approach
Interpretation of gravity data using 2-D continuous wavelet transformation and 3-D inverse modeling
NASA Astrophysics Data System (ADS)
Roshandel Kahoo, Amin; Nejati Kalateh, Ali; Salajegheh, Farshad
2015-10-01
Recently the continuous wavelet transform has been proposed for interpretation of potential field anomalies. In this paper, we introduced a 2D wavelet based method that uses a new mother wavelet for determination of the location and the depth to the top and base of gravity anomaly. The new wavelet is the first horizontal derivatives of gravity anomaly of a buried cube with unit dimensions. The effectiveness of the proposed method is compared with Li and Oldenburg inversion algorithm and is demonstrated with synthetics and real gravity data. The real gravity data is taken over the Mobrun massive sulfide ore body in Noranda, Quebec, Canada. The obtained results of the 2D wavelet based algorithm and Li and Oldenburg inversion on the Mobrun ore body had desired similarities to the drill-hole depth information. In all of the inversion algorithms the model non-uniqueness is the challenging problem. Proposed method is based on a simple theory and there is no model non-uniqueness on it.
Skyrmions in quasi-2D chiral magnets with broken bulk and surface inversion symmetry
NASA Astrophysics Data System (ADS)
Rowland, James; Banerjee, Sumilan; Randeria, Mohit
2015-03-01
Most theoretical studies of skyrmions have focused on chiral magnets with broken bulk inversion symmetry, stabilized by easy-axis anisotropy. Recently, we considered 2D systems with broken surface inversion and showed that skyrmion crystals are more stable than in 3D, pointing out the importance of easy-plane anisotropy. In the present work we investigate quasi-2D systems which break both bulk and surface inversion symmetry. The Landau-Ginzburg free energy functional thus contains two Dzyloshinskii-Moriya terms of strength DD and DR arising from Dresselhaus and Rashba spin-orbit coupling respectively. We trace the evolution of the phase diagram as DD /DR is varied, and find that skyrmions are increasingly destabilized with respect to the cone phase as DD increases relative to DR. We find an evolution from vortex-like skyrmions in the pure Dresselhaus limit to hedgehog-like skyrmions in the pure Rashba limit. We discuss the relevance of these results to existing experiments and the prospects of tuning the ratio of Dresselhaus and Rashba spin-orbit coupling via film thickness and strain. Supported by NSF DMR-1410364 (J.R. and M.R.) and DOE-BES DE-SC0005035 (S.B.)
Takes Electric or Magnetic field data through Inversion process a 2D Distributon
Energy Science and Technology Software Center (ESTSC)
2008-05-01
Program images 2D distributions in electrical conductivity for geophysical applications. The program can treat surface based and cross well measurement geometries, including inductive and grounded source antennas in the quasi-static limit. The algorithm using Krylov iterative methods to solve for the predicted data and model sensitivities. The model update is achieved using a Gauss-newton optimization process for stability. A new line search capability is now included in the algorithm to insure global convergence of themore » inversion iteration.« less
NASA Astrophysics Data System (ADS)
Li, Jinghe; Song, Linping; Liu, Qing Huo
2016-02-01
A simultaneous multiple frequency contrast source inversion (CSI) method is applied to reconstructing hydrocarbon reservoir targets in a complex multilayered medium in two dimensions. It simulates the effects of a salt dome sedimentary formation in the context of reservoir monitoring. In this method, the stabilized biconjugate-gradient fast Fourier transform (BCGS-FFT) algorithm is applied as a fast solver for the 2D volume integral equation for the forward computation. The inversion technique with CSI combines the efficient FFT algorithm to speed up the matrix-vector multiplication and the stable convergence of the simultaneous multiple frequency CSI in the iteration process. As a result, this method is capable of making quantitative conductivity image reconstruction effectively for large-scale electromagnetic oil exploration problems, including the vertical electromagnetic profiling (VEP) survey investigated here. A number of numerical examples have been demonstrated to validate the effectiveness and capacity of the simultaneous multiple frequency CSI method for a limited array view in VEP.
2013-01-01
Background Cardiovascular-MR (CMR) is the gold standard for quantifying myocardial infarction using late gadolinium enhancement (LGE) technique. Both 2D- and 3D-LGE-sequences are used in clinical practise and in clinical and experimental studies for infarct quantification. Therefore the aim of this study was to investigate if image acquisitions with 2D- and 3D-LGE show the same infarct size in patients and ex vivo. Methods Twenty-six patients with previous myocardial infarction who underwent a CMR scan were included. Images were acquired 10-20 minutes after an injection of 0.2 mmol/kg gadolinium-based contrast agent. Two LGE-sequences, 3D-inversion recovery (IR) and 2D-phase-sensitive (PS) IR, were used in all patients to quantify infarction size. Furthermore, six pigs with reperfused infarction in the left anterior descending artery (40 minutes occlusion and 4 hours of reperfusion) were scanned with 2D- and 3D-LGE ex vivo. A high resolution T1-sequence was used as reference for the infarct quantification ex vivo. Spearman’s rank-order correlation, Wilcoxon matched pairs test and bias according to Bland-Altman was used for comparison of infarct size with different LGE-sequences. Results There was no significant difference between the 2D- and 3D-LGE sequence in left ventricular mass (LVM) (2D: 115 ± 25 g; 3D: 117 ± 24 g: p = 0.35). Infarct size in vivo using 2D- and 3D-LGE showed high correlation and low bias for both LGE-sequences both in absolute volume of infarct (r = 0.97, bias 0.47 ± 2.1 ml) and infarct size as part of LVM (r = 0.94, bias 0.16 ± 2.0%). The 2D- and 3D-LGE-sequences ex vivo correlated well (r = 0.93, bias 0.67 ± 2.4%) for infarct size as part of the LVM. The IR LGE-sequences overestimated infarct size as part of the LVM ex vivo compared to the high resolution T1-sequence (bias 6.7 ± 3.0%, 7.3 ± 2.7% for 2D-PSIR and 3D-IR respectively, p < 0.05 for both). Conclusions Infarct quantification with
Calibration of an Ultrasound Tomography System for Medical Imaging with 2D Contrast-Source Inversion
NASA Astrophysics Data System (ADS)
Faucher, Gabriel Paul
This dissertation describes two possible methods for the calibration of an ultrasound tomography system developed at University of Manitoba's Electromagnetic Imaging Laboratory for imaging with the contrast-source inversion algorithm. The calibration techniques are adapted from existing procedures employed for microwave tomography. A theoretical model of these calibration principles is developed in order to provide a rationale for the effectiveness of the proposed procedures. The applicability of such an imaging algorithm and calibration methods in the context of ultrasound are discussed. Also presented are 2D and 3D finite-difference time-domain update equations for the simulation of acoustic wave propagation in inhomogeneous media. Details regarding the application of an absorbing boundary-condition, point-source modelling and the treatment of penetrable objects are included in this document.
Structure-approximating inverse protein folding problem in the 2D HP model.
Gupta, Arvind; Manuch, Ján; Stacho, Ladislav
2005-12-01
The inverse protein folding problem is that of designing an amino acid sequence which has a particular native protein fold. This problem arises in drug design where a particular structure is necessary to ensure proper protein-protein interactions. In this paper, we show that in the 2D HP model of Dill it is possible to solve this problem for a broad class of structures. These structures can be used to closely approximate any given structure. One of the most important properties of a good protein (in drug design) is its stability--the aptitude not to fold simultaneously into other structures. We show that for a number of basic structures, our sequences have a unique fold. PMID:16379538
Implicit adaptive mesh refinement for 2D reduced resistive magnetohydrodynamics
NASA Astrophysics Data System (ADS)
Philip, Bobby; Chacón, Luis; Pernice, Michael
2008-10-01
An implicit structured adaptive mesh refinement (SAMR) solver for 2D reduced magnetohydrodynamics (MHD) is described. The time-implicit discretization is able to step over fast normal modes, while the spatial adaptivity resolves thin, dynamically evolving features. A Jacobian-free Newton-Krylov method is used for the nonlinear solver engine. For preconditioning, we have extended the optimal "physics-based" approach developed in [L. Chacón, D.A. Knoll, J.M. Finn, An implicit, nonlinear reduced resistive MHD solver, J. Comput. Phys. 178 (2002) 15-36] (which employed multigrid solver technology in the preconditioner for scalability) to SAMR grids using the well-known Fast Adaptive Composite grid (FAC) method [S. McCormick, Multilevel Adaptive Methods for Partial Differential Equations, SIAM, Philadelphia, PA, 1989]. A grid convergence study demonstrates that the solver performance is independent of the number of grid levels and only depends on the finest resolution considered, and that it scales well with grid refinement. The study of error generation and propagation in our SAMR implementation demonstrates that high-order (cubic) interpolation during regridding, combined with a robustly damping second-order temporal scheme such as BDF2, is required to minimize impact of grid errors at coarse-fine interfaces on the overall error of the computation for this MHD application. We also demonstrate that our implementation features the desired property that the overall numerical error is dependent only on the finest resolution level considered, and not on the base-grid resolution or on the number of refinement levels present during the simulation. We demonstrate the effectiveness of the tool on several challenging problems.
Electrical resistivity tomography applied to a complex lava dome: 2D and 3D models comparison
NASA Astrophysics Data System (ADS)
Portal, Angélie; Fargier, Yannick; Lénat, Jean-François; Labazuy, Philippe
2015-04-01
The study of volcanic domes growth (e.g. St. Helens, Unzen, Montserrat) shows that it is often characterized by a succession of extrusion phases, dome explosions and collapse events. Lava dome eruptive activity may last from days to decades. Therefore, their internal structure, at the end of the eruption, is complex and includes massive extrusions and lava lobes, talus and pyroclastic deposits as well as hydrothermal alteration. The electrical resistivity tomography (ERT) method, initially developed for environmental and engineering exploration, is now commonly used for volcano structure imaging. Because a large range of resistivity values is often observed in volcanic environments, the method is well suited to study the internal structure of volcanic edifices. We performed an ERT survey on an 11ka years old trachytic lava dome, the Puy de Dôme volcano (French Massif Central). The analysis of a recent high resolution DEM (LiDAR 0.5 m), as well as other geophysical data, strongly suggest that the Puy de Dôme is a composite dome. 11 ERT profiles have been carried out, both at the scale of the entire dome (base diameter of ~2 km and height of 400 m) on the one hand, and at a smaller scale on the summit part on the other hand. Each profile is composed of 64 electrodes. Three different electrode spacing have been used depending on the study area (35 m for the entire dome, 10 m and 5 m for its summit part). Some profiles were performed with half-length roll-along acquisitions, in order to keep a good trade-off between depth of investigation and resolution. Both Wenner-alpha and Wenner-Schlumberger protocols were used. 2-D models of the electrical resistivity distribution were computed using RES2DINV software. In order to constrain inversion models interpretation, the depth of investigation (DOI) method was applied to those results. It aims to compute a sensitivity index on inversion results, illustrating how the data influence the model and constraining models
Even more inversion formulas for the 2D Radon Transform of functions of compact and convex support.
Mohamed, M-S Ould; Mennessier, C; Clackdoyle, R
2007-01-01
In 2004, Clackdoyle and Noo published a class of inversion formulas for the 2D Radon Transform which depends on the known radius of support of the unknown function. In this work, we extend this class of inversion formulas from functions of circular support to functions with any compact and convex support. We point out the potential benefits of these new inversion formulas in the context of reconstruction from truncated projections. A preliminary implementation of these new inversion formulas is also presented. PMID:18002982
A 2D forward and inverse code for streaming potential problems
NASA Astrophysics Data System (ADS)
Soueid Ahmed, A.; Jardani, A.; Revil, A.
2013-12-01
The self-potential method corresponds to the passive measurement of the electrical field in response to the occurrence of natural sources of current in the ground. One of these sources corresponds to the streaming current associated with the flow of the groundwater. We can therefore apply the self- potential method to recover non-intrusively some information regarding the groundwater flow. We first solve the forward problem starting with the solution of the groundwater flow problem, then computing the source current density, and finally solving a Poisson equation for the electrical potential. We use the finite-element method to solve the relevant partial differential equations. In order to reduce the number of (petrophysical) model parameters required to solve the forward problem, we introduced an effective charge density tensor of the pore water, which can be determined directly from the permeability tensor for neutral pore waters. The second aspect of our work concerns the inversion of the self-potential data using Tikhonov regularization with smoothness and weighting depth constraints. This approach accounts for the distribution of the electrical resistivity, which can be independently and approximately determined from electrical resistivity tomography. A numerical code, SP2DINV, has been implemented in Matlab to perform both the forward and inverse modeling. Three synthetic case studies are discussed.
IP4DI: A software for time-lapse 2D/3D DC-resistivity and induced polarization tomography
NASA Astrophysics Data System (ADS)
Karaoulis, M.; Revil, A.; Tsourlos, P.; Werkema, D. D.; Minsley, B. J.
2013-04-01
We propose a 2D/3D forward modelling and inversion package to invert direct current (DC)-resistivity, time-domain induced polarization (TDIP), and frequency-domain induced polarization (FDIP) data. Each cell used for the discretization of the 2D/3D problems is characterized by a DC-resistivity value and a chargeability or complex conductivity for TDIP/FDIP problems, respectively. The governing elliptic partial differential equations are solved with the finite element method, which can be applied for both real and complex numbers. The inversion can be performed either for a single snapshot of data or for a sequence of snapshots in order to monitor a dynamic process such as a salt tracer test. For the time-lapse inversion, we have developed an active time constrained (ATC) approach that is very efficient in filtering out noise in the data that is not correlated over time. The forward algorithm is benchmarked with simple analytical solutions. The inversion package IP4DI is benchmarked with three tests, two including simple geometries. The last one corresponds to a time-lapse resistivity problem for cross-well tomography during enhanced oil recovery. The algorithms are based on MATLAB® code package and a graphical user interface (GUI).
A 2D inverse problem of predicting boiling heat transfer in a long fin
NASA Astrophysics Data System (ADS)
Orzechowski, Tadeusz
2015-12-01
A method for the determination of local values of the heat transfer coefficient on non-isothermal surfaces was analyzed on the example of a long smooth-surfaced fin made of aluminium. On the basis of the experimental data, two cases were taken into consideration: one-dimensional model for Bi < 0.1 and two-dimensional model for thicker elements. In the case when the drop in temperature over the thickness could be omitted, the rejected local values of heat fluxes were calculated from the integral of the equation describing temperature distribution on the fin. The corresponding boiling curve was plotted on the basis of temperature gradient distribution as a function of superheat. For thicker specimens, where Bi > 0.1, the problem was modelled using a 2-D heat conduction equation, for which the boundary conditions were posed on the surface observed with a thermovision camera. The ill-conditioned inverse problem was solved using a method of heat polynomials, which required validation.
An analytical approach to estimate the number of small scatterers in 2D inverse scattering problems
NASA Astrophysics Data System (ADS)
Fazli, Roohallah; Nakhkash, Mansor
2012-07-01
This paper presents an analytical method to estimate the location and number of actual small targets in 2D inverse scattering problems. This method is motivated from the exact maximum likelihood estimation of signal parameters in white Gaussian noise for the linear data model. In the first stage, the method uses the MUSIC algorithm to acquire all possible target locations and in the next stage, it employs an analytical formula that works as a spatial filter to determine which target locations are associated to the actual ones. The ability of the method is examined for both the Born and multiple scattering cases and for the cases of well-resolved and non-resolved targets. Many numerical simulations using both the coincident and non-coincident arrays demonstrate that the proposed method can detect the number of actual targets even in the case of very noisy data and when the targets are closely located. Using the experimental microwave data sets, we further show that this method is successful in specifying the number of small inclusions.
Modeling Coastal Salinity in Quasi 2D and 3D Using a DUALEM-421 and Inversion Software.
Davies, Gareth; Huang, Jingyi; Monteiro Santos, Fernando Acacio; Triantafilis, John
2015-01-01
Rising sea levels, owing to climate change, are a threat to fresh water coastal aquifers. This is because saline intrusions are caused by increases and intensification of medium-large scale influences including sea level rise, wave climate, tidal cycles, and shifts in beach morphology. Methods are therefore required to understand the dynamics of these interactions. While traditional borehole and galvanic contact resistivity (GCR) techniques have been successful they are time-consuming. Alternatively, frequency-domain electromagnetic (FEM) induction is potentially useful as physical contact with the ground is not required. A DUALEM-421 and EM4Soil inversion software package are used to develop a quasi two- (2D) and quasi three-dimensional (3D) electromagnetic conductivity images (EMCI) across Long Reef Beach located north of Sydney Harbour, New South Wales, Australia. The quasi 2D models discern: the dry sand (<10 mS/m) associated with the incipient dune; sand with fresh water (10 to 20 mS/m); mixing of fresh and saline water (20 to 500 mS/m), and; saline sand of varying moisture (more than 500 mS/m). The quasi 3D EMCIs generated for low and high tides suggest that daily tidal cycles do not have a significant effect on local groundwater salinity. Instead, the saline intrusion is most likely influenced by medium-large scale drivers including local wave climate and morphology along this wave-dominated beach. Further research is required to elucidate the influence of spring-neap tidal cycles, contrasting beach morphological states and sea level rise. PMID:25053423
2D numerical simulation of the resistive reconnection layer
D. A. Uzdensky; R. M. Kulsrud
2000-07-21
In this paper the authors present a two-dimensional numerical simulation of a reconnection current layer in incompressible resistive magnetohydrodynamics with uniform resistivity in the limit of very large Lundquist numbers. They use realistic boundary conditions derived consistently from the outside magnetic field, and they also take into account the effect of the backpressure from flow into the separatrix region. They find that within a few Alfven times the system reaches a steady state consistent with the Sweet-Parker model, even if the initial state is Petschek-like.
2D Numerical Simulation of the Resistive Reconnection Layer
Kulsrud, R.M.; Uzdensky, D.A.
1999-03-01
In this paper we present a two-dimensional numerical simulation of a reconnection current layer in incompressible resistive magnetohydrodynamics with uniform resistivity in the limit of very large Lundquist numbers. We use realistic boundary conditions derived consistently from the outside magnetic field, and we also take into account the effect of the back pressure from flow into the separatrix region. We find that within a few Alfvén times the system reaches a steady state consistent with the Sweet-Parker model, even if the initial state is Petschek-like.
2-D resistivity distribution in the seismo-genic zone of the 2011 Tohoku earthquake (M9.0)
NASA Astrophysics Data System (ADS)
Ichihara, H.; Kasaya, T.; Baba, K.
2012-12-01
We discussed 2-D resistivity distribution around the rupture zone of the 2011 Tohoku earthquake (M9.0) in order to clarify structural properties that enabled such a large fault slip. Electro-magnetic time series were collected using ocean bottom electro-magnetometers (OBEMs) along a survey line in the latitude 38N between 2009 and 2011. The MT responses show the out of the quadrant phases in TE mode. The out of quadrant phases were explained by coastal effect inducing strong distortion especially in the magnetic field [Key and Constable, 2011] and indicate strong resistivity contrast between seawater and crustal area. Then we carried out 2-D resistivity inversion after the correction of 3-D bathymetric effect. The inverted resistivity model entirely consists of conductive surface layer (< 3 ohm-m) and subsequence resistive area (> 300 ohm-m). While the main rupture zone shows high resistivity (> 3000 ohm-m), resistivity near the trench and the deeper area along the plate boundary is relatively low. The former conductive zone probably reflects thick sediments containing large amount of fluid. The lower conductive zone possibly indicates connected pore fluids due to dehydration. Thus the main rupture zone seems to be located on a resistive area sandwiched by conductive zones. It implies that heterogeneous distribution of fluid controlled fault activity of the 2011 earthquake.
NASA Astrophysics Data System (ADS)
Portal, Angélie; Fargier, Yannick; Lénat, Jean-François; Labazuy, Philippe
2016-04-01
The electrical resistivity tomography (ERT) method, initially developed for environmental and engineering exploration, is now commonly used for geological structures imaging. Such structures can present complex characteristics that conventional 2D inversion processes cannot perfectly integrate. Here we present a new 3D inversion algorithm named EResI, firstly developed for levee investigation, and presently applied to the study of a complex lava dome (the Puy de Dôme volcano, France). EResI algorithm is based on a conventional regularized Gauss-Newton inversion scheme and a 3D non-structured discretization of the model (double grid method based on tetrahedrons). This discretization allows to accurately model the topography of investigated structure (without a mesh deformation procedure) and also permits a precise location of the electrodes. Moreover, we demonstrate that a complete 3D unstructured discretization limits the number of inversion cells and is better adapted to the resolution capacity of tomography than a structured discretization. This study shows that a 3D inversion with a non-structured parametrization has some advantages compared to classical 2D inversions. The first advantage comes from the fact that a 2D inversion leads to artefacts due to 3D effects (3D topography, 3D internal resistivity). The second advantage comes from the fact that the capacity to experimentally align electrodes along an axis (for 2D surveys) depends on the constrains on the field (topography...). In this case, a 2D assumption induced by 2.5D inversion software prevents its capacity to model electrodes outside this axis leading to artefacts in the inversion result. The last limitation comes from the use of mesh deformation techniques used to accurately model the topography in 2D softwares. This technique used for structured discretization (Res2dinv) is prohibed for strong topography (>60 %) and leads to a small computational errors. A wide geophysical survey was carried out
NASA Astrophysics Data System (ADS)
Kiflu, H. G.; Kruse, S. E.; Harro, D.; Loke, M. H.; Wilkinson, P. B.
2013-12-01
Electrical resistivity tomography is commonly used to identify geologic features associated with sinkhole formation. In covered karst terrain, however, it can be difficult to resolve the depth to top of limestone with this method. This is due to the fact that array lengths, and hence depth of resolution, are often limited by residential or commercial lot dimensions in urban environments. Furthermore, the sediments mantling the limestone are often clay-rich and highly conductive. The resistivity method has limited sensitivity to resistive zones beneath conductive zones. This sensitivity can be improved significantly with electrodes implanted at depth in the cover sediments near the top of limestone. An array of deep electrodes is installed with direct push technology in the karst cover. When combined with a surface array in which each surface electrode is underlain by a deep electrode, the array geometry is similar to a borehole array turned on its side. This method, called the Multi-Electrode Resistivity Implant Technique (MERIT), offers the promise of significantly improved resolution of epikarst and cover collapse development zones in the overlying sediment, the limestone or at the sediment-bedrock interface in heterogeneous karst environments. With a non-traditional array design, the question of optimal array geometries arises. Optimizing array geometries is complicated by the fact that many plausible 4-electrode readings will produce negative apparent resistivity values, even in homogeneous terrain. Negative apparent resistivities cannot be used in inversions based on the logarithm of the apparent resistivity. New algorithms for seeking optimal array geometries have been developed by modifying the 'Compare R' method of Wilkinson and Loke. The optimized arrays show significantly improved resolution over basic arrays adapted from traditional 2D surface geometries. Several MERIT case study surveys have been conducted in covered karst in west-central Florida, with
Burgess, Ian B; Aizenberg, Joanna; Lončar, Marko
2013-12-01
Structural hierarchy and complex 3D architecture are characteristics of biological photonic designs that are challenging to reproduce in synthetic materials. Top-down lithography allows for designer patterning of arbitrary shapes, but is largely restricted to planar 2D structures. Self-assembly techniques facilitate easy fabrication of 3D photonic crystals, but controllable defect-integration is difficult. In this paper we combine the advantages of top-down and bottom-up fabrication, developing two techniques to deposit 2D-lithographically-patterned planar layers on top of or in between inverse-opal 3D photonic crystals and creating hierarchical structures that resemble the architecture of the bright green wing scales of the butterfly, Parides sesostris. These fabrication procedures, combining advantages of both top-down and bottom-up fabrication, may prove useful in the development of omnidirectional coloration elements and 3D-2D photonic crystal devices. PMID:24263010
Time lapse inversion of 2D ERT data for monitoring river water infiltration
NASA Astrophysics Data System (ADS)
Wallin, E. L.; Johnson, T. C.; Greenwood, W. J.
2011-12-01
Uranium transport in the 300 area is driven by both the chemical and physical effects of stage fluctuations in the Columbia River and resulting river water, ground water interaction. Because river water is less conductive than groundwater, it serves as a natural tracer that can be imaged using surface ERT. We've monitored 4 lines for 4 months over the high stage spring runoff interval to identify preferred flowpaths for river water intrusion. The four lines overlay former waste disposal sites including the processing ponds and sanitary leach trenches. We have used this ERT data set to investigate two methods of time lapse inversion, sequential and all-at-once. Each technique is a model difference approach as opposed to a data difference inversion approach. Both use a regularized inversion with model constraints that regularize spatially and temporally. For the sequential inversion, the starting and reference models are taken from the previous inversion. Both starting and reference models are taken from the first inversion for the all-at-once method. In either case an inversion is triggered if data misfit from the starting model exceeds the chi-squared convergence criteria. It was found that starting with a relatively smooth model provided better visualization of temporal conductivity changes when inverting all data sets with the same initial model, while an initial model exhibiting smaller data misfit may be used successfully as the starting point for sequential inversion. Inland conductivity changes within model cells were found to be highly correlated with river stage, and when paired with the characterization model, provide evidence of waste trenches, the processing pond, as well as the existence of a paleo-channel incised into the Ringold Formation and dipping structures on the Hanford-Ringold contact that provide preferred pathways for river water intrusion.
NASA Astrophysics Data System (ADS)
Chuang, Hsun-Jen; Chamlagain, Bhim; Koehler, Michael; Perera, Meeghage Madusanka; Yan, Jiaqiang; Mandrus, David; Tománek, David; Zhou, Zhixian
2016-03-01
We report a new strategy for fabricating 2D/2D low-resistance ohmic contacts for a variety of transition metal dichalcogenides (TMDs) using van der Waals assembly of substitutionally doped TMDs as drain/source contacts and TMDs with no intentional doping as channel materials. We demonstrate that few-layer WSe2 field-effect transistors (FETs) with 2D/2D contacts exhibit low contact resistances of ~ 0.3 k ohm.um, high on/off ratios up to > 109, and high drive currents exceeding 320 uA um-1. These favorable characteristics are combined with a two-terminal field-effect hole mobility ~ 2x102 cm2 V-1 s-1 at room temperature, which increases to >2x103 cm2 V-1 s-1 at cryogenic temperatures. We observe a similar performance also in MoS2 and MoSe2 FETs with 2D/2D drain and source contacts. The 2D/2D low-resistance ohmic contacts presented here represent a new device paradigm that overcomes a significant bottleneck in the performance of TMDs and a wide variety of other 2D materials as the channel materials in post-silicon electronics.
Chuang, Hsun-Jen; Chamlagain, Bhim; Koehler, Michael; Perera, Meeghage Madusanka; Yan, Jiaqiang; Mandrus, David; Tománek, David; Zhou, Zhixian
2016-03-01
We report a new strategy for fabricating 2D/2D low-resistance ohmic contacts for a variety of transition metal dichalcogenides (TMDs) using van der Waals assembly of substitutionally doped TMDs as drain/source contacts and TMDs with no intentional doping as channel materials. We demonstrate that few-layer WSe2 field-effect transistors (FETs) with 2D/2D contacts exhibit low contact resistances of ∼0.3 kΩ μm, high on/off ratios up to >10(9), and high drive currents exceeding 320 μA μm(-1). These favorable characteristics are combined with a two-terminal field-effect hole mobility μFE ≈ 2 × 10(2) cm(2) V(-1) s(-1) at room temperature, which increases to >2 × 10(3) cm(2) V(-1) s(-1) at cryogenic temperatures. We observe a similar performance also in MoS2 and MoSe2 FETs with 2D/2D drain and source contacts. The 2D/2D low-resistance ohmic contacts presented here represent a new device paradigm that overcomes a significant bottleneck in the performance of TMDs and a wide variety of other 2D materials as the channel materials in postsilicon electronics. PMID:26844954
NASA Astrophysics Data System (ADS)
Martin, T. P.; Schultz, A.
2012-12-01
understand the existing baseline subsurface resistivity structure at the Newberry site prior to well stimulation, magnetotelluric (MT) data will be collected in late July 2012 using two long period (1 Hz sampling) Narod Geophysics NIMS MT instruments along with EarthScope MT data aligned in a ~210 km long N-S profile centered on the stimulation zone. A 2-D inverse model will be obtained from the MT data set. The goal of this investigation is to determine the variations in the electrical resistivity in the mid-to-lower crust beneath the western flank of the caldera, providing a deeper view of putative heat sources than existing studies in this
Numerical solution of 2D-vector tomography problem using the method of approximate inverse
NASA Astrophysics Data System (ADS)
Svetov, Ivan; Maltseva, Svetlana; Polyakova, Anna
2016-08-01
We propose a numerical solution of reconstruction problem of a two-dimensional vector field in a unit disk from the known values of the longitudinal and transverse ray transforms. The algorithm is based on the method of approximate inverse. Numerical simulations confirm that the proposed method yields good results of reconstruction of vector fields.
A multiple-scale Pascal polynomial for 2D Stokes and inverse Cauchy-Stokes problems
NASA Astrophysics Data System (ADS)
Liu, Chein-Shan; Young, D. L.
2016-05-01
The polynomial expansion method is a useful tool for solving both the direct and inverse Stokes problems, which together with the pointwise collocation technique is easy to derive the algebraic equations for satisfying the Stokes differential equations and the specified boundary conditions. In this paper we propose two novel numerical algorithms, based on a third-first order system and a third-third order system, to solve the direct and the inverse Cauchy problems in Stokes flows by developing a multiple-scale Pascal polynomial method, of which the scales are determined a priori by the collocation points. To assess the performance through numerical experiments, we find that the multiple-scale Pascal polynomial expansion method (MSPEM) is accurate and stable against large noise.
Three-dimensional induced polarization data inversion for complex resistivity
Commer, M.; Newman, G.A.; Williams, K.H.; Hubbard, S.S.
2011-03-15
The conductive and capacitive material properties of the subsurface can be quantified through the frequency-dependent complex resistivity. However, the routine three-dimensional (3D) interpretation of voluminous induced polarization (IP) data sets still poses a challenge due to large computational demands and solution nonuniqueness. We have developed a flexible methodology for 3D (spectral) IP data inversion. Our inversion algorithm is adapted from a frequency-domain electromagnetic (EM) inversion method primarily developed for large-scale hydrocarbon and geothermal energy exploration purposes. The method has proven to be efficient by implementing the nonlinear conjugate gradient method with hierarchical parallelism and by using an optimal finite-difference forward modeling mesh design scheme. The method allows for a large range of survey scales, providing a tool for both exploration and environmental applications. We experimented with an image focusing technique to improve the poor depth resolution of surface data sets with small survey spreads. The algorithm's underlying forward modeling operator properly accounts for EM coupling effects; thus, traditionally used EM coupling correction procedures are not needed. The methodology was applied to both synthetic and field data. We tested the benefit of directly inverting EM coupling contaminated data using a synthetic large-scale exploration data set. Afterward, we further tested the monitoring capability of our method by inverting time-lapse data from an environmental remediation experiment near Rifle, Colorado. Similar trends observed in both our solution and another 2D inversion were in accordance with previous findings about the IP effects due to subsurface microbial activity.
2D Electrical Resistivity Tomography surveys optimisation of the solutes transports in porous media.
NASA Astrophysics Data System (ADS)
Lekmine, G.; Pessel, M.; Auradou, H.
2009-04-01
Electrical resistivity tomography applied in borehole or cross-borehole is a method often used to follow the invasion process of pollutant [Daily, 1991]. The aim of this work is to test experimentally the electrode arrays and inversion process used to obtain a spatial representation of tracer propagation in porous media. Experiments were conducted in a plexiglas container with glass beads of 166 microns in diameter. The height of the container is 275 mm, its width 85 mm and its thickness 10 mm. 21 electrodes, equally spaced, are placed along each of the lateral sides of the porous medium : these electrodes are used to perform the electrical measurements. The porous medium is lightened from behind and a video camera records the propagation of the fluids. The fluid containing the tracer (i.e the pollutant) is a water solution containing a small amount of dye together with NaCl (0.5g/l up to 2.0g/l). The medium is first saturated by a water solution containing a slight concentration of NaCl so that its density is smaller than the injected fluid. An upward flow is first established, then the denser fluid is injected at the bottom and over the full width of the medium. In this way, the flow is stabilized by gravity avoiding the development of unstable fingers. Still, the fluids are miscible and a mixing front develops during the flow: in the present study, both the determination by optical and electrical imaging of the mean position of the front and its width are of interest. The comparison of the two techniques allows to study the ability of the inversion process to quantify the solute transport. The electrical measurements are acquired by a standard multi electrode system (IRIS Instruments) and the data are inverted with the Res2Dinv software which models the 2D distribution of conductivity contrasts. The obtained bulk conductivity can be related through Archie's law to fluid conductivity by the porosity and the cementation factor which have been experimentally
Impact of the Born approximation on the estimation error in 2D inverse scattering
NASA Astrophysics Data System (ADS)
Diong, M. L.; Roueff, A.; Lasaygues, P.; Litman, A.
2016-06-01
The aim is to quantify the impact of the Born approximation on the estimation error for a simple inverse scattering problem, while taking into account the noise measurement features. The proposed method consists of comparing two estimation errors: the error obtained with the Born approximation and the error obtained without it. The first error is characterized by the mean and variance of the maximum likelihood estimator, which are straightforward to compute with the Born approximation because the corresponding estimator is linear. The second error is evaluated with the Cramer–Rao bound (CRB). The CRB is a lower bound on the variance of unbiased estimators and thus does not depend on the choice of the estimation method. Beyond the conclusions that will be established under the Born approximation, this study lays out a general methodology that can be generalized to any other approximation.
Precision analysis based on Cramer-Rao bound for 2D acoustics and electromagnetic inverse scattering
NASA Astrophysics Data System (ADS)
Diong, M. L.; Roueff, A.; Lasaygues, P.; Litman, A.
2015-07-01
The aim of the present article is to predict the expected precision quantitatively in inverse scattering when one tries to determine the intrinsic properties of a given target from its scattered field. To conduct such a study, we analyze the precision of contrast estimators with the Cramer-Rao bound (CRB) when the target is homogeneous, infinitely-long and with a circular cross-section and with an additive complex circular gaussian noise at the receivers. An unified framework is derived to handle acoustic or electromagnetic imaging configurations equally. Numerical tests enable to quantitatively appraise the variations of the CRB with respect to the considered physical situation parameters: transmission/reflexion, antennas arrangement, weak/strong scatterers, noise level and source frequency. These analyzes are performed with respect to the real and imaginary parts of the contrast.
Inversion Method for Evapotranspiration from 2D Measurements of Water Vapor Concentration
NASA Astrophysics Data System (ADS)
Eichinger, W. E.; Barnhart, B. L.; Plenner, S.; Thompson, A.; Winsky, B.
2011-12-01
A method has been developed to use 2D measurements of water vapor concentration to make spatial estimates of evapotranspiration. The method is based on theoretical solutions to the conservation of water vapor equation (similar to Sutton, Sienfeld) and is related to methods used to generate footprints. The method uses a diffusion coefficient to close the conservation equation and assumes uniformity of the surface with respect to turbulent processes. This solution allows the development of a transfer coefficient between each location at the surface and each of the water vapor measurements. The problem has been constrained to two dimensions for development, but is easily expanded to three. The method assumes water vapor concentration measurements in a plane at arbitrary positions above the field. Evapotranspiration is allowed to vary in increments across the field. This geometry and the transfer coefficients allow the water vapor concentration at each of the measured locations to be expressed as the sum of the contributions from the entire surface. A least squares solution is found that estimates the evapotranspiration across the field as a function of position. To evaluate the accuracy of the approach in solving for evaporation rates two artificial data sets of water vapor concentration (with varying amounts of artificial noise added) for known evaporation rates were generated and used to estimate the surface evaporation. The method is sensitive to the assumptions made near the edges of field. Results and details of the methods used to resolve edge effects will be shown.
HT2DINV: A 2D forward and inverse code for steady-state and transient hydraulic tomography problems
NASA Astrophysics Data System (ADS)
Soueid Ahmed, A.; Jardani, A.; Revil, A.; Dupont, J. P.
2015-12-01
Hydraulic tomography is a technique used to characterize the spatial heterogeneities of storativity and transmissivity fields. The responses of an aquifer to a source of hydraulic stimulations are used to recover the features of the estimated fields using inverse techniques. We developed a 2D free source Matlab package for performing hydraulic tomography analysis in steady state and transient regimes. The package uses the finite elements method to solve the ground water flow equation for simple or complex geometries accounting for the anisotropy of the material properties. The inverse problem is based on implementing the geostatistical quasi-linear approach of Kitanidis combined with the adjoint-state method to compute the required sensitivity matrices. For undetermined inverse problems, the adjoint-state method provides a faster and more accurate approach for the evaluation of sensitivity matrices compared with the finite differences method. Our methodology is organized in a way that permits the end-user to activate parallel computing in order to reduce the computational burden. Three case studies are investigated demonstrating the robustness and efficiency of our approach for inverting hydraulic parameters.
NASA Astrophysics Data System (ADS)
Liu, Shuang; Hu, Xiangyun; Xi, Yufei; Liu, Tianyou
2015-03-01
The regular grid discretization is prevalent in the inverse modeling for gravity and magnetic data. However, this subdivision strategy performs lower precision to represent the rugged observation surface. To deal with this problem, we evaluate a non-structured discretization method in which the subsurface with rolling terrain is divided into numbers of Delaunay triangular cells and each mesh has the uniform physical property distributions. The gravity and magnetic anomalies of a complex-shaped anomalous body are represented as the summaries of the single anomaly produced by each triangle field source. When inverting for the potential field data, we specify a minimization objective function composed of data constraints and then use the preconditioned conjugate gradient algorithm to iteratively solve the matrix minimization equations, where the preconditioner is determined by the distances between triangular cells and surface observers. We test our method using synthetic data; all tests return favorable results. In the case studies involving the gravity and magnetic anomalies of the Mengku and Pobei deposits in Xinjiang, northwest China, the inferred magnetite orebodies and ultrabasic rocks distributions are verified by the additional drilling and geological information. The discretization of constrained Delaunay triangulation provides an useful approach of computing and inverting the potential field data on the situations of undulate topography and complicated objects.
2D multi-parameter elastic seismic imaging by frequency-domain L1-norm full waveform inversion
NASA Astrophysics Data System (ADS)
Brossier, Romain; Operto, Stéphane; Virieux, Jean
2010-05-01
Full waveform inversion (FWI) is becoming a powerful and efficient tool to derive high-resolution quantitative models of the subsurface. In the frequency-domain, computationally efficient FWI algorithms can be designed for wide-aperture acquisition geometries by limiting inversion to few discrete frequencies. However, FWI remains an ill-posed and highly non-linear data-fitting procedure that is sensitive to noise, inaccuracies of the starting model and definition of multiparameter classes. The footprint of the noise in seismic imaging is conventionally mitigated by stacking highly redundant multifold data. However, when the data redundancy is decimated in the framework of efficient frequency-domain FWI, it is essential to assess the sensitivity of the inversion to noise. The impact of the noise in FWI, when applied to decimated data sets, has been marginally illustrated in the past and least-squares minimisation has remained the most popular approach. We investigate in this study the sensitivity of frequency-domain elastic FWI to noise for realistic onshore and offshore synthetic data sets contaminated by ambient random white noise. Four minimisation functionals are assessed in the framework of frequency domain FWI of decimated data: the classical least-square norm (L2), the least-absolute-values norm (L1), and some combinations of both (the Huber and the so-called Hybrid criteria). These functionals are implemented in a massively-parallel, 2D elastic frequency-domain FWI algorithm. A two-level hierarchical algorithm is implemented to mitigate the non-linearity of the inversion in complex environments. The first outer level consists of successive inversions of frequency groups of increasing high-frequency content. This level defines a multi-scale approach while preserving some data redundancy by means of simultaneous inversion of multiple frequencies. The second inner level used complex-valued frequencies for data preconditioning. This preconditioning controls the
An inversion method of 2D NMR relaxation spectra in low fields based on LSQR and L-curve
NASA Astrophysics Data System (ADS)
Su, Guanqun; Zhou, Xiaolong; Wang, Lijia; Wang, Yuanjun; Nie, Shengdong
2016-04-01
The low-field nuclear magnetic resonance (NMR) inversion method based on traditional least-squares QR decomposition (LSQR) always produces some oscillating spectra. Moreover, the solution obtained by traditional LSQR algorithm often cannot reflect the true distribution of all the components. Hence, a good solution requires some manual intervention, for especially low signal-to-noise ratio (SNR) data. An approach based on the LSQR algorithm and L-curve is presented to solve this problem. The L-curve method is applied to obtain an improved initial optimal solution by balancing the residual and the complexity of the solutions instead of manually adjusting the smoothing parameters. First, the traditional LSQR algorithm is used on 2D NMR T1-T2 data to obtain its resultant spectra and corresponding residuals, whose norms are utilized to plot the L-curve. Second, the corner of the L-curve as the initial optimal solution for the non-negative constraint is located. Finally, a 2D map is corrected and calculated iteratively based on the initial optimal solution. The proposed approach is tested on both simulated and measured data. The results show that this algorithm is robust, accurate and promising for the NMR analysis.
High-resistance liquid-crystal lens array for rotatable 2D/3D autostereoscopic display.
Chang, Yu-Cheng; Jen, Tai-Hsiang; Ting, Chih-Hung; Huang, Yi-Pai
2014-02-10
A 2D/3D switchable and rotatable autostereoscopic display using a high-resistance liquid-crystal (Hi-R LC) lens array is investigated in this paper. Using high-resistance layers in an LC cell, a gradient electric-field distribution can be formed, which can provide a better lens-like shape of the refractive-index distribution. The advantages of the Hi-R LC lens array are its 2D/3D switchability, rotatability (in the horizontal and vertical directions), low driving voltage (~2 volts) and fast response (~0.6 second). In addition, the Hi-R LC lens array requires only a very simple fabrication process. PMID:24663563
NASA Astrophysics Data System (ADS)
Bignardi, S.; Mantovani, A.; Abu Zeid, N.
2016-08-01
OpenHVSR is a computer program developed in the Matlab environment, designed for the simultaneous modeling and inversion of large Horizontal-to-Vertical Spectral Ratio (HVSR or H/V) datasets in order to construct 2D/3D subsurface models (topography included). The program is designed to provide a high level of interactive experience to the user and still to be of intuitive use. It implements several effective and established tools already present in the code ModelHVSR by Herak (2008), and many novel features such as: -confidence evaluation on lateral heterogeneity -evaluation of frequency dependent single parameter impact on the misfit function -relaxation of Vp/Vs bounds to allow for water table inclusion -a new cost function formulation which include a slope dependent term for fast matching of peaks, which greatly enhances convergence in case of low quality HVSR curves inversion -capability for the user of editing the subsurface model at any time during the inversion and capability to test the changes before acceptance. In what follows, we shall present many features of the program and we shall show its capabilities on both simulated and real data. We aim to supply a powerful tool to the scientific and professional community capable of handling large sets of HSVR curves, to retrieve the most from their microtremor data within a reduced amount of time and allowing the experienced scientist the necessary flexibility to integrate into the model their own geological knowledge of the sites under investigation. This is especially desirable now that microtremor testing has become routinely used. After testing the code over different datasets, both simulated and real, we finally decided to make it available in an open source format. The program is available by contacting the authors.
3D Inversion of complex resistivity data: Case study on Mineral Exploration Site.
NASA Astrophysics Data System (ADS)
Son, Jeong-Sul; Kim, Jung-ho; Park, Sam-gyu; Park, My-Kyung
2016-04-01
Complex resistivity (CR) method is a frequency domain induced polarization (IP) method. It is also known as Spectral IP (SIP) method, if wider frequencies are used in data acquisition and interpretation. Although it takes more times than conventional time domain IP method, its data quality is more stable because its data acquisition which measures amplitude and phase is done when the source current is being injected. Our research group has been studying the modeling and inversion algorithms of complex resistivity (CR) method since several years ago and recently applied developed algorithms to various real field application. Due to tough terrain in our country, Profile survey and 2D interpretation were generally used. But to get more precise interpretation, three dimensional modeling and inversion algorithm is required. We developed three dimensional inversion algorithm for this purpose. In the inversion, we adopt the method of adaptive lagraingian multiplier which is automatically set based on the size of error misfit and model regularization norm. It was applied on the real data acquired for mineral exploration sites. CR data was acquired with the Zeta system, manufactured by Zonge Co. In the inversion, only the lower frequency data is used considering its quality and developed 3D inversion algorithm was applied to the acquired data set. Its results were compared to those of time domain IP data conducted at the same site. Resistivity image sections of CR and conventional resistivity method were almost identical. Phase anomalies were well matched with chargeability anomalies and the mining history of the test site. Each anomalies were well discriminated in 3D interpretation than those of 2D. From those experiments, we know that CR method was very effective for the mineral exploration.
Monitoring of injected CO2 using the seismic full waveform inversion for 2-D elastic VTI media
NASA Astrophysics Data System (ADS)
Kim, W. K.; Min, D. J.; KIM, S.; Shin, Y.; Moon, S.
2014-12-01
To monitor the injected CO2 in the subsurface, seismic monitoring techniques are extensively applied because of its high resolution. Among the seismic monitoring techniques, seismic full waveform inversion (FWI) has high applicability because it can delineate parameter changes by injected CO2. When seismic FWIs are applied, subsurface media can be generally assumed to be isotropic. However, most subsurface media are not isotropic, and shale is a representative anisotropic medium, particularly vertical transversely isotropic (VTI) medium, which is often encountered as a barrier to injected CO2. Thus, anisotropic properties of subsurface media are important for monitoring of injected CO2. For these issues, we need to consider anisotropy of subsurface media when seismic FWIs are applied as a monitoring tool for CO2 sequestration. In this study, we performed seismic FWI for 2-D elastic VTI media to investigate the effects of anisotropic properties in CO2 monitoring. For this numerical test, we assumed a geological model, which copies after one of CO2 storage prospects in Korea. We also applied seismic FWI algorithm for 2-D elastic isotropic media for comparison. From this comparison, we noticed that we can obtain more reliable results when we apply the anisotropic FWI algorithm. Numerical examples indicate that we should apply the anisotropic FWI algorithm rather than the isotropic FWI algorithm when we interpret seismic monitoring data acquired in anisotropic media to increase the success of monitoring for injected CO2. Our numerical results can also be used as references for real seismic monitoring of the Korea CO2 sequestration projects in the near future. Acknowledgements This work was supported by the Human Resources Development program (No. 20134010200510) of the Korea Institute of Energy Technology Evaluation and Planning (KETEP) grant funded by the Korean government Ministry of Trade, Industry, and Energy and by the "Development of Technology for CO2 Marine
NASA Astrophysics Data System (ADS)
Ali, Nisa'; Saad, Rosli; Muztaza, Nordiana M.; Ismail, Noer E. H.
2013-05-01
The purpose of this study was to locate the geological contact using 2D resistivity method for Light Rail Transit (LRT) track alignment. The resistivity method was conducted on eight survey lines with the length of line 1 was 600m. The length of line 2, 3, 4, 5, 6, and 7 were 200m each while line 8 is 115m. All the survey used minimum electrode spacing of 5m and using Pole-dipole array with minimum current is 2mA and maximum was 20mA. The result obtained from the pseudosection showed that the area generally divided into three main zones, fill materials/residual soil with a resistivity value of <500 Ωm, saturated zone with a resistivity value of 30-100 Ωm and bedrock with a resistivity value of >2000 Ωm. Three fractured zones were detected along line L1 and a lot of boulders were detected at L1, L3, L4, L5 and L6. The geological contact was between the residual soil and granite bedrock.
A Novel Crosstalk Suppression Method of the 2-D Networked Resistive Sensor Array
Wu, Jianfeng; Wang, Lei; Li, Jianqing; Song, Aiguo
2014-01-01
The 2-D resistive sensor array in the row–column fashion suffered from the crosstalk problem for parasitic parallel paths. Firstly, we proposed an Improved Isolated Drive Feedback Circuit with Compensation (IIDFCC) based on the voltage feedback method to suppress the crosstalk. In this method, a compensated resistor was specially used to reduce the crosstalk caused by the column multiplexer resistors and the adjacent row elements. Then, a mathematical equivalent resistance expression of the element being tested (EBT) of this circuit was analytically derived and verified by the circuit simulations. The simulation results show that the measurement method can greatly reduce the influence on the EBT caused by parasitic parallel paths for the multiplexers' channel resistor and the adjacent elements. PMID:25046011
Two-dimensional inversion of resistivity monitoring data from the Cerro Prieto geothermal field
Goldstein, N.E.; Sasaki, Y.; Wilt, M.J.
1985-03-01
Two-dimensional iterative, least-squares inversions were performed on dc resistivity data obtained over the Cerro Prieto geothermal field at five successive times during the 1979-1983 period. The data were taken on a 20-km-long control line centered over the production region. Inversions were performed on the apparent resistivities after they were converted to percent changes in apparent resistivity relative to the base year data of 1979. The resulting solutions gave the percent change in resistivity within each of 47 rectangular blocks representing the reservoir and recharge regions. These changes are compared to and found consistent with hydrogeologic and recharge models proposed by other workers on the basis of geophysical well logs, well cuttings, well production, geochemical and reservoir engineering data. The solutions support the model of a reservoir that is being recharged mainly by cooler, less saline water, causing changes in both pore fluid resistivity and the extent of boiling near the wells. There may be a component of high-temperature recharge from below and to the east, but flow may be impeded by a two-phase zone. Notwithstanding the various sources of error and uncertainty in the data acquisition and 2-D inversions, repetitive, high precision dc resistivity monitoring seems to be a useful method for assessing reservoir conditions when used in conjunction with production and reservoir engineering data and analyses. 17 refs., 6 figs.
NASA Astrophysics Data System (ADS)
Yan, Bo; Li, Yuguo; Liu, Ying
2016-07-01
In this paper, we present an adaptive finite element (FE) algorithm for direct current (DC) resistivity modeling in 2-D generally anisotropic conductivity structures. Our algorithm is implemented on an unstructured triangular mesh that readily accommodates complex structures such as topography and dipping layers and so on. We implement a self-adaptive, goal-oriented grid refinement algorithm in which the finite element analysis is performed on a sequence of refined grids. The grid refinement process is guided by an a posteriori error estimator. The problem is formulated in terms of total potentials where mixed boundary conditions are incorporated. This type of boundary condition is superior to the Dirichlet type of conditions and improves numerical accuracy considerably according to model calculations. We have verified the adaptive finite element algorithm using a two-layered earth with azimuthal anisotropy. The FE algorithm with incorporation of mixed boundary conditions achieves high accuracy. The relative error between the numerical and analytical solutions is less than 1% except in the vicinity of the current source location, where the relative error is up to 2.4%. A 2-D anisotropic model is used to demonstrate the effects of anisotropy upon the apparent resistivity in DC soundings.
Image-guided inversion of electrical resistivity data
NASA Astrophysics Data System (ADS)
Zhou, J.; Revil, A.; Karaoulis, M.; Hale, D.; Doetsch, J.; Cuttler, S.
2014-04-01
Electrical resistivity tomography (ERT) is based on solving a Poisson equation for the electrical potential and is characterized by a good sensitivity only in the vicinity of the electrodes used to gather the data. To provide more information to ERT, we propose an image-guided or structure-constrained inversion of the apparent resistivity data. This approach uses structural information obtained directly from a guiding image. This guiding image can be drawn from a high resolution geophysical method based on the propagation equation (e.g. migrated seismic or ground penetrating radar images) or possibly from a geological cross-section of the subsurface based on some prior geological expertise. The locations and orientations of the structural features can be extracted by image processing methods to determine the structure tensor and the semblances of the guiding image at a set of pixel. Then, we introduce these structural constraints into the inversion of the apparent resistivity data by weighting the four-direction smoothing matrix to smooth along, but not across, structural features. This approach allows preserving both discontinuities and coherences in the inversion of the resistivity data. The image-guided inversion is also combined with an image-guided interpolation approach used to focus a smooth resistivity image. This yields structurally-appealing resistivity tomograms, while the whole process remains computationally efficient. Such a procedure generates a more realistic resistivity distribution (closer to the true ones), which can be, in turn, used quantitatively using appropriate petrophysical transforms, to obtain parameters of interest such as porosity and saturation. We check the validity of this approach using two synthetic case studies as well as two real datasets. For the field data, the image used to guide the inversion of the electrical resistivity data is a GPR section in the first case and a combination of seismic and structural information in the
NASA Astrophysics Data System (ADS)
Brasse, Heinrich; Schäfer, Anja; Díaz, Daniel; Alvarado, Guillermo E.; Muñoz, Angélica; Mütschard, Lutz
2015-11-01
A long-period magnetotelluric (MT) experiment was conducted in early 2009 in western Nicaragua to study the electrical resistivity and thus fluid/melt distribution at the Central American continental margin where the Cocos plate subducts beneath the Caribbean plate. Strike analysis yields a preference direction perpendicular to the profile, with moderate deviation from two-dimensionality, however. Two-dimensional modeling maps the sediments of the Nicaraguan Depression and a high-conductivity zone in the mid-crust, slightly offset from the arc. Further conductors are modeled in the backarc. However, these features are probably artifacts when a 2-D program is applied to data which show moderate 3-D characteristics. 3-D inversion clarifies the situation, and the major remaining conductive structure is now quasi directly beneath the volcanic chain and interpreted as a deep-seated magma deposit. Conductivity in the backarc is also relatively high and may either be caused by still existing partial melts beneath the Paleocene to Miocene volcanic arcs or by related metallic deposits in the aureoles of hydrothermal alteration.
NASA Astrophysics Data System (ADS)
Gao, J.; Zhang, H.
2015-12-01
Near surface geophysical exploration for the purpose of engineering design or construction For this reason, geophysical imaging demands a higher resolution and a better quantitative interpretation. Seismic travel time tomography and direct current resistivity tomography are two main methods for the near surface survey. Because of the limited coverage of observation system and the complex physical relationship between physical parameters and observations, individual geophysical method suffers issues of non-uniqueness and resolution limitation to some degree. We have developed a joint inversion method to combine seismic travel time tomography and full channel resistivity tomography. For the full channel resistivity survey, it uses two electrodes for power supply and all the other electrodes for recording. Compared with the traditional resistivity method, it collects more data and has a better model converge. Our joint inversion strategy relies on the structure constraint enforced through minimizing cross gradients between seismic velocity and resistivity models (Gallardo, 2003). For resistivity tomography, sensitivity kernels are obtained through the adjoint method by solving the electrostatic field equation with the finite-difference method. For seismic travel time tomography, ray paths and travel times are calculated using the fast marching method. We have tested our joint inversion method for a 2D cross-hole problem where two small zones with high and low velocity/resistivity anomalies. Seismic/electrical sources/receivers are installed in two boreholes. For separate seismic inversion, the smearing effect is evident and two anomaly zones are distorted and misplaced. For separate electric resistivity inversion, although two anomaly zones are positioned correctly their values are not accurate. By joint inversion, two velocity anomaly zones are clearly imaged and the smearing effect is greatly reduced. In comparison, for the resistivity model, the two anomaly zones
NASA Astrophysics Data System (ADS)
Kutter, S.; Chen, X.; Weckmann, U.
2011-12-01
Magnetotelluric (MT) measurements in areas with electromagnetic (EM) noise sources such as electric fences, power and railway lines pose severe challenges to the standard processing procedures. In order to significantly improve the data quality advanced filtering and processing techniques need to be applied. The presented 5-component MT data set from two field campaigns in 2009 and 2010 in the Barberton/Badplaas area, South Africa, was acquired within the framework of the German-South African geo-scientific research initiative Inkaba yeAfrica. Approximately 200 MT sites aligned along six profiles provide a good areal coverage of the southern part of the Barberton Greenstone Belt (BGB). Since it is one of the few remaining well-preserved geological formations from the Archean, it presents an ideal area to study the tectonic evolution and the role of plate tectonics on Early Earth. Comparing the electric properties, the surrounding high and low grade metamorphic rocks are characteristically resistive whereas mineralized shear zones are possible areas of higher electrical conductivity. Mapping their depth extension is a crucial step towards understanding the formation and the evolution of the BGB. Unfortunately, in the measurement area numerous noise sources were active, producing severe spikes and steps in the EM fields. These disturbances mainly affect long periods which are needed for resolving the deepest structures. The Remote Reference technique as well as two filtering techniques are applied to improve the data in different period ranges. Adjusting their parameters for each site is necessary to obtain the best possible results. The improved data set is used for two-dimensional inversion studies for the six profiles applying the RLM2DI algorithm by Rodi and Mackie (2001, implemented in WinGlink). In the models, areas with higher conductivity can be traced beneath known faults throughout the entire array along different profiles. Resistive zones seem to correlate
NASA Astrophysics Data System (ADS)
Schaa, R.; Gross, L.; du Plessis, J.
2016-04-01
We present a general finite-element solver, escript, tailored to solve geophysical forward and inverse modeling problems in terms of partial differential equations (PDEs) with suitable boundary conditions. Escript’s abstract interface allows geoscientists to focus on solving the actual problem without being experts in numerical modeling. General-purpose finite element solvers have found wide use especially in engineering fields and find increasing application in the geophysical disciplines as these offer a single interface to tackle different geophysical problems. These solvers are useful for data interpretation and for research, but can also be a useful tool in educational settings. This paper serves as an introduction into PDE-based modeling with escript where we demonstrate in detail how escript is used to solve two different forward modeling problems from applied geophysics (3D DC resistivity and 2D magnetotellurics). Based on these two different cases, other geophysical modeling work can easily be realized. The escript package is implemented as a Python library and allows the solution of coupled, linear or non-linear, time-dependent PDEs. Parallel execution for both shared and distributed memory architectures is supported and can be used without modifications to the scripts.
Modeling and inversion Matlab algorithms for resistivity, induced polarization and seismic data
NASA Astrophysics Data System (ADS)
Karaoulis, M.; Revil, A.; Minsley, B. J.; Werkema, D. D.
2011-12-01
M. Karaoulis (1), D.D. Werkema (3), A. Revil (1,2), A., B. Minsley (4), (1) Colorado School of Mines, Dept. of Geophysics, Golden, CO, USA. (2) ISTerre, CNRS, UMR 5559, Université de Savoie, Equipe Volcan, Le Bourget du Lac, France. (3) U.S. EPA, ORD, NERL, ESD, CMB, Las Vegas, Nevada, USA . (4) USGS, Federal Center, Lakewood, 10, 80225-0046, CO. Abstract We propose 2D and 3D forward modeling and inversion package for DC resistivity, time domain induced polarization (IP), frequency-domain IP, and seismic refraction data. For the resistivity and IP case, discretization is based on rectangular cells, where each cell has as unknown resistivity in the case of DC modelling, resistivity and chargeability in the time domain IP modelling, and complex resistivity in the spectral IP modelling. The governing partial-differential equations are solved with the finite element method, which can be applied to both real and complex variables that are solved for. For the seismic case, forward modeling is based on solving the eikonal equation using a second-order fast marching method. The wavepaths are materialized by Fresnel volumes rather than by conventional rays. This approach accounts for complicated velocity models and is advantageous because it considers frequency effects on the velocity resolution. The inversion can accommodate data at a single time step, or as a time-lapse dataset if the geophysical data are gathered for monitoring purposes. The aim of time-lapse inversion is to find the change in the velocities or resistivities of each model cell as a function of time. Different time-lapse algorithms can be applied such as independent inversion, difference inversion, 4D inversion, and 4D active time constraint inversion. The forward algorithms are benchmarked against analytical solutions and inversion results are compared with existing ones. The algorithms are packaged as Matlab codes with a simple Graphical User Interface. Although the code is parallelized for multi
Arnal, B; Pinton, G; Garapon, P; Pernot, M; Fink, M; Tanter, M
2013-10-01
Shear wave imaging (SWI) maps soft tissue elasticity by measuring shear wave propagation with ultrafast ultrasound acquisitions (10 000 frames s(-1)). This spatiotemporal data can be used as an input for an inverse problem that determines a shear modulus map. Common inversion methods are local: the shear modulus at each point is calculated based on the values of its neighbour (e.g. time-of-flight, wave equation inversion). However, these approaches are sensitive to the information loss such as noise or the lack of the backscattered signal. In this paper, we evaluate the benefits of a global approach for elasticity inversion using a least-squares formulation, which is derived from full waveform inversion in geophysics known as the adjoint method. We simulate an acoustic waveform in a medium with a soft and a hard lesion. For this initial application, full elastic propagation and viscosity are ignored. We demonstrate that the reconstruction of the shear modulus map is robust with a non-uniform background or in the presence of noise with regularization. Compared to regular local inversions, the global approach leads to an increase of contrast (∼+3 dB) and a decrease of the quantification error (∼+2%). We demonstrate that the inversion is reliable in the case when there is no signal measured within the inclusions like hypoechoic lesions which could have an impact on medical diagnosis. PMID:24018867
Optimized arrays for 2-D resistivity survey lines with a large number of electrodes
NASA Astrophysics Data System (ADS)
Loke, M. H.; Wilkinson, P. B.; Chambers, J. E.; Uhlemann, S. S.; Sorensen, J. P. R.
2015-01-01
Previous studies show that optimized arrays generated using the 'Compare R' method have significantly better resolution than conventional arrays. This method determines the optimum set of arrays by selecting those that give the maximum model resolution. The number of possible arrays (the comprehensive data set) increases with the fourth power of the number of electrodes. The optimization method faces practical limitations for 2-D survey lines with more than 60 electrodes where the number of possible arrays exceeds a million. Several techniques are proposed to reduce the calculation time for such survey lines. A single-precision version of the 'Compare R' algorithm using a new ranking function reduces the calculation time by two to eight times while providing results similar to the double-precision version. Recent improvements in computer GPU technology can reduce the calculation time by about seven times. The calculation time is reduced by half by using the fact that arrays that are symmetrical about the center of the line produce identical changes in the model resolution values. It is further reduced by more than thirty times by calculating the Sherman-Morrison update for all the possible two-electrode combinations, which are then used to calculate the model resolution values for the four-electrode arrays. The calculation time is reduced by more then ten times by using a subset of the comprehensive data set consisting of only symmetrical arrays. Tests with a synthetic model and field data set show that optimized arrays derived from this subset produce inversion models with differences of less than 10% from those derived using the full comprehensive data set. The optimized data sets produced models that are more accurate than the Wenner-Schlumberger array data sets in all the tests.
NASA Astrophysics Data System (ADS)
Blacic, Tanya M.; Jun, Hyunggu; Rosado, Hayley; Shin, Changsoo
2016-02-01
In seismic oceanography, processed images highlight small temperature changes, but inversion is needed to obtain absolute temperatures. Local search-based full waveform inversion has a lower computational cost than global search but requires accurate starting models. Unfortunately, most marine seismic data have little associated hydrographic data and the band-limited nature of seismic data makes extracting the long wavelength sound speed trend directly from seismic data inherently challenging. Laplace and Laplace-Fourier domain inversion (LDI) can use rudimentary starting models without prior information about the medium. Data are transformed to the Laplace domain, and a smooth sound speed model is extracted by examining the zero and low frequency components of the damped wavefield. We applied LDI to five synthetic data sets based on oceanographic features and recovered smoothed versions of our synthetic models, showing the viability of LDI for creating starting models suitable for more detailed inversions.
Altered NK Cell Development and Enhanced NK Cell-Mediated Resistance to MCMV in NKG2D-Deficient Mice
Zafirova, Biljana; Mandarić, Sanja; Antulov, Ronald; Krmpotić, Astrid; Jonsson, Helena; Yokoyama, Wayne M.; Jonjić, Stipan; Polić, Bojan
2009-01-01
Summary NKG2D is a potent activating receptor on NK cells which acts as a molecular sensor for stressed cells expressing NKG2D ligands such as infected or tumor transformed cells. Although NKG2D is expressed on NK cell precursors, its role in NK cell development is still not known. We have generated NKG2D-deficient mice by targeting the Klrk1 locus. Here we provide evidence for an important regulatory role of NKG2D in the development of NK cells. The absence of NKG2D causes faster division of NK cells, perturbation in size of some NK cell subpopulations and their augmented sensitivity to apoptosis. As expected, NKG2D−/− NK cells are less responsive to tumor targets expressing NKG2D ligands. NKG2D−/− mice, however, show an enhanced NK cell-mediated resistance to MCMV infection as a consequence of NK cell dysregulation. Altogether, these findings provide evidence for yet unknown regulatory function of NKG2D in NK cell physiology. PMID:19631564
Residual resistance of 2D and 3D structures and Joule heat release.
Gurevich, V L; Kozub, V I
2011-06-22
We consider a residual resistance and Joule heat release in 2D nanostructures as well as in ordinary 3D conductors. We assume that elastic scattering of conduction electrons by lattice defects is predominant. Within a rather intricate situation in such systems we discuss in detail two cases. (1) The elastic scattering alone (i.e. without regard of inelastic mechanisms of scattering) leads to a transition of the mechanical energy (stored by the electrons under the action of an electric field) into heat in a traditional way. This process can be described by the Boltzmann equation where it is possible to do the configuration averaging over defect positions in the electron-impurity collision term. The corresponding conditions are usually met in metals. (2) The elastic scattering can be considered with the help of the standard electron-impurity collision integral only in combination with some additional averaging procedure (possibly including inelastic scattering or some mechanisms of electron wavefunction phase destruction). This situation is typical for degenerate semiconductors with a high concentration of dopants and conduction electrons. Quite often, heat release can be observed via transfer of heat to the lattice, i.e. via inelastic processes of electron-phonon collisions and can take place at distances much larger than the size of the device. However, a direct heating of the electron system can be registered too by, for instance, local measurements of the current noise or direct measurement of an electron distribution function. PMID:21628783
2D-resistivity surveys of deteriorating historic stonework in Oxford, UK
NASA Astrophysics Data System (ADS)
Sass, O.; Viles, H. A.
2009-04-01
Historic stonework deteriorates in often very complex ways and despite many years of research on the topic, we are still a long way from being able to predict its occurrence and severity. As most deterioration processes involve water, techniques which provide a better picture of the moisture contents and fluctuations within stonework are very valuable in attempts to improve understanding. 2D resistivity methods can provide useful information about moisture distributions within porous historic stonework. We report on a series of experiments on historic walls within the centre of Oxford, UK, which illustrate varying degrees of deterioration including catastrophic decay. Using medical electrodes we have been able to carry out non-invasive and non-destructive 2D resistivity surveys to study the distribution and amount of water stored in deteriorating limestone walls. Fifteen vertical profiles, each 2-2.5 m in length, have been monitored at five sites. Furthermore, simulated driving rain experiments have been carried out at two sites. The data indicate the diversity and complexity of moisture distributions within these walls. Replacement stone patches show consistently higher moisture conditions than the surrounding stone. Some profiles show wetter sections towards the base of the wall, usually where a plinth is absent. Conversely, hard stone plinths obviously reduce capillary rise from ground water. However, at several sites we noticed a wetter zone immediately above the top of the plinth which often correlates with the occurrence of catastrophic decay - indicating that the plinth may encourage concentration of decay. Most profiles indicate the presence of wetter patches 5-10cm behind the wall face under blackened crusts. Such patches of heightened absolute moisture contents could play a very important role in encouraging catastrophic decay. Severely decayed sections of profiles often exhibit wetter near-surface conditions than surrounding stonework, whilst areas with
NASA Astrophysics Data System (ADS)
Yang, X.; Zhu, P.; Gu, Y.; Xu, Z.
2015-12-01
Small scale heterogeneities of subsurface medium can be characterized conveniently and effectively using a few simple random medium parameters (RMP), such as autocorrelation length, angle and roughness factor, etc. The estimation of these parameters is significant in both oil reservoir prediction and metallic mine exploration. Poor accuracy and low stability existed in current estimation approaches limit the application of random medium theory in seismic exploration. This study focuses on improving the accuracy and stability of RMP estimation from post-stacked seismic data and its application in the seismic inversion. Experiment and theory analysis indicate that, although the autocorrelation of random medium is related to those of corresponding post-stacked seismic data, the relationship is obviously affected by the seismic dominant frequency, the autocorrelation length, roughness factor and so on. Also the error of calculation of autocorrelation in the case of finite and discrete model decreases the accuracy. In order to improve the precision of estimation of RMP, we design two improved approaches. Firstly, we apply region growing algorithm, which often used in image processing, to reduce the influence of noise in the autocorrelation calculated by the power spectrum method. Secondly, the orientation of autocorrelation is used as a new constraint in the estimation algorithm. The numerical experiments proved that it is feasible. In addition, in post-stack seismic inversion of random medium, the estimated RMP may be used to constrain inverse procedure and to construct the initial model. The experiment results indicate that taking inversed model as random medium and using relatively accurate estimated RMP to construct initial model can get better inversion result, which contained more details conformed to the actual underground medium.
Lin, Lin; Yang, Chao; Lu, Jiangfeng; Ying, Lexing; E, Weinan
2009-09-25
We present an efficient parallel algorithm and its implementation for computing the diagonal of $H^-1$ where $H$ is a 2D Kohn-Sham Hamiltonian discretized on a rectangular domain using a standard second order finite difference scheme. This type of calculation can be used to obtain an accurate approximation to the diagonal of a Fermi-Dirac function of $H$ through a recently developed pole-expansion technique \\cite{LinLuYingE2009}. The diagonal elements are needed in electronic structure calculations for quantum mechanical systems \\citeHohenbergKohn1964, KohnSham 1965,DreizlerGross1990. We show how elimination tree is used to organize the parallel computation and how synchronization overhead is reduced by passing data level by level along this tree using the technique of local buffers and relative indices. We analyze the performance of our implementation by examining its load balance and communication overhead. We show that our implementation exhibits an excellent weak scaling on a large-scale high performance distributed parallel machine. When compared with standard approach for evaluating the diagonal a Fermi-Dirac function of a Kohn-Sham Hamiltonian associated a 2D electron quantum dot, the new pole-expansion technique that uses our algorithm to compute the diagonal of $(H-z_i I)^-1$ for a small number of poles $z_i$ is much faster, especially when the quantum dot contains many electrons.
NASA Astrophysics Data System (ADS)
Key, K.
2013-12-01
This work announces the public release of an open-source inversion code named MARE2DEM (Modeling with Adaptively Refined Elements for 2D Electromagnetics). Although initially designed for the rapid inversion of marine electromagnetic data, MARE2DEM now supports a wide variety of acquisition configurations for both offshore and onshore surveys that utilize electric and magnetic dipole transmitters or magnetotelluric plane waves. The model domain is flexibly parameterized using a grid of arbitrarily shaped polygonal regions, allowing for complicated structures such as topography or seismically imaged horizons to be easily assimilated. MARE2DEM efficiently solves the forward problem in parallel by dividing the input data parameters into smaller subsets using a parallel data decomposition algorithm. The data subsets are then solved in parallel using an automatic adaptive finite element method that iterative solves the forward problem on successively refined finite element meshes until a specified accuracy tolerance is met, thus freeing the end user from the burden of designing an accurate numerical modeling grid. Regularized non-linear inversion for isotropic or anisotropic conductivity is accomplished with a new implementation of Occam's method referred to as fast-Occam, which is able to minimize the objective function in much fewer forward evaluations than the required by the original method. This presentation will review the theoretical considerations behind MARE2DEM and use a few recent offshore EM data sets to demonstrate its capabilities and to showcase the software interface tools that streamline model building and data inversion.
NASA Astrophysics Data System (ADS)
Pérez-Corona, M.; García, J. A.; Taller, G.; Polgár, D.; Bustos, E.; Plank, Z.
2016-02-01
The purpose of geophysical electrical surveys is to determine the subsurface resistivity distribution by making measurements on the ground surface. From these measurements, the true resistivity of the subsurface can be estimated. The ground resistivity is related to various geological parameters, such as the mineral and fluid content, porosity and degree of water saturation in the rock. Electrical resistivity surveys have been used for many decades in hydrogeological, mining and geotechnical investigations. More recently, they have been used for environmental surveys. To obtain a more accurate subsurface model than is possible with a simple 1-D model, a more complex model must be used. In a 2-D model, the resistivity values are allowed to vary in one horizontal direction (usually referred to as the x direction) but are assumed to be constant in the other horizontal (the y) direction. A more realistic model would be a fully 3-D model where the resistivity values are allowed to change in all three directions. In this research, a simulation of the cone penetration test and 2D imaging resistivity are used as tools to simulate the distribution of hydrocarbons in soil.
NASA Astrophysics Data System (ADS)
Kumar, Dewashish; Rao, V. Ananda; Sarma, V. S.
2014-04-01
Electrical resistivity method is a versatile and economical technique for groundwater prospecting in different geological settings due to wide spectrum of resistivity compared to other geophysical parameters. Exploration and exploitation of groundwater, a vital and precious resource, is a challenging task in hard rock, which exhibits inherent heterogeneity. In the present study, two-dimensional Electrical Resistivity Tomography (2D-ERT) technique using two different arrays, viz., pole-dipole and pole-pole, were deployed to look into high signal strength data in a tectonically disturbed hard rock ridge region for groundwater. Four selected sites were investigated. 2D subsurface resistivity tomography data were collected using Syscal Pro Switch-10 channel system and covered a 2 km long profile in a tough terrain. The hydrogeological interpretation based on resistivity models reveal the water horizons trap within the clayey sand and weathered/fractured quartzite formations. Aquifer resistivity lies between ˜3-35 and 100-200 Ωm. The results of the resistivity models decipher potential aquifer lying between 40 and 88 m depth, nevertheless, it corroborates with the static water level measurements in the area of study. The advantage of using pole-pole in conjunction with the pole-dipole array is well appreciated and proved worth which gives clear insight of the aquifer extent, variability and their dimension from shallow to deeper strata from the hydrogeological perspective in the present geological context.
NASA Astrophysics Data System (ADS)
Schmoldt, J.; Jones, A. G.; Muller, M. R.; Kiyan, D.; Hogg, C.; Rosell, O.
2010-12-01
Two-dimensional (2D) inversions of magnetotelluric (MT) data are presently far more commonly used than three-dimensional (3D) inversions as they still significantly outperform 3D inversions in terms of speed, thus allowing for much better resolution of the subsurface through a larger feasible number of grid cells. The suitability of 2D inversion needs thereby to be tested for cases where the electric resistivity structure of the subsurface is potentially 3D to some extent. One particular case of a 3D subsurface structure consists of lateral interfaces with varying orientations at crustal and mantle depths. Such a case might emerge, for instance, where crustal faulting, originating from present day tectonics, is situated above a mantle where structures are dominated by earlier or current plate tectonic processes. Those plate tectonic processes could comprise continental collision from an oblique direction, or lattice preferred orientation in the lithosphere-asthenosphere transition zone due to an oblique relative motion between lithosphere and asthenosphere. Whereas recovery of crustal structures can usually be achieved in a straightforward manner by confining the modelled frequency range to the crustal depths, deriving mantle structures is more challenging. Different approaches for this case have been investigated here using synthetic model studies as well as inversion of a real MT dataset collected in southern Spain as part of the PICASSO fieldwork campaign. The PICASSO project intends to enhance knowledge about the geological setting of the Alboran Domain beneath the western Mediterranean Sea and its surrounding regions, and through that knowledge to understand processes related to continent-continent collision. The Iberian Peninsula is the focus of the first phase of DIAS’s PICASSO efforts, and comprised a magnetotelluric profile crossing the Tajo Basin and Betics Cordillera. Analyses of MT responses and seismic tomography data indicate varying geologic strike
NASA Astrophysics Data System (ADS)
Azwin, I. N.; Saad, Rosli; Saidin, Mokhtar; Nordiana, M. M.; Anderson Bery, Andy; Hidayah, I. N. E.
2015-01-01
Interest in studying impact crater on earth has increased tremendously due to its importance in geologic events, earth inhabitant history as well as economic value. The existences of few shock metamorphism and crater morphology evidences are discovered in Bukit Bunuh, Malaysia thus detailed studies are performed using geophysical and geotechnical methods to verify the type of the crater and characteristics accordingly. This paper presents the combined analysis of 2-D electrical resistivity, seismic refraction, geotechnical SPT N value, moisture content and RQD within the study area. Three stages of data acquisition are made starting with regional study followed by detailed study on West side and East side. Bulk resistivity and p-wave seismic velocity were digitized from 2-D resistivity and seismic sections at specific distance and depth for corresponding boreholes and samples taken. Generally, Bukit Bunuh shows the complex crater characteristics. Standard table of bulk resistivity and p-wave seismic velocity against SPT N value, moisture content and RQD are produce according to geological classifications of impact crater; inside crater, rim/slumped terrace and outside crater.
NASA Astrophysics Data System (ADS)
Mashayekhi, Parisa; Ghorbani-Dashtaki, Shoja; Mosaddeghi, Mohammad Reza; Shirani, Hossein; Nodoushan, Ali Reza Mohammadi
2016-04-01
In this study, HYDRUS-2D/3D was used to simulate ponded infiltration through double-ring infiltrometers into a hypothetical loamy soil profile. Twelve scenarios of inverse modelling (divided into three groups) were considered for estimation of Mualem-van Genuchten hydraulic parameters. In the first group, simulation was carried out solely using cumulative infiltration data. In the second group, cumulative infiltration data plus water content at h = -330 cm (field capacity) were used as inputs. In the third group, cumulative infiltration data plus water contents at h = -330 cm (field capacity) and h = -15 000 cm (permanent wilting point) were used simultaneously as predictors. The results showed that numerical inverse modelling of the double-ring infiltrometer data provided a reliable alternative method for determining soil hydraulic parameters. The results also indicated that by reducing the number of hydraulic parameters involved in the optimization process, the simulation error is reduced. The best one in infiltration simulation which parameters α, n, and Ks were optimized using the infiltration data and field capacity as inputs. Including field capacity as additional data was important for better optimization/definition of soil hydraulic functions, but using field capacity and permanent wilting point simultaneously as additional data increased the simulation error.
Baker, Theresa; Nerle, Sujata; Pritchard, Justin; Zhao, Boyang; Rivera, Victor M.
2015-01-01
Although targeted therapies have revolutionized cancer treatment, overcoming acquired resistance remains a major clinical challenge. EZH2 inhibitors (EZH2i), EPZ-6438 and GSK126, are currently in the early stages of clinical evaluation and the first encouraging signs of efficacy have recently emerged in the clinic. To anticipate mechanisms of resistance to EZH2i, we used a forward genetic platform combining a mutagenesis screen with next generation sequencing technology and identified a hotspot of secondary mutations in the EZH2 D1 domain (Y111 and I109). Y111D mutation within the WT or A677G EZH2 allele conferred robust resistance to both EPZ-6438 and GSK126, but it only drove a partial resistance within the Y641F allele. EZH2 mutants required histone methyltransferase (HMT) catalytic activity and the polycomb repressive complex 2 (PRC2) components, SUZ12 and EED, to drive drug resistance. Furthermore, D1 domain mutations not only blocked the ability of EZH2i to bind to WT and A677G mutant, but also abrogated drug binding to the Y641F mutant. These data provide the first cellular validation of the mechanistic model underpinning the oncogenic function of WT and mutant EZH2. Importantly, our findings suggest that acquired-resistance to EZH2i may arise in WT and mutant EZH2 patients through a single mutation that remains targetable by second generation EZH2i. PMID:26360609
Hangel, Gilbert; Strasser, Bernhard; Považan, Michal; Gruber, Stephan; Chmelík, Marek; Gajdošík, Martin; Trattnig, Siegfried; Bogner, Wolfgang
2015-11-01
This work presents a new approach for high-resolution MRSI of the brain at 7 T in clinically feasible measurement times. Two major problems of MRSI are the long scan times for large matrix sizes and the possible spectral contamination by the transcranial lipid signal. We propose a combination of free induction decay (FID)-MRSI with a short acquisition delay and acceleration via in-plane two-dimensional generalised autocalibrating partially parallel acquisition (2D-GRAPPA) with adiabatic double inversion recovery (IR)-based lipid suppression to allow robust high-resolution MRSI. We performed Bloch simulations to evaluate the magnetisation pathways of lipids and metabolites, and compared the results with phantom measurements. Acceleration factors in the range 2-25 were tested in a phantom. Five volunteers were scanned to verify the value of our MRSI method in vivo. GRAPPA artefacts that cause fold-in of transcranial lipids were suppressed via double IR, with a non-selective symmetric frequency sweep. The use of long, low-power inversion pulses (100 ms) reduced specific absorption rate requirements. The symmetric frequency sweep over both pulses provided good lipid suppression (>90%), in addition to a reduced loss in metabolite signal-to-noise ratio (SNR), compared with conventional IR suppression (52-70%). The metabolic mapping over the whole brain slice was not limited to a rectangular region of interest. 2D-GRAPPA provided acceleration up to a factor of nine for in vivo FID-MRSI without a substantial increase in g-factors (<1.1). A 64 × 64 matrix can be acquired with a common repetition time of ~1.3 s in only 8 min without lipid artefacts caused by acceleration. Overall, we present a fast and robust MRSI method, using combined double IR fat suppression and 2D-GRAPPA acceleration, which may be used in (pre)clinical studies of the brain at 7 T. PMID:26370781
Body edge delineation in 2D DC resistivity imaging using differential method
NASA Astrophysics Data System (ADS)
Susanto, Kusnahadi; Fitrah Bahari, Mohammad
2016-01-01
DC resistivity is widely used to identify the kind of rock and the lithology contact. However, the image resulting from resistivity processing is shown in a contour image. There is be a problem to interpret where the edge of body location is. This study uses differential method to delineate the edge of body in DC resistivity contour. This method was applied to the boundary between gravel and underlying clay layer. The first and the second order differential method is applied to the delineation of lithology contact. The profiling curve has to be sliced and extracted from the resistivity contour before the differential method can be used. The spectral analysis shows the frequency and wavenumber of the profiling curve used to make gridding. The slicing process was conducted horizontally and vertically in order to get the mesh size which will be used in the differential method. The second order differential, the Laplace operator, is able to show the edge of body more clearly than the first order differential and shows the contact between gravel and clay.
Drug Abuse Resistance Education Act of 1990. House of Representatives, 101st Congress, 2d Session.
ERIC Educational Resources Information Center
Congress of the U.S., Washington, DC. House Committee on Education and Labor.
This document is a congressional report on the Drug Abuse Resistance Education Act of 1990, House of Representatives (H.R.) 5064 which requires the Secretary of Education to reserve $15 million from the appropriations made for the Drug Free Schools and Communities Act for grants to consortia of local educational agencies and other agencies to…
Resistively detected high-order magnetoplasmons in a high-quality 2D electron gas
NASA Astrophysics Data System (ADS)
Zudov, M. A.; Shi, Q.; Pfeiffer, L. N.; West, K. W.; Watson, J. D.; Manfra, M. J.
We report on high-order magnetoplasmon resonances detected in photoresistance in high-mobility GaAs quantum wells. These resonances manifest themselves as a series of resistance extrema in the regime of Shubnikov-de Haas oscillations. Extending to orders above 20, the extrema exhibit alternating strength, being less (more) pronounced at even (odd) order magnetoplasmon modes. The lower magnetoplasmon modes reveal the importance of retardation effects.
Tuning and simulating a 193-nm resist for 2D applications
NASA Astrophysics Data System (ADS)
Howard, William B.; Wiaux, Vincent; Ercken, Monique; Bui, Bang; Byers, Jeff D.; Pochkowski, Mike
2002-07-01
For some applications, the usefulness of lithography simulation results depends strongly on the matching between experimental conditions and the simulation input parameters. If this matching is optimized and other sources of error are minimized, then the lithography model can be used to explain printed wafer experimental results. Further, simulation can be useful in predicting the results or in choosing the correct set of experiments. In this paper, PROLITH and ProDATA AutoTune were used to systematically vary simulation input parameters to match measured results on printed wafers used in a 193 nm process. The validity of the simulation parameters was then checked using 3D simulation compared to 2D top-down SEM images. The quality of matching was evaluated using the 1D metrics of average gate CD and Line End Shortening (LES). To ensure the most accurate simulation, a new approach was taken to create a compound mask from GDSII contextual information surrounding an accurate SEM image of the reticle region of interest. Corrections were made to account for all metrology offsets.
NASA Astrophysics Data System (ADS)
Kamaruddin, Nur Aminuda; Saad, Rosli; Nordiana, M. M.; Azwin, I. N.
2015-04-01
The Great Sumatra Fault system was split into two sub-parallel lines or segments at the Northern Sumatra. This event is one of the impacts of powerful earthquakes that hit Sumatra Island especially one that occurred in 2004. These two sub-parallel segments known as Aceh and Seulimeum fault. The study is focused on the Seulimeum fault and two geophysical methods chosen aimed to compare and verified the result obtained respectively. 2-D resistivity method is a common geophysical method used in determination of near surface structures such as faults, cavities, voids and sinkholes. Meanwhile, the magnetic method often chosen to delineate subsurface structures, determine depth of magnetic source bodies and possibly sediment thickness. Three survey lines of resistivity method and randomly magnetic stations were carried out covering Krueng district. The resistivity data processed using Res2Dinv and result presented using Surfer software. The fault identified by the contrast of low and high resistivity value. Meanwhile, the magnetic data were presented in magnetic residual contour map and the extended fault system is suspected represent by the contrast value of the magnetic anomalies. Within suspected fault zone, the results of resistivity are tally with magnetic result.
A 2D vector map watermarking algorithm resistant to simplication attack
NASA Astrophysics Data System (ADS)
Wang, Chuanjian; Liang, Bin; Zhao, Qingzhan; Qiu, Zuqi; Peng, Yuwei; Yu, Liang
2009-12-01
Vector maps are valuable asset of data producers. How to protect copyright of vector maps effectively using digital watermarking is a hot research issue. In this paper, we propose a new robust and blind watermarking algorithm resilient to simplification attack. We proof that spatial topological relation between map objects bears an important property of approximate simplification invariance. We choose spatial topological relations as watermark feature domain and embed watermarks by slightly modifying spatial topological relation between map objects. Experiment shows that our algorithm has good performance to resist simplification attack and tradeoff of the robustness and data fidelity is acquired.
2D Resistive Magnetohydrodynamics Calculations with an Arbitrary Lagrange Eulerian Code
NASA Astrophysics Data System (ADS)
Rousculp, C. L.; Gianakon, T. A.; Lipnikov, K. N.; Nelson, E. M.
2015-11-01
Single fluid resistive MHD is useful for modeling Z-pinch configurations in cylindrical geometry. One such example is thin walled liners for shock physics or HEDP experiments driven by capacitor banks such as the LANL's PHELIX or Sandia-Z. MHD is also useful for modeling high-explosive-driven flux compression generators (FCGs) and their high-current switches. The resistive MHD in our arbitrary Lagrange Eulerian (ALE) code operates in one and two dimensions in both Cartesian and cylindrical geometry. It is implemented as a time-step split operator, which consists of, ideal MHD connected to the explicit hydro momentum and energy equations and a second order mimetic discretization solver for implicit solution of the magnetic diffusion equation. In a staggered grid scheme, a single-component of cell-centered magnetic flux is conserved in the Lagrangian frame exactly, while magnetic forces are accumulated at the nodes. Total energy is conserved to round off. Total flux is conserved under the ALE relaxation and remap. The diffusion solver consistently computes Ohmic heating. Both Neumann and Dirichlet boundary conditions are available with coupling to external circuit models. Example calculations will be shown.
NASA Astrophysics Data System (ADS)
Sourbier, F.; Operto, S.; Virieux, J.
2006-12-01
We present a distributed-memory parallel algorithm for 2D visco-acoustic full-waveform inversion of wide-angle seismic data. Our code is written in fortran90 and use MPI for parallelism. The algorithm was applied to real wide-angle data set recorded by 100 OBSs with a 1-km spacing in the eastern-Nankai trough (Japan) to image the deep structure of the subduction zone. Full-waveform inversion is applied sequentially to discrete frequencies by proceeding from the low to the high frequencies. The inverse problem is solved with a classic gradient method. Full-waveform modeling is performed with a frequency-domain finite-difference method. In the frequency-domain, solving the wave equation requires resolution of a large unsymmetric system of linear equations. We use the massively parallel direct solver MUMPS (http://www.enseeiht.fr/irit/apo/MUMPS) for distributed-memory computer to solve this system. The MUMPS solver is based on a multifrontal method for the parallel factorization. The MUMPS algorithm is subdivided in 3 main steps: a symbolic analysis step that performs re-ordering of the matrix coefficients to minimize the fill-in of the matrix during the subsequent factorization and an estimation of the assembly tree of the matrix. Second, the factorization is performed with dynamic scheduling to accomodate numerical pivoting and provides the LU factors distributed over all the processors. Third, the resolution is performed for multiple sources. To compute the gradient of the cost function, 2 simulations per shot are required (one to compute the forward wavefield and one to back-propagate residuals). The multi-source resolutions can be performed in parallel with MUMPS. In the end, each processor stores in core a sub-domain of all the solutions. These distributed solutions can be exploited to compute in parallel the gradient of the cost function. Since the gradient of the cost function is a weighted stack of the shot and residual solutions of MUMPS, each processor
Luo, Y.; Xia, J.; Liu, J.; Xu, Y.; Liu, Q.
2008-01-01
Multichannel Analysis of Surface Waves utilizes a multichannel recording system to estimate near-surface shear (S)-wave velocities from high-frequency Rayleigh waves. A pseudo-2D S-wave velocity (vS) section is constructed by aligning 1D models at the midpoint of each receiver spread and using a spatial interpolation scheme. The horizontal resolution of the section is therefore most influenced by the receiver spread length and the source interval. The receiver spread length sets the theoretical lower limit and any vS structure with its lateral dimension smaller than this length will not be properly resolved in the final vS section. A source interval smaller than the spread length will not improve the horizontal resolution because spatial smearing has already been introduced by the receiver spread. In this paper, we first analyze the horizontal resolution of a pair of synthetic traces. Resolution analysis shows that (1) a pair of traces with a smaller receiver spacing achieves higher horizontal resolution of inverted S-wave velocities but results in a larger relative error; (2) the relative error of the phase velocity at a high frequency is smaller than at a low frequency; and (3) a relative error of the inverted S-wave velocity is affected by the signal-to-noise ratio of data. These results provide us with a guideline to balance the trade-off between receiver spacing (horizontal resolution) and accuracy of the inverted S-wave velocity. We then present a scheme to generate a pseudo-2D S-wave velocity section with high horizontal resolution using multichannel records by inverting high-frequency surface-wave dispersion curves calculated through cross-correlation combined with a phase-shift scanning method. This method chooses only a pair of consecutive traces within a shot gather to calculate a dispersion curve. We finally invert surface-wave dispersion curves of synthetic and real-world data. Inversion results of both synthetic and real-world data demonstrate that
2D Time-lapse Resistivity Monitoring of an Organic Produced Gas Plume in a Landfill using ERT.
NASA Astrophysics Data System (ADS)
Amaral, N. D.; Mendonça, C. A.; Doherty, R.
2014-12-01
This project has the objective to study a landfill located on the margins of Tietê River, in São Paulo, Brazil, using the electroresistivity tomography method (ERT). Due to huge organic matter concentrations in the São Paulo Basin quaternary sediments, there is subsurface depth related biogas accumulation (CH4 and CO2), induced by anaerobic degradation of the organic matter. 2D resistivity sections were obtained from a test area since March 2012, a total of 7 databases, being the last one dated from October 2013. The studied line has the length of 56m, the electrode interval is of 2m. In addition, there are two boreholes along the line (one with 3 electrodes and the other one with 2) in order to improve data quality and precision. The boreholes also have a multi-level sampling system that indicates the fluid (gas or water) presence in relation to depth. With our results it was possible to map the gas plume position and its area of extension in the sections as it is a positive resistivity anomaly, with the gas level having approximately 5m depth. With the time-lapse analysis (Matlab script) between the obtained 2D resistivity sections from the site, it was possible to map how the biogas volume and position change in the landfill in relation to time. Our preliminary results show a preferential gas pathway through the subsurface studied area. A consistent relation between the gas depth and obtained microbiological data from archea and bacteria population was also observed.
Gállego, Isaac; Oncins, Gerard; Sisquella, Xavier; Fernàndez-Busquets, Xavier; Daban, Joan-Ramon
2010-01-01
In a previous study, we found that metaphase chromosomes are formed by thin plates, and here we have applied atomic force microscopy (AFM) and friction force measurements at the nanoscale (nanotribology) to analyze the properties of these planar structures in aqueous media at room temperature. Our results show that high concentrations of NaCl and EDTA and extensive digestion with protease and nuclease enzymes cause plate denaturation. Nanotribology studies show that native plates under structuring conditions (5 mM Mg2+) have a relatively high friction coefficient (μ ≈ 0.3), which is markedly reduced when high concentrations of NaCl or EDTA are added (μ ≈ 0.1). This lubricant effect can be interpreted considering the electrostatic repulsion between DNA phosphate groups and the AFM tip. Protease digestion increases the friction coefficient (μ ≈ 0.5), but the highest friction is observed when DNA is cleaved by micrococcal nuclease (μ ≈ 0.9), indicating that DNA is the main structural element of plates. Whereas nuclease-digested plates are irreversibly damaged after the friction measurement, native plates can absorb kinetic energy from the AFM tip without suffering any damage. These results suggest that plates are formed by a flexible and mechanically resistant two-dimensional network which allows the safe storage of DNA during mitosis. PMID:21156137
NASA Astrophysics Data System (ADS)
Jinmin, M.; Saad, R.; Saidin, M.; Ismail, N. A.
2015-03-01
A 2-D resistivity imaging (2-DERI) study was conducted at Bukit Bunuh, Lenggong, Perak. Archaeological Global Research Centre, Universiti Sains Malaysia shows the field evidence of shock metamorphisms (suevite breccia) and crater morphology at Bukit Bunuh. A regional 2-DERI study focusing at Bukit Bunuh to identify the features of subsurface and detail study was then executed to verify boundary of the crater with the rebound effects at Bukit Bunuh which covered approximately 132.25 km2. 2-DERI survey used resistivity equipment by ABEM SAS4000 Terrameter and ES10-64C electrode slector with pole-dipole array. The survey lines were carried out using `roll-along' technique. The data were processed and analysed using RES2DINV, Excel and Surfer software to obtain resistivity results for qualitative interpretations. Bedrock depths were digitized from section by sections obtained. 2-DERI results gives both regional and detail study shows that the study area was divided into two main zones, overburden consists of alluvium mix with boulders embedded with resistivity value of 10-800 Ωm and granitic bedrock with resistivity value of >1500 Ωm and depth 5-50 m. The low level bedrock was circulated by high level bedrock (crater rim) was formed at the same area with few spots of high level bedrock which appeared at the centre of the rim which suspected as rebound zones (R). Assimilations of 2-DERI with boreholes are successful give valid and reliable results. The results of the study indicates geophysical method are capable to retrieve evidence of meteorite impact subsurface of the studied area.
ERIC Educational Resources Information Center
Brown, Malcolm
2009-01-01
Inversions are fascinating phenomena. They are reversals of the normal or expected order. They occur across a wide variety of contexts. What do inversions have to do with learning spaces? The author suggests that they are a useful metaphor for the process that is unfolding in higher education with respect to education. On the basis of…
NASA Astrophysics Data System (ADS)
Escalona, A.; Sena, A.; Mann, P.
2003-12-01
We have reprocessed and reinterpreted more than 10,000 km of "Gulfrex" multi-channel 2D seismic reflection lines collected by Gulf Oil Corporation in 1972 along the northern margin of South America (offshore Venezuela and Trinidad). These digital data were donated to the University of Texas Institute for Geophysics and represent the largest single, digital reflection survey of the region. Reprocessing of these data included: format correction, filtering, post-stack multiple suppression, and fk migration. Reprocessed data were loaded and interpreted on a workstation. The data straddle a 2,000,000 km2 zone of Paleocene-Recent, time-transgressive, oblique collision between the Caribbean arc system and the passive continental margin of northern South America. Free-air, satellite gravity data shows the remarkable 1000-km-scale continuity of four basement ridges between the uncollided part of the Caribbean arc system (NS-trending Lesser Antilles arc) and the EW-trending collisional area north of Venezuela. The basement ridges involved in the Venezuelan collisional zone include: 1) Aruba-Bonaire-Curacao ridge that can be traced as a continuous feature to the Aves ridge remnant arc of the Lesser Antilles; 2) the partially inverted Blanquilla-Bonaire basin that can be traced into the Grenada back-arc basin; 3) Margarita-Los Testigos platform that can be traced to the Lesser Antilles volcanic arc; and 4) foreland basins and fold-thrust belts of eastern Venezuela (Serrania del Interior and Maturin basin) that can be traced to the Tobago forearc basin and Barbados accretionary prism. Gulfrex data document the progressive change of basinal fault systems from NS-striking normal faults formed in extensional, Lesser Antilles intra-arc settings to rotated and inverted, NE and EW-striking normal faults deformed in the collisional area north of Venezuela. Age of initial shortening of basinal areas and inversion of normal faults setting does not follow the simple, expected pattern of
NASA Astrophysics Data System (ADS)
Auer, L.; Greenhalgh, S. A.; Maurer, H. R.; Marelli, S.; Nuber, A.
2012-04-01
Seismic full waveform inversion is often based on forward modeling in the computationally attractive 2-D domain. Any solution of the 2-D cartesian wave equation inherently carries the implicit assumption of a line source extended in the out-of-plane medium invariant direction. This implies that the source energy in homogeneous media spreads over the surface of an approximately expanding cylinder, such that the wavefield amplitudes (at least in the far field) scale inversely with the square-root of distance. However, realistic point sources like explosives or airguns, fired in a 3-D medium, generate amplitudes that decay inversely with the first power of distance, since the wavefield expands quasi-spherically in all three dimensions. Usually, practitioners correct for this amplitude difference and the associated phase shift of π/4 by transforming the recorded 3-D field data to the approximate 2-D situation by using simplistic, asymptotic filter algorithms. Such filters operate on a square root of time-sample convolutional basis and implicitly assume straight ray paths and a constant velocity medium. The unsubstantiated usage of these asymptotic filters is in contradiction to their well known limitations. In this study, we present an extensive quantitative appraisal of 3D-to-2D data transformation procedures. Our analysis relies on a simple numerical modeling study, based on propagating 3-D and 2-D wavefields through 2-D media and comparing the true 2-D and the filtered 3-D synthetic data. It is shown that the filtering errors are moderate in purely acoustic situations but become substantial in complex media when arrivals overlap each other or ray paths deviate strongly from straight lines. Normalized root-mean-square deviations up to 5% and maximum relative time domain errors of up to 40% were found in high contrast media, when full elastic treatment was considered. In order to examine if this error translates into a deficient model reconstruction in full waveform
Joint inversion of acoustic and resistivity data for the estimation of gas hydrate concentration
Lee, Myung W.
2002-01-01
Downhole log measurements, such as acoustic or electrical resistivity logs, are frequently used to estimate in situ gas hydrate concentrations in the pore space of sedimentary rocks. Usually the gas hydrate concentration is estimated separately based on each log measurement. However, measurements are related to each other through the gas hydrate concentration, so the gas hydrate concentrations can be estimated by jointly inverting available logs. Because the magnitude of slowness of acoustic and resistivity values differs by more than an order of magnitude, a least-squares method, weighted by the inverse of the observed values, is attempted. Estimating the resistivity of connate water and gas hydrate concentration simultaneously is problematic, because the resistivity of connate water is independent of acoustics. In order to overcome this problem, a coupling constant is introduced in the Jacobian matrix. In the use of different logs to estimate gas hydrate concentration, a joint inversion of different measurements is preferred to the averaging of each inversion result.
2.5D complex resistivity modeling and inversion using unstructured grids
NASA Astrophysics Data System (ADS)
Xu, Kaijun; Sun, Jie
2016-04-01
The characteristic of complex resistivity on rock and ore has been recognized by people for a long time. Generally we have used the Cole-Cole Model(CCM) to describe complex resistivity. It has been proved that the electrical anomaly of geologic body can be quantitative estimated by CCM parameters such as direct resistivity(ρ0), chargeability(m), time constant(τ) and frequency dependence(c). Thus it is very important to obtain the complex parameters of geologic body. It is difficult to approximate complex structures and terrain using traditional rectangular grid. In order to enhance the numerical accuracy and rationality of modeling and inversion, we use an adaptive finite-element algorithm for forward modeling of the frequency-domain 2.5D complex resistivity and implement the conjugate gradient algorithm in the inversion of 2.5D complex resistivity. An adaptive finite element method is applied for solving the 2.5D complex resistivity forward modeling of horizontal electric dipole source. First of all, the CCM is introduced into the Maxwell's equations to calculate the complex resistivity electromagnetic fields. Next, the pseudo delta function is used to distribute electric dipole source. Then the electromagnetic fields can be expressed in terms of the primary fields caused by layered structure and the secondary fields caused by inhomogeneities anomalous conductivity. At last, we calculated the electromagnetic fields response of complex geoelectric structures such as anticline, syncline, fault. The modeling results show that adaptive finite-element methods can automatically improve mesh generation and simulate complex geoelectric models using unstructured grids. The 2.5D complex resistivity invertion is implemented based the conjugate gradient algorithm.The conjugate gradient algorithm doesn't need to compute the sensitivity matrix but directly computes the sensitivity matrix or its transpose multiplying vector. In addition, the inversion target zones are
NASA Astrophysics Data System (ADS)
Nenna, Vanessa; Pidlisecky, Adam; Knight, Rosemary
2011-10-01
We apply an extended Kalman filter (EKF) approach to inversion of time-lapse electrical resistivity imaging (ERI) field data. The EKF is a method of time series signal processing that incorporates both a state evolution model, describing changes in the physical system, and an observation model, incorporating the physics of the electrical resistivity measurement. We test the feasibility of using an EKF approach to inverting ERI data collected with 2-D surface array geometries. As a first test, we invert synthetic data generated using a simulated recharge event and water saturation distributions converted to electrical conductivity values using an Archie's law relationship. In the synthetic example we demonstrate the impact that the noise structure of the state evolution and the regularization weight have on EKF-estimated model parameters and errors. We then apply the method to inversion of field data collected to monitor changes in electrical conductivity beneath a recharge pond that is part of an aquifer storage and recovery project in northern California. Using lines of electrodes buried at a depth of 0.25 m when the base of the pond is dry, we monitor the wetting front associated with the diversion of stormflow runoff to the pond. Using field data, we demonstrate that by oversampling in time, we are able to apply the so-called random walk model for the state evolution and to build the model of observation noise directly from collected data. EKF-estimated values track changes in conductivity associated with both increasing water content in subsurface sediments and changes in the properties of the pore water, showing the method is a feasible approach for inversion of time-lapse ERI field data.
Ottaviani, R A; Ashton-Miller, J A; Kothari, S U; Wojtys, E M
1995-01-01
To determine if the height of a basketball shoe alters the maximal inversion and eversion moment that can be actively resisted by the ankle in the frontal plane, we tested 20 healthy, young adult men with no recent ankle injuries. Subjects underwent unipedal functional ankle strength testing under weightbearing conditions at 0 degrees, 16 degrees, and 32 degrees of ankle plantar flexion using a specially designed testing apparatus. Testing was performed with the subject wearing either a low- or a three quarter-top basketball shoe. Shoe height did not significantly affect an individual's ability to actively resist an eversion moment at any angle of ankle plantar flexion. However, tests at 0 degrees of ankle plantar flexion demonstrated that the three quarter-top basketball shoe we tested significantly increased the maximal resistance to an inversion moment by 29.4%. At 16 degrees of ankle plantar flexion, inversion resistance was also significantly improved by 20.4%. These results show that athletic shoe height can significantly increase the active resistance to an inversion moment in moderate ankle plantar flexion. The findings apply to a neutral foot position in the frontal plane, an orientation equivalent to the early phase of a potential ankle sprain. PMID:7573650
Stochastic Inversion of Electrical Resistivity Changes Using a Markov Chain, Monte Carlo Approach
Ramirez, A; Nitao, J; Hanley, W; Aines, R; Glaser, R; Sengupta, S; Dyer, K; Hickling, T; Daily, W
2004-09-21
We describe a stochastic inversion method for mapping subsurface regions where the electrical resistivity is changing. The technique combines prior information, electrical resistance data and forward models to produce subsurface resistivity models that are most consistent with all available data. Bayesian inference and a Metropolis simulation algorithm form the basis for this approach. Attractive features include its ability to: (1) provide quantitative measures of the uncertainty of a generated estimate and, (2) allow alternative model estimates to be identified, compared and ranked. Methods that monitor convergence and summarize important trends of the posterior distribution are introduced. Results from a physical model test and a field experiment were used to assess performance. The stochastic inversions presented provide useful estimates of the most probable location, shape, and volume of the changing region, and the most likely resistivity change. The proposed method is computationally expensive, requiring the use of extensive computational resources to make its application practical.
Nonlinear inversion of electrical resistivity imaging using pruning Bayesian neural networks
NASA Astrophysics Data System (ADS)
Jiang, Fei-Bo; Dai, Qian-Wei; Dong, Li
2016-06-01
Conventional artificial neural networks used to solve electrical resistivity imaging (ERI) inversion problem suffer from overfitting and local minima. To solve these problems, we propose to use a pruning Bayesian neural network (PBNN) nonlinear inversion method and a sample design method based on the K-medoids clustering algorithm. In the sample design method, the training samples of the neural network are designed according to the prior information provided by the K-medoids clustering results; thus, the training process of the neural network is well guided. The proposed PBNN, based on Bayesian regularization, is used to select the hidden layer structure by assessing the effect of each hidden neuron to the inversion results. Then, the hyperparameter α k , which is based on the generalized mean, is chosen to guide the pruning process according to the prior distribution of the training samples under the small-sample condition. The proposed algorithm is more efficient than other common adaptive regularization methods in geophysics. The inversion of synthetic data and field data suggests that the proposed method suppresses the noise in the neural network training stage and enhances the generalization. The inversion results with the proposed method are better than those of the BPNN, RBFNN, and RRBFNN inversion methods as well as the conventional least squares inversion.
NASA Astrophysics Data System (ADS)
Zhai, Cuili; Zhang, Ting
2016-09-01
In this article, we consider the global existence and uniqueness of the solution to the 2D incompressible non-resistive MHD system with non-equilibrium background magnetic field. Our result implies that a strong enough non-equilibrium background magnetic field will guarantee the stability of the nonlinear MHD system. Beside the classical energy method, the interpolation inequalities and the algebraic structure of the equations coming from the incompressibility of the fluid are crucial in our arguments.
Parallelized Three-Dimensional Resistivity Inversion Using Finite Elements And Adjoint State Methods
NASA Astrophysics Data System (ADS)
Schaa, Ralf; Gross, Lutz; Du Plessis, Jaco
2015-04-01
The resistivity method is one of the oldest geophysical exploration methods, which employs one pair of electrodes to inject current into the ground and one or more pairs of electrodes to measure the electrical potential difference. The potential difference is a non-linear function of the subsurface resistivity distribution described by an elliptic partial differential equation (PDE) of the Poisson type. Inversion of measured potentials solves for the subsurface resistivity represented by PDE coefficients. With increasing advances in multichannel resistivity acquisition systems (systems with more than 60 channels and full waveform recording are now emerging), inversion software require efficient storage and solver algorithms. We developed the finite element solver Escript, which provides a user-friendly programming environment in Python to solve large-scale PDE-based problems (see https://launchpad.net/escript-finley). Using finite elements, highly irregular shaped geology and topography can readily be taken into account. For the 3D resistivity problem, we have implemented the secondary potential approach, where the PDE is decomposed into a primary potential caused by the source current and the secondary potential caused by changes in subsurface resistivity. The primary potential is calculated analytically, and the boundary value problem for the secondary potential is solved using nodal finite elements. This approach removes the singularity caused by the source currents and provides more accurate 3D resistivity models. To solve the inversion problem we apply a 'first optimize then discretize' approach using the quasi-Newton scheme in form of the limited-memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) method (see Gross & Kemp 2013). The evaluation of the cost function requires the solution of the secondary potential PDE for each source current and the solution of the corresponding adjoint-state PDE for the cost function gradients with respect to the subsurface
NASA Astrophysics Data System (ADS)
Knight, R.; Weiss, C. J.
2009-05-01
Electrical resistivity imaging (ERI) can be used to obtain information about subsurface structure, properties, and processes for a wide range of near-surface applications. A critical step in the use of ERI is the inversion of the acquired data to obtain an image that displays the magnitude of the electrical resistivity throughout the subsurface region of interest. In order to obtain this image, a number of critical choices need to be made - a choice is made of an inversion algorithm, and further choices are made in terms of its implementation and the incorporation of prior geologic knowledge and constraints. These choices can significantly affect the obtained resistivity image in ways that are not often documented or well understood. It is important that the near-surface geophysics community, and other users of ERI data, engage in an ongoing discussion of how we develop, use, and share inversion algorithms. Two 2D surface electrical resistivity field data sets have been made available to the community to be used as the starting point for this discussion; they are available through the NS website. One data set was acquired in Mack Creek in the H. J. Andrews Experimental Forest to map out the interface between the alluvial sediments and the underlying bedrock. Information about the thickness of the sediments is needed to constrain hydrogeological models and has been difficult to obtain due to the remote location of the site and the inability to install boreholes. An important issue is the need to quantify the uncertainty in the ERI-derived location of the sediment-bedrock interface; a measure of uncertainty in the location of this interface could be included in the hydrogeologic modeling. The second data set was acquired during the monitoring of an infiltration test in the Mojave Desert. Measurements of the dynamically changing distribution of infiltrated water by standard means, using single-point probes installed in the subsurface, cannot adequately characterize
Three dimensional modeling and inversion of Borehole-surface Electrical Resistivity Data
NASA Astrophysics Data System (ADS)
Zhang, Y.; Liu, D.; Liu, Y.; Qin, M.
2013-12-01
After a long time of exploration, many oil fields have stepped into the high water-cut period. It is sorely needed to determining the oil-water distribution and water flooding front. Borehole-surface electrical resistivity tomography (BSERT) system is a low-cost measurement with wide measuring scope and small influence on the reservoir. So it is gaining more and more application in detecting water flooding areas and evaluating residual oil distribution in oil fields. In BSERT system, current is connected with the steel casing of the observation well. The current flows along the long casing and transmits to the surface through inhomogeneous layers. Then received electric potential difference data on the surface can be used to inverse the deep subsurface resistivity distribution. This study presents the 3D modeling and inversion method of electrical resistivity data. In an extensive literature, the steel casing is treated as a transmission line current source with infinite small radius and constant current density. However, in practical multi-layered formations with different resistivity, the current density along the casing is not constant. In this study, the steel casing is modeled by a 2.5e-7 ohm-m physical volume that the casing occupies in the finite element mesh. Radius of the casing can be set to a little bigger than the true radius, and this helps reduce the element number and computation time. The current supply point is set on the center of the top surface of the physical volume. The homogeneous formation modeling result shows the same precision as the transmission line current source model. The multi-layered formation modeling result shows that the current density along the casing is high in the low-resistivity layer, and low in the high-resistivity layer. These results are more reasonable. Moreover, the deviated and horizontal well can be simulated as simple as the vertical well using this modeling method. Based on this forward modeling method, the
Improved resistivity imaging of groundwater solute plumes using POD-based inversion
NASA Astrophysics Data System (ADS)
Oware, E. K.; Moysey, S. M.; Khan, T.
2012-12-01
We propose a new approach for enforcing physics-based regularization in electrical resistivity imaging (ERI) problems. The approach utilizes a basis-constrained inversion where an optimal set of basis vectors is extracted from training data by Proper Orthogonal Decomposition (POD). The key aspect of the approach is that Monte Carlo simulation of flow and transport is used to generate a training dataset, thereby intrinsically capturing the physics of the underlying flow and transport models in a non-parametric form. POD allows for these training data to be projected onto a subspace of the original domain, resulting in the extraction of a basis for the inversion that captures characteristics of the groundwater flow and transport system, while simultaneously allowing for dimensionality reduction of the original problem in the projected space We use two different synthetic transport scenarios in heterogeneous media to illustrate how the POD-based inversion compares with standard Tikhonov and coupled inversion. The first scenario had a single source zone leading to a unimodal solute plume (synthetic #1), whereas, the second scenario had two source zones that produced a bimodal plume (synthetic #2). For both coupled inversion and the POD approach, the conceptual flow and transport model used considered only a single source zone for both scenarios. Results were compared based on multiple metrics (concentration root-mean square error (RMSE), peak concentration, and total solute mass). In addition, results for POD inversion based on 3 different data densities (120, 300, and 560 data points) and varying number of selected basis images (100, 300, and 500) were compared. For synthetic #1, we found that all three methods provided qualitatively reasonable reproduction of the true plume. Quantitatively, the POD inversion performed best overall for each metric considered. Moreover, since synthetic #1 was consistent with the conceptual transport model, a small number of basis
NASA Astrophysics Data System (ADS)
Xiang, S.; Zhang, H.
2015-12-01
2D wide-angle seismic reflection/refraction survey has been widely used to investigate crustal structure and Moho topography. Similarly gravity survey is also very important in the study of local and regional earth features. Seismic survey is sensitive to the seismic velocity parameters and interface variations. For gravity survey, it is sensitive to density parameters of the medium but the resolution along the vertical direction is relatively poor. In this study, we have developed a strategy to jointly invert for seismic velocity model, density model and interface positions using the gravity observations and seismic arrival times from different phases. For the joint inversion of seismic and gravity data, it often relies on the empirical relationship between seismic velocity and density. In comparison, our joint inversion strategy also includes the cross-gradient based structure constraint for seismic velocity and density models in addition to the empirical relationship between them. The objective function for the joint inversion includes data misfit terms for seismic travel times and gravity observations, the cross-gradient constraint, the smoothness terms for two models, and the data misfit term between predicted gravity data based on density model converted from velocity model using the empirical relationship. Each term has its respective weight. We have applied the new joint inversion method to the Riwoqe-Yushu-Maduo profile in northwest China. The profile crosses through the Qiangtang block and Bayan Har block from southwest to northeast, respectively. The 2010 Ms 7.1 Yushu earthquake is located on the profile, around the Ganzi-Yushu fault zone. The joint inversion produces the velocity and density models that are similar in structure and at the same time fit their respective data sets well. Compared to separate seismic inversion using seismic travel times, the joint inversion with gravity data gives a velocity model that better delineates the fault zones. Low
Extraction of semiconductor dopant profiles from spreading resistance data: An inverse problem
NASA Astrophysics Data System (ADS)
Choo, S. C.; Leong, M. S.; Liem, Chandra B. T.; Kong, K. C.
1990-06-01
The traditional method of solving, on a layer-by-layer basis, the inverse problem of extracting resistivity values from spreading resistance measurements is found to produce wildly oscillatory, physically unacceptable resistivity profiles in the case of p-type silicon structures, where a resistivity-dependent probe contact radius is used in conjunction with the probe calibration data. These oscillations are manifestations of the fact that the inverse problem has non-unique solutions; they occur because the problem is inherently ill-posed. The well-known Tikhonov regularisation technique, which converts the present set of highly non-linear integral equations to an equivalent variational problem, is applied to stabilise the solution. Tests are performed on a variety of simulated profiles, and they reveal the existence of an optimum value for the regularisation parameter that is to be used with a second difference expression for the stabiliser of the cost function. When applied to measured spreading resistance data, the technique is found to produce results of reconstruction that are stable and physically reasonable.
NASA Astrophysics Data System (ADS)
Sava, P.; Revil, A.; Karaoulis, M.
2014-08-01
We propose a new, simple and efficient method to image electrical resistivity between a set of wells. Our procedure consists of two steps: first, we map the interfaces between various subsurface formations using seismoelectric conversions; second, we derive the formation resistivity using image-guided cross-well electric tomography. In the first step, we focus seismic energy at a set of points located on a regular grid between wells, which enables us to map the geological formations in terms of heterogeneities in electrical, hydraulic and/or seismic properties. The density of the scanning points (i.e. the seismoelectric image resolution) is related to the wavelength of the seismic impulse used to scan the formations. Each time the seismic energy is focused at a point, the resulting electrical potential burst (equivalent to the one generated by a volumetric seismic source) is recorded remotely at a set of electrodes positioned in wells (the reference electrode can be located on the ground surface or far enough to be considered at infinity). We construct a high-resolution `seismoelectric' image by assigning the electrical potential simulated at these fixed electrodes to the location of the seismic focus. In a follow-up step, the structure of this image is used in image-guided inversion to improve electrical resistivity tomography between the two wells. The structural information from the seismoelectric image is used to impose constraints on the model covariance matrix used in the inversion of the electrical resistivity data. This approach offers new perspectives in recovering fine structure of resistivity (high definition resistivity tomography) between the wells, which cannot be resolved through conventional cross-well resistivity or from seismic tomography alone.
Shetty, N J; Hariprasad, T P N; Sanil, D; Zin, T
2013-11-01
Polytene chromosomes were prepared from the ovarian nurse cells of semi-gravid females of ten insecticide-resistant strains of Anopheles stephensi. Altogether, 16 heterozygous paracentric inversions, namely b/+ (11D-16C) in alphamethrin; i/+ (14B-18A) and h/+ (27B-28A) in DDT; j/+ (14A-16B) in chlorpyrifos; k/+ (11D-16B) in cyfluthrin; l/+ (11A-16C) in deltamethrin; m/+ (14B-15C) and e/+ (32A-33B) in bifenthrin; n/+ (12D-14B), f/+ (33A-36A) and g/+ (33C-34A) in propoxur; o/+ (11A-12D), h/+ (37A-37C) and i/+ (31C-32C) in temephos; d/+ (33D-35C) in carbofuran and a/+ (41C-43B) in neem strains, were reported. No inversions were observed in X chromosome so far. The frequency of inversions in different insecticides was found to be highest in the 2R arm, followed by the 3R arm. Such inversions were not reported in the corresponding susceptible strains or in the parental stocks. PMID:23982309
Three-dimensional anisotropic inversion of resistivity tomography data in an abandoned mine area
NASA Astrophysics Data System (ADS)
Yi, Myeong-Jong; Kim, Jung-Ho; Son, Jeong-Sul
2011-02-01
We have developed an inversion code for three-dimensional (3D) resistivity tomography including the anisotropy effect. The algorithm is based on the finite element approximations for the forward modelling and Active Constraint Balancing method is adopted to enhance the resolving power of the smoothness constraint least-squares inversion. Using numerical experiments, we have shown that anisotropic inversion is viable to get an accurate image of the subsurface when the subsurface shows strong electrical anisotropy. Moreover, anisotropy can be used as additional information in the interpretation of subsurface. This algorithm was also applied to the field dataset acquired in the abandoned old mine area, where a high-rise apartment block has been built up over a mining tunnel. The main purpose of the investigation was to evaluate the safety analysis of the building due to old mining activities. Strong electrical anisotropy has been observed and it was proven to be caused by geological setting of the site. To handle the anisotropy problem, field data were inverted by a 3D anisotropic tomography algorithm and we could obtain 3D subsurface images, which matches well with geology mapping observations. The inversion results have been used to provide the subsurface model for the safety analysis in rock engineering and we could assure the residents that the apartment has no problem in its safety after the completion of investigation works.
NASA Astrophysics Data System (ADS)
Downs, C. M.; Nowicki, R. S.; Rains, M. C.; Kruse, S.
2015-12-01
In west-central Florida, wetland and lake distribution is strongly controlled by karst landforms. Sandhill wetlands and lakes are sand-filled upland basins whose water levels are groundwater driven. Lake dimensions only reach wetland edges during extreme precipitation events. Current wetland classification schemes are inappropriate for identifying sandhill wetlands due to their unique hydrologic regime and ecologic expression. As a result, it is difficult to determine whether or not a wetland is impacted by groundwater pumping, development, and climate change. A better understanding of subsurface structures and how they control the hydrologic regime is necessary for development of an identification and monitoring protocol. Long-term studies record vegetation diversity and distribution, shallow ground water levels and surface water levels. The overall goals are to determine the hydrologic controls (groundwater, seepage, surface water inputs). Most recently a series of geophysical surveys was conducted at select sites in Hernando and Pasco County, Florida. Electrical resistivity and ground penetrating radar were employed to image sand-filled basins and the top of the limestone bedrock and stratigraphy of wetland slopes, respectively. The deepest extent of these sand-filled basins is generally reflected in topography as shallow depressions. Resistivity along inundated wetlands suggests the pools are surface expressions of the surficial aquifer. However, possible breaches in confining clay layers beneath topographic highs between depressions are seen in resistivity profiles as conductive anomalies and in GPR as interruptions in otherwise continuous horizons. These data occur at sites where unconfined and confined water levels are in agreement, suggesting communication between shallow and deep groundwater. Wetland plants are observed outside the historic wetland boundary at many sites, GPR profiles show near-surface layers dipping towards the wetlands at a shallower
NASA Astrophysics Data System (ADS)
Wilkinson, Paul; Chambers, Jonathan; Uhlemann, Sebastian; Meldrum, Philip; Smith, Alister; Dixon, Neil; Loke, Meng Heng
2016-02-01
Reliable tomographic inversion of geoelectrical monitoring data from unstable slopes relies critically on knowing the electrode positions, which may move over time. We develop and present an innovative inverse method to recover movements in both surface directions from geoelectrical measurements made on a grid of monitoring electrodes. For the first time, we demonstrate this method using field data from an active landslide to recover sequences of movement over timescales of days to years. Comparison with GPS measurements demonstrated an accuracy of within 10% of the electrode spacing, sufficient to correct the majority of artifacts that would occur in subsequent image reconstructions if incorrect positions are used. Over short timescales where the corresponding subsurface resistivity changes were smaller, the constraints could be relaxed and an order-of-magnitude better accuracy was achievable. This enabled the onset and acceleration of landslide activity to be detected with a temporal resolution of a few days.
Control of wire heating with resistively guided fast electrons through an inverse conical taper
Robinson, A. P. L. Schmitz, H.; Green, J. S.; Booth, N.; Ridgers, C. P.; Pasley, J.
2015-04-15
The heating of a solid wire embedded in a solid substrate (of lower Z material) with relativistic electrons generated by ultra-intense laser irradiation is considered. Previously, it has been noted that the initial angular distribution of the fast electrons is a highly important factor in the efficacy of the heating [Robinson et al., Phys. Plasmas 20, 122701 (2013)]. We show that, using 3D numerical simulations, the addition of an inverse conical taper at the front of wire can considerably improve the heating of the wire due to the reduction of angular spread of the fast electrons which is caused by transport through the inverse conical taper [Robinson et al., “Guiding of laser-generated fast electrons by exploiting the resistivity-gradients around a conical guide element,” Plasma Phys. Controlled Fusion (to be published)].
NASA Astrophysics Data System (ADS)
Chang, S.; Chang, P.; Wu, Y.; Chen, J.; Huang, C.; Wang, Y.; Chen, W.
2011-12-01
In the study we try to characteristize the hydraulic boundary that separates the Pleistocene Toukeshan Formation from the Holocene sediments with Electrical Resistivity Imaging (ERI) method near the Bizetou Pass of the Choushui River. Comparing 2D Electrical resistivity imaging (ERI) with water level logs and core data from observation wells, we attempted to map the distribution of the shallow groundwater surfaces and the composition variations of the shallow unconfined aquifer from upper to lower fan across the Bizetou Pass. We found that the shallow groundwater surface drops sudenly from about 110m to 70 m after passing the Bizetou Pass from observation wells at the east and west side of the Pass and with about 3 km apart. The inverted resistivity images also show that the hydraulic gradient estimated from the resistivity images is getting larger to about 7% near the Bizetou Pass and then becoming less than 3% in the west the Pass. In addition, we found a significant change from about 1500 ohm-m to 500 ohm-m in the unconfined aquifer after passing the Bizetou Pass from the upstream side. The high resistivity value (about 1500 ohm-m) in the east of Bizetou Pass may represent the compacted Pleistocene Tokeshan sand and gravel layers, and the relatively low resistivity value (about 500 ohm-m) in the west of Bizetou Pass was the loose Holocene alluvium sediments. Since the previous study shows that the Changhua fault is not outcropped at the ground surface, our findings imply that the hydraulic boundary may be due to the less permeable Toukoshan Formation. To confirm this, the future work will extend the study area in North-South direction between Changhua County and Yunlin County.
NASA Astrophysics Data System (ADS)
Kellerer-Pirklbauer, Andreas
2015-04-01
Mountain permafrost covers some 2000 km² of the Austrian Alps which is less than 2.5% of the national territory. Delineating the altitudinal lower limit of permafrost in the mountains of Austria is difficult due the complex topography, the rather sparseness of field verification data and the lack of long-term permafrost monitoring data. Such monitoring data should cover different slope aspects, different elevations, different substrates and different mountain regions of Austria. In this study it was attempted to delineate the lower limit of permafrost at two study sites in the Tauern Range, Austria, applying two-dimensional geoelectrical resistivity tomography (ERT). In addition, multi-annual ground temperature data collected by miniature temperature datalogger (MDT) were used to validate the results. At the study site Hochreichart (maximum elevation 2416 m asl), located in the Seckauer Tauern Range, 14 ERT profiles (lengths 48-196 m; electrode spacing 2, 2.5 or 4 m) were measured at elevations between 1805 and 2416 m asl. Measurements were carried out at two cirques (Reichart, Schöneben) and at the summit plateau of Hochreichart. Results at this site indicate that permafrost lenses are detectable at elevations down to c.1900 m asl at radiation-sheltered sites. Furthermore, at the summit plateau permafrost only occurs as rather small lenses. The ERT-based permafrost pattern is generally confirmed by the MTD data with negative mean annual ground temperature values at only a few monitoring sites. However, the possibility of air-filled cavities causing higher resistive zones faking permafrost existence cannot be excluded because coarse-grained sediments (i.e. relict rock glaciers and autochthonous block fields) are widespread at this study site. At the second study site Kögele Cirque (maximum elevation 3030 m asl) located in the Schober Mountains 12 ERT profiles (lengths 48 m; electrode spacing 2 m) were measured at elevations between 2631 and 2740 m asl. Spatially
Kura, Nura Umar; Ramli, Mohammad Firuz; Ibrahim, Shaharin; Sulaiman, Wan Nor Azmin; Zaudi, Muhammad Amar; Aris, Ahmad Zaharin
2014-01-01
The existing knowledge regarding seawater intrusion and particularly upconing, in which both problems are linked to pumping, entirely relies on theoretical assumptions. Therefore, in this paper, an attempt is made to capture the effects of pumping on seawater intrusion and upconing using 2D resistivity measurement. For this work, two positions, one perpendicular and the other parallel to the sea, were chosen as profile line for resistivity measurement in the coastal area near the pumping wells of Kapas Island, Malaysia. Subsequently, water was pumped out of two pumping wells simultaneously for about five straight hours. Then, immediately after the pumping stopped, resistivity measurements were taken along the two stationed profile lines. This was followed by additional measurements after four and eight hours. The results showed an upconing with low resistivity of about 1-10 Ωm just beneath the pumping well along the first profile line that was taken just after the pumping stopped. The resistivity image also shows an intrusion of saline water (water enriched with diluted salt) from the sea coming towards the pumping well with resistivity values ranging between 10 and 25 Ωm. The subsequent measurements show the recovery of freshwater in the aquifer and how the saline water is gradually diluted or pushed out of the aquifer. Similarly the line parallel to the sea (L2) reveals almost the same result as the first line. However, in the second and third measurements, there were some significant variations which were contrary to the expectation that the freshwater may completely flush out the saline water from the aquifer. These two time series lines show that as the areas with the lowest resistivity (1 Ωm) shrink with time, the low resistivity (10 Ωm) tends to take over almost the entire area implying that the freshwater-saltwater equilibrium zone has already been altered. These results have clearly enhanced our current understanding and add more scientific
Ramli, Mohammad Firuz; Ibrahim, Shaharin; Sulaiman, Wan Nor Azmin; Aris, Ahmad Zaharin
2014-01-01
The existing knowledge regarding seawater intrusion and particularly upconing, in which both problems are linked to pumping, entirely relies on theoretical assumptions. Therefore, in this paper, an attempt is made to capture the effects of pumping on seawater intrusion and upconing using 2D resistivity measurement. For this work, two positions, one perpendicular and the other parallel to the sea, were chosen as profile line for resistivity measurement in the coastal area near the pumping wells of Kapas Island, Malaysia. Subsequently, water was pumped out of two pumping wells simultaneously for about five straight hours. Then, immediately after the pumping stopped, resistivity measurements were taken along the two stationed profile lines. This was followed by additional measurements after four and eight hours. The results showed an upconing with low resistivity of about 1–10 Ωm just beneath the pumping well along the first profile line that was taken just after the pumping stopped. The resistivity image also shows an intrusion of saline water (water enriched with diluted salt) from the sea coming towards the pumping well with resistivity values ranging between 10 and 25 Ωm. The subsequent measurements show the recovery of freshwater in the aquifer and how the saline water is gradually diluted or pushed out of the aquifer. Similarly the line parallel to the sea (L2) reveals almost the same result as the first line. However, in the second and third measurements, there were some significant variations which were contrary to the expectation that the freshwater may completely flush out the saline water from the aquifer. These two time series lines show that as the areas with the lowest resistivity (1 Ωm) shrink with time, the low resistivity (10 Ωm) tends to take over almost the entire area implying that the freshwater-saltwater equilibrium zone has already been altered. These results have clearly enhanced our current understanding and add more scientific
Stefan, Christopher P.; Koehler, Jeffrey W.; Minogue, Timothy D.
2016-01-01
Antibiotic resistance (AR) is an epidemic of increasing magnitude requiring rapid identification and profiling for appropriate and timely therapeutic measures and containment strategies. In this context, ciprofloxacin is part of the first-line of countermeasures against numerous high consequence bacteria. Significant resistance can occur via single nucleotide polymorphisms (SNP) and deletions within ciprofloxacin targeted genes. Ideally, use of ciprofloxacin would be prefaced with AR determination to avoid overuse or misuse of the antibiotic. Here, we describe the development and evaluation of a panel of 44 single-stranded molecular inversion probes (MIPs) coupled to next-generation sequencing (NGS) for the detection of genetic variants known to confer ciprofloxacin resistance in Bacillus anthracis, Yersinia pestis, and Francisella tularensis. Sequencing results demonstrate MIPs capture and amplify targeted regions of interest at significant levels of coverage. Depending on the genetic variant, limits of detection (LOD) for high-throughput pooled sequencing ranged from approximately 300–1800 input genome copies. LODs increased 10-fold in the presence of contaminating human genome DNA. In addition, we show that MIPs can be used as an enrichment step with high resolution melt (HRM) real-time PCR which is a sensitive assay with a rapid time-to-answer. Overall, this technology is a multiplexable upfront enrichment applicable with multiple downstream molecular assays for the detection of targeted genetic regions. PMID:27174456
Stefan, Christopher P; Koehler, Jeffrey W; Minogue, Timothy D
2016-01-01
Antibiotic resistance (AR) is an epidemic of increasing magnitude requiring rapid identification and profiling for appropriate and timely therapeutic measures and containment strategies. In this context, ciprofloxacin is part of the first-line of countermeasures against numerous high consequence bacteria. Significant resistance can occur via single nucleotide polymorphisms (SNP) and deletions within ciprofloxacin targeted genes. Ideally, use of ciprofloxacin would be prefaced with AR determination to avoid overuse or misuse of the antibiotic. Here, we describe the development and evaluation of a panel of 44 single-stranded molecular inversion probes (MIPs) coupled to next-generation sequencing (NGS) for the detection of genetic variants known to confer ciprofloxacin resistance in Bacillus anthracis, Yersinia pestis, and Francisella tularensis. Sequencing results demonstrate MIPs capture and amplify targeted regions of interest at significant levels of coverage. Depending on the genetic variant, limits of detection (LOD) for high-throughput pooled sequencing ranged from approximately 300-1800 input genome copies. LODs increased 10-fold in the presence of contaminating human genome DNA. In addition, we show that MIPs can be used as an enrichment step with high resolution melt (HRM) real-time PCR which is a sensitive assay with a rapid time-to-answer. Overall, this technology is a multiplexable upfront enrichment applicable with multiple downstream molecular assays for the detection of targeted genetic regions. PMID:27174456
Papinutto, Nico; Schlaeger, Regina; Panara, Valentina; Zhu, Alyssa H.; Caverzasi, Eduardo; Stern, William A.; Hauser, Stephen L.; Henry, Roland G.
2015-01-01
The source of inter-subject variability and the influence of age and gender on morphometric characteristics of the spinal cord, such as the total cross-sectional area (TCA), the gray matter (GM) and white matter (WM) areas, currently remain under investigation. Understanding the effect of covariates such as age, gender, brain volumes, and skull- and vertebra-derived metrics on cervical and thoracic spinal cord TCA and GM areas in healthy subjects would be fundamental for exploring compartment specific changes in neurological diseases affecting the spinal cord. Using Magnetic Resonance Imaging at 3T we investigated 32 healthy subjects using a 2D phase sensitive inversion recovery sequence and we measured TCA, GM and WM areas at 4 cervical and thoracic levels of the spinal cord. We assessed age and gender relationships of cord measures and explored associations between cord measures and a) brain volumes and b) skull- and vertebra-derived metrics. Age and gender had a significant effect on TCA, WM and GM areas (with women and elderly having smaller values than men and younger people respectively), but not on the GM area/TCA ratio. The total intracranial volume and C3 vertebra dimensions showed the highest correlations with cord measures. When used in multi-regression models, they reduced cord areas group variability by approximately a third. Age and gender influences on cord measures and normalization strategies here presented might be of use in the study of compartment specific changes in various neurological diseases affecting the spinal cord. PMID:25781178
Lehikoinen, A.; Huttunen, J.M.J.; Finsterle, S.; Kowalsky, M.B.; Kaipio, J.P.
2009-08-01
We propose an approach for imaging the dynamics of complex hydrological processes. The evolution of electrically conductive fluids in porous media is imaged using time-lapse electrical resistance tomography. The related dynamic inversion problem is solved using Bayesian filtering techniques, that is, it is formulated as a sequential state estimation problem in which the target is an evolving posterior probability density of the system state. The dynamical inversion framework is based on the state space representation of the system, which involves the construction of a stochastic evolution model and an observation model. The observation model used in this paper consists of the complete electrode model for ERT, with Archie's law relating saturations to electrical conductivity. The evolution model is an approximate model for simulating flow through partially saturated porous media. Unavoidable modeling and approximation errors in both the observation and evolution models are considered by computing approximate statistics for these errors. These models are then included in the construction of the posterior probability density of the estimated system state. This approximation error method allows the use of approximate - and therefore computationally efficient - observation and evolution models in the Bayesian filtering. We consider a synthetic example and show that the incorporation of an explicit model for the model uncertainties in the state space representation can yield better estimates than a frame-by-frame imaging approach.
Wallin, Erin L.; Johnson, Timothy C.; Greenwood, William J.; Zachara, John M.
2013-03-29
The Hanford 300 Area is located adjacent to the Columbia River in south-central Washington State, USA, and was a former site for nuclear fuel processing operations. Waste disposal practices resulted in persistent unsaturated zone and groundwater contamination, the primary contaminant of concern being uranium. Uranium behavior at the site is intimately linked with river stage driven groundwater-river water exchange such that understanding the nature of river water intrusion into the 300 Area is critical for predicting uranium desorption and transport. In this paper we use time-lapse electrical resistivity tomography (ERT) to image the inland intrusion of river during high stage conditions. We demonstrate a modified time-lapse inversion approach, whereby the transient water table elevation is explicitly modeled by removing regularization constraints across the water table boundary. This implementation was critical for producing meaningful imaging results. We inverted approximately 1200 data sets (400 per line over 3 lines) using high performance computing resources to produce a time-lapse sequence of changes in bulk conductivity caused by river water intrusion during the 2011 spring runoff cycle over approximately 125 days. The resulting time series for each mesh element was then analyzed using common time series analysis to reveal the timing and location of river water intrusion beneath each line. The results reveal non-uniform flows characterized by preferred flow zones where river water enters and exits quickly with stage increase and decrease, and low permeability zones with broader bulk conductivity ‘break through’ curves and longer river water residence times. The time-lapse ERT inversion approach removes the deleterious effects of changing water table elevation and enables remote and spatial continuous groundwater-river water exchange monitoring using surface based ERT arrays under conditions where groundwater and river water conductivity are in contrast.
NASA Astrophysics Data System (ADS)
Tran, A. P.; Dafflon, B.; Hubbard, S. S.; Kowalsky, M. B.; Tokunaga, T. K.; Faybishenko, B.; Long, P.
2014-12-01
Evaluation of spatiotemporal dynamics of heat transport and water flow in terrestrial environments is essential for understanding hydrological and biogeochemical processes. Electrical resistance tomography has been increasingly well used for monitoring subsurface hydrological processes and estimating soil hydraulic properties through coupled hydrogeophysical inversion. However, electrical resistivity depends on a variety of factors such as temperature, which may limit the accuracy of hydrogeophysical inversion. The main objective of this study is to develop a hydrogeophysical inversion framework to enable the incorporation of nonisothermal processes into the hydrogeophysical inversion procedure, and use of this procedure to investigate the effect of hydrological controls on biogeochemical cycles in terrestrial environments. We developed the coupled hydro-thermal-geophysical inversion approach, using the iTOUGH2 framework. In this framework, the heat transport and water flow are simultaneously modeled with TOUGH2 code, which effectively accounts for the multiphase, multi-component and nonisothermal flow in porous media. A flexible approach is used to incorporate petrophysical relationships and uncertainty to link soil moisture and temperature with the electrical resistivity. The developed approach was applied to both synthetic and field case studies. At the DOE subsurface biogeochemistry field site located near Rifle CO, seasonal snowmelt delivers a hydrological pulse to the system, which in turn influences the cycles of nitrogen, carbon and other critical elements. Using the new approach, we carried out numerical inversion of electrical resistance data collected along a 100 m transect at the Rifle site, and compared the results with field investigations of the soil, vadose zone, including the capillary fringe, and groundwater, as well as temperature and tensiometer measurements. Preliminary results show the importance of accounting for nonisothermal conditions to
2D superconductivity by ionic gating
NASA Astrophysics Data System (ADS)
Iwasa, Yoshi
2D superconductivity is attracting a renewed interest due to the discoveries of new highly crystalline 2D superconductors in the past decade. Superconductivity at the oxide interfaces triggered by LaAlO3/SrTiO3 has become one of the promising routes for creation of new 2D superconductors. Also, the MBE grown metallic monolayers including FeSe are also offering a new platform of 2D superconductors. In the last two years, there appear a variety of monolayer/bilayer superconductors fabricated by CVD or mechanical exfoliation. Among these, electric field induced superconductivity by electric double layer transistor (EDLT) is a unique platform of 2D superconductivity, because of its ability of high density charge accumulation, and also because of the versatility in terms of materials, stemming from oxides to organics and layered chalcogenides. In this presentation, the following issues of electric filed induced superconductivity will be addressed; (1) Tunable carrier density, (2) Weak pinning, (3) Absence of inversion symmetry. (1) Since the sheet carrier density is quasi-continuously tunable from 0 to the order of 1014 cm-2, one is able to establish an electronic phase diagram of superconductivity, which will be compared with that of bulk superconductors. (2) The thickness of superconductivity can be estimated as 2 - 10 nm, dependent on materials, and is much smaller than the in-plane coherence length. Such a thin but low resistance at normal state results in extremely weak pinning beyond the dirty Boson model in the amorphous metallic films. (3) Due to the electric filed, the inversion symmetry is inherently broken in EDLT. This feature appears in the enhancement of Pauli limit of the upper critical field for the in-plane magnetic fields. In transition metal dichalcogenide with a substantial spin-orbit interactions, we were able to confirm the stabilization of Cooper pair due to its spin-valley locking. This work has been supported by Grant-in-Aid for Specially
Wang, Chong; Liu, Chang-Ming; Wei, Li-Liang; Shi, Li-Ying; Pan, Zhi-Fen; Mao, Lian-Gen; Wan, Xiao-Chen; Ping, Ze-Peng; Jiang, Ting-Ting; Chen, Zhong-Liang; Li, Zhong-Jie; Li, Ji-Cheng
2016-01-01
The epidemic of pulmonary tuberculosis (TB), especially multidrug-resistance tuberculosis (MDR-TB) presented a major challenge for TB treatment today. We performed iTRAQ labeling coupled with two-dimensional liquid chromatography-tandem mass spectrometry (2D LC-MS/MS) and Solexa sequencing among MDR-TB patients, drug-sensitive tuberculosis (DS-TB) patients, and healthy controls. A total of 50 differentially expressed proteins and 43 differentially expressed miRNAs (fold change >1.50 or <0.60, P<0.05) were identified in the MDR-TB patients compared to both DS-TB patients and healthy controls. We found that 22.00% of differentially expressed proteins and 32.56% of differentially expressed miRNAs were related, and could construct a network mainly in complement and coagulation cascades. Significant differences in CD44 antigen (CD44), coagulation factor XI (F11), kininogen-1 (KNG1), miR-4433b-5p, miR-424-5p, and miR-199b-5p were found among MDR-TB patients, DS-TB patients and healthy controls (P<0.05) by enzyme-linked immunosorbent assay (ELISA) and SYBR green qRT-PCR validation. A strong negative correlation, consistent with the target gene prediction, was found between miR-199b-5p and KNG1 (r=-0.232, P=0.017). Moreover, we established the MDR-TB diagnostic model based on five biomarkers (CD44, KNG1, miR-4433b-5p, miR-424-5p, and miR-199b-5p). Our study proposes potential biomarkers for MDR-TB diagnosis, and also provides a new experimental basis to understand the pathogenesis of MDR-TB. PMID:26884721
NASA Astrophysics Data System (ADS)
Zhang, Ying-Ying; Liu, De-Jun; Ai, Qing-Hui; Qin, Min-Jun
2014-10-01
Electrical resistivity tomography using a steel cased borehole as a long electrode is an advanced technique for geoelectrical survey based on the conventional mise-à-la-masse measurement. In most previous works, the steel casing is simplified as a transmission line current source with an infinitely small radius and constant current density. However, in practical stratified formations with different resistivity values, the current density along the casing cannot be constant. In this study, the steel casing is modeled by a conductive physical volume that the casing occupies in the finite element mesh. The current supply point is set on the center of the top surface of the physical volume. Synthetic modeling, using both a homogenous and layered formation, demonstrates reasonability of the forward modeling method proposed herein. Based on this forward modeling method, the inversion procedure can be implemented by using a freeware R3t (Lancaster University, UK). Inversion results of synthetic modeling data match fairly well with the defined target location and validate that the method works on the inversion of the casing-surface electrical resistivity data. Finally, a field example of Changqing oil field in China is carried out using the inversion method to image water flooding results and to discover wells with great potential to enhance residual oil recovery.
Re-Inversion of Surface Electrical Resistivity Tomography Data from the Hanford Site B-Complex
Johnson, Timothy C.; Wellman, Dawn M.
2013-05-01
This report documents the three-dimensional (3D) inversion results of surface electrical resistivity tomography (ERT) data collected over the Hanford Site B-Complex. The data were collected in order to image the subsurface distribution of electrically conductive vadose zone contamination resulting from both planned releases of contamination into subsurface infiltration galleries (cribs, trenches, and tile fields), as well as unplanned releases from the B, BX, and BY tank farms and/or associated facilities. Electrically conductive contaminants are those which increase the ionic strength of pore fluids compared to native conditions, which comprise most types of solutes released into the subsurface B-Complex. The ERT data were collected and originally inverted as described in detail in report RPP-34690 Rev 0., 2007, which readers should refer to for a detailed description of data collection and waste disposal history. Although the ERT imaging results presented in that report successfully delineated the footprint of vadose zone contamination in areas outside of the tank farms, imaging resolution was not optimized due to the inability of available inversion codes to optimally process the massive ERT data set collected at the site. Recognizing these limitations and the potential for enhanced ERT characterization and time-lapse imaging at contaminated sites, a joint effort was initiated in 2007 by the U.S. Department of Energy – Office of Science (DOE-SC), with later support by the Office of Environmental Management (DOE-EM), and the U.S. Department of Defense (DOD), to develop a high-performance distributed memory parallel 3D ERT inversion code capable of optimally processing large ERT data sets. The culmination of this effort was the development of E4D (Johnson et al., 2010,2012) In 2012, under the Deep Vadose Zone Applied Field Research Initiative (DVZ-AFRI), the U.S. Department of Energy – Richland Operations Office (DOE-RL) and CH2M Hill Plateau Remediation
NASA Astrophysics Data System (ADS)
Ahmad, J.; Ogunsuyi, O.; Schmitt, D. R.; Rokosh, C. D.; Pawlowicz, J. G.
2008-12-01
Buried valleys are common in the regions of the Northern Hemisphere covered by ice sheets during the last glaciation. The valleys are filled by a variety of glacio-lacustrine and glacio-fluvial sedimentation. These valleys are important sources of fresh water, aggregates, and even shallow methane deposits. The surface expressions of the buried valleys, however, are often not apparent and they are often only serendipitously found during drilling for water or petroleum. Geophysical investigations of one deep (~ 350 m) buried valley in northern Alberta that had been located on the basis of geophysical logs were carried out. This buried valley was somewhat unique in that shallow (~ 30 m) methane deposits, providing a significant hazard to drilling, exist. A 10 km high resolution seismic profile (vibrator sweep 14 Hz to 250 Hz, 40 Hz geophone singles at 4 m spacing) was obtained. Significant differences in the raw shot records were apparent across the array due to the lateral differences in compressional wave velocity between the untouched bedrock and the valley fill sediments. Travel time inversion of first arrivals and deeper reflections further quantifies this lateral variation showing that the valley fill and bedrock velocities differ by more than 50% ranging from about 1700 m/s to nearly 3000 m/s, respectively. The reflection seismic image agrees well with the refraction inversion. The gross structure of the steep-sided valley is apparent. The internal architecture, however, shows a variety of clino-form dipping reflectors at the edge of the valley that are possibly related to subglacial sedimentation, a strong dipping reflector that is unconformable with the others and may be representative of recurrent discharge events, and numerous flat lying reflectors that are likely related to lacustrine sedimentation. Co-incident electrical resistivity tomography, too, is largely in agreement with the gross structure. The clay rich bedrock shales are substantially more
NASA Astrophysics Data System (ADS)
Hermans, Thomas; Nguyen, Frédéric
2015-04-01
Electrical resistivity tomography (ERT) has become a popular imaging methodology in a broad range of applications given its large sensitivity to subsurface parameters and its relative simplicity to implement. More particularly, time-lapse ERT is now increasingly used for monitoring purposes in many contexts such as water content, permafrost, landslide, seawater intrusion, solute transport or heat transport experiments. Specific inversion schemes have been developed for time-lapse data sets. However, in contrast with static inversions for which many techniques including geostatistical, minimum support or structural inversion are commonly applied, most of the methodologies for time-lapse inversion still rely on non-physically based spatial and/or temporal smoothing of the parameters or parameter changes. In this work, we propose a time-lapse ERT inversion scheme based on the difference inversion scheme. We replace the standard smoothness-constraint regularization operator by the parameter change covariance matrix. The objective function can be expressed as ψdiff(Δm ) = ||Wd [d - d0 + f(m0)- f(m )]||2 + λ ||||C-Δ0m.5Δm ||||2 where Wd is the data weighting matrix, d and d0 are the data sets corresponding to the considered time-step and to the background, f() is the forward operator, m and m0 are the models corresponding to the considered time-step and to the background, Δm is the parameter change (resistivity), CΔm is the parameter change covariance matrix and λ the regularization parameter. This operator takes into account the correlation between changes in resistivity at different locations through a variogram computed using independent data (e.g., electromagnetic logs). It may vary for subsequent time-steps if the correlation length is time-dependent. The methodology is first validated and compared to the standard smoothness-constraint inversion using a synthetic benchmark simulating the injection of a conductive tracer into a homogeneous aquifer inducing
NASA Astrophysics Data System (ADS)
Nielsen, O. F.; Ploug, C.; Mendoza, J. A.; Martínez, K.
2009-05-01
-soundings. The SMCI here combines the capability of the DC- resistivity method to map resistive areas with the capability of the EM method to map conductive zones in order to produce models that are more in agreement with the expected geological formations. Furthermore, the results from the SMCI inversion are compared with traditional inversion methods. The comparison is made from resistivity maps where the resolution of the method can be evaluated with known lithology, water levels, etc. The results from SMCI inversion demonstrate significant improvement over those derived from standard 2D constrained inversion. The SMCI is applicable in general geological mapping, groundwater surveys, mineral exploration or other surveys where geophysical information of different types can be combined or where geological information is available.
Timothy C. Johnson; Roeof J. Versteeg; Andy Ward; Frederick D. Day-Lewis; Andre Revil
2010-09-01
Electrical geophysical methods have found wide use in the growing discipline of hydrogeophysics, both for characterizing the electrical properties of the subsurface, and for monitoring subsurface processes in terms of the spatiotemporal changes in subsurface conductivity, chargeability, and source currents they govern. Current multichannel and multielectrode data collections systems are able to collect large amounts of data in relatively short periods of time. However, practitioners are often unable to fully utilize these large data sets and the information they contain due to the processing limitations of standard desktop computers. This limitation can be addressed by utilizing the storage and processing capabilities of high-performance parallel computing environments. We present a parallel distributed-memory forward and inverse modeling algorithm for analyzing resistivity and time-domain induced polarization data. The primary components of the parallel computations include distributed computation of the pole solutions in forward mode, distributed storage and computation of the Jacobian matrix in inverse mode, and parallel execution of the inverse equation solver. We demonstrate the corresponding parallel code for three efforts: (1) resistivity characterization of the Hanford 300 Area Integrated Field Research Challenge site in Hanford, WA; (2) resistivity characterization of a volcanic island in the southern Tyrrhenian Sea in Italy; and (3) resistivity and IP monitoring of biostimulation at a superfund site in Brandywine, MD. Inverse analysis of each of these data sets would be limited (or impossible) in a standard serial computing environment which underscores the need for high-performance computing to fully utilize the potential of electrical geophysical methods in hydrogeophysical applications.
Tripp, A.C.; Killpack, T.J.
1981-01-01
IPINV.REV1 is a batch program that is capable of inverting resistivity data to two-dimensional models for a dipole-dipole array. The forward problem is computed using a transmission surface analogy. IPINV.REV1 is capable of inverting resistivity data to two-dimensional models of arbitrary complexity. The two-dimensional forward modeling routine is based on the transmission surface analogy (Madden, 1972). The inversion algorithm is a linearized least-squares technique. Step-size stabilization is provided by either the Box-Kanemaesu (1972) method or by using a Marquardt step. The program uses log derivatives to increase the rate of convergence.
NASA Astrophysics Data System (ADS)
Fernández-López, Sheila; Carrera, Jesús; Ledo, Juanjo; Queralt, Pilar; Luquot, Linda; Martínez, Laura; Bellmunt, Fabián
2016-04-01
Seawater intrusion in aquifers is a complex phenomenon that can be characterized with the help of electric resistivity tomography (ERT) because of the low resistivity of seawater, which underlies the freshwater floating on top. The problem is complex because of the need for joint inversion of electrical and hydraulic (density dependent flow) data. Here we present an adjoint-state algorithm to treat electrical data. This method is a common technique to obtain derivatives of an objective function, depending on potentials with respect to model parameters. The main advantages of it are its simplicity in stationary problems and the reduction of computational cost respect others methodologies. The relationship between the concentration of chlorides and the resistivity values of the field is well known. Also, these resistivities are related to the values of potentials measured using ERT. Taking this into account, it will be possible to define the different resistivities zones from the field data of potential distribution using the basis of inverse problem. In this case, the studied zone is situated in Argentona (Baix Maresme, Catalonia), where the values of chlorides obtained in some wells of the zone are too high. The adjoint-state method will be used to invert the measured data using a new finite element code in C ++ language developed in an open-source framework called Kratos. Finally, the information obtained numerically with our code will be checked with the information obtained with other codes.
He, Weiqi; Kuang, Yongqin; Xing, Xuemin; Simpson, Richard J; Huang, Haidong; Yang, Tao; Chen, Jingmin; Yang, Libin; Liu, Enyu; He, Weifeng; Gu, Jianwen
2014-05-01
Three-dimensional cell culture techniques can better reflect the in vivo characteristics of tumor cells compared with traditional monolayer cultures. Compared with their 2D counterparts, 3D-cultured tumor cells showed enhanced resistance to the cytotoxic T cell-mediated immune response. However, it remains unclear whether 3D-cultured tumor cells have an enhanced resistance to NK cell cytotoxicity. In this study, a total of 363 differentially expressed proteins were identified between the 2D- and 3D-cultured U251 cells by comparative proteomics, and an immune-associated protein-protein interaction (PPI) network based on these differential proteins was constructed by bioinformatics. Within the network, HLA-E, as a molecule for inhibiting NK cell activation, was significantly up-regulated in the 3D-cultured tumor cells. Then, we found that the 3D-cultured U251 cells exhibited potent resistance to NK cell cytotoxicity in vitro and were prone to tumor formation in vivo. The resistance of the 3D-cultured tumor cells to NK cell lysis was mediated by the HLA-E/NKG2A interaction because the administration of antibodies that block either HLA-E or NKG2A completely eliminated this resistance and significantly decreased tumor formation. Taken together, our findings indicate that HLA-E up-regulation in 3D-cultured cells may result in enhanced tumor resistance to NK cell-mediated immune response. PMID:24742303
NASA Astrophysics Data System (ADS)
Liu, Y.; Li, T.; Zhu, C.; Zhang, R.; Wu, Y.
2015-12-01
Three-dimensional (3-D) electromagnetic (EM) forward modelling and inversion continues to be an important issue for the correct interpretation of EM data.To this end,approximate solutions have been developed that allow the construction of relatively fast forward modelling and inversion schemes.We have developed an improved quasi-linear approximation which is more appropriate in solving the linear equation for greatly shortening calculation time.We achieved this by using green's function properties.Then we introduced the improved quasi-linear approximation to spectral induced polarization (SIP) to tackle the problem of the resolution and the efficiency.The localized quasi-linear (LQL) approximation theory is appropriate for multisource array-type surveys assuming that the normal field is slowly varying within the inhomogeneity domain.However,the normal field of attenuates severely which dose not satisfy the assumption of the LQL approximation.As a consenquence,the imaginary part is not accurate when LQL approximation is adopted for the simulation.The improved quasi-linear approximation provide a new approach with the same resolution of QL approximation and much less calculation time.We have also constructed three-dimensional SIP forward modeling based on improved quasi-linear approximation method.It only takes 0.8s for forward modeling when inhomogeneity domain is divided into 2000 blocks.Beyond that, we have introduced the Cole-Cole model to the algorithm and complete the three-dimensional complex resistivity conjugate gradient inversion with parameter restraint.The model trial results show that this method can obtain good inversion results in physical parameters such as zero frequency resistivity, polarization.The results demonstrate the stability and the efficiency of the improved quasi-linear approximation and the method may be a practical solution for3-D EM forward modelling and inversion of SIP.
Technology Transfer Automated Retrieval System (TEKTRAN)
The wheat stem rust resistance gene Sr6, present in several wheat cultivars, confers a high level of resistance against a wide range of races of Puccinia graminis f. sp. tritici. The inheritance and expression of the Sr6 gene is complex. Resistance conferred by Sr6 is influenced by temperature, lig...
NASA Astrophysics Data System (ADS)
Liu, Renqiang; Duan, Yonggang; Tan, Fengqi; Wang, Guochang; Qin, Jianhua; Neupane, Bhupati
2015-10-01
An accurate inversion of original reservoir resistivity is an important problem for waterflood development in oilfields in the middle-late development period. This paper describes the theoretical model of original resistivity recovery for a conglomerate reservoir established by petrophysical models, based on the stratigraphic model of reservoir vertical invasion of the conglomerate reservoir of an oilfield. Likewise two influencing factors of the resistivity change with a water-flooded reservoir were analyzed. The first one is the clay volume decrease due to an injected water wash argillaceous particle and the reservoir resistivity changes are influenced by it, and the other is to inject water to displace crude oil in the pore space leading to the increase of the water-bearing volume. Moreover the conductive ions of the injected water and the original formation water exchange and balance because of their salinity difference, and the reservoir resistivity changes are also influenced by them. Through the analysis of the above influential factors based on the fine identification of conglomerate lithologies the inversion models of three variables, including changes in the amount of clay, the resistivity of the irreducible water and the increase of the water bearing volume, were established by core analysis data, production performance and well logging curves information, and accurately recovered the original reservoir resistivity of the conglomerate. The original oil saturation of the reservoir was calculated according to multiple linear regression models. Finally, the produced index is defined as the difference of the original oil saturation and current oil saturation to the original oil saturation ratio, and it eliminates the effects of conglomerate lithologies and heterogeneity for the quantitative evaluation of flooded layers by the use of the principle of relative value. Compared with traditional flooding sensitive parameters which are oil saturation and water
Takamuku, Shinya; Gomi, Hiroaki
2015-07-22
How our central nervous system (CNS) learns and exploits relationships between force and motion is a fundamental issue in computational neuroscience. While several lines of evidence have suggested that the CNS predicts motion states and signals from motor commands for control and perception (forward dynamics), it remains controversial whether it also performs the 'inverse' computation, i.e. the estimation of force from motion (inverse dynamics). Here, we show that the resistive sensation we experience while moving a delayed cursor, perceived purely from the change in visual motion, provides evidence of the inverse computation. To clearly specify the computational process underlying the sensation, we systematically varied the visual feedback and examined its effect on the strength of the sensation. In contrast to the prevailing theory that sensory prediction errors modulate our perception, the sensation did not correlate with errors in cursor motion due to the delay. Instead, it correlated with the amount of exposure to the forward acceleration of the cursor. This indicates that the delayed cursor is interpreted as a mechanical load, and the sensation represents its visually implied reaction force. Namely, the CNS automatically computes inverse dynamics, using visually detected motions, to monitor the dynamic forces involved in our actions. PMID:26156766
Chen, Yingying; Stabryla, Lisa; Wei, Na
2016-01-01
Development of acetic acid-resistant Saccharomyces cerevisiae is important for economically viable production of biofuels from lignocellulosic biomass, but the goal remains a critical challenge due to limited information on effective genetic perturbation targets for improving acetic acid resistance in the yeast. This study employed a genomic-library-based inverse metabolic engineering approach to successfully identify a novel gene target, WHI2 (encoding a cytoplasmatic globular scaffold protein), which elicited improved acetic acid resistance in S. cerevisiae. Overexpression of WHI2 significantly improved glucose and/or xylose fermentation under acetic acid stress in engineered yeast. The WHI2-overexpressing strain had 5-times-higher specific ethanol productivity than the control in glucose fermentation with acetic acid. Analysis of the expression of WHI2 gene products (including protein and transcript) determined that acetic acid induced endogenous expression of Whi2 in S. cerevisiae. Meanwhile, the whi2Δ mutant strain had substantially higher susceptibility to acetic acid than the wild type, suggesting the important role of Whi2 in the acetic acid response in S. cerevisiae. Additionally, overexpression of WHI2 and of a cognate phosphatase gene, PSR1, had a synergistic effect in improving acetic acid resistance, suggesting that Whi2 might function in combination with Psr1 to elicit the acetic acid resistance mechanism. These results improve our understanding of the yeast response to acetic acid stress and provide a new strategy to breed acetic acid-resistant yeast strains for renewable biofuel production. PMID:26826231
NASA Astrophysics Data System (ADS)
Kirczenow, George
2015-09-01
Valley currents and nonlocal resistances of graphene nanostructures with broken inversion symmetry are considered theoretically in the linear response regime. Scattering state wave functions of electrons entering the nanostructure from the contacts represented by groups of ideal leads are calculated by solving the Lippmann-Schwinger equation and are projected onto the valley state subspaces to obtain the valley velocity fields and total valley currents in the nanostructures. In the tunneling regime when the Fermi energy is in the spectral gap around the Dirac point energy, inversion symmetry breaking is found to result in strong enhancement of the nonlocal four-terminal Büttiker-Landauer resistance and in valley currents several times stronger than the conventional electric current. These strong valley currents are the direct result of the injection of electrons from a contact into the graphene in the tunneling regime. They are chiral and occur near contacts from which electrons are injected into the nanostructure whether or not a net electric current flows through the contact. It is also pointed out that enhanced nonlocal resistances in the linear response regime are not a signature of valley currents arising from the combined effect of the electric field and Berry curvature on the velocities of electrons.
Takamuku, Shinya; Gomi, Hiroaki
2015-01-01
How our central nervous system (CNS) learns and exploits relationships between force and motion is a fundamental issue in computational neuroscience. While several lines of evidence have suggested that the CNS predicts motion states and signals from motor commands for control and perception (forward dynamics), it remains controversial whether it also performs the ‘inverse’ computation, i.e. the estimation of force from motion (inverse dynamics). Here, we show that the resistive sensation we experience while moving a delayed cursor, perceived purely from the change in visual motion, provides evidence of the inverse computation. To clearly specify the computational process underlying the sensation, we systematically varied the visual feedback and examined its effect on the strength of the sensation. In contrast to the prevailing theory that sensory prediction errors modulate our perception, the sensation did not correlate with errors in cursor motion due to the delay. Instead, it correlated with the amount of exposure to the forward acceleration of the cursor. This indicates that the delayed cursor is interpreted as a mechanical load, and the sensation represents its visually implied reaction force. Namely, the CNS automatically computes inverse dynamics, using visually detected motions, to monitor the dynamic forces involved in our actions. PMID:26156766
NASA Astrophysics Data System (ADS)
Maiti, Saumen; Erram, Vinit C.; Gupta, Gautam; Tiwari, Ram Krishna
2012-09-01
SummaryInversion of vertical electrical sounding (VES) data, especially from the crystalline hard rock area, assumes a special significance for groundwater exploration. Here we used a newly developed algorithm based on the Bayesian neural network (BNN) theory combined with Hybrid Monte Carlo (HMC)/Markov Chain Monte Carlo (MCMC) simulation scheme to invert the Direct Current (DC) VES measurements obtained from 30-locations around Tenduli-Vengurla, Sindhudurg district, Maharashtra, India. The inversion results suggest that the top layer is mostly comprised of laterites followed by mixture of clay/clayey sand and garnulites/granite as basement rocks. The source of groundwater appears to be accessible in weathered/semi-weathered layer of laterite/clayey sand that exists within the depth of 10-15 m from the surface. The NW-SE trending major lineaments and its criss-crosses are also identified from the apparent and true resistivity surface map. The pseudo-section at different depths in the western part of the area, near Nivti, shows extensive influence of saltwater intrusion and its impact reaching up to the depth of 30 m from the surface along the coastal area. Our results also show that intrusion of saline water decreases from the western part to the eastern part of the region. Two dimensional modeling of four resistivity profiles from the study region identified two potential groundwater reservoirs; one lying between Path-Tenduli and another in between Mat and Zaraph. The deduced true electrical resistivity section against depth correlates well with available borehole lithology in the area. The results presented here would be useful for interpreting the geological signatures like fractures, major joints and lineaments, which in turn will be helpful for identifying groundwater reservoirs and drainage pattern in the crystalline hard rock area. The newly developed HMC-based BNN method is robust and would provide insights for constraining the geophysical models and
Four-dimensional inversion of resistivity monitoring data through Lp norm minimizations
NASA Astrophysics Data System (ADS)
Kim, Jung-Ho; Supper, Robert; Tsourlos, Panagiotis; Yi, Myeong-Jong
2013-12-01
A new 4-D inversion algorithm is developed so that any data misfits and model roughness in the space and time domains can be selectively minimized, in terms of either the L1 norm or the L2 norm. This study is motivated by the experience that a 4-D inversion adopting full L2 norm minimization may sometimes result in a model that is too smoothly varying with time. It is further encouraged by the realization that a particular criterion of either L1 or L2 norm cannot be universally optimal for accurately reconstructing the subsurface condition. In addition, we try to overcome difficulties of jointly choosing two optimal regularization parameters in space and time domains. To achieve this, we devise automatic determination methods, not only of the Lagrangian multipliers for the space-domain smoothness constraint, but also of the regularization parameter for penalizing the model roughness along the time axis. Both kinds of regularization parameters are actively updated at each iteration, according to variations in data misfit and model roughness. We conducted inversion experiments using synthetic and field monitoring data to test the proposed algorithms, and further to compare the performance of L1 and L2 norm minimizations. Both the synthetic and field data experiments demonstrated that the proposed automatic determination method produced ground changes that were more similar to the true changes than those of approaches using pre-determined parameter values. Inversion experiments showed that L1 norm minimization of the time-domain roughness could reduce the problem of overly smooth model changes when the subsurface changes are locally confined, but that the L2 norm approach would be more reasonable when the changes are expected to be widespread.
NASA Astrophysics Data System (ADS)
Represas, Patricia; Monteiro Santos, F. A.; Ribeiro, Joana A.; Andrade Afonso, A.; Ribeiro, José; Moreira, Mário; Mendes-Victor, L. A.
2015-10-01
The Chaves basin is a pull-apart tectonic depression implanted on granites, schists, and graywackes, and filled with a sedimentary sequence of variable thickness. It is a rather complex structure, as it includes an intricate network of faults and hydrogeological systems. The topography of the basement of the Chaves basin still remains unclear, as no drill hole has ever intersected the bottom of the sediments, and resistivity surveys suffer from severe equivalence issues resulting from the geological setting. In this work, a joint inversion approach of 1D resistivity and gravity data designed for layered environments is used to combine the consistent spatial distribution of the gravity data with the depth sensitivity of the resistivity data. A comparison between the results from the inversion of each data set individually and the results from the joint inversion show that although the joint inversion has more difficulty adjusting to the observed data, it provides more realistic and geologically meaningful models than the ones calculated by the inversion of each data set individually. This work provides a contribution for a better understanding of the Chaves basin, while using the opportunity to study further both the advantages and difficulties comprising the application of the method of joint inversion of gravity and resistivity data.
Energy Science and Technology Software Center (ESTSC)
2005-07-01
Aniso2d is a two-dimensional seismic forward modeling code. The earth is parameterized by an X-Z plane in which the seismic properties Can have monoclinic with x-z plane symmetry. The program uses a user define time-domain wavelet to produce synthetic seismograms anrwhere within the two-dimensional media.
Ward, Stanley H.
1989-01-01
Multiple arrays of electric or magnetic transmitters and receivers are used in a borehole geophysical procedure to obtain a multiplicity of redundant data suitable for processing into a resistivity or induced polarization model of a subsurface region of the earth.
NASA Astrophysics Data System (ADS)
Meqbel, Naser M.; Egbert, Gary D.; Wannamaker, Philip E.; Kelbert, Anna; Schultz, Adam
2014-09-01
Long period (10-20,000 s) magnetotelluric (MT) data are being acquired across the continental USA on a quasi-regular grid of ˜70 km spacing as an electromagnetic component of the National Science Foundation EarthScope/USArray Program. These data are sensitive to fluids, melts, and other orogenic indicators, and thus provide a valuable complement to other components of EarthScope. We present and interpret results of 3-D MT data inversion from 325 sites acquired from 2006-2011 to provide a regional scale view of electrical resistivity from the middle crust to nearly the mantle transition zone, covering an area from NW Washington to NW Colorado. Beneath the active extensional subprovinces in the south-central region, on average we see a resistive upper crust, and then extensive areas of low resistivity in the lower crust and uppermost mantle. Further below, much of the upper half of the upper mantle appears moderately resistive, then subsequently the lower upper mantle becomes moderately conductive. This column suggests a dynamic process of moderately hydrated and fertile deeper upper mantle upwelling during extension, intersection of that material with the damp solidus causing dehydration and melting, and upward exodus of generated mafic melts to pond and exsolve saline fluids near Moho levels. Lithosphere here is very thin. To the east and northeast, thick sections of resistive lithosphere are imaged under the Wyoming and Medicine Hat Cratons. These are punctuated with numerous electrically conductive sutures presumably containing graphitic or sulfide-bearing meta-sediments deeply underthrust and emplaced during ancient collisions. Below Cascadia, the subducting Juan de Fuca and Gorda lithosphere appears highly resistive. Suspected oceanic lithosphere relicts in the central NW part of the model domain also are resistive, including the accreted “Siletzia” terrane beneath the Coast Ranges and Columbia Embayment, and the seismically fast “slab curtain” beneath
Kang, Dongwei; Fang, Zengjun; Li, Zhenyu; Huang, Boshi; Zhang, Heng; Lu, Xueyi; Xu, Haoran; Zhou, Zhongxia; Ding, Xiao; Daelemans, Dirk; De Clercq, Erik; Pannecouque, Christophe; Zhan, Peng; Liu, Xinyong
2016-09-01
We designed and synthesized a series of human immunodeficiency virus type 1 (HIV-1) non-nucleoside reverse transcriptase inhibitors (NNRTIs) with a piperidine-substituted thiophene[3,2-d]pyrimidine scaffold, employing a strategy of structure-based molecular hybridization and substituent decorating. Most of the synthesized compounds exhibited broad-spectrum activity with low (single-digit) nanomolar EC50 values toward a panel of wild-type (WT), single-mutant, and double-mutant HIV-1 strains. Compound 27 was the most potent; compared with ETV, its antiviral efficacy was 3-fold greater against WT, 5-7-fold greater against Y181C, Y188L, E138K, and F227L+V106A, and nearly equipotent against L100I and K103N, though somewhat weaker against K103N+Y181C. Importantly, 27 has lower cytotoxicity (CC50 > 227 μM) and a huge selectivity index (SI) value (ratio of CC50/EC50) of >159101. 27 also showed favorable, drug-like pharmacokinetic and safety properties in rats in vivo. Molecular docking studies and the structure-activity relationships provide important clues for further molecular elaboration. PMID:27541578
Greg Flach, Frank Smith
2011-12-31
Mesh2d is a Fortran90 program designed to generate two-dimensional structured grids of the form [x(i),y(i,j)] where [x,y] are grid coordinates identified by indices (i,j). The x(i) coordinates alone can be used to specify a one-dimensional grid. Because the x-coordinates vary only with the i index, a two-dimensional grid is composed in part of straight vertical lines. However, the nominally horizontal y(i,j0) coordinates along index i are permitted to undulate or otherwise vary. Mesh2d also assigns an integer material type to each grid cell, mtyp(i,j), in a user-specified manner. The complete grid is specified through three separate input files defining the x(i), y(i,j), and mtyp(i,j) variations.
Energy Science and Technology Software Center (ESTSC)
2011-12-31
Mesh2d is a Fortran90 program designed to generate two-dimensional structured grids of the form [x(i),y(i,j)] where [x,y] are grid coordinates identified by indices (i,j). The x(i) coordinates alone can be used to specify a one-dimensional grid. Because the x-coordinates vary only with the i index, a two-dimensional grid is composed in part of straight vertical lines. However, the nominally horizontal y(i,j0) coordinates along index i are permitted to undulate or otherwise vary. Mesh2d also assignsmore » an integer material type to each grid cell, mtyp(i,j), in a user-specified manner. The complete grid is specified through three separate input files defining the x(i), y(i,j), and mtyp(i,j) variations.« less
NASA Astrophysics Data System (ADS)
Lotsch, Bettina V.
2015-07-01
Graphene's legacy has become an integral part of today's condensed matter science and has equipped a whole generation of scientists with an armory of concepts and techniques that open up new perspectives for the postgraphene area. In particular, the judicious combination of 2D building blocks into vertical heterostructures has recently been identified as a promising route to rationally engineer complex multilayer systems and artificial solids with intriguing properties. The present review highlights recent developments in the rapidly emerging field of 2D nanoarchitectonics from a materials chemistry perspective, with a focus on the types of heterostructures available, their assembly strategies, and their emerging properties. This overview is intended to bridge the gap between two major—yet largely disjunct—developments in 2D heterostructures, which are firmly rooted in solid-state chemistry or physics. Although the underlying types of heterostructures differ with respect to their dimensions, layer alignment, and interfacial quality, there is common ground, and future synergies between the various assembly strategies are to be expected.
Ward, S.H.
1989-10-17
Multiple arrays of electric or magnetic transmitters and receivers are used in a borehole geophysical procedure to obtain a multiplicity of redundant data suitable for processing into a resistivity or induced polarization model of a subsurface region of the earth. 30 figs.
Rapid approximate inversion of airborne TEM
NASA Astrophysics Data System (ADS)
Fullagar, Peter K.; Pears, Glenn A.; Reid, James E.; Schaa, Ralf
2015-11-01
Rapid interpretation of large airborne transient electromagnetic (ATEM) datasets is highly desirable for timely decision-making in exploration. Full solution 3D inversion of entire airborne electromagnetic (AEM) surveys is often still not feasible on current day PCs. Therefore, two algorithms to perform rapid approximate 3D interpretation of AEM have been developed. The loss of rigour may be of little consequence if the objective of the AEM survey is regional reconnaissance. Data coverage is often quasi-2D rather than truly 3D in such cases, belying the need for `exact' 3D inversion. Incorporation of geological constraints reduces the non-uniqueness of 3D AEM inversion. Integrated interpretation can be achieved most readily when inversion is applied to a geological model, attributed with lithology as well as conductivity. Geological models also offer several practical advantages over pure property models during inversion. In particular, they permit adjustment of geological boundaries. In addition, optimal conductivities can be determined for homogeneous units. Both algorithms described here can operate on geological models; however, they can also perform `unconstrained' inversion if the geological context is unknown. VPem1D performs 1D inversion at each ATEM data location above a 3D model. Interpretation of cover thickness is a natural application; this is illustrated via application to Spectrem data from central Australia. VPem3D performs 3D inversion on time-integrated (resistive limit) data. Conversion to resistive limits delivers a massive increase in speed since the TEM inverse problem reduces to a quasi-magnetic problem. The time evolution of the decay is lost during the conversion, but the information can be largely recovered by constructing a starting model from conductivity depth images (CDIs) or 1D inversions combined with geological constraints if available. The efficacy of the approach is demonstrated on Spectrem data from Brazil. Both separately and in
Altman, Michael D; Ali, Akbar; Reddy, G S Kiran Kumar; Nalam, Madhavi N L; Anjum, Saima Ghafoor; Cao, Hong; Chellappan, Sripriya; Kairys, Visvaldas; Fernandes, Miguel X; Gilson, Michael K; Schiffer, Celia A; Rana, Tariq M; Tidor, Bruce
2008-05-14
The acquisition of drug-resistant mutations by infectious pathogens remains a pressing health concern, and the development of strategies to combat this threat is a priority. Here we have applied a general strategy, inverse design using the substrate envelope, to develop inhibitors of HIV-1 protease. Structure-based computation was used to design inhibitors predicted to stay within a consensus substrate volume in the binding site. Two rounds of design, synthesis, experimental testing, and structural analysis were carried out, resulting in a total of 51 compounds. Improvements in design methodology led to a roughly 1000-fold affinity enhancement to a wild-type protease for the best binders, from a Ki of 30-50 nM in round one to below 100 pM in round two. Crystal structures of a subset of complexes revealed a binding mode similar to each design that respected the substrate envelope in nearly all cases. All four best binders from round one exhibited broad specificity against a clinically relevant panel of drug-resistant HIV-1 protease variants, losing no more than 6-13-fold affinity relative to wild type. Testing a subset of second-round compounds against the panel of resistant variants revealed three classes of inhibitors: robust binders (maximum affinity loss of 14-16-fold), moderate binders (35-80-fold), and susceptible binders (greater than 100-fold). Although for especially high-affinity inhibitors additional factors may also be important, overall, these results suggest that designing inhibitors using the substrate envelope may be a useful strategy in the development of therapeutics with low susceptibility to resistance. PMID:18412349
Altman, Michael D.; Ali, Akbar; Reddy, G. S. Kiran Kumar; Nalam, Madhavi N. L.; Anjum, Saima Ghafoor; Cao, Hong; Chellappan, Sripriya; Kairys, Visvaldas; Fernandes, Miguel X.; Gilson, Michael K.; Schiffer, Celia A.; Rana, Tariq M.; Tidor, Bruce
2010-01-01
The acquisition of drug-resistance mutations by infectious pathogens remains a pressing health concern, and the development of strategies to combat this threat is a priority. Here we have applied a general strategy, inverse design using the substrate envelope, to develop inhibitors of HIV-1 protease. Structure-based computation was used to design inhibitors predicted to stay within a consensus substrate volume in the binding site. Two rounds of design, synthesis, experimental testing, and structural analysis were carried out, resulting in a total of 51 compounds. Improvements in design methodology led to a roughly 1000-fold affinity enhancement to a wild-type protease for the best binders, from Ki of 30–50 nM in round one to below 100 pM in round two. Crystal structures of a subset of complexes revealed a binding mode similar to each design that respected the substrate envelope in nearly all cases. All four best binders from round one exhibited broad specificity against a clinically relevant panel of drug-resistant HIV-1 protease variants, losing no more than 6–13 fold affinity relative to wild type. Testing a subset of second-round compounds against the panel of resistant variants revealed three classes of inhibitors — robust binders (maximum affinity loss of 14–16 fold), moderate binders (35–80 fold), and susceptible binders (greater than 100 fold). Although for especially high-affinity inhibitors additional factors may also be important, overall, these results suggest that designing inhibitors using the substrate envelope may be a useful strategy in the development of therapeutics with low susceptibility to resistance. PMID:18412349
NASA Astrophysics Data System (ADS)
Nguyen, Frédéric; Hermans, Thomas
2015-04-01
Inversion of time-lapse resistivity data allows obtaining 'snapshots' of changes occurring in monitored systems for applications such as aquifer storage, geothermal heat exchange, site remediation or tracer tests. Based on these snapshots, one can infer qualitative information on the location and morphology of changes occurring in the subsurface but also quantitative estimates on the degree of changes in certain property such as temperature or total dissolved solid content. Analysis of these changes can provide direct insight into flow and transport and associated processes and controlling parameters. However, the reliability of the analysis is dependent on survey geometry, measurement schemes, data error, and regularization. Survey design parameters may be optimized prior to the monitoring survey. Regularization, on the other hand, may be chosen depending on available information collected during the monitoring. Common approaches consider smoothing model changes both in space and time but it is often needed to obtain a sharp temporal anomaly, for example in fractured aquifers. We here propose to use the alternative regularization approach based on minimum gradient support (MGS) (Zhdanov, 2002) for time-lapse surveys which will focus the changes in tomograms snapshots. MGS will limit the occurrences of changes in electrical resistivity but will also restrict the variations of these changes inside the different zones. A commonly encountered difficulty by practitioners in this type of regularization is the choice of an additional parameter, the so-called β, required to define the MGS functional. To the best of our knowledge, there is no commonly accepted or standard methodology to optimize the MGS parameter β. The inversion algorithm used in this study is CRTomo (Kemna 2000). It uses a Gauss-Newton scheme to iteratively minimize an objective function which consists of a data misfit functional and a model constraint functional. A univariate line search is performed
NASA Astrophysics Data System (ADS)
Johnson, T. C.; Versteeg, R. J.; Ward, A. L.; Strickland, C. E.; Greenwood, J.
2009-12-01
The Hanford 300 Area Integrated Field Research Challenge is a DOE funded multi-institution field study designed to better understand the field scale subsurface processes governing uranium transport in the Hanford 300 Area and other DOE sites. The primary focus area includes 28 monitoring wells spaced approximately 10 meters apart and extending downward 18 meters through the unconfined aquifer to a lower fine-grained confining unit. Each well is equipped with 15 permanent electrodes in the vadose zone and 15 removable electrodes in the saturated zone for a total 840 electrodes. To characterize the electrical conductivity structure of the site, a set of 32 overlapping 3D DC resistivity surveys were conducted where each survey included 4 adjacent wells. The complete set of surveys comprises approximately 150,000 resistivity measurements. We demonstrate the characterization by inverting the complete data set using a parallel resistivity modeling/inversion code developed at the Idaho National Laboratory. We compare two approaches; 1) an Occam minimum-structure type inversion, and 2) an inversion using geostatistical and borehole conductivity constraints. The first inversion shows the larger scale conductivity features of the site as resolved by the resistivity data only. The second inversion provides solutions which honor the resistivity data, the directional and zonal semivariograms derived from borehole conductivity logs, and the corresponding borehole conductivities. This approach allows an ensemble of solutions to be generated which incorporate uncertainty in both the semivariograms and the resistivity data. Assuming electrical conductivity is well correlated with hydrogeologic properties at the site, this ensemble of models can be used, along with appropriate petrophyiscal transforms, to generate an ensemble of hydrogeologic models which can ultimately be used to investigate uncertainty in flow and transport processes.
Buitrago-Lopez, Adriana; van den Hooven, Edith H; Rueda-Clausen, Christian F; Serrano, Norma; Ruiz, Alvaro J; Pereira, Mark A; Mueller, Noel T
2015-01-01
Background Low socioeconomic status (SES) has been associated with higher risk of cardiometabolic diseases in developed societies, but investigation of SES and cardiometabolic risk in children in less economically developed populations is sparse. We aimed to examine associations among SES and cardiometabolic risk factors in Colombian children. Methods We used data from a population-based study of 1282 children aged 6–10 years from Bucaramanga, Colombia. SES was classified according to household wealth, living conditions and access to public utilities. Anthropometric and biochemical parameters were measured at a clinic visit. Cardiometabolic risk factors were analysed continuously using linear regression and as binary outcomes—according to established paediatric cut points—using logistic regression to calculate OR and 95% CIs. Results Mean age of the children was 8.4 (SD 1.4) and 51.1% of the sample were boys. Odds of overweight/obesity, abdominal obesity and insulin resistance were greater among higher SES. Compared with the lowest SES stratum, children in the highest SES had higher odds of overweight/obesity (OR=3.25, 95% CI 1.89 to 5.57), abdominal obesity (OR=2.74, 95% CI 1.41 to 5.31) and insulin resistance (OR=2.60, 95% CI 1.81 to 3.71). In contrast, children in the highest SES had lower odds of hypertriglyceridaemia (triglycerides ≥90th centile; OR=0.28, 95% CI 0.14 to 0.54) and low (≤10th centile) high-density lipoprotein (HDL) cholesterol (OR=0.35, 95% CI 0.15 to 0.78). Conclusions In Colombian children, SES is directly associated with obesity and insulin resistance, but inversely associated with dyslipidaemia (hypertriglyceridaemia and low HDL cholesterol). Our findings highlight the need to analyse cardiometabolic risk factors separately in children and to carefully consider a population's level of economic development when studying their social determinants of cardiometabolic disease. PMID:25691273
Commer, Michael; Helwig, Stefan, L.; Hordt, Andreas; Scholl,Carsten; Tezkan, Bulent
2006-06-14
Three long-offset transient electromagnetic (LOTEM) surveyswerecarried out at the active volcano Merapi in Central Java (Indonesia)during the years 1998, 2000, and 2001. The measurements focused on thegeneral resistivity structure of the volcanic edifice at depths of 0.5-2km and the further investigation of a southside anomaly. The measurementswere insufficient for a full 3D inversion scheme, which could enable theimaging of finely discretized resistivity distributions. Therefore, astable, damped least-squares joint-inversion approach is used to optimize3D models with a limited number of parameters. The mode ls feature therealistic simulation of topography, a layered background structure, andadditional coarse 3D blocks representing conductivity anomalies.Twenty-eight LOTEM transients, comprising both horizontal and verticalcomponents of the magnetic induction time derivative, were analyzed. Inview of the few unknowns, we were able to achieve reasonable data fits.The inversion results indicate an upwelling conductor below the summit,suggesting hydrothermal activity in the central volcanic complex. Ashallow conductor due to a magma-filled chamber, at depths down to 1 kmbelow the summit, suggested by earlier seismic studies, is not indicatedby the inversion results. In conjunction with an anomalous-density model,derived from arecent gravity study, our inversion results provideinformation about the southern geological structure resulting from amajor sector collapse during the Middle Merapi period. The density modelallows to assess a porosity range andthus an estimated vertical salinityprofile to explain the high conductivities on a larger scale, extendingbeyond the foothills of Merapi.
NASA Astrophysics Data System (ADS)
Silva, Pablo G.
2010-05-01
The ancient roman city of Baelo Claudia has been subject of several papers on earthquake environmental effects (EEE) and well as earthquake archaeological effects (EAE). During the field training course on archaeoseismology and palaeoseismology conducted in September 2009 (INQUA-IGCP567 Workshop) held at Baelo Claudia, four Electric Resistivity Tomography (ERT) profiles were carried out, by the teams of the Salamanca University (Spain), RWTH Aachen University (Germany) and the Geological Survey of Spain (IGME). ERT surveys were developed in the eastern side of the ancient roman Forum across the unexcavated sector of the archaeological site heading on the 1st Century AD Isis Temple. Each ERT profile was constituted by a 48 multielectrode array with spacing of 2 m resulting in a total length of investigation of around 384 m. ERT lines were separated 10 m each other resulting in a total research area of 3840 m2 to a mean investigation depth of 16 m. The selected survey configurations were Pole-Dipole and Wenner in order to get detailed information about lateral resistivity contrasts, but with a reasonable depth of investigation. The resulting 2D resistivity pseudosections clearly display deformations of the buried roman pavements which propagated in depth within the pre-roman clayey substratum of the Bolonia Bay area.. 3D modelling of the 2D pseudosections indicates that the observed deformations are related to near-surface landsliding, being possible to calculate the minimum volume of mobilized material. ERT 3D imaging allow to refine previous GPR surveys conducted at this same area and to get a subsurface picture of ground deformations caused by repeated earthquakes during the 1st and 3rd Centuries AD. Preliminary calculated volume for the mobilized materials affecting the foundations of the Isis Temple and Forum clearly points to a minimum ESI-07 VIII Intensity validating previous research in the zone. This study has been supported by the Spanish Research Projects
Johnson, Timothy C.; Wellman, Dawn M.
2015-06-26
Electrical resistivity tomography (ERT) has been widely used in environmental applications to study processes associated with subsurface contaminants and contaminant remediation. Anthropogenic alterations in subsurface electrical conductivity associated with contamination often originate from highly industrialized areas with significant amounts of buried metallic infrastructure. The deleterious influence of such infrastructure on imaging results generally limits the utility of ERT where it might otherwise prove useful for subsurface investigation and monitoring. In this manuscript we present a method of accurately modeling the effects of buried conductive infrastructure within the forward modeling algorithm, thereby removing them from the inversion results. The method is implemented in parallel using immersed interface boundary conditions, whereby the global solution is reconstructed from a series of well-conditioned partial solutions. Forward modeling accuracy is demonstrated by comparison with analytic solutions. Synthetic imaging examples are used to investigate imaging capabilities within a subsurface containing electrically conductive buried tanks, transfer piping, and well casing, using both well casings and vertical electrode arrays as current sources and potential measurement electrodes. Results show that, although accurate infrastructure modeling removes the dominating influence of buried metallic features, the presence of metallic infrastructure degrades imaging resolution compared to standard ERT imaging. However, accurate imaging results may be obtained if electrodes are appropriately located.
NASA Astrophysics Data System (ADS)
Wang, Jin; Ma, Jianyong; Zhou, Changhe
2014-11-01
A 3×3 high divergent 2D-grating with period of 3.842μm at wavelength of 850nm under normal incidence is designed and fabricated in this paper. This high divergent 2D-grating is designed by the vector theory. The Rigorous Coupled Wave Analysis (RCWA) in association with the simulated annealing (SA) is adopted to calculate and optimize this 2D-grating.The properties of this grating are also investigated by the RCWA. The diffraction angles are more than 10 degrees in the whole wavelength band, which are bigger than the traditional 2D-grating. In addition, the small period of grating increases the difficulties of fabrication. So we fabricate the 2D-gratings by direct laser writing (DLW) instead of traditional manufacturing method. Then the method of ICP etching is used to obtain the high divergent 2D-grating.
Time-Lapse Joint Inversion of Cross-Well DC Resistivity and Seismic Data: A Numerical Investigation
Time-lapse joint inversion of geophysical data is required to image the evolution of oil reservoirs during production and enhanced oil recovery, CO2 sequestration, geothermal fields during production, and to monitor the evolution of contaminant plumes. Joint inversion schemes red...
Tomosynthesis imaging with 2D scanning trajectories
NASA Astrophysics Data System (ADS)
Khare, Kedar; Claus, Bernhard E. H.; Eberhard, Jeffrey W.
2011-03-01
Tomosynthesis imaging in chest radiography provides volumetric information with the potential for improved diagnostic value when compared to the standard AP or LAT projections. In this paper we explore the image quality benefits of 2D scanning trajectories when coupled with advanced image reconstruction approaches. It is intuitively clear that 2D trajectories provide projection data that is more complete in terms of Radon space filling, when compared with conventional tomosynthesis using a linearly scanned source. Incorporating this additional information for obtaining improved image quality is, however, not a straightforward problem. The typical tomosynthesis reconstruction algorithms are based on direct inversion methods e.g. Filtered Backprojection (FBP) or iterative algorithms that are variants of the Algebraic Reconstruction Technique (ART). The FBP approach is fast and provides high frequency details in the image but at the same time introduces streaking artifacts degrading the image quality. The iterative methods can reduce the image artifacts by using image priors but suffer from a slow convergence rate, thereby producing images lacking high frequency details. In this paper we propose using a fast converging optimal gradient iterative scheme that has advantages of both the FBP and iterative methods in that it produces images with high frequency details while reducing the image artifacts. We show that using favorable 2D scanning trajectories along with the proposed reconstruction method has the advantage of providing improved depth information for structures such as the spine and potentially producing images with more isotropic resolution.
Unitary quantum lattice gas representation of 2D quantum turbulence
NASA Astrophysics Data System (ADS)
Zhang, Bo; Vahala, George; Vahala, Linda; Soe, Min
2011-05-01
Quantum vortex structures and energy cascades are examined for two dimensional quantum turbulence (2D QT) using a special unitary evolution algorithm. The qubit lattice gas (QLG) algorithm, is employed to simulate the weakly-coupled Bose-Einstein condensate (BEC) governed by the Gross-Pitaevskii (GP) equation. A parameter regime is uncovered in which, as in 3D QT, there is a very short Poincare recurrence time. This short recurrence time is destroyed as the nonlinear interaction energy is increased. Energy cascades for 2D QT are considered to examine whether 2D QT exhibits the inverse cascades of 2D classical turbulence. In the parameter regime considered, the spectra analysis reveals no such dual cascades---dual cascades being a hallmark of 2D classical turbulence.
Energy Science and Technology Software Center (ESTSC)
2004-08-01
AnisWave2D is a 2D finite-difference code for a simulating seismic wave propagation in fully anisotropic materials. The code is implemented to run in parallel over multiple processors and is fully portable. A mesh refinement algorithm has been utilized to allow the grid-spacing to be tailored to the velocity model, avoiding the over-sampling of high-velocity materials that usually occurs in fixed-grid schemes.
NASA Astrophysics Data System (ADS)
Tønning, Erik; Polders, Daniel; Callaghan, Paul T.; Engelsen, Søren B.
2007-09-01
This paper demonstrates how the multi-linear PARAFAC model can with advantage be used to decompose 2D diffusion-relaxation correlation NMR spectra prior to 2D-Laplace inversion to the T2- D domain. The decomposition is advantageous for better interpretation of the complex correlation maps as well as for the quantification of extracted T2- D components. To demonstrate the new method seventeen mixtures of wheat flour, starch, gluten, oil and water were prepared and measured with a 300 MHz nuclear magnetic resonance (NMR) spectrometer using a pulsed gradient stimulated echo (PGSTE) pulse sequence followed by a Carr-Purcell-Meiboom-Gill (CPMG) pulse echo train. By varying the gradient strength, 2D diffusion-relaxation data were recorded for each sample. From these double exponentially decaying relaxation data the PARAFAC algorithm extracted two unique diffusion-relaxation components, explaining 99.8% of the variation in the data set. These two components were subsequently transformed to the T2- D domain using 2D-inverse Laplace transformation and quantitatively assigned to the oil and water components of the samples. The oil component was one distinct distribution with peak intensity at D = 3 × 10 -12 m 2 s -1 and T2 = 180 ms. The water component consisted of two broad populations of water molecules with diffusion coefficients and relaxation times centered around correlation pairs: D = 10 -9 m 2 s -1, T2 = 10 ms and D = 3 × 10 -13 m 2 s -1, T2 = 13 ms. Small spurious peaks observed in the inverse Laplace transformation of original complex data were effectively filtered by the PARAFAC decomposition and thus considered artefacts from the complex Laplace transformation. The oil-to-water ratio determined by PARAFAC followed by 2D-Laplace inversion was perfectly correlated with known oil-to-water ratio of the samples. The new method of using PARAFAC prior to the 2D-Laplace inversion proved to have superior potential in analysis of diffusion-relaxation spectra, as it
In situ fluid typing and quantification with 1D and 2D NMR logging.
Sun, Boqin
2007-05-01
In situ nuclear magnetic resonance (NMR) fluid typing has recently gained momentum due to data acquisition and inversion algorithm enhancement of NMR logging tools. T(2) distributions derived from NMR logging contain information on bulk fluids and pore size distributions. However, the accuracy of fluid typing is greatly overshadowed by the overlap between T(2) peaks arising from different fluids with similar apparent T(2) relaxation times. Nevertheless, the shapes of T(2) distributions from different fluid components are often different and can be predetermined. Inversion with predetermined T(2) distributions allows us to perform fluid component decomposition to yield individual fluid volume ratios. Another effective method for in situ fluid typing is two-dimensional (2D) NMR logging, which results in proton population distribution as a function of T(2) relaxation time and fluid diffusion coefficient (or T(1) relaxation time). Since diffusion coefficients (or T(1) relaxation time) for different fluid components can be very different, it is relatively easy to separate oil (especially heavy oil) from water signal in a 2D NMR map and to perform accurate fluid typing. Combining NMR logging with resistivity and/or neutron/density logs provides a third method for in situ fluid typing. We shall describe these techniques with field examples. PMID:17466778
NASA Astrophysics Data System (ADS)
Mayor, Louise
2016-05-01
Graphene might be the most famous example, but there are other 2D materials and compounds too. Louise Mayor explains how these atomically thin sheets can be layered together to create flexible “van der Waals heterostructures”, which could lead to a range of novel applications.
Graphene suspensions for 2D printing
NASA Astrophysics Data System (ADS)
Soots, R. A.; Yakimchuk, E. A.; Nebogatikova, N. A.; Kotin, I. A.; Antonova, I. V.
2016-04-01
It is shown that, by processing a graphite suspension in ethanol or water by ultrasound and centrifuging, it is possible to obtain particles with thicknesses within 1-6 nm and, in the most interesting cases, 1-1.5 nm. Analogous treatment of a graphite suspension in organic solvent yields eventually thicker particles (up to 6-10 nm thick) even upon long-term treatment. Using the proposed ink based on graphene and aqueous ethanol with ethylcellulose and terpineol additives for 2D printing, thin (~5 nm thick) films with sheet resistance upon annealing ~30 MΩ/□ were obtained. With the ink based on aqueous graphene suspension, the sheet resistance was ~5-12 kΩ/□ for 6- to 15-nm-thick layers with a carrier mobility of ~30-50 cm2/(V s).
Radiative heat transfer in 2D Dirac materials
Rodriguez-López, Pablo; Tse, Wang -Kong; Dalvit, Diego A. R.
2015-05-12
We compute the radiative heat transfer between two sheets of 2D Dirac materials, including topological Chern insulators and graphene, within the framework of the local approximation for the optical response of these materials. In this approximation, which neglects spatial dispersion, we derive both numerically and analytically the short-distance asymptotic of the near-field heat transfer in these systems, and show that it scales as the inverse of the distance between the two sheets. In conclusion, we discuss the limitations to the validity of this scaling law imposed by spatial dispersion in 2D Dirac materials.
NASA Astrophysics Data System (ADS)
Fernández Martínez, Juan L.; García Gonzalo, Esperanza; Fernández Álvarez, José P.; Kuzma, Heidi A.; Menéndez Pérez, César O.
2010-05-01
PSO is an optimization technique inspired by the social behavior of individuals in nature (swarms) that has been successfully used in many different engineering fields. In addition, the PSO algorithm can be physically interpreted as a stochastic damped mass-spring system. This analogy has served to introduce the PSO continuous model and to deduce a whole family of PSO algorithms using different finite-differences schemes. These algorithms are characterized in terms of convergence by their respective first and second order stability regions. The performance of these new algorithms is first checked using synthetic functions showing a degree of ill-posedness similar to that found in many geophysical inverse problems having their global minimum located on a very narrow flat valley or surrounded by multiple local minima. Finally we present the application of these PSO algorithms to the analysis and solution of a VES inverse problem associated with a seawater intrusion in a coastal aquifer in southern Spain. PSO family members are successfully compared to other well known global optimization algorithms (binary genetic algorithms and simulated annealing) in terms of their respective convergence curves and the sea water intrusion depth posterior histograms.
Energy Science and Technology Software Center (ESTSC)
2001-01-31
This software reduces the data from two-dimensional kSA MOS program, k-Space Associates, Ann Arbor, MI. Initial MOS data is recorded without headers in 38 columns, with one row of data per acquisition per lase beam tracked. The final MOSS 2d data file is reduced, graphed, and saved in a tab-delimited column format with headers that can be plotted in any graphing software.
NASA Astrophysics Data System (ADS)
Zhang, Zhi-Yong; Tan, Han-Dong; Wang, Kun-Peng; Lin, Chang-Hong; Zhang, Bin; Xie, Mao-Bi
2016-03-01
Traditional two-dimensional (2D) complex resistivity forward modeling is based on Poisson's equation but spectral induced polarization (SIP) data are the coproducts of the induced polarization (IP) and the electromagnetic induction (EMI) effects. This is especially true under high frequencies, where the EMI effect can exceed the IP effect. 2D inversion that only considers the IP effect reduces the reliability of the inversion data. In this paper, we derive differential equations using Maxwell's equations. With the introduction of the Cole-Cole model, we use the finite-element method to conduct 2D SIP forward modeling that considers the EMI and IP effects simultaneously. The data-space Occam method, in which different constraints to the model smoothness and parametric boundaries are introduced, is then used to simultaneously obtain the four parameters of the Cole—Cole model using multi-array electric field data. This approach not only improves the stability of the inversion but also significantly reduces the solution ambiguity. To improve the computational efficiency, message passing interface programming was used to accelerate the 2D SIP forward modeling and inversion. Synthetic datasets were tested using both serial and parallel algorithms, and the tests suggest that the proposed parallel algorithm is robust and efficient.
Rodríguez, Ana; Shimada, Takehiko; Cervera, Magdalena; Redondo, Ana; Alquézar, Berta; Rodrigo, María Jesús; Zacarías, Lorenzo; Palou, Lluís; López, María M; Peña, Leandro
2015-01-01
Volatile organic compounds (VOCs) are secondary metabolites acting as a language for the communication of plants with the environment. In orange fruits, the monoterpene D-limonene accumulates at very high levels in oil glands from the peel. Drastic down-regulation of D-limonene synthase gene expression in the peel of transgenic oranges harboring a D-limonene synthase transgene in antisense (AS) configuration altered the monoterpene profile in oil glands, mainly resulting in reduced accumulation of D-limonene. This led to fruit resistance against Penicillium digitatum (Pd), Xanthomonas citri subsp. citri (Xcc) and other specialized pathogens. Here, we analyze resistance to pathogens in independent AS and empty vector (EV) lines, which have low, medium or high D-limonene concentrations and show that the level of resistance is inversely related to the accumulation of D-limonene in orange peels, thus explaining the need of high D-limonene accumulation in mature oranges in nature for the efficient attraction of specialized microorganism frugivores. PMID:26023857
Rodríguez, Ana; Shimada, Takehiko; Cervera, Magdalena; Redondo, Ana; Alquézar, Berta; Rodrigo, María Jesús; Zacarías, Lorenzo; Palou, Lluís; López, María M; Peña, Leandro
2015-01-01
Volatile organic compounds (VOCs) are secondary metabolites acting as a language for the communication of plants with the environment. In orange fruits, the monoterpene D-limonene accumulates at very high levels in oil glands from the peel. Drastic down-regulation of D-limonene synthase gene expression in the peel of transgenic oranges harboring a D-limonene synthase transgene in antisense (AS) configuration altered the monoterpene profile in oil glands, mainly resulting in reduced accumulation of D-limonene. This led to fruit resistance against Penicillium digitatum (Pd), Xanthomonas citri subsp. citri (Xcc) and other specialized pathogens. Here, we analyze resistance to pathogens in independent AS and empty vector (EV) lines, which have low, medium or high D-limonene concentrations and show that the level of resistance is inversely related to the accumulation of D-limonene in orange peels, thus explaining the need of high D-limonene accumulation in mature oranges in nature for the efficient attraction of specialized microorganism frugivores. PMID:26023857
Nanoimprint lithography: 2D or not 2D? A review
NASA Astrophysics Data System (ADS)
Schift, Helmut
2015-11-01
Nanoimprint lithography (NIL) is more than a planar high-end technology for the patterning of wafer-like substrates. It is essentially a 3D process, because it replicates various stamp topographies by 3D displacement of material and takes advantage of the bending of stamps while the mold cavities are filled. But at the same time, it keeps all assets of a 2D technique being able to pattern thin masking layers like in photon- and electron-based traditional lithography. This review reports about 20 years of development of replication techniques at Paul Scherrer Institut, with a focus on 3D aspects of molding, which enable NIL to stay 2D, but at the same time enable 3D applications which are "more than Moore." As an example, the manufacturing of a demonstrator for backlighting applications based on thermally activated selective topography equilibration will be presented. This technique allows generating almost arbitrary sloped, convex and concave profiles in the same polymer film with dimensions in micro- and nanometer scale.
NASA Astrophysics Data System (ADS)
Revil, A.; Karaoulis, M.; Jardani, A.
2013-12-01
I have developed a simple yet powerful model named POLARIS to explain the electrical conductivity of porous materials including the effect of polarization of the grains. I show that surface conductivity cannot be neglected even for field applications in clean (clay-free) materials but that surface conductivity can be accurately predicted from the quadrature conductivity. The model can also be used to predict reasonably well the permeability and can be used in joint inversion strategies to help hydraulic tomography. I will present the theory, laboratory data (with a variety of clayey and non-clayey rocks and soils), and field data. I will also discuss how electrical resistivity, induced polarization, and self-potential data can be used in concert to infer more information on preferential ground water flow paths. Here again several field examples will be discussed to a variety of case studies relevant to groundwater evaluation and management problems.
Alkerwi, Ala'a; Sauvageot, Nicolas; Crichton, Georgina E; Elias, Merrill F; Stranges, Saverio
2016-05-01
This study examined the association of chocolate consumption with insulin resistance and serum liver enzymes in a national sample of adults in Luxembourg. A random sample of 1153 individuals, aged 18-69 years, was recruited to participate in the cross-sectional Observation of Cardiovascular Risk Factors in Luxembourg study. Chocolate consumption (g/d) was obtained from a semi-quantitative FFQ. Blood glucose and insulin levels were used for the homoeostasis model assessment of insulin resistance (HOMA-IR). Hepatic biomarkers such as serum γ-glutamyl-transpeptidase (γ-GT), serum aspartate transaminase and serum alanine transaminase (ALT) (mg/l) were assessed using standard laboratory assays. Chocolate consumers (81·8 %) were more likely to be younger, physically active, affluent people with higher education levels and fewer chronic co-morbidities. After excluding subjects taking antidiabetic medications, higher chocolate consumption was associated with lower HOMA-IR (β=-0·16, P=0·004), serum insulin levels (β=-0·16, P=0·003) and γ-GT (β=-0·12, P=0·009) and ALT (β=-0·09, P=0·004), after adjustment for age, sex, education, lifestyle and dietary confounding factors, including intakes of fruits and vegetables, alcohol, polyphenol-rich coffee and tea. This study reports an independent inverse relationship between daily chocolate consumption and levels of insulin, HOMA-IR and liver enzymes in adults, suggesting that chocolate consumption may improve liver enzymes and protect against insulin resistance, a well-established risk factor for cardiometabolic disorders. Further observational prospective research and well-designed randomised-controlled studies are needed to confirm this cross-sectional relationship and to comprehend the role and mechanisms that different types of chocolate may play in insulin resistance and cardiometabolic disorders. PMID:26983749
2D Turbulence with Complicated Boundaries
NASA Astrophysics Data System (ADS)
Roullet, G.; McWilliams, J. C.
2014-12-01
We examine the consequences of lateral viscous boundary layers on the 2D turbulence that arises in domains with complicated boundaries (headlands, bays etc). The study is carried out numerically with LES. The numerics are carefully designed to ensure all global conservation laws, proper boundary conditions and a minimal range of dissipation scales. The turbulence dramatically differs from the classical bi-periodic case. Boundary layer separations lead to creation of many small vortices and act as a continuing energy source exciting the inverse cascade of energy throughout the domain. The detachments are very intermittent in time. In free decay, the final state depends on the effective numerical resolution: laminar with a single dominant vortex for low Re and turbulent with many vortices for large enough Re. After very long time, the turbulent end-state exhibits a striking tendency for the emergence of shielded vortices which then interact almost elastically. In the forced case, the boundary layers allow the turbulence to reach a statistical steady state without any artificial hypo-viscosity or other large-scale dissipation. Implications are discussed for the oceanic mesoscale and submesoscale turbulence.
Sparse radar imaging using 2D compressed sensing
NASA Astrophysics Data System (ADS)
Hou, Qingkai; Liu, Yang; Chen, Zengping; Su, Shaoying
2014-10-01
Radar imaging is an ill-posed linear inverse problem and compressed sensing (CS) has been proved to have tremendous potential in this field. This paper surveys the theory of radar imaging and a conclusion is drawn that the processing of ISAR imaging can be denoted mathematically as a problem of 2D sparse decomposition. Based on CS, we propose a novel measuring strategy for ISAR imaging radar and utilize random sub-sampling in both range and azimuth dimensions, which will reduce the amount of sampling data tremendously. In order to handle 2D reconstructing problem, the ordinary solution is converting the 2D problem into 1D by Kronecker product, which will increase the size of dictionary and computational cost sharply. In this paper, we introduce the 2D-SL0 algorithm into the reconstruction of imaging. It is proved that 2D-SL0 can achieve equivalent result as other 1D reconstructing methods, but the computational complexity and memory usage is reduced significantly. Moreover, we will state the results of simulating experiments and prove the effectiveness and feasibility of our method.
Corona-Meraz, Fernanda-Isadora; Navarro-Hernández, Rosa-Elena; Ruíz-Quezada, Sandra-Luz; Madrigal-Ruíz, Perla-Monserrat; Castro-Albarrán, Jorge; Chavarría-Ávila, Efraín; Guzmán-Ornelas, Milton-Omar; Gómez-Bañuelos, Eduardo; Petri, Marcelo-Herón; Ramírez-Cedano, Joel-Isidro; Aguilar-Aldrete, María-Elena; Ríos-Ibarra, Clara; Vázquez-Del Mercado, Mónica
2016-01-01
Background. In obesity there is a subclinical chronic low-grade inflammatory response where insulin resistance (IR) may develop. Chemerin is secreted in white adipose tissue and promotes low-grade inflammatory process, where it expressed CMKLR1 receptor. The role of chemerin and CMKLR1 in inflammatory process secondary to obesity is not defined yet. Methods. Cross-sectional study with 134 individuals classified as with and without obesity by body mass index (BMI) and IR. Body fat storage measurements and metabolic and inflammatory markers were measured by routine methods. Soluble chemerin and basal levels of insulin by ELISA and relative expression of CMKLR1 were evaluated with qPCR and 2−ΔΔCT method. Results. Differences (P < 0.05) were observed between obesity and lean individuals in body fat storage measurements and metabolic-inflammatory markers. Both CMKLR1 expression and chemerin levels were increased in obesity without IR. Soluble chemerin levels correlate with adiposity and metabolic markers (r = 8.8% to 38.5%), P < 0.05. Conclusion. The increment of CMKLR1 expression was associated with insulin production. Increased serum levels of chemerin in obesity were observed, favoring a dysmetabolic response. The results observed in this study suggest that both chemerin and CMKLR1 have opposite expression in the context of low-grade inflammatory response manifested in the development of IR. PMID:27239101
Metrology for graphene and 2D materials
NASA Astrophysics Data System (ADS)
Pollard, Andrew J.
2016-09-01
The application of graphene, a one atom-thick honeycomb lattice of carbon atoms with superlative properties, such as electrical conductivity, thermal conductivity and strength, has already shown that it can be used to benefit metrology itself as a new quantum standard for resistance. However, there are many application areas where graphene and other 2D materials, such as molybdenum disulphide (MoS2) and hexagonal boron nitride (h-BN), may be disruptive, areas such as flexible electronics, nanocomposites, sensing and energy storage. Applying metrology to the area of graphene is now critical to enable the new, emerging global graphene commercial world and bridge the gap between academia and industry. Measurement capabilities and expertise in a wide range of scientific areas are required to address this challenge. The combined and complementary approach of varied characterisation methods for structural, chemical, electrical and other properties, will allow the real-world issues of commercialising graphene and other 2D materials to be addressed. Here, examples of metrology challenges that have been overcome through a multi-technique or new approach are discussed. Firstly, the structural characterisation of defects in both graphene and MoS2 via Raman spectroscopy is described, and how nanoscale mapping of vacancy defects in graphene is also possible using tip-enhanced Raman spectroscopy (TERS). Furthermore, the chemical characterisation and removal of polymer residue on chemical vapour deposition (CVD) grown graphene via secondary ion mass spectrometry (SIMS) is detailed, as well as the chemical characterisation of iron films used to grow large domain single-layer h-BN through CVD growth, revealing how contamination of the substrate itself plays a role in the resulting h-BN layer. In addition, the role of international standardisation in this area is described, outlining the current work ongoing in both the International Organization of Standardization (ISO) and the
NKG2D ligands as therapeutic targets
Spear, Paul; Wu, Ming-Ru; Sentman, Marie-Louise; Sentman, Charles L.
2013-01-01
The Natural Killer Group 2D (NKG2D) receptor plays an important role in protecting the host from infections and cancer. By recognizing ligands induced on infected or tumor cells, NKG2D modulates lymphocyte activation and promotes immunity to eliminate ligand-expressing cells. Because these ligands are not widely expressed on healthy adult tissue, NKG2D ligands may present a useful target for immunotherapeutic approaches in cancer. Novel therapies targeting NKG2D ligands for the treatment of cancer have shown preclinical success and are poised to enter into clinical trials. In this review, the NKG2D receptor and its ligands are discussed in the context of cancer, infection, and autoimmunity. In addition, therapies targeting NKG2D ligands in cancer are also reviewed. PMID:23833565
Development of 2D SIP Data Processing Software for a Metallic Mineral Deposit Exploration
NASA Astrophysics Data System (ADS)
PARK, M.; Kim, K. S.; Seo, H. K.; Son, J.; Park, S.; Kim, C.; Kim, J. H.
2015-12-01
algorithm, we could get more reliable inversion result for the time-domain IP data. Therefore, it is expected the SIP method should be more popular if the developed 2D SIP processing software is widely used.
Interpretation of Magnetic Phase Anomalies over 2D Tabular Bodies
NASA Astrophysics Data System (ADS)
Subrahmanyam, M.
2016-05-01
In this study, phase angle (inverse tangent of the ratio of the horizontal to vertical gradients of magnetic anomalies) profile over two-dimensional tabular bodies has been subjected to detailed analysis for determining the source parameters. Distances between certain characteristic positions on this phase curve are related to the parameters of two-dimensional tabular magnetic sources. In this paper, I have derived the mathematical expressions for these relations. It has been demonstrated here that for locating the origin of the 2D tabular source, knowledge on the type of the model (contact, sheet, dyke, and fault) is not necessary. A procedure is evolved to determine the location, depth, width and magnetization angle of the 2D sources from the mathematical expressions. The method is tested on real field data. The effect of the overlapping bodies is also discussed with two synthetic examples. The interpretation technique is developed for contact, sheet, dike and inclined fault bodies.
Development of the Borehole 2-D Seismic Tomography Software Using MATLAB
NASA Astrophysics Data System (ADS)
Nugraha, A. D.; Syahputra, A.; Fatkhan, F.; Sule, R.; Hendriyana, A.
2011-12-01
We developed 2-D borehole seismic tomography software that we called "EARTHMAX-2D TOMOGRAPHY" to image subsurface physical properties including P-wave and S-wave velocities between two boreholes. We used Graphic User Interface (GUI) facilities of MATLAB programming language to create the software. In this software, we used travel time of seismic waves from source to receiver by using pseudo bending ray tracing method as input for tomography inversion. We can also set up a model parameterization, initial velocity model, ray tracing processes, conduct borehole seismic tomography inversion, and finally visualize the inversion results. The LSQR method was applied to solve of tomography inversion solution. We provided the Checkerboard Test Resolution (CTR) to evaluate the model resolution of the tomography inversion. As validation of this developed software, we tested it for geotechnical purposes. We then conducted data acquisition in the "ITB X-field" that is located on ITB campus. We used two boreholes that have a depth of 39 meters. Seismic wave sources were generated by impulse generator and sparker and then they were recorded by borehole hydrophone string type 3. Later on, we analyzed and picked seismic arrival time as input for tomography inversion. As results, we can image the estimated weathering layer, sediment layer, and basement rock in the field depicted by seismic wave structures. More detailed information about the developed software will be presented. Keywords: borehole, tomography, earthmax-2D, inversion
Perspectives for spintronics in 2D materials
NASA Astrophysics Data System (ADS)
Han, Wei
2016-03-01
The past decade has been especially creative for spintronics since the (re)discovery of various two dimensional (2D) materials. Due to the unusual physical characteristics, 2D materials have provided new platforms to probe the spin interaction with other degrees of freedom for electrons, as well as to be used for novel spintronics applications. This review briefly presents the most important recent and ongoing research for spintronics in 2D materials.
A Geophysical Inversion Model Enhancement Technique Based on the Blind Deconvolution
NASA Astrophysics Data System (ADS)
Zuo, B.; Hu, X.; Li, H.
2011-12-01
A model-enhancement technique is proposed to enhance the geophysical inversion model edges and details without introducing any additional information. Firstly, the theoretic correctness of the proposed geophysical inversion model-enhancement technique is discussed. An inversion MRM (model resolution matrix) convolution approximating PSF (Point Spread Function) method is designed to demonstrate the correctness of the deconvolution model enhancement method. Then, a total-variation regularization blind deconvolution geophysical inversion model-enhancement algorithm is proposed. In previous research, Oldenburg et al. demonstrate the connection between the PSF and the geophysical inverse solution. Alumbaugh et al. propose that more information could be provided by the PSF if we return to the idea of it behaving as an averaging or low pass filter. We consider the PSF as a low pass filter to enhance the inversion model basis on the theory of the PSF convolution approximation. Both the 1D linear and the 2D magnetotelluric inversion examples are used to analyze the validity of the theory and the algorithm. To prove the proposed PSF convolution approximation theory, the 1D linear inversion problem is considered. It shows the ratio of convolution approximation error is only 0.15%. The 2D synthetic model enhancement experiment is presented. After the deconvolution enhancement, the edges of the conductive prism and the resistive host become sharper, and the enhancement result is closer to the actual model than the original inversion model according the numerical statistic analysis. Moreover, the artifacts in the inversion model are suppressed. The overall precision of model increases 75%. All of the experiments show that the structure details and the numerical precision of inversion model are significantly improved, especially in the anomalous region. The correlation coefficient between the enhanced inversion model and the actual model are shown in Fig. 1. The figure
NASA Astrophysics Data System (ADS)
Kasprzak, Marek; Traczyk, Andrzej
2014-06-01
In urbanized areas, particularly in lowland terrains and floors of large river valleys, the natural land configuration is often hard to recognize due to a long history of human activity. Accordingly, archaeological works in cities, which supply knowledge on settlement conditions, are usually accompanied by geological and geomophological research. Lately, data from light detection and ranging (LiDAR) have become a valuable source of information on urban land configuration. Geophysical methods are also becoming increasingly popular in background studies. The paper presents a method of using and linking these sources of spatial information about landforms in such areas. The main aim is to identify to what extent these complementary sources of data and the proposed method can be used in such a specific environment to reconstruct natural, buried terrain morphology. The city of Wrocław in Central Europe serves as an example. To this end geomorphometric studies were conducted with the use of digital elevation models (DEMs) based on LiDAR scanning and derivated land-surface parameters—SAGA Wetness Index, Channel Network Base Level and Altitude above Channel Network. The study also involved determining morphological edges and measurements of the meanders of the Odra, as well as expanding information on the spatial distribution of alluvia and the structure of slope breaks. To this end, geophysical measurements were conducted using the Two-Dimensional Electrical Resistivity Tomography method. Additionally, five typical sequences of man-made ground present within the perimeter of the city were distinguished. As a result, a map of the main landforms of Wrocław is presented. Finally, we argue that although high resolution DEM and derivate land-surface parameters are very useful in terrain analysis, places with thick man-made ground or strongly levelled areas must be recognized by geoarchaeological excavations or geological bore holes. The geophysical survey is useful to
Numerical simulation of ( T 2, T 1) 2D NMR and fluid responses
NASA Astrophysics Data System (ADS)
Tan, Mao-Jin; Zou, You-Long; Zhang, Jin-Yan; Zhao, Xin
2012-12-01
One-dimensional nuclear magnetic resonance (1D NMR) logging technology is limited for fluid typing, while two-dimensional nuclear magnetic resonance (2D NMR) logging can provide more parameters including longitudinal relaxation time ( T 1) and transverse relaxation time ( T 2) relative to fluid types in porous media. Based on the 2D NMR relaxation mechanism in a gradient magnetic field, echo train simulation and 2D NMR inversion are discussed in detail. For 2D NMR inversion, a hybrid inversion method is proposed based on the damping least squares method (LSQR) and an improved truncated singular value decomposition (TSVD) algorithm. A series of spin echoes are first simulated with multiple waiting times ( T W s) in a gradient magnetic field for given fluid models and these synthesized echo trains are inverted by the hybrid method. The inversion results are consistent with given models. Moreover, the numerical simulation of various fluid models such as the gas-water, light oil-water, and vicious oil-water models were carried out with different echo spacings ( T E s) and T W s by this hybrid method. Finally, the influences of different signal-to-noise ratios (SNRs) on inversion results in various fluid models are studied. The numerical simulations show that the hybrid method and optimized observation parameters are applicable to fluid typing of gas-water and oil-water models.
Annotated Bibliography of EDGE2D Use
J.D. Strachan and G. Corrigan
2005-06-24
This annotated bibliography is intended to help EDGE2D users, and particularly new users, find existing published literature that has used EDGE2D. Our idea is that a person can find existing studies which may relate to his intended use, as well as gain ideas about other possible applications by scanning the attached tables.
Staring 2-D hadamard transform spectral imager
Gentry, Stephen M.; Wehlburg, Christine M.; Wehlburg, Joseph C.; Smith, Mark W.; Smith, Jody L.
2006-02-07
A staring imaging system inputs a 2D spatial image containing multi-frequency spectral information. This image is encoded in one dimension of the image with a cyclic Hadamarid S-matrix. The resulting image is detecting with a spatial 2D detector; and a computer applies a Hadamard transform to recover the encoded image.
Anisotropic resistivity tomography: A model study for characterization of fractured rocks
Sasaki, Yutaka
1994-12-31
Since fractured rocks often exhibit anisotropy with respect to hydraulic conductivity, it is expected that anisotropy may play an important factor in describing their electrical properties. Based upon this observation, numerical experiments have been carried out to determine whether anisotropic resistivity tomography can be used for characterization of fractured rocks. In fractured rock masses, the conventional 2-D inversion in which anisotropy is ignored produces peculiar distortions of the resistivity distribution. In contrast, the inversion accounting for anisotropy reconstructs anisotropic background media, as well as conductive anomalies associated with the zones of concentrated fracturing. It is also found that the fracture planes inclined with respect to the strike direction may have no significant effects on 2-D inversion if the angle is within about 20{degree}.
2d PDE Linear Symmetric Matrix Solver
Energy Science and Technology Software Center (ESTSC)
1983-10-01
ICCG2 (Incomplete Cholesky factorized Conjugate Gradient algorithm for 2d symmetric problems) was developed to solve a linear symmetric matrix system arising from a 9-point discretization of two-dimensional elliptic and parabolic partial differential equations found in plasma physics applications, such as resistive MHD, spatial diffusive transport, and phase space transport (Fokker-Planck equation) problems. These problems share the common feature of being stiff and requiring implicit solution techniques. When these parabolic or elliptic PDE''s are discretized withmore » finite-difference or finite-element methods,the resulting matrix system is frequently of block-tridiagonal form. To use ICCG2, the discretization of the two-dimensional partial differential equation and its boundary conditions must result in a block-tridiagonal supermatrix composed of elementary tridiagonal matrices. The incomplete Cholesky conjugate gradient algorithm is used to solve the linear symmetric matrix equation. Loops are arranged to vectorize on the Cray1 with the CFT compiler, wherever possible. Recursive loops, which cannot be vectorized, are written for optimum scalar speed. For matrices lacking symmetry, ILUCG2 should be used. Similar methods in three dimensions are available in ICCG3 and ILUCG3. A general source containing extensions and macros, which must be processed by a pre-compiler to obtain the standard FORTRAN source, is provided along with the standard FORTRAN source because it is believed to be more readable. The pre-compiler is not included, but pre-compilation may be performed by a text editor as described in the UCRL-88746 Preprint.« less
Light field morphing using 2D features.
Wang, Lifeng; Lin, Stephen; Lee, Seungyong; Guo, Baining; Shum, Heung-Yeung
2005-01-01
We present a 2D feature-based technique for morphing 3D objects represented by light fields. Existing light field morphing methods require the user to specify corresponding 3D feature elements to guide morph computation. Since slight errors in 3D specification can lead to significant morphing artifacts, we propose a scheme based on 2D feature elements that is less sensitive to imprecise marking of features. First, 2D features are specified by the user in a number of key views in the source and target light fields. Then the two light fields are warped view by view as guided by the corresponding 2D features. Finally, the two warped light fields are blended together to yield the desired light field morph. Two key issues in light field morphing are feature specification and warping of light field rays. For feature specification, we introduce a user interface for delineating 2D features in key views of a light field, which are automatically interpolated to other views. For ray warping, we describe a 2D technique that accounts for visibility changes and present a comparison to the ideal morphing of light fields. Light field morphing based on 2D features makes it simple to incorporate previous image morphing techniques such as nonuniform blending, as well as to morph between an image and a light field. PMID:15631126
2D materials for nanophotonic devices
NASA Astrophysics Data System (ADS)
Xu, Renjing; Yang, Jiong; Zhang, Shuang; Pei, Jiajie; Lu, Yuerui
2015-12-01
Two-dimensional (2D) materials have become very important building blocks for electronic, photonic, and phononic devices. The 2D material family has four key members, including the metallic graphene, transition metal dichalcogenide (TMD) layered semiconductors, semiconducting black phosphorous, and the insulating h-BN. Owing to the strong quantum confinements and defect-free surfaces, these atomically thin layers have offered us perfect platforms to investigate the interactions among photons, electrons and phonons. The unique interactions in these 2D materials are very important for both scientific research and application engineering. In this talk, I would like to briefly summarize and highlight the key findings, opportunities and challenges in this field. Next, I will introduce/highlight our recent achievements. We demonstrated atomically thin micro-lens and gratings using 2D MoS2, which is the thinnest optical component around the world. These devices are based on our discovery that the elastic light-matter interactions in highindex 2D materials is very strong. Also, I would like to introduce a new two-dimensional material phosphorene. Phosphorene has strongly anisotropic optical response, which creates 1D excitons in a 2D system. The strong confinement in phosphorene also enables the ultra-high trion (charged exciton) binding energies, which have been successfully measured in our experiments. Finally, I will briefly talk about the potential applications of 2D materials in energy harvesting.
Inertial solvation in femtosecond 2D spectra
NASA Astrophysics Data System (ADS)
Hybl, John; Albrecht Ferro, Allison; Farrow, Darcie; Jonas, David
2001-03-01
We have used 2D Fourier transform spectroscopy to investigate polar solvation. 2D spectroscopy can reveal molecular lineshapes beneath ensemble averaged spectra and freeze molecular motions to give an undistorted picture of the microscopic dynamics of polar solvation. The transition from "inhomogeneous" to "homogeneous" 2D spectra is governed by both vibrational relaxation and solvent motion. Therefore, the time dependence of the 2D spectrum directly reflects the total response of the solvent-solute system. IR144, a cyanine dye with a dipole moment change upon electronic excitation, was used to probe inertial solvation in methanol and propylene carbonate. Since the static Stokes' shift of IR144 in each of these solvents is similar, differences in the 2D spectra result from solvation dynamics. Initial results indicate that the larger propylene carbonate responds more slowly than methanol, but appear to be inconsistent with rotational estimates of the inertial response. To disentangle intra-molecular vibrations from solvent motion, the 2D spectra of IR144 will be compared to the time-dependent 2D spectra of the structurally related nonpolar cyanine dye HDITCP.
Internal Photoemission Spectroscopy of 2-D Materials
NASA Astrophysics Data System (ADS)
Nguyen, Nhan; Li, Mingda; Vishwanath, Suresh; Yan, Rusen; Xiao, Shudong; Xing, Huili; Cheng, Guangjun; Hight Walker, Angela; Zhang, Qin
Recent research has shown the great benefits of using 2-D materials in the tunnel field-effect transistor (TFET), which is considered a promising candidate for the beyond-CMOS technology. The on-state current of TFET can be enhanced by engineering the band alignment of different 2D-2D or 2D-3D heterostructures. Here we present the internal photoemission spectroscopy (IPE) approach to determine the band alignments of various 2-D materials, in particular SnSe2 and WSe2, which have been proposed for new TFET designs. The metal-oxide-2-D semiconductor test structures are fabricated and characterized by IPE, where the band offsets from the 2-D semiconductor to the oxide conduction band minimum are determined by the threshold of the cube root of IPE yields as a function of photon energy. In particular, we find that SnSe2 has a larger electron affinity than most semiconductors and can be combined with other semiconductors to form near broken-gap heterojunctions with low barrier heights which can produce a higher on-state current. The details of data analysis of IPE and the results from Raman spectroscopy and spectroscopic ellipsometry measurements will also be presented and discussed.
Towards 2D Bayesian Tomography of Receiver Functions
NASA Astrophysics Data System (ADS)
Ray, A.; Bodin, T.; Key, K.
2014-12-01
Receiver function analysis is a powerful tool widely used to isolate and interpret receiver-side structure effects in teleseismic records. The idea is to deconvolve the vertical component from the horizontal components to produce a time series, thus eliminating the influence of the source and distant path effects. The receiver function is usually migrated and directly interpreted by visual inspection. However, deconvolution is a numerically unstable procedure that needs to be stabilized, and the solution depends on the choice of regularization parameters (e.g. water level and the width of a low pass filter). Since the solution is blurred with multiple reflections from the subsurface that produce apparent discontinuities, qualitative interpretation of receiver functions is subjective. Alternatively, waveforms can be directly inverted for a 1D S-wave velocity model beneath the receiver. An inversion procedure is more quantitative, as a forward model will take into account all possible reflections and conversions. If cast in a Bayesian framework, an inversion also enables one to assess model uncertainties and quantify parameter trade-offs. However, seismologists have preferred migration techniques as they are easier to implement, computationally cheaper, and allow construction of 2D or 3D sections. Inversions have been limited thus far to the 1D case. In this work we present a method for inversion of converted waveforms measured at a number of aligned stations. The unknown model is a 2D vertical cross section parameterized with a variable number of discontinuities, although the forward model used to compute synthetics under individual stations is 1D. Body waves are inverted jointly with surface wave dispersion measurements to reduce the range of possible solutions. The problem is solved with a fully non linear Bayesian inversion scheme where the posterior velocity distribution is sampled with a Markov Chain Monte Carlo Algorithm. Our approach uses the 'trans
Nutter, C.
1980-11-01
GRAV2D is an interactive computer program used for modeling 2-1/2 dimensional gravity data. A forward algorithm is used to give the theoretical attraction of gravity intensity at a station due to a perturbing body given by the initial model. The resultant model can then be adjusted for a better fit by a combination of manual adjustment, one-dimensional automatic search, and Marquardt inversion. GRAV2D has an interactive data management system for data manipulation and display built around subroutines to do a forward problem, a one-dimensional direct search and an inversion. This is a user's guide and documentation for GRAV2D.
NASA Astrophysics Data System (ADS)
Sergienko, Olga
2013-04-01
Since Doug MacAyeal's pioneering studies of the ice-stream basal traction optimizations by control methods, inversions for unknown parameters (e.g., basal traction, accumulation patterns, etc) have become a hallmark of the present-day ice-sheet modeling. The common feature of such inversion exercises is a direct relationship between optimized parameters and observations used in the optimization procedure. For instance, in the standard optimization for basal traction by the control method, ice-stream surface velocities constitute the control data. The optimized basal traction parameters explicitly appear in the momentum equations for the ice-stream velocities (compared to the control data). The inversion for basal traction is carried out by minimization of the cost (or objective, misfit) function that includes the momentum equations facilitated by the Lagrange multipliers. Here, we build upon this idea, and demonstrate how to optimize for parameters indirectly related to observed data using a suite of nested constraints (like Russian dolls) with additional sets of Lagrange multipliers in the cost function. This method opens the opportunity to use data from a variety of sources and types (e.g., velocities, radar layers, surface elevation changes, etc.) in the same optimization process.
Brittle damage models in DYNA2D
Faux, D.R.
1997-09-01
DYNA2D is an explicit Lagrangian finite element code used to model dynamic events where stress wave interactions influence the overall response of the system. DYNA2D is often used to model penetration problems involving ductile-to-ductile impacts; however, with the advent of the use of ceramics in the armor-anti-armor community and the need to model damage to laser optics components, good brittle damage models are now needed in DYNA2D. This report will detail the implementation of four brittle damage models in DYNA2D, three scalar damage models and one tensor damage model. These new brittle damage models are then used to predict experimental results from three distinctly different glass damage problems.
Ginsparg, P.
1991-01-01
These are introductory lectures for a general audience that give an overview of the subject of matrix models and their application to random surfaces, 2d gravity, and string theory. They are intentionally 1.5 years out of date.
Ginsparg, P.
1991-12-31
These are introductory lectures for a general audience that give an overview of the subject of matrix models and their application to random surfaces, 2d gravity, and string theory. They are intentionally 1.5 years out of date.
2D electronic materials for army applications
NASA Astrophysics Data System (ADS)
O'Regan, Terrance; Perconti, Philip
2015-05-01
The record electronic properties achieved in monolayer graphene and related 2D materials such as molybdenum disulfide and hexagonal boron nitride show promise for revolutionary high-speed and low-power electronic devices. Heterogeneous 2D-stacked materials may create enabling technology for future communication and computation applications to meet soldier requirements. For instance, transparent, flexible and even wearable systems may become feasible. With soldier and squad level electronic power demands increasing, the Army is committed to developing and harnessing graphene-like 2D materials for compact low size-weight-and-power-cost (SWAP-C) systems. This paper will review developments in 2D electronic materials at the Army Research Laboratory over the last five years and discuss directions for future army applications.
2-d Finite Element Code Postprocessor
Energy Science and Technology Software Center (ESTSC)
1996-07-15
ORION is an interactive program that serves as a postprocessor for the analysis programs NIKE2D, DYNA2D, TOPAZ2D, and CHEMICAL TOPAZ2D. ORION reads binary plot files generated by the two-dimensional finite element codes currently used by the Methods Development Group at LLNL. Contour and color fringe plots of a large number of quantities may be displayed on meshes consisting of triangular and quadrilateral elements. ORION can compute strain measures, interface pressures along slide lines, reaction forcesmore » along constrained boundaries, and momentum. ORION has been applied to study the response of two-dimensional solids and structures undergoing finite deformations under a wide variety of large deformation transient dynamic and static problems and heat transfer analyses.« less
Chemical Approaches to 2D Materials.
Samorì, Paolo; Palermo, Vincenzo; Feng, Xinliang
2016-08-01
Chemistry plays an ever-increasing role in the production, functionalization, processing and applications of graphene and other 2D materials. This special issue highlights a selection of enlightening chemical approaches to 2D materials, which nicely reflect the breadth of the field and convey the excitement of the individuals involved in it, who are trying to translate graphene and related materials from the laboratory into a real, high-impact technology. PMID:27478083
2D Time-lapse Seismic Tomography Using An Active Time Constraint (ATC) Approach
We propose a 2D seismic time-lapse inversion approach to image the evolution of seismic velocities over time and space. The forward modeling is based on solving the eikonal equation using a second-order fast marching method. The wave-paths are represented by Fresnel volumes rathe...
Extended 2D generalized dilaton gravity theories
NASA Astrophysics Data System (ADS)
de Mello, R. O.
2008-09-01
We show that an anomaly-free description of matter in (1+1) dimensions requires a deformation of the 2D relativity principle, which introduces a non-trivial centre in the 2D Poincaré algebra. Then we work out the reduced phase space of the anomaly-free 2D relativistic particle, in order to show that it lives in a noncommutative 2D Minkowski space. Moreover, we build a Gaussian wave packet to show that a Planck length is well defined in two dimensions. In order to provide a gravitational interpretation for this noncommutativity, we propose to extend the usual 2D generalized dilaton gravity models by a specific Maxwell component, which guages the extra symmetry associated with the centre of the 2D Poincaré algebra. In addition, we show that this extension is a high energy correction to the unextended dilaton theories that can affect the topology of spacetime. Further, we couple a test particle to the general extended dilaton models with the purpose of showing that they predict a noncommutativity in curved spacetime, which is locally described by a Moyal star product in the low energy limit. We also conjecture a probable generalization of this result, which provides strong evidence that the noncommutativity is described by a certain star product which is not of the Moyal type at high energies. Finally, we prove that the extended dilaton theories can be formulated as Poisson Sigma models based on a nonlinear deformation of the extended Poincaré algebra.
2D-2D tunneling field-effect transistors using WSe2/SnSe2 heterostructures
NASA Astrophysics Data System (ADS)
Roy, Tania; Tosun, Mahmut; Hettick, Mark; Ahn, Geun Ho; Hu, Chenming; Javey, Ali
2016-02-01
Two-dimensional materials present a versatile platform for developing steep transistors due to their uniform thickness and sharp band edges. We demonstrate 2D-2D tunneling in a WSe2/SnSe2 van der Waals vertical heterojunction device, where WSe2 is used as the gate controlled p-layer and SnSe2 is the degenerately n-type layer. The van der Waals gap facilitates the regulation of band alignment at the heterojunction, without the necessity of a tunneling barrier. ZrO2 is used as the gate dielectric, allowing the scaling of gate oxide to improve device subthreshold swing. Efficient gate control and clean interfaces yield a subthreshold swing of ˜100 mV/dec for >2 decades of drain current at room temperature, hitherto unobserved in 2D-2D tunneling devices. The subthreshold swing is independent of temperature, which is a clear signature of band-to-band tunneling at the heterojunction. A maximum switching ratio ION/IOFF of 107 is obtained. Negative differential resistance in the forward bias characteristics is observed at 77 K. This work bodes well for the possibilities of two-dimensional materials for the realization of energy-efficient future-generation electronics.
2D FEM Heat Transfer & E&M Field Code
Energy Science and Technology Software Center (ESTSC)
1992-04-02
TOPAZ and TOPAZ2D are two-dimensional implicit finite element computer codes for heat transfer analysis. TOPAZ2D can also be used to solve electrostatic and magnetostatic problems. The programs solve for the steady-state or transient temperature or electrostatic and magnetostatic potential field on two-dimensional planar or axisymmetric geometries. Material properties may be temperature or potential-dependent and either isotropic or orthotropic. A variety of time and temperature-dependent boundary conditions can be specified including temperature, flux, convection, and radiation.more » By implementing the user subroutine feature, users can model chemical reaction kinetics and allow for any type of functional representation of boundary conditions and internal heat generation. The programs can solve problems of diffuse and specular band radiation in an enclosure coupled with conduction in the material surrounding the enclosure. Additional features include thermal contact resistance across an interface, bulk fluids, phase change, and energy balances.« less
2D FEM Heat Transfer & E&M Field Code
1992-04-02
TOPAZ and TOPAZ2D are two-dimensional implicit finite element computer codes for heat transfer analysis. TOPAZ2D can also be used to solve electrostatic and magnetostatic problems. The programs solve for the steady-state or transient temperature or electrostatic and magnetostatic potential field on two-dimensional planar or axisymmetric geometries. Material properties may be temperature or potential-dependent and either isotropic or orthotropic. A variety of time and temperature-dependent boundary conditions can be specified including temperature, flux, convection, and radiation. By implementing the user subroutine feature, users can model chemical reaction kinetics and allow for any type of functional representation of boundary conditions and internal heat generation. The programs can solve problems of diffuse and specular band radiation in an enclosure coupled with conduction in the material surrounding the enclosure. Additional features include thermal contact resistance across an interface, bulk fluids, phase change, and energy balances.
NASA Astrophysics Data System (ADS)
Revil, A.
2015-12-01
Geological expertise and petrophysical relationships can be brought together to provide prior information while inverting multiple geophysical datasets. The merging of such information can result in more realistic solution in the distribution of the model parameters, reducing ipse facto the non-uniqueness of the inverse problem. We consider two level of heterogeneities: facies, described by facies boundaries and heteroegenities inside each facies determined by a correlogram. In this presentation, we pose the geophysical inverse problem in terms of Gaussian random fields with mean functions controlled by petrophysical relationships and covariance functions controlled by a prior geological cross-section, including the definition of spatial boundaries for the geological facies. The petrophysical relationship problem is formulated as a regression problem upon each facies. The inversion of the geophysical data is performed in a Bayesian framework. We demonstrate the usefulness of this strategy using a first synthetic case for which we perform a joint inversion of gravity and galvanometric resistivity data with the stations located at the ground surface. The joint inversion is used to recover the density and resistivity distributions of the subsurface. In a second step, we consider the possibility that the facies boundaries are deformable and their shapes are inverted as well. We use the level set approach to perform such deformation preserving prior topological properties of the facies throughout the inversion. With the help of prior facies petrophysical relationships and topological characteristic of each facies, we make posterior inference about multiple geophysical tomograms based on their corresponding geophysical data misfits. The method is applied to a second synthetic case showing that we can recover the heterogeneities inside the facies, the mean values for the petrophysical properties, and, to some extent, the facies boundaries using the 2D joint inversion of
Calculating tissue shear modulus and pressure by 2D Log-Elastographic methods
McLaughlin, Joyce R; Zhang, Ning; Manduca, Armando
2010-01-01
Shear modulus imaging, often called elastography, enables detection and characterization of tissue abnormalities. In this paper the data is two displacement components obtained from successive MR or ultrasound data sets acquired while the tissue is excited mechanically. A 2D plane strain elastic model is assumed to govern the 2D displacement, u. The shear modulus, μ, is unknown and whether or not the first Lamé parameter, λ, is known the pressure p = λ∇ · u which is present in the plane strain model cannot be measured and is unreliably computed from measured data and can be shown to be an order one quantity in the units kPa. So here we present a 2D Log-Elastographic inverse algorithm that: (1) simultaneously reconstructs the shear modulus, μ, and p, which together satisfy a first order partial differential equation system, with the goal of imaging μ; (2) controls potential exponential growth in the numerical error; and (3) reliably reconstructs the quantity p in the inverse algorithm as compared to the same quantity computed with a forward algorithm. This work generalizes the Log-Elastographic algorithm in [20] which uses one displacement component, is derived assuming the component satisfies the wave equation, and is tested on synthetic data computed with the wave equation model. The 2D Log-Elastographic algorithm is tested on 2D synthetic data and 2D in-vivo data from Mayo Clinic. We also exhibit examples to show that the 2D Log-Elastographic algorithm improves the quality of the recovered images as compared to the Log-Elastographic and Direct Inversion algorithms. PMID:21822349
Optical modulators with 2D layered materials
NASA Astrophysics Data System (ADS)
Sun, Zhipei; Martinez, Amos; Wang, Feng
2016-04-01
Light modulation is an essential operation in photonics and optoelectronics. With existing and emerging technologies increasingly demanding compact, efficient, fast and broadband optical modulators, high-performance light modulation solutions are becoming indispensable. The recent realization that 2D layered materials could modulate light with superior performance has prompted intense research and significant advances, paving the way for realistic applications. In this Review, we cover the state of the art of optical modulators based on 2D materials, including graphene, transition metal dichalcogenides and black phosphorus. We discuss recent advances employing hybrid structures, such as 2D heterostructures, plasmonic structures, and silicon and fibre integrated structures. We also take a look at the future perspectives and discuss the potential of yet relatively unexplored mechanisms, such as magneto-optic and acousto-optic modulation.
Large Area Synthesis of 2D Materials
NASA Astrophysics Data System (ADS)
Vogel, Eric
Transition metal dichalcogenides (TMDs) have generated significant interest for numerous applications including sensors, flexible electronics, heterostructures and optoelectronics due to their interesting, thickness-dependent properties. Despite recent progress, the synthesis of high-quality and highly uniform TMDs on a large scale is still a challenge. In this talk, synthesis routes for WSe2 and MoS2 that achieve monolayer thickness uniformity across large area substrates with electrical properties equivalent to geological crystals will be described. Controlled doping of 2D semiconductors is also critically required. However, methods established for conventional semiconductors, such as ion implantation, are not easily applicable to 2D materials because of their atomically thin structure. Redox-active molecular dopants will be demonstrated which provide large changes in carrier density and workfunction through the choice of dopant, treatment time, and the solution concentration. Finally, several applications of these large-area, uniform 2D materials will be described including heterostructures, biosensors and strain sensors.
2D microwave imaging reflectometer electronics
Spear, A. G.; Domier, C. W. Hu, X.; Muscatello, C. M.; Ren, X.; Luhmann, N. C.; Tobias, B. J.
2014-11-15
A 2D microwave imaging reflectometer system has been developed to visualize electron density fluctuations on the DIII-D tokamak. Simultaneously illuminated at four probe frequencies, large aperture optics image reflections from four density-dependent cutoff surfaces in the plasma over an extended region of the DIII-D plasma. Localized density fluctuations in the vicinity of the plasma cutoff surfaces modulate the plasma reflections, yielding a 2D image of electron density fluctuations. Details are presented of the receiver down conversion electronics that generate the in-phase (I) and quadrature (Q) reflectometer signals from which 2D density fluctuation data are obtained. Also presented are details on the control system and backplane used to manage the electronics as well as an introduction to the computer based control program.
2D microwave imaging reflectometer electronics
NASA Astrophysics Data System (ADS)
Spear, A. G.; Domier, C. W.; Hu, X.; Muscatello, C. M.; Ren, X.; Tobias, B. J.; Luhmann, N. C.
2014-11-01
A 2D microwave imaging reflectometer system has been developed to visualize electron density fluctuations on the DIII-D tokamak. Simultaneously illuminated at four probe frequencies, large aperture optics image reflections from four density-dependent cutoff surfaces in the plasma over an extended region of the DIII-D plasma. Localized density fluctuations in the vicinity of the plasma cutoff surfaces modulate the plasma reflections, yielding a 2D image of electron density fluctuations. Details are presented of the receiver down conversion electronics that generate the in-phase (I) and quadrature (Q) reflectometer signals from which 2D density fluctuation data are obtained. Also presented are details on the control system and backplane used to manage the electronics as well as an introduction to the computer based control program.
2D microwave imaging reflectometer electronics.
Spear, A G; Domier, C W; Hu, X; Muscatello, C M; Ren, X; Tobias, B J; Luhmann, N C
2014-11-01
A 2D microwave imaging reflectometer system has been developed to visualize electron density fluctuations on the DIII-D tokamak. Simultaneously illuminated at four probe frequencies, large aperture optics image reflections from four density-dependent cutoff surfaces in the plasma over an extended region of the DIII-D plasma. Localized density fluctuations in the vicinity of the plasma cutoff surfaces modulate the plasma reflections, yielding a 2D image of electron density fluctuations. Details are presented of the receiver down conversion electronics that generate the in-phase (I) and quadrature (Q) reflectometer signals from which 2D density fluctuation data are obtained. Also presented are details on the control system and backplane used to manage the electronics as well as an introduction to the computer based control program. PMID:25430247
2D-Crystal-Based Functional Inks.
Bonaccorso, Francesco; Bartolotta, Antonino; Coleman, Jonathan N; Backes, Claudia
2016-08-01
The possibility to produce and process graphene, related 2D crystals, and heterostructures in the liquid phase makes them promising materials for an ever-growing class of applications as composite materials, sensors, in flexible optoelectronics, and energy storage and conversion. In particular, the ability to formulate functional inks with on-demand rheological and morphological properties, i.e., lateral size and thickness of the dispersed 2D crystals, is a step forward toward the development of industrial-scale, reliable, inexpensive printing/coating processes, a boost for the full exploitation of such nanomaterials. Here, the exfoliation strategies of graphite and other layered crystals are reviewed, along with the advances in the sorting of lateral size and thickness of the exfoliated sheets together with the formulation of functional inks and the current development of printing/coating processes of interest for the realization of 2D-crystal-based devices. PMID:27273554
Experimental validation of 2D profile photoresist shrinkage model
NASA Astrophysics Data System (ADS)
Bunday, Benjamin; Cordes, Aaron; Self, Andy; Ferry, Lorena; Danilevsky, Alex
2011-03-01
For many years, lithographic resolution has been the main obstacle in allowing the pace of transistor densification to meet Moore's Law. For the 32 nm node and beyond, new lithography techniques will be used, including immersion ArF (iArF) lithography and extreme ultraviolet lithography (EUVL). As in the past, these techniques will use new types of photoresists with the capability to print smaller feature widths and pitches. These smaller feature sizes will also require the use of thinner layers of photoresists, such as under 100 nm. In previous papers, we focused on ArF and iArF photoresist shrinkage. We evaluated the magnitude of shrinkage for both R&D and mature resists as a function of chemical formulation, lithographic sensitivity, scanning electron microscope (SEM) beam condition, and feature size. Shrinkage results were determined by the well accepted methodology described in SEMATECH's CD-SEM Unified Specification. In other associated works, we first developed a 1-D model for resist shrinkage for the bottom linewidth and then a 2-D profile model that accounted for shrinkage of all aspects of a trapezoidal profile along a given linescan. A fundamental understanding of the phenomenology of the shrinkage trends was achieved, including how the shrinkage behaves differently for different sized and shaped features. In the 1-D case, calibration of the parameters to describe the photoresist material and the electron beam was all that was required to fit the models to real shrinkage data, as long as the photoresist was thick enough that the beam could not penetrate the entire layer of resist. The later 2-D model included improvements for solving the CD shrinkage in thin photoresists, which is now of great interest for upcoming realistic lithographic processing to explore the change in resist profile with electron dose and to predict the influence of initial resist profile on shrinkage characteristics. The 2-D model also included shrinkage due to both the primary
The 2D lingual appliance system.
Cacciafesta, Vittorio
2013-09-01
The two-dimensional (2D) lingual bracket system represents a valuable treatment option for adult patients seeking a completely invisible orthodontic appliance. The ease of direct or simplified indirect bonding of 2D lingual brackets in combination with low friction mechanics makes it possible to achieve a good functional and aesthetic occlusion, even in the presence of a severe malocclusion. The use of a self-ligating bracket significantly reduces chair-side time for the orthodontist, and the low-profile bracket design greatly improves patient comfort. PMID:24005953
Inkjet printing of 2D layered materials.
Li, Jiantong; Lemme, Max C; Östling, Mikael
2014-11-10
Inkjet printing of 2D layered materials, such as graphene and MoS2, has attracted great interests for emerging electronics. However, incompatible rheology, low concentration, severe aggregation and toxicity of solvents constitute critical challenges which hamper the manufacturing efficiency and product quality. Here, we introduce a simple and general technology concept (distillation-assisted solvent exchange) to efficiently overcome these challenges. By implementing the concept, we have demonstrated excellent jetting performance, ideal printing patterns and a variety of promising applications for inkjet printing of 2D layered materials. PMID:25169938
Measurement of 2D birefringence distribution
NASA Astrophysics Data System (ADS)
Noguchi, Masato; Ishikawa, Tsuyoshi; Ohno, Masahiro; Tachihara, Satoru
1992-10-01
A new measuring method of 2-D birefringence distribution has been developed. It has not been an easy job to get a birefringence distribution in an optical element with conventional ellipsometry because of its lack of scanning means. Finding an analogy between the rotating analyzer method in ellipsometry and the phase-shifting method in recently developed digital interferometry, we have applied the phase-shifting algorithm to ellipsometry, and have developed a new method that makes the measurement of 2-D birefringence distribution easy and possible. The system contains few moving parts, assuring reliability, and measures a large area of a sample at one time, making the measuring time very short.
Role of defects in frictional properties of 2-D materials
NASA Astrophysics Data System (ADS)
Kavalur, Aditya; Kim, Woo Kyun
Graphene and other 2-D materials have provided a promising prospect to improve the tribological properties of small length scale devices such as MEMS/NEMS due to their low friction coefficient and excellent wear resistance. Several recent research efforts have been devoted to unveiling the physical origin of the superior tribological properties of these 2-D materials from both experimental and theoretical standpoints, however, many of them still remain far from clearly understood. Recently, it was shown that lamellar materials do not conform to the predictions of the Prandtl-Tomlinson model due to additional friction mechanisms of delamination and visco-elastic ploughing. These mechanisms are critical as they explain the low and negative coefficients of friction observed in recent AFM experiments. However, thus far, most simulation and theoretical studies about these novel friction mechanisms have focused on only pristine graphene whereas real graphene sheets prepared by CVD and other conventional techniques possess various forms of defects such as vacancies and non-hexagonal rings. In this study we examine the role of these defects in frictional properties of 2-D materials in relation to delamination and visco-elastic ploughing.
Intersections, ideals, and inversion
Vasco, D.W.
1998-10-01
Techniques from computational algebra provide a framework for treating large classes of inverse problems. In particular, the discretization of many types of integral equations and of partial differential equations with undetermined coefficients lead to systems of polynomial equations. The structure of the solution set of such equations may be examined using algebraic techniques.. For example, the existence and dimensionality of the solution set may be determined. Furthermore, it is possible to bound the total number of solutions. The approach is illustrated by a numerical application to the inverse problem associated with the Helmholtz equation. The algebraic methods are used in the inversion of a set of transverse electric (TE) mode magnetotelluric data from Antarctica. The existence of solutions is demonstrated and the number of solutions is found to be finite, bounded from above at 50. The best fitting structure is dominantly onedimensional with a low crustal resistivity of about 2 ohm-m. Such a low value is compatible with studies suggesting lower surface wave velocities than found in typical stable cratons.
Parallel stitching of 2D materials
Ling, Xi; Wu, Lijun; Lin, Yuxuan; Ma, Qiong; Wang, Ziqiang; Song, Yi; Yu, Lili; Huang, Shengxi; Fang, Wenjing; Zhang, Xu; et al
2016-01-27
Diverse parallel stitched 2D heterostructures, including metal–semiconductor, semiconductor–semiconductor, and insulator–semiconductor, are synthesized directly through selective “sowing” of aromatic molecules as the seeds in the chemical vapor deposition (CVD) method. Lastly, the methodology enables the large-scale fabrication of lateral heterostructures, which offers tremendous potential for its application in integrated circuits.
Parallel Stitching of 2D Materials.
Ling, Xi; Lin, Yuxuan; Ma, Qiong; Wang, Ziqiang; Song, Yi; Yu, Lili; Huang, Shengxi; Fang, Wenjing; Zhang, Xu; Hsu, Allen L; Bie, Yaqing; Lee, Yi-Hsien; Zhu, Yimei; Wu, Lijun; Li, Ju; Jarillo-Herrero, Pablo; Dresselhaus, Mildred; Palacios, Tomás; Kong, Jing
2016-03-01
Diverse parallel stitched 2D heterostructures, including metal-semiconductor, semiconductor-semiconductor, and insulator-semiconductor, are synthesized directly through selective "sowing" of aromatic molecules as the seeds in the chemical vapor deposition (CVD) method. The methodology enables the large-scale fabrication of lateral heterostructures, which offers tremendous potential for its application in integrated circuits. PMID:26813882
Baby universes in 2d quantum gravity
NASA Astrophysics Data System (ADS)
Ambjørn, Jan; Jain, Sanjay; Thorleifsson, Gudmar
1993-06-01
We investigate the fractal structure of 2d quantum gravity, both for pure gravity and for gravity coupled to multiple gaussian fields and for gravity coupled to Ising spins. The roughness of the surfaces is described in terms of baby universes and using numerical simulations we measure their distribution which is related to the string susceptibility exponent γstring.
TOPAZ2D heat transfer code users manual and thermal property data base
NASA Astrophysics Data System (ADS)
Shapiro, A. B.; Edwards, A. L.
1990-05-01
TOPAZ2D is a two dimensional implicit finite element computer code for heat transfer analysis. This user's manual provides information on the structure of a TOPAZ2D input file. Also included is a material thermal property data base. This manual is supplemented with The TOPAZ2D Theoretical Manual and the TOPAZ2D Verification Manual. TOPAZ2D has been implemented on the CRAY, SUN, and VAX computers. TOPAZ2D can be used to solve for the steady state or transient temperature field on two dimensional planar or axisymmetric geometries. Material properties may be temperature dependent and either isotropic or orthotropic. A variety of time and temperature dependent boundary conditions can be specified including temperature, flux, convection, and radiation. Time or temperature dependent internal heat generation can be defined locally be element or globally by material. TOPAZ2D can solve problems of diffuse and specular band radiation in an enclosure coupled with conduction in material surrounding the enclosure. Additional features include thermally controlled reactive chemical mixtures, thermal contact resistance across an interface, bulk fluid flow, phase change, and energy balances. Thermal stresses can be calculated using the solid mechanics code NIKE2D which reads the temperature state data calculated by TOPAZ2D. A three dimensional version of the code, TOPAZ3D is available.
Screening and transport in 2D semiconductor systems at low temperatures
Das Sarma, S.; Hwang, E. H.
2015-01-01
Low temperature carrier transport properties in 2D semiconductor systems can be theoretically well-understood within RPA-Boltzmann theory as being limited by scattering from screened Coulomb disorder arising from random quenched charged impurities in the environment. In this work, we derive a number of analytical formula, supported by realistic numerical calculations, for the relevant density, mobility, and temperature range where 2D transport should manifest strong intrinsic (i.e., arising purely from electronic effects) metallic temperature dependence in different semiconductor materials arising entirely from the 2D screening properties, thus providing an explanation for why the strong temperature dependence of the 2D resistivity can only be observed in high-quality and low-disorder 2D samples and also why some high-quality 2D materials manifest much weaker metallicity than other materials. We also discuss effects of interaction and disorder on the 2D screening properties in this context as well as compare 2D and 3D screening functions to comment why such a strong intrinsic temperature dependence arising from screening cannot occur in 3D metallic carrier transport. Experimentally verifiable predictions are made about the quantitative magnitude of the maximum possible low-temperature metallicity in 2D systems and the scaling behavior of the temperature scale controlling the quantum to classical crossover. PMID:26572738
Microwave Imaging with Infrared 2-D Lock-in Amplifier
NASA Astrophysics Data System (ADS)
Chiyo, Noritaka; Arai, Mizuki; Tanaka, Yasuhiro; Nishikata, Atsuhiro; Maeno, Takashi
We have developed a 3-D electromagnetic field measurement system using 2-D lock-in amplifier. This system uses an amplitude modulated electromagnetic wave source to heat a resistive screen. A very small change of temperature on a screen illuminated with the modulated electromagnetic wave is measured using an infrared thermograph camera. In this paper, we attempted to apply our system to microwave imaging. By placing conductor patches in front of the resistive screen and illuminating with microwave, the shape of each conductor was clearly observed as the temperature difference image of the screen. In this way, the conductor pattern inside the non-contact type IC card could be visualized. Moreover, we could observe the temperature difference image reflecting the shape of a Konnyaku (a gelatinous food made from devil's-tonge starch) or a dried fishbone, both as non-conducting material resembling human body. These results proved that our method is applicable to microwave see-through imaging.
NASA Astrophysics Data System (ADS)
Nath, Saurabh; Mukherjee, Anish; Chatterjee, Souvick; Ganguly, Ranjan; Sen, Swarnendu; Mukhopadhyay, Achintya; Boreyko, Jonathan
2014-11-01
We have observed that capillarity forces may cause floatation in a few non-intuitive configurations. These may be divided into 2 categories: i) floatation of heavier liquid droplets on lighter immiscible ones and ii) fully submerged floatation of lighter liquid droplets in a heavier immiscible medium. We call these counter-intuitive because of the inverse floatation configuration. For case (i) we have identified and studied in detail the several factors affecting the shape and maximum volume of the floating drop. We used water and vegetable oil combinations as test fluids and established the relation between Bond Number and maximum volume contained in a floating drop (in the order of μL). For case (ii), we injected vegetable oil drop-wise into a pool of water. The fully submerged configuration of the drop is not stable and a slight perturbation to the system causes the droplet to burst and float in partially submerged condition. Temporal variation of a characteristic length of the droplet is analyzed using MATLAB image processing. The constraint of small Bond Number establishes the assumption of lubrication regime in the thin gap. A brief theoretical formulation also shows the temporal variation of the gap thickness. Jadavpur University, Jagadis Bose Centre of Excellence, Virginia Tech.
Application of 2D Non-Graphene Materials and 2D Oxide Nanostructures for Biosensing Technology
Shavanova, Kateryna; Bakakina, Yulia; Burkova, Inna; Shtepliuk, Ivan; Viter, Roman; Ubelis, Arnolds; Beni, Valerio; Starodub, Nickolaj; Yakimova, Rositsa; Khranovskyy, Volodymyr
2016-01-01
The discovery of graphene and its unique properties has inspired researchers to try to invent other two-dimensional (2D) materials. After considerable research effort, a distinct “beyond graphene” domain has been established, comprising the library of non-graphene 2D materials. It is significant that some 2D non-graphene materials possess solid advantages over their predecessor, such as having a direct band gap, and therefore are highly promising for a number of applications. These applications are not limited to nano- and opto-electronics, but have a strong potential in biosensing technologies, as one example. However, since most of the 2D non-graphene materials have been newly discovered, most of the research efforts are concentrated on material synthesis and the investigation of the properties of the material. Applications of 2D non-graphene materials are still at the embryonic stage, and the integration of 2D non-graphene materials into devices is scarcely reported. However, in recent years, numerous reports have blossomed about 2D material-based biosensors, evidencing the growing potential of 2D non-graphene materials for biosensing applications. This review highlights the recent progress in research on the potential of using 2D non-graphene materials and similar oxide nanostructures for different types of biosensors (optical and electrochemical). A wide range of biological targets, such as glucose, dopamine, cortisol, DNA, IgG, bisphenol, ascorbic acid, cytochrome and estradiol, has been reported to be successfully detected by biosensors with transducers made of 2D non-graphene materials. PMID:26861346