Wavenumber Analysis Research Articles

Automatized localization of induced geothermal seismicity using robust time-domain array processing

The surveillance of geothermal seismicity is typically conducted using seismic networks, deployed around the power plants and subject to noise conditions in often highly urbanized areas. In contrast, seismic arrays can be situated at greater distances and allow monitoring of different power plants from one central location, less affected by noise interference. However, the effectiveness of arrays to monitor geothermal reservoirs is not well investigated and the increased distance to the source coincides with a decreased accuracy of the earthquake localizations. It is therefore essential to establish robust data processing and to obtain precise estimates of the location uncertainties. Here, we use time-domain array data processing and solve for the full 3-D slowness vector using robust linear regression. The approach implements a Biweight M-estimator, which yields stable parameter estimates and is well suited for real-time applications. We compare its performance to conventional least squares regression and frequency wavenumber analysis. Additionally, we implement a statistical approach based on changepoint analysis to automatically identify P- and S-wave arrivals within the recorded waveforms. The method can be seen as a simplification of autoregressive prediction. The estimated onsets facilitate reliable calculations of epicentral distances. We assess the performance of our methodology by comparison to network localizations for 77 induced earthquakes from the Landau and Insheim deep-geothermal reservoirs, situated in Rhineland-Palatinate, Germany. Our results demonstrate that we can differentiate earthquakes originating from both reservoirs and successfully localize the majority of events within the magnitude range of ML -0.2 to ML 1.3. The discrepancy between the two localization methods is mostly less than 1 km, which falls within the statistical errors. However, a few localizations deviate significantly, which can be attributed to poor observations during the winter of 2021/2022.

Source localization based on steered frequency-wavenumber analysis for sparse array.

When using a sparse array, locating the target signal of a high-frequency component is difficult. Although forecasting the direction in a sparse situation is challenging, the frequency-wavenumber (f-k) spectrum can simultaneously determine the direction and frequency of the analyzed signal. The striation of the f-k spectrum shifts along the wavenumber axis in a sparse situation, which reduces the spatial resolution required to determine the target's direction using the f-k spectrum. In this study, f-k spectra of a high-frequency signal were used for near-field source localization. Snapping shrimp sounds (5-24 kHz) from SAVEX15 (a shallow-water acoustic variability experiment conducted in May 2015) were used as the data source, and a simulation was used to evaluate the proposed method. Beam steering was performed before creating the f-k spectrum to improve spatial resolution. We found that the spatial resolution was improved, and the location of the sound source could be determined when a signal with beam steering was utilized. The shrimp sound from SAVEX15, a near-field broadband signal, was used to determine the shrimp's location (range, 38 m; depth, 100 m) and the tilt of the vertical line array. These results suggest that the proposed analysis helps to accurately estimate the location of sound source.

Characterization of Gravity Waves in Three Dimensions in the Daytime Thermosphere Using Combined Optical and Radio Measurements and Estimation of Horizontal Neutral Winds

AbstractGravity waves, which are considered omnipresent in the Earth's upper atmosphere, are generally investigated by monitoring the fluctuations in different atmospheric parameters. Here, we report the propagation characteristics of thermospheric gravity waves both in horizontal and vertical directions obtained using collocated optical and radio measurements from Ahmedabad, India for February 2021. The measurements of OI 630.0 nm dayglow emission rates over zenith are used to derive time periods of gravity waves. Wave number analyses of variations in the emission over a large field‐of‐view have been performed to derive gravity wave scale sizes and propagation characteristics in the horizontal direction. Time periods, horizontal scale sizes, and propagation directions are found to be in the range of 31–125 min, 78–243 km, and 203°–248° from east, respectively. Vertical wavelengths of the gravity waves are obtained from radio measurements and are in the range of 26–247 km. As the gravity wave characteristics are influenced by the ambient neutral winds, the measured gravity wave propagation characteristics in three dimensions have been used as inputs into the gravity wave dispersion relation to estimate the magnitudes of thermospheric horizontal neutral winds. These estimated daytime winds in the direction of wave propagation are found to be in the range of 1–105 ms−1, and they compare well with those measured independently from MIGHTI added HWM14 model‐derived winds. The daytime winds estimated by this approach are possibly the first of their kind as obtained from ground‐based measurements.

Impact of the Low Wavenumber Structure in the Initial Vortex Wind Analyses on the Prediction of the Intensification of Hurricane Patricia (2015)

