skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Inversion of tracer test data using tomographic constraints

Authors:
; ;
Publication Date:
Research Org.:
Subsurface Biogeochemical Research (SBR)
Sponsoring Org.:
USDOE Office of Science (SC), Biological and Environmental Research (BER)
OSTI Identifier:
1154404
Resource Type:
Journal Article
Resource Relation:
Journal Name: Water Resources Research; Journal Volume: 42; Journal Issue: 4
Country of Publication:
United States
Language:
English

Citation Formats

Niklas,Linde, Stefan,Finsterle, and Susan,Hubbard. Inversion of tracer test data using tomographic constraints. United States: N. p., 2006. Web. doi:10.1029/2004WR003806.
Niklas,Linde, Stefan,Finsterle, & Susan,Hubbard. Inversion of tracer test data using tomographic constraints. United States. doi:10.1029/2004WR003806.
Niklas,Linde, Stefan,Finsterle, and Susan,Hubbard. Sat . "Inversion of tracer test data using tomographic constraints". United States. doi:10.1029/2004WR003806.
@article{osti_1154404,
title = {Inversion of tracer test data using tomographic constraints},
author = {Niklas,Linde and Stefan,Finsterle and Susan,Hubbard},
abstractNote = {},
doi = {10.1029/2004WR003806},
journal = {Water Resources Research},
number = 4,
volume = 42,
place = {United States},
year = {Sat Apr 01 00:00:00 EST 2006},
month = {Sat Apr 01 00:00:00 EST 2006}
}
  • A reasonable description of the hydraulic conductivity structure is a prerequisite for modeling contaminant transport. However, formulations of hydrogeological inverse problems utilizing hydrogeological data only often fail to reliably resolve features at a resolution required for accurately predicting transport. Incorporation of geophysical data into the inverse problem offers the potential to increase this resolution. In this study, we invert hydrological tracer test data using the shape and relative magnitude variations derived from geophysical tomographic data to regionalize a hydrogeological inverse problem in order to estimate the hydraulic conductivity structure. Our approach does not require that the petrophysical relationship be knownmore » a-priori, but that it is linear and stationary within each geophysical anomaly. However, tomograms are imperfect models of geophysical properties and geophysical properties are not necessarily strongly linked to hydraulic conductivity. Therefore, we focus on synthetic examples where the correlation between radar velocity and hydraulic conductivity, as well as the geophysical data acquisition errors, are varied in order to assess what aspects of the hydraulic conductivity structure we can expect to resolve under different conditions. The results indicate that regularization of the tracer inversion procedure using geophysical data improves estimates of hydraulic conductivity. We find that even under conditions of corrupted geophysical data, we can accurately estimate the effective hydraulic conductivity and areas of high and low hydraulic conductivity. However, given imperfect geophysical data, our results suggest that we cannot expect accurate estimates of the variability of the hydraulic conductivity structure.« less
  • Accurate tomographic inversion is important for diagnostic systems on stellarators and tokamaks which rely on measurements of line integrated emission spectra. A tomographic inversion technique based on spline optimization with enforcement of constraints is described that can produce unique and physically relevant inversions even in situations with noisy or incomplete input data. This inversion technique is routinely used in the analysis of data from the x-ray imaging crystal spectrometer (XICS) installed at the Large Helical Device. The XICS diagnostic records a 1D image of line integrated emission spectra from impurities in the plasma. Through the use of Doppler spectroscopy andmore » tomographic inversion, XICS can provide profile measurements of the local emissivity, temperature, and plasma flow. Tomographic inversion requires the assumption that these measured quantities are flux surface functions, and that a known plasma equilibrium reconstruction is available. In the case of low signal levels or partial spatial coverage of the plasma cross-section, standard inversion techniques utilizing matrix inversion and linear-regularization often cannot produce unique and physically relevant solutions. The addition of physical constraints, such as parameter ranges, derivative directions, and boundary conditions, allow for unique solutions to be reliably found. The constrained inversion technique described here utilizes a modified Levenberg-Marquardt optimization scheme, which introduces a condition avoidance mechanism by selective reduction of search directions. The constrained inversion technique also allows for the addition of more complicated parameter dependencies, for example, geometrical dependence of the emissivity due to asymmetries in the plasma density arising from fast rotation. The accuracy of this constrained inversion technique is discussed, with an emphasis on its applicability to systems with limited plasma coverage.« less
  • Accurate tomographic inversion is important for diagnostic systems on stellarators and tokamaks which rely on measurements of line integrated emission spectra. A tomographic inversion technique based on spline optimization with enforcement of constraints is described that can produce unique and physically relevant inversions even in situations with noisy or incomplete input data. This inversion technique is routinely used in the analysis of data from the x-ray imaging crystal spectrometer (XICS) installed at LHD. The XICS diagnostic records a 1D image of line integrated emission spectra from impurities in the plasma. Through the use of Doppler spectroscopy and tomographic inversion, XICSmore » can provide pro file measurements of the local emissivity, temperature and plasma flow. Tomographic inversion requires the assumption that these measured quantities are flux surface functions, and that a known plasma equilibrium reconstruction is available. In the case of low signal levels or partial spatial coverage of the plasma cross-section, standard inversion techniques utilizing matrix inversion and linear-regularization often cannot produce unique and physically relevant solutions. The addition of physical constraints, such as parameter ranges, derivative directions, and boundary conditions, allow for unique solutions to be reliably found. The constrained inversion technique described here utilizes a modifi ed Levenberg-Marquardt optimization scheme, which introduces a condition avoidance mechanism by selective reduction of search directions. The constrained inversion technique also allows for the addition of more complicated parameter dependencies, for example geometrical dependence of the emissivity due to asymmetries in the plasma density arising from fast rotation. The accuracy of this constrained inversion technique is discussed, with an emphasis on its applicability to systems with limited plasma coverage.« less
  • No abstract prepared.
  • The realization of local earthquake tomography is usually conducted by removing distant events outside the study region, because these events may increase errors. In this study, tomographic inversion has been conducted using the travel time data of local and regional events in order to improve the structural resolution, especially for deep structures. We used the local MERapi Amphibious EXperiments (MERAMEX) data catalog that consists of 292 events from May to October 2004. The additional new data of regional events in the Java region were taken from the Meteorological, Climatological, and Geophysical Agency (MCGA) of Indonesia, which consist of 882 events,more » having at least 10 recording phases at each seismographic station from April 2009 to February 2011. We have conducted joint inversions of the combined data sets using double-difference tomography to invert for velocity structures and to conduct hypocenter relocation simultaneously. The checkerboard test results of Vp and Vs structures demonstrate a significantly improved spatial resolution from the shallow crust down to a depth of 165 km. Our tomographic inversions reveal a low velocity anomaly beneath the Lawu - Merapi zone, which is consistent with the results from previous studies. A strong velocity anomaly zone with low Vp, low Vs and low Vp/Vs is also identified between Cilacap and Banyumas. We interpret this anomaly as a fluid content material with large aspect ratio or sediment layer. This anomaly zone is in a good agreement with the existence of a large dome containing sediment in this area as proposed by previous geological studies. A low velocity anomaly zone is also detected in Kebumen, where it may be related to the extensional oceanic basin toward the land.« less