Wavefield Data Research Articles

Impact induced damage assessment by means of Lamb wave image processing

Abstract The aim of this research is an analysis of full wavefield Lamb wave interaction with impact-induced damage at various impact energies in order to find out the limitation of the wavenumber adaptive image filtering method. In other words, the relation between impact energy and damage detectability will be shown. A numerical model based on the time domain spectral element method is used for modeling of Lamb wave propagation and interaction with barely visible impact damage in a carbon-epoxy laminate. Numerical studies are followed by experimental research on the same material with an impact damage induced by various energy and also a Teflon insert simulating delamination. Wavenumber adaptive image filtering and signal processing are used for damage visualization and assessment for both numerical and experimental full wavefield data. It is shown that it is possible to visualize and assess the impact damage location, size and to some extent severity by using the proposed technique.

Mathematical design of a novel input/instruction device using a moving acoustic emitter

This paper is concerned with the mathematical design of a novel input/instruction device using a moving emitter. The emitter acts as a point source and can be installed on a digital pen or worn on the finger of the human being who desires to interact/communicate with the computer. The input/instruction can be recognized by identifying the moving trajectory of the emitter performed by the human being from the collected wave field data. The identification process is modelled as an inverse source problem where one intends to identify the trajectory of a moving point source. There are several salient features of our study which distinguish our result from the existing ones in the literature. First, the point source is moving in an inhomogeneous background medium, which models the human body. Second, the dynamical wave field data are collected in a limited aperture. Third, the reconstruction method is independent of the background medium, and it is totally direct without any matrix inversion. Hence, it is efficient and robust with respect to the measurement noise. Both theoretical justifications and computational experiments are presented to verify our novel findings.

Detecting buried wave-penetrable scatterers in a two-layered medium

In a two-layered medium, we prove that a buried inhomogeneous scatterer is uniquely determined from the wave field data measured in the upper half-space with respect to many incident point sources. Moreover, we extend the multilevel sampling method in Liu and Zou (2013) to numerically reconstruct the buried scatterer applying only few incident fields and partial scattered data. The extended recovery scheme only involves matrix–vector operations and does not need to solve any large-scale ill-posed linear systems or any optimization process. It is feasible to deal with the scatterers of different features and easy to implement, highly tolerant to noise and computationally quite cheap. We can regard it as an effective yet simple computational method to provide a reliable initial guess for the implementation in existing more accurate and refined optimization-type reconstruction algorithms.

Statistical partial wavefield imaging using Lamb wave signals

This article presents a baseline-free, model-driven, statistical damage detection and imaging framework for guided waves measured from partial (i.e. non-dense) wavefield scans. Wavefield analysis is an effective non-contact technique for non-destructive evaluation. Yet, there are several limitations to practically implement wavefield methods. These limitations include slow data acquisition and a lack of statistical reliability. Our approach addresses both of these challenges. We use sparse wavenumber analysis, sparse wavenumber synthesis, and data-fitting optimization to accurately model damage-free wavefield data. We then combine this model with matched field processing to image damage from a small number of partial wavefield measurements. We further derive a hypothesis test based on extreme value theory to statistically detect damage. We test our framework with Lamb wave measurements from a steel plate. With 70 experimental wavefield measurements, we achieve an empirical probability of damage detection of more than 98%, an empirical probability of false alarm of less than 0.17%, and an accurate image of the damage.

Wave Energy Resource Assessment off the Coast of China around the Zhoushan Islands

Based on ERA-Interim reanalysis wave field data for the 36 years from 1979 to 2014, the temporal and spatial distributions and development potential of wave energy are studied in detail in the offshore and relatively nearshore waters adjacent to the Zhoushan Islands. The results show that areas of relatively high wave energy are located in the offshore waters to the east and southeast of the Zhoushan Islands. The potential wave energy in relatively nearshore waters (water depths from 10 m to 65 m) is relatively higher and especially in the nearshore waters to the southeast of the city of Taizhou, Zhejiang Province, which are suitable locations for wave energy development. The conclusions provide scientific guidance for wave energy development in the sea areas adjacent to the Zhoushan Islands.