AbstractEnsemble clustering analysis was performed to explore the role of the initial hurricane vortex‐scale wind structure in the prediction of the intensification of Hurricane Patricia (2015). Convection‐allowing ensemble forecasts were classified into spin‐down (SPD) and spin‐up (SPU) groups. Specifically, 10 members with an intensification rate >0 m/s and 10 members with an intensification rate <0 m/s for the first 6 hr were defined as the SPD and SPU members. The result showed that the tangential winds outside the inner‐core region in the SPD members were weaker compared to the SPU members. Additionally, the SPD members had a weaker inflow near the surface and a weaker outflow between the heights of 8 and 12 km than the SPU members. The SPU members showed more significant azimuthal asymmetry than the SPD members in the surface, tangential and radial winds. Wavenumber analysis showed that the low wavenumber components dominated the differences between the SPD and SPU members. Numerical experiments were conducted to test the hypothesis generated by the clustering analysis. It was found that the storm's maximum wind speed (MWS) intensified during the first 6 hr of the model forecast if only the low wavenumber structure in the SPU members was included in the initial conditions, whereas it decayed during the first 6 hr if only the low wavenumber structure in the SPD members was included. This result confirms that the low wavenumber structure of the initial wind analyses was important in predicting the intensity changes of Hurricane Patricia (2015).

On the occurrence of sandbars

In channelized alluvial rivers, alternating sandbars are formed by alternating scouring and sedimentation in the downstream direction. The conditions for the occurrence of sandbars are obtained in the stability analysis when the ratio of channel width to depth, β, exceeds a critical value due to innate instability of bottom. However, neither the mechanism nor the conditions for the occurrence of sandbars are well understood because it is difficult to quantify bars before or early in the occurrence of that using conventional measurement methods. We conducted experiments in a laboratory flume with slope as a variable to obtain dense data on the geometry of the water surface, bottom surface, and flow depth. We then performed a wavenumber analysis of the measured data to verify the validity of the results estimated by the stability analysis. We calculated the temporal variation of the correlation coefficients between the depth–water surface rwl and depth–bottom surface rbl along the development stage of the sandbar since flow state could be estimated from the relationship between the flow depth and the water surface/bottom surface. The temporal variation of the correlation suggests that the development process of sandbars can be divided into three phases, and the flow state is different for each phase. We also found that the occurrence time of sandbars would be defined based on the intersection of the temporal variation of rwl and rbl. Furthermore, we show that the formation process of the sandbar is independent of Fr.

Modified high-resolution wavenumber analysis for detection of pulse wave velocity using coefficient of variation of arterial wall acceleration waveforms

Modified high-resolution wavenumber analysis for detection of pulse wave velocity using coefficient of variation of arterial wall acceleration waveforms

On the potential of offshore sensors and array processing for improving seismic event detection and locations in the North Sea

SUMMARY The capability to detect and locate seismicity in the North Sea is necessary to ensure the safety of future CO2 storage operations. Currently, the event location capabilities are hampered by a one-sided network configuration, as only sensors deployed on the Norwegian mainland are providing good signal-to-noise ratios. However, applying array processing techniques to groups of offshore sensors deployed at oil and gas reservoirs can improve the signal-to-noise ratio at such sites, which generally suffer from high noise contamination originating from platforms or seismic shooting. In this study we test the feasibility of array processing to sensors deployed at the Grane oil field to improve earthquake event locations. Out of the 3400 nodes installed, we have access to two data subsets: (i) continuous data from 10 sensors spread out over the whole field and (ii) short segments of data from 30 sensors selected with the purpose to apply array methods. As the average interstation distance between the 10 sensors is 6 km, the recorded wavefield is not coherent and traditional array processing methods are not applicable. To optimize the usage of the 10-sensor offshore data for event location purposes, we develop a new method for incoherent array processing. We compute a characteristic function based on the kurtosis of the seismic traces prior to a frequency–wavenumber analysis. This method successfully determines slowness and backazimuth for almost all the evaluated earthquakes in this study. Using common coherent array processing, the 30 sensors show superior detection capabilities over the 10 sensors and will therefore decrease the detection threshold for seismic events in the region significantly if continuous data are available. We conclude that offshore sensors at Grane can be included as a part of a system for passive seismic monitoring in the North Sea. We recommend using the 30-sensors data for this purpose. However, as we only have access to continuous data from the 10 sensors in this study, we found that they are a good substitute when the 30-sensor data are not available.

