Geology topics

J. Ivanov

Publications and source records attributed to J. Ivanov.

11 recordsLinked to original sources

USGS2011

Multi-channel analysis of surface waves MASW of models with high shear-wave velocity contrast

We use the multi-channel analysis of surface waves MASW method to analyze synthetic seismic data calculated using models with high shear-wave velocity Vs contrast. The MASW dispersion-curve images of the Rayleigh wave are obtained using various sets of source-offset and spread-size configurations from the synthetic seismic data and compared with the theoretically calculated fundamental- and higher-mode dispersion-curves. Such tests showed that most of the dispersion-curve images are dominated by higher-mode energy at the low frequencies, especially when analyzing data from long receiver offsets and thus significantly divert from numerically expected dispersion-curve trends, which can lead to significant Vs overestimation. Further analysis showed that using data with relatively short spread lengths and source offsets can image the desired fundamental-mode of the Rayleigh wave that matches the numerically expected dispersion-curve pattern. As a result, it was concluded that it might be possible to avoid higher-mode contamination at low frequencies at sites with high Vs contrast by appropriate selection of spread size and seismic source offset. ?? 2011 Society of Exploration Geophysicists.

SEG Technical Program Expanded Abstracts↗

USGS2011

Near-surface shear-wave velocity measurements in unlithified sediment

S-wave velocity can be directly correlated to material stiffness and lithology making it a valuable physical property that has found uses in construction, engineering, and environmental projects. This study compares different methods for measuring S-wave velocities, investigating and identifying the differences among the methods' results, and prioritizing the different methods for optimal S-wave use at the U. S. Army's Yuma Proving Grounds YPG. Multichannel Analysis of Surface Waves MASW and S-wave tomography were used to generate S-wave velocity profiles. Each method has advantages and disadvantages. A strong signal-to-noise ratio at the study site gives the MASW method promising resolution. S-wave first arrivals are picked on impulsive sledgehammer data which were then used for the tomography process. Three-component downhole seismic data were collected in-line with a locking geophone, providing ground truth to compare the data and to draw conclusions about the validity of each data set. Results from these S-wave measurement techniques are compared with borehole seismic data and with lithology data from continuous samples to help ascertain the accuracy, and therefore applicability, of each method. This study helps to select the best methods for obtaining S-wave velocities for media much like those found in unconsolidated sediments at YPG. ?? 2011 Society of Exploration Geophysicists.

SEG Technical Program Expanded Abstracts↗

USGS2009

High-frequency Rayleigh-wave method

High-frequency (≥2 Hz) Rayleigh-wave data acquired with a multichannel recording system have been utilized to determine shear (S)-wave velocities in near-surface geophysics since the early 1980s. This overview article discusses the main research results of high-frequency surface-wave techniques achieved by research groups at the Kansas Geological Survey and China University of Geosciences in the last 15 years. The multichannel analysis of surface wave (MASW) method is a non-invasive acoustic approach to estimate near-surface S-wave velocity. The differences between MASW results and direct borehole measurements are approximately 15% or less and random. Studies show that simultaneous inversion with higher modes and the fundamental mode can increase model resolution and an investigation depth. The other important seismic property, quality factor ( Q ), can also be estimated with the MASW method by inverting attenuation coefficients of Rayleigh waves. An inverted model (S-wave velocity or Q ) obtained using a damped least-squares method can be assessed by an optimal damping vector in a vicinity of the inverted model determined by an objective function, which is the trace of a weighted sum of model-resolution and model-covariance matrices. Current developments include modeling high-frequency Rayleigh-waves in near-surface media, which builds a foundation for shallow seismic or Rayleigh-wave inversion in the time-offset domain; imaging dispersive energy with high resolution in the frequency-velocity domain and possibly with data in an arbitrary acquisition geometry, which opens a door for 3D surface-wave techniques; and successfully separating surface-wave modes, which provides a valuable tool to perform S-wave velocity profiling with high-horizontal resolution.

Earth Science↗

USGS2008

Refraction tomography mapping of near-surface dipping layers using landstreamer data at East Canyon Dam, Utah

We apply the P-wave refraction-tomography method to seismic data collected with a landstreamer. Refraction-tomography inversion solutions were determined using regularization parameters that provided the most realistic near-surface solutions that best matched the dipping layer structure of nearby outcrops. A reasonably well matched solution was obtained using an unusual set of optimal regularization parameters. In comparison, the use of conventional regularization parameters did not provide as realistic results. Thus, we consider that even if there is only qualitative a-priori information about a site (i.e., visual) - in the case of the East Canyon Dam, Utah - it might be possible to minimize the refraction nonuniqueness by estimating the most appropriate regularization parameters.

Conference Paper↗

USGS2007

Sensitivity of high-frequency Rayleigh-wave data revisited