Optimizing Acquisition Geometry in Shallow Gas Cloud Using Particle Swarm Optimization Approach

Many hydrocarbon explorations in mature fields have been severely affected by complex and overburdening issues, such as shallow gas accumulation, gas pockets, and gas seepage. In this work, a new forward modelling technique is proposed in evaluating the potential survey design for fields affected by shallow gas cloud. In recent years, the implementation of innovative acquisition layouts has been producing significantly better seismic images, especially in the low illumination subsurface area. However, the uncertainty of the effectiveness in new acquisition design subsurface coverage always become a major stumbling block. To overcome this constraint, an optimization approach is suggested through the smart source and receiver location arrangement on the surface, with significant alignment to the conventional source and receiver arrangement approach. The particle swarm optimization (PSO) method is used to find the source-receiver configuration with maximum subsurface illumination coverage for the gas affected field situated in Malaysia Basin. Implementation of the PSO algorithm requires both a velocity model building process and wave field extrapolation from a target reflector to the surface level. The wave field data then was used to simulate receiver optimization outputs which eventually determined the subsurface illumination coverage. The results from the new optimization method for both synthetic model and Malaysia Basin data, offer a greater understanding of the consequences of obstacles caused by shallow anomalies with respect to seismic acquisition, data processing, and interpretation.

The derivative-free Fourier shell identity for photoacoustics.

In X-ray tomography, the Fourier slice theorem provides a relationship between the Fourier components of the object being imaged and the measured projection data. The Fourier slice theorem is the basis for X-ray Fourier-based tomographic inversion techniques. A similar relationship, referred to as the ‘Fourier shell identity’ has been previously derived for photoacoustic applications. However, this identity relates the pressure wavefield data function and its normal derivative measured on an arbitrary enclosing aperture to the three-dimensional Fourier transform of the enclosed object evaluated on a sphere. Since the normal derivative of pressure is not normally measured, the applicability of the formulation is limited in this form. In this paper, alternative derivations of the Fourier shell identity in 1D, 2D polar and 3D spherical polar coordinates are presented. The presented formulations do not require the normal derivative of pressure, thereby lending the formulas directly adaptable for Fourier based absorber reconstructions.

Multifrequency magnetic resonance elastography of the brain reveals tissue degeneration in neuromyelitis optica spectrum disorder.

Application of multifrequency magnetic resonance elastography (MMRE) of the brain parenchyma in patients with neuromyelitis optica spectrum disorder (NMOSD) compared to age matched healthy controls (HC). 15 NMOSD patients and 17 age- and gender-matched HC were examined using MMRE. Two three-dimensional viscoelastic parameter maps, the magnitude |G*| and phase angle φ of the complex shear modulus were reconstructed by simultaneous inversion of full wave-field data in 1.9-mm isotropic resolution at 7 harmonic drive frequencies from 30 to 60Hz. In NMOSD patients, a significant reduction of |G*| was observed within the white matter fraction (p = 0.017), predominantly within the thalamic regions (p = 0.003), compared to HC. These parameters exceeded the reduction in brain volume measured in patients versus HC (p = 0.02 whole-brain volume reduction). Volumetric differences in white matter fraction and the thalami were not detectable between patients and HC. However, phase angle φ was decreased in patients within the white matter (p = 0.03) and both thalamic regions (p = 0.044). MMRE reveals global tissue degeneration with accelerated softening of the brain parenchyma in patients with NMOSD. The predominant reduction of stiffness is found within the thalamic region and related white matter tracts, presumably reflecting Wallerian degeneration. • Magnetic resonance elastography reveals diffuse cerebral tissue changes in patients with NMOSD. • Premature tissue softening in NMOSD patients indicates tissue degeneration. • Hypothesis of a widespread cerebral neurodegeneration in form of diffuse tissue alteration.

Robust contrast source inversion method with automatic choice rule of regularization parameters for ultrasound waveform tomography