Two-Dimensional Wavenumber Analysis Implemented in Ultrasonic Vector Doppler Method with Focused Transmit Beams

The multi-angle Doppler method was introduced for the estimation of velocity vectors by measuring axial velocities from multiple directions. We have recently reported that the autocorrelation-based velocity vector estimation could be ameliorated significantly by estimating the wavenumbers in two dimensions. Since two-dimensional wavenumber estimation requires a snapshot of an ultrasonic field, the method was first implemented in plane wave imaging. Although plane wave imaging is predominantly useful for examining blood flows at an extremely high temporal resolution, it was reported that the contrast in a B-mode image obtained with a few plane wave emissions was lower than that obtained with focused beams. In this study, the two-dimensional wavenumber analysis was first implemented in a framework with focused transmit beams. The simulations showed that the proposed method achieved an accuracy in velocity estimation comparable to that of the method with plane wave imaging. Furthermore, the performances of the methods implemented in focused beam and plane wave imaging were compared by measuring human common carotid arteries in vivo. Image contrasts were analyzed in normal and clutter-filtered B-mode images. The method with focused beam imaging achieved a better contrast in normal B-mode imaging, and similar velocity magnitudes and angles were obtained by both the methods with focused beam and plane wave imaging. In contrast, the method with plane wave imaging gave a better contrast in a clutter-filtered B-mode image and smaller variances in velocity magnitudes than those with focused beams.

RBF-FD discretization of the Navier-Stokes equations on scattered but staggered nodes

A semi-implicit fractional-step method that uses a staggered node layout and radial basis function-finite differences (RBF-FD) to solve the incompressible Navier-Stokes equations is developed. Polyharmonic splines (PHS) with polynomial augmentation (PHS+poly) are used to construct the global differentiation matrices. A systematic parameter study identifies a combination of stencil size, PHS exponent, and polynomial degree that minimizes the truncation error for a wave-like test function on scattered nodes. Classical modified wavenumber analysis is extended to RBF-FDs on heterogeneous node distributions and used to confirm that the accuracy of the selected 28-point stencil is comparable to that of spectral-like, 6th-order Padé-type finite differences. The Navier-Stokes solver is demonstrated on two benchmark problems, internal flow in a lid-driven cavity in the Reynolds number regime 102≤Re≤104, and open flow around a cylinder at Re=100 and 200. The combination of grid staggering and careful parameter selection facilitates accurate and stable simulations at significantly lower resolutions than previously reported, using more compact RBF-FD stencils, without special treatment near solid walls, and without the need for hyperviscosity or other means of regularization.

Attenuation band splitting in a finite plate strip with two-dimensional acoustic black holes

Acoustic black holes (ABH) lattice is an effective way to realize the bandgaps (BGs) through lightweight design. The predicted BGs using infinitely large periodic structures are commonly believed to enable corresponding attenuation bands (ABs) even the structure only contains a limited number of unit cells. This paper reports an unusual splitting phenomenon in the ABs of a finite plate strip with two-dimensional (2D) ABHs. A wide AB predicted by BG is split into two ABs in the finite strip separated by a high resonance peak, which challenges the aforementioned common belief. Results show that the increasing number of periodic unit cells has a negligible effect on the frequency and amplitude of the resonance peak. Its underlying mechanism is clarified by investigating the relationship between the BGs and ABs, obtained from infinite and finite periodic plate strips with the same 2D ABHs respectively. Complex wavenumber analysis results show the existence of two flexural wave modes in the considered frequency range due to the finite width of the strip. Two types of BGs, in which one satisfies the Bragg condition, are found when both wave modes fail to propagate. Analysis of structural modes further uncover the physical reason for the emerging AB splitting: a certain natural frequency of the boundary portion of the 2D ABHs strip falling into the theoretical BG. Due to the boundary reflection and energy concentration, this structural portion can loosely be regarded as a clamped-free plate strip. It is demonstrated that the resonance peak splitting the AB becomes tunable via disrupting the structural details of finite ABH strips. In particular, extending the terminal end of a finite ABH strip allows the change in boundary portion, which is conducive to eliminating the splitting peak. The reported study is helpful for the design of finite periodic structures with broadband attenuation bandgap.

Can oceanic flows be heard? Abyssal melodies