Rayleigh-wave phase velocity of a layered earth model is a function of frequency and four groups of earth properties: P-wave velocity, S-wave velocity (Vs), density, and thickness of layers. Analysis of the Jacobian matrix (or the difference method) provides a measure of dispersion curve sensitivity to earth properties. Vs is the dominant influence for the fundamental mode (Xia et al., 1999) and higher modes (Xia et al., 2003) of dispersion curves in a high frequency range (>2 Hz) followed by layer thickness. These characteristics are the foundation of determining S-wave velocities by inversion of Rayleigh-wave data. More applications of surface-wave techniques show an anomalous velocity layer such as a high-velocity layer (HVL) or a low-velocity layer (LVL) commonly exists in near-surface materials. Spatial location (depth) of an anomalous layer is usually the most important information that surface-wave techniques are asked to provide. Understanding and correctly defining the sensitivity of high-frequency Rayleigh-wave data due to depth of an anomalous velocity layer are crucial in applying surface-wave techniques to obtain a Vs profile and/or determine the depth of an anomalous layer. Because depth is not a direct earth property of a layered model, changes in depth will result in changes in other properties. Modeling results show that sensitivity at a given depth calculated by the difference method is dependent on the Vs difference (contrast) between an anomalous layer and surrounding layers. The larger the contrast is, the higher the sensitivity due to depth of the layer. Therefore, the Vs contrast is a dominant contributor to sensitivity of Rayleigh-wave data due to depth of an anomalous layer. Modeling results also suggest that the most sensitive depth for an HVL is at about the middle of the depth to the half-space, but for an LVL it is near the ground surface. ?? 2007 Society of Exploration Geophysicists.

Conference Paper↗

USGS2007

Applications of the JARS method to study levee sites in southern Texas and southern New Mexico

We apply the joint analysis of refractions with surface waves (JARS) method to several sites and compare its results to traditional refraction-tomography methods in efforts of finding a more realistic solution to the inverse refraction-traveltime problem. The JARS method uses a reference model, derived from surface-wave shear-wave velocity estimates, as a constraint. In all of the cases JARS estimates appear more realistic than those from the conventional refraction-tomography methods. As a result, we consider, the JARS algorithm as the preferred method for finding solutions to the inverse refraction-tomography problems. ?? 2007 Society of Exploration Geophysicists.

Conference Paper↗

USGS2007

Multichannel analysis of surface waves (MASW) - Active and passive methods

The conventional seismic approaches for near-surface investigation have usually been either high-resolution reflection or refraction surveys that deal with a depth range of a few tens to hundreds meters. Seismic signals from these surveys consist of wavelets with frequencies higher than 50 Hz. The multichannel analysis of surface waves (MASW) method deals with surface waves in the lower frequencies (e.g., 1-30 Hz) and uses a much shallower depth range of investigation (e.g., a few to a few tens of meters). ?? 2007 Society of Exploration Geophysicists.

Leading Edge (Tulsa, OK)↗

USGS2006

Time-lapse seismic study of levees in southern New Mexico

The primary objective of this work was to measure changes in compressional- (Vp) and shear-wave (Vs) velocities in an earthen levee during a ponding experiment designed to simulate flood conditions on the Rio Grande in south New Mexico. Although similar to such experiment, performed an year earlier on the Rio Grande in south Texas, the levee seismic response results are different. This work was similar to previous Preliminary testing at three levee sites, all within a 1 km radius and each with unique physical, EM, and core characteristics, was completed and a single low-conductivity, highly fractured site was selected for investigation. Several different types of seismic data were recorded. Seismic data analysis techniques appraised included P-refraction tomography and Rayleigh surface-wave analysis using multichannel analysis of surface waves (MASW). P-wave velocity change (decrease) was rapid and isolated to one section within the pool confines, which already had anomalously high velocity most likely related to burrowing animals modification of the levee structure. S-wave velocity change was gradual and could be observed along the whole width of the pond within and below the levee. The results within the levee sand core were consistent with the observations of sand S-wave velocity changed due to saturation. ?? 2005 Society of Exploration Geophysicists.

SEG Technical Program Expanded Abstracts↗

USGS2005

The inverse problem of refraction travel times, part I: Types of Geophysical Nonuniqueness through Minimization