We consider ultrasound waveform tomography using an ultrasound prototype equipped with the ring-array transducers. For this purpose, we use robust contrast source inversion (robust CSI), viz extended contrast source inversion, to reconstruct the sound-speed image from the wave-field data. The robust CSI method is implemented by the alternating minimization method. An automatic choice rule is employed into the alternating minimization method in order to heuristically determine a suitable regularization parameter while iterating. We prove the convergence of this algorithm. The numerical examples show that the robust CSI method with the automatic choice rule improves the spatial resolution of medical images and enhances the robustness, even when the wave-field data of a wavelength of 6.16 mm contaminated by 5% noise are used. The numerical results also show that the images reconstructed by the proposed method yield a spatial resolution of approximately half the wavelength that may be adequate for imaging a breast tumor at Stage I.

Complete wavefield seismic structural, stratigraphic, geomorphic and elastic interpretation in the offshore Browse Basin

The first quantitative interpretation project in the Browse Basin that used complete wavefields from dual-sensor streamer data to characterise a variety of overburden and reservoir seismic structural, sedimentary and geomorphic features is discussed in this extended abstract. Having established that complete wavefield pre-stack angle data are amplitude versus angle (AVA) compliant, the authors analyse the elastic response from thin gas sands in a location where well control predicts a weak gas response in terms of the Vp/Vs ratio and negligible response in terms of the acoustic impedance. Various seismic expressions in the near-surface are better resolved using complete wavefield imaging, including deltaic deposits, potentially prospective patch reefs and intersecting channel systems, and shallow hazards associated with historic drilling failures appear to be resolvable with greater confidence. On the larger scale, the overall data interpretability is substantially enhanced in texture and character by broadband imaging of primary reflections, and then by subsequent reservoir characterisation. Initial (deeper) interpretation focused around mapping of the prospective Early Jurassic Plover Formation and the Early Cretaceous Upper Vulcan Formation, in which evidence of channelisation can be demonstrated on conventional migration volumes and reinforced using complete wavefield migration. Ultra-low frequency signal penetration to the deeper Triassic level has been limited historically, but the first complete seismic wavefield data analysis presented demonstrates that a combination of volume attributes provide substantially improved geologic information from the seafloor to the deepest targets.

Isolation of ultrasonic scattering by wavefield baseline subtraction

Wavefield imaging generally refers to the measurement of signals over a two-dimensional rectilinear grid that originate from a spatially fixed source. Subtraction of such wavefields is investigated as a means of separating scattered signals from the total wavefield; that is, baseline wavefield data acquired from a defect-free specimen are subtracted from analogous data acquired after introduction of a defect. The wavefields considered here are generated by a 5MHz angle-beam probe and measured over an area of the accessible specimen surface using a laser vibrometer. The primary challenge in isolating the scattered waves is imperfect temporal and spatial alignment of the two wavefields. Two methods for aligning the wavefields in space and time prior to performing baseline subtraction are presented and their efficacy is evaluated using data acquired before and after introducing notches that originate from a through-hole. Although perfect baseline subtraction is not achieved, the improvement in performance after alignment using either method allows for scattered waves from small defects to be separated and visualized, even when their amplitudes are much smaller than those of the incident waves.

Exploitable wave energy assessment based on ERA-Interim reanalysis data—A case study in the East China Sea and the South China Sea