Fluid flows generate an acoustic noise field. In principle, oceanic flows on varying time and length scales produce a sound field and its detectability is considered here. A fragile lower bound analysis is made of the acoustic signature, using the Lighthill theory, of a simple train of boundary vortices generated by baroclinic tidal flows. Subject to numerous assumptions, the accompanying sound should be detectable within the hum band of seismo-acoustic pressure fields, and more generally, across the entire oceanic spectrum-likely through wave number analyses of spatially coherent acoustic array data.

Subsurface characterization using passive seismic in the urban area of Dublin City, Ireland

AbstractWe investigate the shallow geological features and previously unknown depths of the main geological interfaces in the highly urbanized Dublin City area by using a set of complementary ambient noise seismic methods. A sparse seismic array, composed of 19 broadband sensors, was installed in a 5 × 5 km area across the city centre. The dataset was acquired over a 5‐month‐long deployment and used to perform horizontal/vertical spectral ratios and frequency–wave number analysis cross‐correlation interferometry. From the horizontal/vertical high‐frequency resonance peaks, we estimate the depth to bedrock. Then we use a subset of six seismic stations to obtain the frequency–wavenumber Rayleigh wave fundamental mode in the 0.8–3 Hz frequency band. Next, the ambient noise dataset is cross‐correlated in order to extract the empirical Green's functions before measuring surface‐wave phase velocities by undertaking a dispersion analysis in the 0.5–9 Hz frequency band. Then, a Monte Carlo global optimization algorithm is used to invert the phase velocity dispersion measurements. We generate a reference one‐dimensional depth S velocity profile along with a set of localized one‐dimensional S velocity profiles. Finally, a smooth three‐dimensional shear wave velocity model is derived for the top 1000 m in the sedimentary Dublin Basin. From the one‐dimensional velocity profiles and three‐dimensional shear velocity model, with some sensitivity down to approximately 1.2 km, we observe velocity changes with depth that highlight the presence of three consistent interfaces. We discuss interpretative scenarios that correlate the velocity features to potential stratigraphic boundaries occurring within the sedimentary Dublin Basin, suggesting that basement rocks could be located at a depth significantly greater than the top kilometre. The results of this study may indicate that future geothermal studies should be directed at structures deeper than 1 km, towards the bottom of the sedimentary basin beneath the city.

Investigations of equatorial plasma bubbles as observed in the OI 630 nm nightglow emissions over off-equatorial and low-latitudinal locations over Indian longitudes

Equatorial plasma bubbles (EPBs), generated over the dip-equator during the sunset time, get mapped over to the low-latitudes through the geomagnetic field lines and contribute to structures of different scales in the nighttime thermosphere-ionosphere system. In the present study, the EPBs are observed in the OI 630 nm nightglow emissions over Mt. Abu and Kolhapur using a High Throughput Imaging Echelle spectrograph (HiTIES) and All-Sky Imager (ASI), respectively. Similar periodicities obtained in these measurements consistently corroborate the presence of the EPBs at both these locations. The strength of the Pre-Reversal Enhancement (PRE) in the zonal electric field has been investigated using the ionosonde data over Trivandrum, a dip-equatorial location. The strength of PRE decides the latitudinal extent of the EPBs. The eastward movement of the EPBs has been estimated to be in the range of 190–90 ms−1, decreasing in magnitude from post-sunset to midnight. The wave number analysis carried out using the observed OI 630 nm emissions in the zonal direction has resulted in a contrasting behaviour during the presence/absence of the EPBs. Based on the analysis of around 1300 images of data, it is revealed that the scale sizes in the range of 250–300 km are omnipresent, whereas, shorter scale sizes (50–250 km) are present only during the presence of EPBs. It is inferred that these shorter scale size gravity waves played a significant role in the seeding the perturbation of the EPBs.