In a set of two papers we study the inverse problem of refraction travel times. The purpose of this work is to use the study as a basis for development of more sophisticated methods for finding more reliable solutions to the inverse problem of refraction travel times, which is known to be nonunique. The first paper, "Types of Geophysical Nonuniqueness through Minimization," emphasizes the existence of different forms of nonuniqueness in the realm of inverse geophysical problems. Each type of nonuniqueness requires a different type and amount of a priori information to acquire a reliable solution. Based on such coupling, a nonuniqueness classification is designed. Therefore, since most inverse geophysical problems are nonunique, each inverse problem must be studied to define what type of nonuniqueness it belongs to and thus determine what type of a priori information is necessary to find a realistic solution. The second paper, "Quantifying Refraction Nonuniqueness Using a Three-layer Model," serves as an example of such an approach. However, its main purpose is to provide a better understanding of the inverse refraction problem by studying the type of nonuniqueness it possesses. An approach for obtaining a realistic solution to the inverse refraction problem is planned to be offered in a third paper that is in preparation. The main goal of this paper is to redefine the existing generalized notion of nonuniqueness and a priori information by offering a classified, discriminate structure. Nonuniqueness is often encountered when trying to solve inverse problems. However, possible nonuniqueness diversity is typically neglected and nonuniqueness is regarded as a whole, as an unpleasant "black box" and is approached in the same manner by applying smoothing constraints, damping constraints with respect to the solution increment and, rarely, damping constraints with respect to some sparse reference information about the true parameters. In practice, when solving geophysical problems different types of nonuniqueness exist, and thus there are different ways to solve the problems. Nonuniqueness is usually regarded as due to data error, assuming the true geology is acceptably approximated by simple mathematical models. Compounding the nonlinear problems, geophysical applications routinely exhibit exact-data nonuniqueness even for models with very few parameters adding to the nonuniqueness due to data error. While nonuniqueness variations have been defined earlier, they have not been linked to specific use of a priori information necessary to resolve each case. Four types of nonuniqueness, typical for minimization problems are defined with the corresponding methods for inclusion of a priori information to find a realistic solution without resorting to a non-discriminative approach. The above-developed stand-alone classification is expected to be helpful when solving any geophysical inverse problems. ?? Birkha??user Verlag, Basel, 2005.

Pure and Applied Geophysics↗

USGS2005

The inverse problem of refraction travel times, part II: Quantifying refraction nonuniqueness using a three-layer model

This paper is the second of a set of two papers in which we study the inverse refraction problem. The first paper, "Types of Geophysical Nonuniqueness through Minimization," studies and classifies the types of nonuniqueness that exist when solving inverse problems depending on the participation of a priori information required to obtain reliable solutions of inverse geophysical problems. In view of the classification developed, in this paper we study the type of nonuniqueness associated with the inverse refraction problem. An approach for obtaining a realistic solution to the inverse refraction problem is offered in a third paper that is in preparation. The nonuniqueness of the inverse refraction problem is examined by using a simple three-layer model. Like many other inverse geophysical problems, the inverse refraction problem does not have a unique solution. Conventionally, nonuniqueness is considered to be a result of insufficient data and/or error in the data, for any fixed number of model parameters. This study illustrates that even for overdetermined and error free data, nonlinear inverse refraction problems exhibit exact-data nonuniqueness, which further complicates the problem of nonuniqueness. By evaluating the nonuniqueness of the inverse refraction problem, this paper targets the improvement of refraction inversion algorithms, and as a result, the achievement of more realistic solutions. The nonuniqueness of the inverse refraction problem is examined initially by using a simple three-layer model. The observations and conclusions of the three-layer model nonuniqueness study are used to evaluate the nonuniqueness of more complicated n-layer models and multi-parameter cell models such as in refraction tomography. For any fixed number of model parameters, the inverse refraction problem exhibits continuous ranges of exact-data nonuniqueness. Such an unfavorable type of nonuniqueness can be uniquely solved only by providing abundant a priori information. Insufficient a priori information during the inversion is the reason why refraction methods often may not produce desired results or even fail. This work also demonstrates that the application of the smoothing constraints, typical when solving ill-posed inverse problems, has a dual and contradictory role when applied to the ill-posed inverse problem of refraction travel times. This observation indicates that smoothing constraints may play such a two-fold role when applied to other inverse problems. Other factors that contribute to inverse-refraction-problem nonuniqueness are also considered, including indeterminacy, statistical data-error distribution, numerical error and instability, finite data, and model parameters. ?? Birkha??user Verlag, Basel, 2005.

Pure and Applied Geophysics↗

USGS2002

Comparing shear-wave velocity profiles inverted from multichannel surface wave with borehole measurements

Recent field tests illustrate the accuracy and consistency of calculating near-surface shear (S)-wave velocities using multichannel analysis of surface waves (MASW). S-wave velocity profiles (S-wave velocity vs. depth) derived from MASW compared favorably to direct borehole measurements at sites in Kansas, British Columbia, and Wyoming. Effects of changing the total number of recording channels, sampling interval, source offset, and receiver spacing on the inverted S-wave velocity were studied at a test site in Lawrence, Kansas. On the average, the difference between MASW calculated Vs and borehole measured Vs in eight wells along the Fraser River in Vancouver, Canada was less than 15%. One of the eight wells was a blind test well with the calculated overall difference between MASW and borehole measurements less than 9%. No systematic differences were observed in derived Vs values from any of the eight test sites. Surface wave analysis performed on surface data from Wyoming provided S-wave velocities in near-surface materials. Velocity profiles from MASW were confirmed by measurements based on suspension log analysis. ?? 2002 Elsevier Science Ltd. All rights reserved.

Soil Dynamics and Earthquake Engineering↗