Wave energy resources assessment is a very important process before the exploitation and utilization of the wave energy. At present, the existing wave energy assessment is focused on theoretical wave energy conditions for interesting areas. While the evaluation for exploitable wave energy conditions is scarcely ever performed. Generally speaking, the wave energy are non-exploitable under a high sea state and a lower sea state which must be ignored when assessing wave energy. Aiming at this situation, a case study of the East China Sea and the South China Sea is performed. First, a division basis between the theoretical wave energy and the exploitable wave energy is studied. Next, based on recent 20 a ERA-Interim wave field data, some indexes including the spatial and temporal distribution of wave power density, a wave energy exploitable ratio, a wave energy level, a wave energy stability, a total wave energy density, the seasonal variation of the total wave energy and a high sea condition frequency are calculated. And then the theoretical wave energy and the exploitable wave energy are compared each other; the distributions of the exploitable wave energy are assessed and a regional division for exploitable wave energy resources is carried out; the influence of the high sea state is evaluated. The results show that considering collapsing force of the high sea state and the utilization efficiency for wave energy, it is determined that the energy by wave with a significant wave height being not less 1 m or not greater than 4 m is the exploitable wave energy. Compared with the theoretical wave energy, the average wave power density, energy level, total wave energy density and total wave energy of the exploitable wave energy decrease obviously and the stability enhances somewhat. Pronounced differences between the theoretical wave energy and the exploitable wave energy are present. In the East China Sea and the South China Sea, the areas of an abundant and stable exploitable wave energy are primarily located in the north-central part of the South China Sea, the Luzon Strait, east of Taiwan, China and north of Ryukyu Islands; annual average exploitable wave power density values in these areas are approximately 10–15 kW/m; the exploitable coefficient of variation (COV) and seasonal variation (SV) values in these areas are less than 1.2 and 1, respectively. Some coastal areas of the Beibu Gulf, the Changjiang Estuary, the Hangzhou Bay and the Zhujiang Estuary are the poor areas of the wave energy. The areas of the high wave energy exploitable ratio is primarily in nearshore waters. The influence of the high sea state for the wave energy in nearshore waters is less than that in offshore waters. In the areas of the abundant wave energy, the influence of the high sea state for the wave energy is prominent and the utilization of wave energy is relatively difficult. The developed evaluation method may give some references for an exploitable wave energy assessment and is valuable for practical applications.

Compressive sensing of full wave field data for structural health monitoring applications.

Numerous nondestructive evaluations and structural health monitoring approaches based on guide waves rely on analysis of wave fields recorded through scanning laser Doppler vibrometers (SLDVs) or ultrasonic scanners. The informative content which can be extracted from these inspections is relevant; however, the acquisition process is generally time-consuming, posing a limit in the applicability of such approaches. To reduce the acquisition time, we use a random sampling scheme based on compressive sensing (CS) to minimize the number of points at which the field is measured. The CS reconstruction performance is mostly influenced by the choice of a proper decomposition basis to exploit the sparsity of the acquired signal. Here, different bases have been tested to recover the guided waves wave field acquired on both an aluminum and a composite plate. Experimental results show that the proposed approach allows a reduction of the measurement locations required for accurate signal recovery to less than 34% of the original sampling grid.

Crack imaging and quantification in aluminum plates with guided wave wavenumber analysis methods

Guided wavefield analysis methods for detection and quantification of crack damage in an aluminum plate are presented in this paper. New wavenumber components created by abrupt wave changes at the structural discontinuity are identified in the frequency–wavenumber spectra. It is shown that the new wavenumbers can be used to detect and characterize the crack dimensions. Two imaging based approaches, filter reconstructed imaging and spatial wavenumber imaging, are used to demonstrate how the cracks can be evaluated with wavenumber analysis. The filter reconstructed imaging is shown to be a rapid method to map the plate and any existing damage, but with less precision in estimating crack dimensions; while the spatial wavenumber imaging provides an intensity image of spatial wavenumber values with enhanced resolution of crack dimensions. These techniques are applied to simulated wavefield data, and the simulation based studies show that spatial wavenumber imaging method is able to distinguish cracks of different severities. Laboratory experimental validation is performed for a single crack case to confirm the methods’ capabilities for imaging cracks in plates.

Full waveform inversion in the frequency domain using direct iterative T-matrix methods