Synthesis, spectral, structural and antimicrobial activities of Ethyl-4-｛[-(1-(2-(4-nitrobenzoyl)hydrazono)ethyl]｝-3,5-dimethyl-1H-pyrrole-2-carboxylate

A new ethyl-4-{[-(1-(2-(4-nitrobenzoyl)hydrazono)ethyl]}-3,5-dimethyl-1H-pyrrole-2-carboxylate (ENBHEDPC) has been synthesized, characterized (FT-IR, 1H and 13C NMR, UV-Visible, DART-Mass spectroscopy, elemental analysis) and evaluated for antituberculosis and antibacterial activity. Quantum chemical calculations have been performed by DFT level of theory using the B3LYP functional and 6–311++G(d,p) as basis set. The 1H and 13C NMR chemical shifts and electronic transitions are calculated using gage including atomic orbitals (GIAO) and time dependent density functional theory (TD-DFT) approach in DMSO as solvent. Combined theoretical and experimental wavenumber analyses confirm the existence of dimer. Hydrogen bonding and other multiple interactions have been analyzed by Bader's ‘Atoms in molecules’ AIM theory in detail. The intermolecular hydrogen bond energy of dimer is calculated as -12.28 kcal/mol. The calculated thermodynamic parameters show that reaction is exothermic and non-spontaneous at room temperature. The local reactivity descriptors analyses determine the reactive sites within molecule by nucleophilic attack. Nonlinear optical (NLO) behavior of title compound is investigated by the computed value of the first hyperpolarizability (β0). The synthesized compound show good antimicrobial and antituberculosis activity.

Rbf-Fd Discretization of the Navier-Stokes Equations Using Staggered Nodes

A semi-implicit fractional-step method that uses a staggered node layout and radial basis function-finite differences (RBF-FD) to solve the incompressible Navier-Stokes equations is developed. Polyharmonic splines (PHS) with polynomial augmentation (PHS+poly) are used to construct the global differentiation matrices. A systematic parameter study identifies a combination of stencil size, PHS exponent, and polynomial degree that minimizes the truncation error for a wave-like test function on scattered nodes. Classical modified wavenumber analysis is extended to RBF-FDs on heterogeneous node distributions and used to confirm that the accuracy of the selected 28-point stencil is comparable to that of spectral-like, 6th-order Pad\'e-type finite differences. The Navier-Stokes solver is demonstrated on two benchmark problems, internal flow in a lid-driven cavity in the Reynolds number regime $10^2 \leq \mathrm{Re}\leq 10^4$, and open flow around a cylinder at $\mathrm{Re}= 100$ and $200$. The combination of grid staggering and careful parameter selection facilitates accurate and stable simulations at significantly lower resolutions than previously reported, using more compact RBF-FD stencils, without special treatment near solid walls, and without the need for hyperviscosity or other means of regularization.

Wavefield Analysis Tools for Wavenumber and Velocities Extraction in Polar Coordinates.

Experimental characterization of Lamb waves in plate-like structures overcomes the intrinsic limits of a priori semianalytical finite element simulations, where material inaccuracies and nonidealities cannot be easily considered. Unfortunately, the experimental extraction of guided wave dispersion curves, and especially their polar representation along different directions of propagation at a given frequency, is not trivial. In nonisotropic materials, such analysis is a key aspect for a reliable and robust characterization of the behavior of waves. In this work, by exploiting scanning laser Doppler vibrometer measurements with narrowband excitation, two different signal processing methods for the extraction of the wavenumber polar representation at the excitation frequency are investigated and characterized. The first method is based on a distance regularized level set (DRLSE) algorithm, widely used in image processing and computer vision but, to the best of the author's knowledge, never used in the Lamb waves' field. The second method is based on the 2-D sparse wavenumber analysis which exploits the wavefield sparse representation in the wavenumber domain. With a precise and reliable extraction of the wavenumber characteristic in the k -space, the polar representations at the excitation frequency of phase and group velocities can be estimated. The former, by exploiting the well-known wavenumber-frequency relation, the latter, instead, by computing numerical derivative among wavenumbers at multiple frequencies. The methodology has been validated on three different composite plates with different degrees of nonisotropy properties. The results show the effectiveness of the two methods, highlighting the advantages and disadvantages of both.

Effects of non-axisymmetric internal structures on vibro-acoustic characteristics of a submerged cylindrical shell using wavenumber analysis

Non-axisymmetric structures (foundations, equipment etc.) have been extensively arranged in axisymmetric shells and proved to change the mechanism of vibration and sound radiation of the shell. To study the effect of non-axisymmetric internal structures, a hybrid analytic-numerical method is developed for vibro-acoustic characteristics analyses of a submerged cylindrical shell. The wave based method (WBM) and finite element method (FEM) are respectively employed to model the shell and internal structures. A precise coupling is carried out to assemble the shell and internal structures, which fully considers all six degree of freedoms (DOFs) at junctions between the two substructures. The acoustic pressure is described by the Helmholtz integral equation. According to the axisymmetric property of the shell, the surface integral is reduced into the curve integral. Subsequently, the final vibro-acoustic model is obtained by introducing the acoustic pressure into the dynamic equation of the coupled structure. The high accuracy of the present method is firstly demonstrated by FEM and boundary element method (BEM). Then by comparing vibro-acoustic responses of the shell with and without internal structures in wave domain, the effects of non-axisymmetric internal structures on vibro-acoustic characteristics of the shell are discussed.

Estimation of Shear-Wave Velocity Structures in Taichung, Taiwan, Using Array Measurements of Microtremors

Near-surface S-wave velocity structures (VS) are crucial in site-effect studies and ground-motion simulations or predictions. We explored S-wave velocity structures in Taichung, the second-largest city in Taiwan by population, by employing array measurements of microtremors at a total of 53 sites. First, the fundamental-mode dispersion curves of Rayleigh waves were estimated by adopting the frequency–wavenumber analysis method. Second, the surface-wave inversion technique was used to calculate the S-wave velocity structures of the area. At many sites, observed phase velocities were almost flat, with a phase velocity of approximately 800–1300 m/s in the frequency range of 0.6–2 Hz. A high-velocity zone (VS of 900–1500 m/s) with a convex shape was observed at the shallow S-wave structures of these sites (depths of 50–500 m). On the basis of the inversion results, we constructed two-dimensional and three-dimensional contour maps to elucidate the variations of VS structures in Taichung. According to VS-contour maps at different depths, lowest S-wave velocities are found at the western coastal plain, whereas highest S-wave velocities appear on the eastern side. The S-wave velocity gradually decreases from east to west. Moreover, the S-wave velocity of the Tertiary bedrock is assumed to be 1500 m/s in the area. According to the depth-contour map (VS = 1500 m/s), the depths of the bedrock range from 250 m (the eastern part) to 1550 m (the western part). The thicknesses of the alluvium gradually decrease from west to east. Our results are consistent with the geology of the Taichung area.

Shallow S-Wave Velocity Structure in the Middle-Chelif Basin, Algeria, Using Ambient Vibration Single-Station and Array Measurements

In order to better assess the seismic hazard in the northern region of Algeria, the shear-wave velocity structure in the Middle-Chelif Basin is estimated using ambient vibration single-station and array measurements. The Middle-Chelif Basin is located in the central part of the Chelif Basin, the largest of the Neogene sedimentary basins in northern Algeria. This basin hosts the El-Asnam fault, one of the most important active faults in the Mediterranean area. In this seismically active region, most towns and villages are built on large unconsolidated sedimentary covers. Application of the horizontal-to-vertical spectral ratio (HVSR) technique at 164 sites, and frequency–wavenumber (F–K) analysis at 7 other sites, allowed for the estimation of the ground resonance frequencies, shear-wave velocity profiles, and sedimentary cover thicknesses. The electrical resistivity tomography method was used at some sites to further constrain the thickness of the superficial sedimentary layers. The soil resonance frequencies range from 0.75 Hz to 12 Hz and the maximum frequency peak amplitude is 6.2. The structure of the estimated shear-wave velocities is presented in some places as 2D profiles to help interpret the existing faults. The ambient vibration data allowed us to estimate the maximum depth in the Middle-Chelif Basin, which is 760 m near the city of El-Abadia.

Model-based autonomous plate defects visualization method for quantitative wall-thinning estimation

The visualization of the wall thinning of plate-like structures using scanning laser Doppler vibrometry (SLDV) is a promising method in nondestructive evaluation using laser ultrasonics. In particular, the Lamb-wave-based SLDV method that uses continuous excitation exhibits excellent performance for the estimation of the wall thinning of plates. Currently, plate thickness is quantitatively evaluated based on wavenumber analysis using measured signals. However, it is difficult to estimate plate thickness automatically below the product of frequency and thickness of 6 MHz·mm without knowing the wavenumber sensitivity and minimum wavenumber distance from reference mode owing to the lack of the physical understanding of Lamb waves. In this study, a model-based autonomous plate defects visualization method is proposed for the quantitative imaging of the wall thinning of plates so that inspectors can use scanning laser Doppler vibrometry (SLDV) without any knowledge of Lamb waves and its signal processing. Interdigital-transducer-based SLDV is utilized to validate the proposed method, and a 6-mm-thick carbon steel plate with 1–8% wall thinning, and a 2-mm-thick aluminum plate with Y-shaped wall thinning are used. Experiments demonstrate that the capability of the proposed method for detecting wall thinning in plates is equivalent to that of manual plate defects visualization method.