We present two direct iterative solutions to the nonlinear seismic waveform inversion problem that are based on volume integral equation methods for seismic forward modelling in the acoustic approximation. The solutions are presented in the frequency domain, where accurate inversion results can often be obtained using a relatively low number of frequency components. Our inverse scattering approach effectively replaces an ill-posed nonlinear inverse problem with a series of linear ill-posed inverse problems, for which there already exist efficient (regularized) solution methods. Both these solutions update the wavefield within the scattering domain after each iteration. The main difference is that the background medium Green functions are kept fixed in the first solution, but updated after each iteration in the second solution. This means that our solutions are very similar to the Born iterative (BI) and the distorted Born iterative (DBI) methods that are commonly used in acoustic and electromagnetic inverse scattering. However, we have eliminated the need to perform a full forward simulation (or to invert a huge matrix) at each iteration via the use of an iterative T-matrix method for fixed background media for the BI method and a variational T-matrix method for dynamic background media for the DBI method. The T-matrix (variation) is linearly related with the seismic wavefield data (residuals), but related with the unknown scattering potential model parameter (updates) in a non-linear manner, which is independent of the source-receiver configuration. This mathematical structure, which allows one to peel off the effects of the source-receiver configuration, is very attractive when dealing with multiple (simultaneous) sources, and is also compatible with the (future) use of renormalization methods for dealing with local minima problems. To illustrate the performance and potential of the two direct iterative methods for FWI, we performed a series of numerical experiments on synthetic seismic waveform data associated with a simple 2D model and the more complicated Marmousi model. The results of these numerical experiments suggest that the use of a fixed (e.g. smooth and ray-tracing friendly) background medium may be adequate for some applications with moderately large velocity contrasts, but the solution based on a dynamic (non-smooth and constantly updated) background medium will normally provide superiour inversion results; also in the case of low signal-to-noise ratios.

Efficient elastic reverse-time migration for the decomposed P-wavefield using stress tensor in the time domain

Because complex mixed waves are typically generated in elastic media, wavefield decomposition is required for such media to obtain migration images accurately. In isotropic media, this is achieved according to the Helmholtz decomposition theorem; in particular, the divergence operator is commonly applied to P-wavefield decomposition. In this study, two types of elastic reverse-time migration algorithms are proposed for decomposition of the P-wavefield without requiring the divergence operator. The first algorithm involves formulation of the stress tensor by spatially differentiated displacement according to the stress–strain relationship and is utilized to construct an imaging condition for the decomposed P-wavefield. We demonstrate this approach through numerical testing. The second algorithm allows us to obtain emphasized interfaces through the application of the absolute value function to decomposed wavefield in imaging condition.Because reverse-time migration can be defined by a zero-lag cross-correlation relationship between the partial-derivative wavefield and the observed wavefield data, we derive the virtual source to construct the partial-derivative wavefield based on a 2D staggered-grid finite-difference modeling method in the time domain. The explicitly computed partial-derivative wavefield from virtual sources with the stress tensor is in agreement with the partial-derivative wavefield directly computed from residual by between with and without a perturbation point in the subsurface. Moreover, the back-propagation technique is used to enhance the computational efficiency.To validate our two types of imaging conditions, numerical tests are conducted. The migration images created according to our imaging conditions can represent the subsurface structure accurately. Thus, we can confirm the feasibility of obtaining migration images of the decomposed P-wavefield without requiring the application of the divergence operator.

A wave energy resource assessment in the China’s seas based on multi-satellite merged radar altimeter data

Wave energy resources are abundant in both offshore and nearshore areas of the China’s seas. A reliable assessment of the wave energy resources must be performed before they can be exploited. First, for a water depth in offshore waters of China, a parameterized wave power density model that considers the effects of the water depth is introduced to improve the calculating accuracy of the wave power density. Second, wave heights and wind speeds on the surface of the China’s seas are retrieved from an AVISO multi-satellite altimeter data set for the period from 2009 to 2013. Three mean wave period inversion models are developed and used to calculate the wave energy period. Third, a practical application value for developing the wave energy is analyzed based on buoy data. Finally, the wave power density is then calculated using the wave field data. Using the distribution of wave power density, the energy level frequency, the time variability indexes, the total wave energy and the distribution of total wave energy density according to a wave state, the offshore wave energy in the China’s seas is assessed. The results show that the areas of abundant and stable wave energy are primarily located in the north-central part of the South China Sea, the Luzon Strait, southeast of Taiwan in the China’s seas; the wave power density values in these areas are approximately 14.0–18.5 kW/m. The wave energy in the China’s seas presents obvious seasonal variations and optimal seasons for a wave energy utilization are in winter and autumn. Except for very coastal waters, in other sea areas in the China’s seas, the energy is primarily from the wave state with 0.5 m ⩽ H s ⩽ 4 m, 4 s ⩽ T e ⩽ 10 s where H s is a significant wave height and T e is an energy period; within this wave state, the wave energy accounts for 80% above of the total wave energy. This characteristic is advantageous to designing wave energy convertors (WECs). The practical application value of the wave energy is higher which can be as an effective supplement for an energy consumption in some areas. The above results are consistent with the wave model which indicates fully that this new microwave remote sensing method altimeter is effective and feasible for the wave energy assessment.

항만공진주기와 선박동요량을 고려한 항만가동율 산정

This study proposed an estimation method of allowable wave height for loading and unloading of the ship and evaluation of effective working days considering moored ship motion that is affected by sip sizes, mooring conditions, wave periods and directions. The method was examined validity by comparison with wave field data at pier 8th in Pohang New Harbor. The wave field data obtained with wave height of 0.10~0.75 m and wave period of 7~13 s in ship sizes of 800~35,000 ton when a downtimes have occurred. On the other hand, the results of allowable wave height for loading and unloading of the ship in this method have obtained with wave heights of 0.19~0.50 m and wave periods of 8~12 s for ship sizes of 5,000, 10,000 and 30,000 ton. Thus this method well reproduced the field data respond to various a ship sizes and wave periods. And the results of this in Korea are didn’t respond to various the ship sizes and wave periods, and we h method tended to decrease in 16~62 percent when have considered long wave, and it is decreased in 0~46 percent when didn’t consider long wave than design standards in case of the ship sizes of 5,000~30,000 ton, wave period of 12 s and wave angle of 75o. The allowable wave heights for loading and unloading of the ship proposed by design standards in Korea have found that overestimated on smaller than 10,000 ton. On the other hand, the rate of effective working days considering ship motion at pier 8th in Pohang New Harbor reduced in 6.5 percent when compare with the results without considering ship motion. 핵심용어: Allowable wave height for loading and unloading of ship, Oscillation quantity of moored ship, ship model for analysis of moored ship motion, effective working days in a harbor, Pohang New harbor, í•˜ì—­í•œê³„íŒŒê³ , ì„ ë°•ë™ìš”ëŸ‰, ì„ ë°•ë™ìš”í•´ì„ëª¨í˜•, 항만가동율, í¬í•­ì‹ í•­

Case study of guided wave propagation in a one-side water-immersed steel plate

This paper presents a case study of guided waves in a steel plate with one side immersed in water. A hybrid sensing system that uses PZT as the guided wave actuator and a scanning laser vibrometer as the guided wave receiver is employed to acquire the time–space wavefield data. By using the two dimensional Fourier transform, the time–space wavefield is transformed into the frequency–wavenumber domain where the wave modes and the guided wave dispersion can be determined. The study confirms the existence of quasi-Scholte wave in the one-sidedly water-immersed plate in addition to fundamental guided waves. The results also show the quasi-Scholte wave can be directly generated and measured in the immersed plate at low frequencies using the present sensing system. Through pitch–catch sensing tests, the influence of water on guided wave propagation in the one-side water-immersed plate is investigated. It is seen that the water level affects the wave propagation time linearly and can be potentially used for estimation of water level in a container.

Instantaneous and local wavenumber estimations for damage quantification in composites

The continued and expanded use of composite materials in aerospace applications necessitates structural health monitoring and/or nondestructive evaluation techniques that can provide quantitative and detailed damage information for layered plate-like components (such as composite laminates). Guided wavefield methods are at the basis of a number of promising techniques for the detection and the characterization of damage in plate-like structures. Among the processing techniques that have been proposed for guided wavefield analysis, the estimation of instantaneous and local wavenumbers can lead to effective metrics that quantify the size and the depth of delaminations in composite laminates. This article reports the application of both instantaneous and local wavenumber damage quantification techniques to guided wavefield data for delaminated composite laminates. The techniques are applied to experimental data for a simple single delamination case and to simulated data for a more complex multi-ply delamination case. The two techniques are compared in terms of accuracy in damage characterization and computational demand. The proposed methodologies can be considered as steps toward a hybrid structural health monitoring/nondestructive evaluation approach for damage assessment in composites.