Enhancing microseismic/acoustic emission source localization accuracy with an outlier-robust kernel density estimation approach

Monitoring sensors in complex engineering environments often record abnormal data,leading to significant positioning errors.To reduce the influence of abnormal arrival times,we introduce an innovative,outlier-robust localization method that integrates kernel density estimation(KDE)with damping linear correction to enhance the precision of microseismic/acoustic emission(MS/AE)source positioning.Our approach systematically addresses abnormal arrival times through a three-step process:initial location by 4-arrival combinations,elimination of outliers based on three-dimensional KDE,and refinement using a linear correction with an adaptive damping factor.We validate our method through lead-breaking experiments,demonstrating over a 23%improvement in positioning accuracy with a maximum error of 9.12 mm(relative error of 15.80%)—outperforming 4 existing methods.Simulations under various system errors,outlier scales,and ratios substantiate our method’s superior performance.Field blasting experiments also confirm the practical applicability,with an average positioning error of 11.71 m(relative error of 7.59%),compared to 23.56,66.09,16.95,and 28.52 m for other methods.This research is significant as it enhances the robustness of MS/AE source localization when confronted with data anomalies.It also provides a practical solution for real-world engineering and safety monitoring applications.

Piecewise Acoustic Source Imaging with Unknown Speed of Sound Using a Level-Set Method

We investigate the following inverse problem:starting from the acoustic wave equation,reconstruct a piecewise constant passive acoustic source from a single boundary temporal measurement without knowing the speed of sound.When the amplitudes of the source are known a priori,we prove a unique determination result of the shape and propose a level set algorithm to reconstruct the singularities.When the singularities of the source are known a priori,we show unique determination of the source amplitudes and propose a least-squares fitting algorithm to recover the source amplitudes.The analysis bridges the low-frequency source inversion problem and the inverse problem of gravimetry.The proposed algorithms are validated and quantitatively evaluated with numerical experiments in 2D and 3D.

3D near-surface P-wave velocity structure imaging with Distributed Acoustic Sensing and electric hammer source

Distributed Acoustic Sensing(DAS) is an emerging technique for ultra-dense seismic observation, which provides a new method for high-resolution sub-surface seismic imaging. Recently a large number of linear DAS arrays have been used for two-dimensional S-wave near-surface imaging in urban areas. In order to explore the feasibility of three-dimensional(3D) structure imaging using a DAS array, we carried out an active source experiment at the Beijing National Earth Observatory. We deployed a 1 km optical cable in a rectangular shape, and the optical cable was recast into 250 sensors with a channel spacing of 4 m. The DAS array clearly recorded the P, S and surface waves generated by a hammer source. The first-arrival P wave travel times were first picked with a ShortTerm Average/Long-Term Average(STA/LTA) method and further manually checked. The P-wave signals recorded by the DAS are consistent with those recorded by the horizontal components of short-period seismometers. At shorter source-receiver distances, the picked P-wave arrivals from the DAS recording are consistent with vertical component recordings of seismometers, but they clearly lag behind the latter at greater distances.This is likely due to a combination of the signal-to-noise ratio and the polarization of the incoming wave. Then,we used the Tomo DD software to invert the 3D P-wave velocity structure for the uppermost 50 m with a resolution of 10 m. The inverted P-wave velocity structures agree well with the S-wave velocity structure previously obtained through ambient noise tomography. Our study indicates the feasibility of 3D near-surface imaging with the active source and DAS array. However, the inverted absolute velocity values at large depths may be biased due to potential time shifts between the DAS recording and seismometer at large source-receiver distances.

Study on Joint Method of 3D Acoustic Emission Source Localization Simplex and Grid Search Scanning

Acoustic emission(AE)source localization is a fundamental element of rock fracture damage imaging.To improve the efficiency and accuracy of AE source localization,this paper proposes a joint method comprising a three-dimensional(3D)AE source localization simplex method and grid search scanning.Using the concept of the geometry of simplexes,tetrahedral iterations were first conducted to narrow down the suspected source region.This is followed by a process of meshing the region and node searching to scan for optimal solutions,until the source location is determined.The resulting algorithm was tested using the artificial excitation source localization and uniaxial compression tests,after which the localization results were compared with the simplex and exhaustive methods.The results revealed that the localization obtained using the proposed method is more stable and can be effectively avoided compared with the simplex localization method.Furthermore,compared with the global scanning method,the proposed method is more efficient,with an average time of 10%–20%of the global scanning localization algorithm.Thus,the proposed algorithm is of great significance for laboratory research focused on locating rupture damages sustained by large-sized rock masses or test blocks.

Distribution pattern of acoustic and streaming field during multi-source ultrasonic melt treatment process

The ultrasonic melt treatment(UMT)is widely used in the fields of casting and metallurgy.However,there are certain drawbacks associated with the conventional process of single-source ultrasonic(SSU)treatment,such as the fast attenuation of energy and limited range of effectiveness.In this study,the propagation models of SSU and four-source ultrasonic(FSU)in Al melt were respectively established,and the distribution patterns of acoustic and streaming field during the ultrasonic treatment process were investigated by numerical simulation and physical experiments.The simulated results show that the effective cavitation zone is mainly located in a small spherical region surrounding the end of ultrasonic horn during the SSU treatment process.When the FSU is applied,the effective cavitation zone is obviously expanded in the melt.It increases at first and then decreases with increasing the vibration-source spacing(Lv)from 30 mm to 100 mm.Especially,when the Lv is 80 mm,the area of effective cavitation zone reaches the largest,indicating the best effect of cavitation.Moreover,the acoustic streaming level and flow pattern in the melt also change with the increase of Lv.When the Lv is 80 mm,both the average flow rate and maximum flow rate of the melt reach the highest,and the flow structure is more stable and uniform,with the typical morphological characteristics of angular vortex,thus significantly expanding the range of acoustic streaming.The accuracy of the simulation results was verified by physical experiments of glycerol aqueous solution and tracer particles.

Three-dimensional analytical solution of acoustic emission source location for cuboid monitoring network without pre-measured wave velocity

To find analytical solutions of nonlinear systems for locating the acoustic emission/microseismic（AE/MS） source without knowing the wave velocity of structures, the sensor location coordinates were simplified as a cuboid monitoring network. Different locations of sensors on upper and lower surfaces were considered and used to establish nonlinear equations. Based on the proposed functions of time difference of arrivals, the analytical solutions were obtained using five sensors under three networks. The proposed analytical solutions were validated using authentic data of numerical tests and experiments. The results show that located results are consistent with authentic data, and the outstanding characteristics of the new solution are that the solved process is not influenced by the wave velocity knowledge and iterated algorithms.

Effects of sonic speed on location accuracy of acoustic emission source in rocks

To quantitatively study the location errors induced by deviation of sonic speed, the line and plane location tests were carried out. A broken pencil was simulated as acoustic emission source in the rocks. The line and plane location tests were carried out in the granite rod using two sensors and the cube of marble using four sensors, respectively. To compare the position accuracy between line and plane positions, the line poison test was also carried out on the marble surface. The results show that for line positioning, the maximum error of absolute distance is about 0.8 cm. With the speed difference of 200 m/s, the average value of absolute difference from the position error is about 0.4 cm. For the plane positioning, in the case of the sensor array of 30 cm, the absolute positioning distance is up to 8.7 cm. It can be seen that the sonic speed seriously impacts on the plane positioning accuracy. The plane positioning error is lager than the line positioning error, which means that when the line position can satisfy the need in practical engineering, it is better to use the line position instead of the plane location. The plane positioning error with the diagonal speed is the minimum one.

Three-dimensional analytical solution of acoustic emission or microseismic source location under cube monitoring network

To find the analytical solution of the acoustic emission/microseismic(AE/MS) source location coordinates, the sensor location coordinates were optimized and simplified. A cube monitoring network of sensor location was selected, and the AE/MS source localization equations were established. A location method with P-wave velocity by analytical solutions (P-VAS) was obtained with these equations. The virtual location tests show that the relocation results of analytical method are fully consistent with the actual coordinates for events both inside and outside the monitoring network; whereas the location error of traditional time difference method is between 0.01 and 0.03 m for events inside the sensor array, and the location errors are larger, which is up to 1080986 m for events outside the sensor array. The broken pencil location tests were carried out in the cross section of 100 mm×98 mm, 350 mm-length granite rock specimen using five AE sensors. Five AE sources were relocated with the conventional method and the P-VAS method. For the four events outside monitoring network, the positioning accuracy by P-VAS method is higher than that by the traditional method, and the location accuracy of the larger one can be increased by 17.61 mm. The results of both virtual and broken pencil location tests show that the proposed analytical solution is effective to improve the positioning accuracy. It can locate the coordinates of AE/MS source only using simple four arithmetic operations, without determining the fitting initial value and iterative calculation, which can be solved by a conventional calculator or Microsoft Excel.

Acoustic emission source identification based on harmonic wavelet packet and support vector machine

In order to solve the fatigue damage identification problem of helicopter moving components, a new approach for acoustic emission （AE） source type identification based on the harmonic wavelet packet （HWPT） feature extraction and the hierarchy support vector machine （H-SVM） classifier is proposed. After a four-level decomposition of the HWPT, the energy feature of AE signals in different frequency bands is extracted, which overcomes the shortcomings of the traditional wavelet packet including energy leakage, and inflexible frequency band selection and different frequency resolutions on different levels. The H-SVM classifier is trained with a subset of the experimental data for known AE source types and tested using the remaining set of data. The results of pressure-off experiments on the specimens of carbon fiber materials indicate that the proposed approach can effectively implement the AE source type identification, and has a better performance in terms of computational efficiency and identification accuracy than the wavelet packet （WPT） feature extraction.

Mechanical properties and acoustic emission characteristics of soft rock with different water contents under dynamic disturbance

Uniaxial compression tests and cyclic loading acoustic emission tests were conducted on 20%,40%,60%,80%,dry and saturated muddy sandstone by using a creep impact loading system to investigate the mechanical properties and acoustic emission characteristics of soft rocks with different water contents under dynamic disturbance.The mechanical properties and acoustic emission characteristics of muddy sandstones at different water contents were analysed.Results of experimental studies show that water is a key factor in the mechanical properties of rocks,softening them,increasing their porosity,reducing their brittleness and increasing their plasticity.Under uniaxial compression,the macroscopic damage characteristics of the muddy sandstone change from mono-bevel shear damage and‘X’type conjugate bevel shear damage to a roadway bottom-drum type damage as the water content increases.Dynamic perturbation has a strengthening effect on the mechanical properties of samples with 60%and less water content,and a weakening effect on samples with 80%and more water content,but the weakening effect is not obvious.Macroscopic damage characteristics of dry samples remain unchanged,water samples from shear damage and tensile–shear composite damage gradually transformed into cleavage damage,until saturation transformation monoclinic shear damage.The evolution of acoustic emission energy and event number is mainly divided into four stages:loading stage(Ⅰ),dynamic loading stage(Ⅱ),yield failure stage(Ⅲ),and post-peak stage(Ⅳ),the acoustic emission characteristics of the stages were different for different water contents.The characteristic value of acoustic emission key point frequency gradually decreases,and the damage degree of the specimen increases,corresponding to low water content—high main frequency—low damage and high water content—low main frequency—high damage.

ROBUST ACOUSTIC SOURCE LOCALIZATION FOR DIGITAL HEARING AIDS IN NOISE AND REVERBERANT ENVIRONMENT

A new method in digital hearing aids to adaptively localize the speech source in noise and reverberant environment is proposed. Based on the room reverberant model and the multichannel adaptive eigenvalue decomposition （MCAED） algorithm, the proposed method can iteratively estimate impulse response coefficients between the speech source and microphones by the adaptive subgradient projection method. Then, it acquires the time delays of microphone pairs, and calculates the source position by the geometric method. Compared with the traditional normal least mean square （NLMS） algorithm, the adaptive subgradient projection method achieves faster and more accurate convergence in a low signal-to-noise ratio （SNR） environment. Simulations for glasses digital hearing aids with four-component square array demonstrate the robust performance of the proposed method.

A Review on Sources,Extractions and Analysis Methods of a Sustainable Biomaterial:Tannins

Condensed and hydrolysable tannins are non-toxic natural polyphenols that are a commercial commodity industrialized for tanning hides to obtain leather and for a growing number of other industrial applications mainly to substitute petroleum-based products.They are a definite class of sustainable materials of the forestry industry.They have been in operation for hundreds of years to manufacture leather and now for a growing number of applications in a variety of other industries,such as wood adhesives,metal coating,pharmaceutical/medical applications and several others.This review presents the main sources,either already or potentially commercial of this forestry by-materials,their industrial and laboratory extraction systems,their systems of analysis with their advantages and drawbacks,be these methods so simple to even appear primitive but nonetheless of proven effectiveness,or very modern and instrumental.It constitutes a basic but essential summary of what is necessary to know of these sustainable materials.In doing so,the review highlights some of the main challenges that remain to be addressed to deliver the quality and economics of tannin supply necessary to fulfill the industrial production requirements for some materials-based uses.

Time-lagged Effects of the Spring Atmospheric Heat Source over the Tibetan Plateau on Summer Precipitation in Northeast China during 1961–2020:Role of Soil Moisture

The spring atmospheric heat source(AHS)over the Tibetan Plateau(TP)has been suggested to affect the Asian summer monsoon and summer precipitation over South China.However,its influence on the summer precipitation in Northeast China(NEC)remains unknown.The connection between spring TP AHS and subsequent summer precipitation over NEC from 1961 to 2020 is analyzed in this study.Results illustrate that stronger spring TP AHS can enhance subsequent summer NEC precipitation,and higher soil moisture in the Yellow River Valley-North China region(YRVNC)acts as a bridge.During spring,the strong TP AHS could strengthen the transportation of water vapor to East China and lead to excessive rainfall in the YRVNC.Thus,soil moisture increases,which regulates local thermal conditions by decreasing local surface skin temperature and sensible heat.Owing to the memory of soil moisture,the lower spring sensible heat over the YRVNC can last until mid-summer,decrease the land–sea thermal contrast,and weaken the southerly winds over the East Asia–western Pacific region and convective activities over the South China Sea and tropical western Pacific.This modulates the East Asia–Pacific teleconnection pattern,which leads to a cyclonic anomaly and excessive summer precipitation over NEC.

Quantitative effect of kerogen type on the hydrocarbon generation potential of Paleogene lacustrine source rocks,Liaohe Western Depression,China

Kerogen types exert a decisive effect on the onset and capacity of hydrocarbon generation of source rocks.Lacustrine source rocks in the Liaohe Western Depression are characterized by thick deposition,high total organic carbon(TOC)content,various kerogen types,and a wide range of thermal maturity.Consequently,their hydrocarbon generation potential and resource estimation can be misinterpreted.In this study,geochemical tests,numerical analysis,hydrocarbon generation kinetics,and basin modeling were integrated to investigate the differential effects of kerogen types on the hydrocarbon generation potential of lacustrine source rocks.Optimized hydrocarbon generation and expulsion(HGE)models of different kerogen types were established quantitatively upon abundant Rock-Eval/TOC/vitrinite reflectance(R_(o))datasets.Three sets of good-excellent source rocks deposited in the fourth(Es4),third(Es3),and first(Es1)members of Paleogene Shahejie Formation,are predominantly types I-II_(1),II_(1)-II_(2),and II-III,respectively.The activation energy of types I-II_(2)kerogen is concentrated(180-230 kcal/mol),whereas that of type III kerogen is widely distributed(150-280 kcal/mol).The original hydrocarbon generation potentials of types I,II_(1),II_(2),and III kerogens are 790,510,270,and 85 mg/g TOC,respectively.The Ro values of the hydrocarbon generation threshold for type I-III source rocks gradually increase from 0.42%to 0.74%,and Ro values of the hydrocarbon expulsion threshold increase from 0.49%to 0.87%.Types I and II_(1)source rocks are characterized by earlier hydrocarbon generation,more rapid hydrocarbon expulsion,and narrower hydrocarbon generation windows than types II_(2)and III source rocks.The kerogen types also affect the HGE history and resource potential.Three types(conventional,tight,and shale oil/gas)and three levels(realistic,expected,and prospective)of hydrocarbon resources of different members in the Liaohe Western Depression are evaluated.Findings suggest that the Es3 member has considerable conventional and unconventional hydrocarbon resources.This study can quantitatively characterize the hydrocarbon generation potential of source rocks with different kerogen types,and facilitate a quick and accurate assessment of hydrocarbon resources,providing strategies for future oil and gas exploration.

The coupling control of biological precursors and environmental factors onβ-carotane enrichment in alkaline lacustrine source rocks:A case study from the Fengcheng formation in the western Junggar Basin,NW China

The organic-rich mudstones and dolostones of the Permian Fengcheng Formation(Fm.)are typically alkaline lacustrine source rocks,which are typified by impressively abundantβ-carotane.Abundant β-carotane has been well acknowledged as an effective indicator of biological sources or depositional environments.However,the specific biological sources of β-carotane and the coupling control of biological sources and environmental factors on the enrichment of β-carotane in the Fengcheng Fm.remains obscure.Based on a comprehensive investigation of the bulk,molecular geochemistry,and organic petrology of sedimentary rocks and the biochemistry of phytoplankton in modern alkaline lakes,we proposed a new understanding of the biological precursors of β-carotane and elucidated the enrichment mechanism of β-carotane in the Fengcheng Fm.The results show that the biological precursors crucially control the enrichment of β-carotane in the Fengcheng Fm.The haloalkaliphilic cyanobacteria are the primary biological sources of β-carotane,which is suggested by a good positive correlation between the 2-methylhopane index,7-+8-methyl heptadecanes/C_(max),C_(29%),and β-carotane/C_(max)in sedimentary rocks and the predominance of cyanobacteria with abundantβ-carotene in modern alkaline lakes.The enrichment of β-carotane requires the reducing condition,and the paleoredox state that affects the enrichment of β-carotane appears to have a threshold.The paleoclimate conditions do not considerably impact the enrichment of β-carotane,but they have some influence on the water's paleosalinity by affecting evaporation and precipitation.While it does not directly affect the enrichment of β-carotane in the Fengcheng Fm.,paleosalinity does have an impact on the cyanobacterial precursor supply and the preservation conditions.

Reservoir heterogeneity analysis using multi-directional textural attributes from deep learning-based enhanced acoustic impedance inversion:A study from Poseidon,NW shelf Australia

Reservoir heterogeneities play a crucial role in governing reservoir performance and management.Traditionally,detailed and inter-well heterogeneity analyses are commonly performed by mapping seismic facies change in the seismic data,which is a time-intensive task.Many researchers have utilized a robust Grey-level co-occurrence matrix(GLCM)-based texture attributes to map reservoir heterogeneity.However,these attributes take seismic data as input and might not be sensitive to lateral lithology variation.To incorporate the lithology information,we have developed an innovative impedance-based texture approach using GLCM workflow by integrating 3D acoustic impedance volume(a rock propertybased attribute)obtained from a deep convolution network-based impedance inversion.Our proposed workflow is anticipated to be more sensitive toward mapping lateral changes than the conventional amplitude-based texture approach,wherein seismic data is used as input.To evaluate the improvement,we applied the proposed workflow to the full-stack 3D seismic data from the Poseidon field,NW-shelf,Australia.This study demonstrates that a better demarcation of reservoir gas sands with improved lateral continuity is achievable with the presented approach compared to the conventional approach.In addition,we assess the implication of multi-stage faulting on facies distribution for effective reservoir characterization.This study also suggests a well-bounded potential reservoir facies distribution along the parallel fault lines.Thus,the proposed approach provides an efficient strategy by integrating the impedance information with texture attributes to improve the inference on reservoir heterogeneity,which can serve as a promising tool for identifying potential reservoir zones for both production benefits and fluid storage.

A robust triaxial localization method of AE source using refraction path

Acoustic emission(AE)localization algorithms based on homogeneous media or single-velocity are less accurate when applied to the triaxial localization experiments.To the end,a robust triaxial localization method of AE source using refraction path is proposed.Firstly,the control equation of the refraction path is established according to the sensor coordinates and arrival times.Secondly,considering the influence of time-difference-of-arrival(TDOA)errors,the residual of the governing equation is calculated to estimate the equation weight.Thirdly,the refraction points in different directions are solved using Snell’s law and orthogonal constraints.Finally,the source coordinates are iteratively solved by weighted correction terms.The feasibility and accuracy of the proposed method are verified by pencil-lead breaking experiments.The simulation results show that the new method is almost unaffected by the refraction ratio,and always holds more stable and accurate positioning performance than the traditional method under different ratios and scales of TDOA outliers.

Acoustic radiation force on a cylindrical composite particle with an elastic thin shell and an internal eccentric liquid column in a plane ultrasonic wave field

A model with three-layer structure is introduced to explore the acoustic radiation force(ARF)on composite particles with an elastic thin shell.Combing acoustic scattering of cylinder and the thin-shell theorem,the ARF expression was derived,and the longitudinal and transverse components of the force and axial torque for an eccentric liquid-filled composite particle was obtained.It was found that many factors,such as medium properties,acoustic parameters,eccentricity,and radius ratio of the inner liquid column,affect the acoustic scattering field of the particle,which in turn changes the forces and torque.The acoustic response varies with the particle structures,so the resonance peaks of the force function and torque shift with the eccentricity and radii ratio of particle.The acoustic response of the particle is enhanced and exhibits higher force values due to the presence of the elastic thin shell and the coupling effect with the eccentricity of the internal liquid column.The decrease of the inner liquid density may suppress the high-order resonance peaks,and internal fluid column has less effects on the change in force on composite particle at ka>3,while limited differences exist at ka<3.The axial torque on particles due to geometric asymmetry is closely related to ka and the eccentricity.The distribution of positive and negative force and torque along the axis ka exhibits that composite particle can be manipulated or separated by ultrasound.Our theoretical analysis can provide support for the acoustic manipulation,sorting,and targeting of inhomogeneous particles.

The viscous strip approach to simplify the calculation of the surface acoustic wave generated streaming

In recent decades,the importance of surface acoustic waves,as a biocompatible tool to integrate with microfluidics,has been proven in various medical and biological applications.The numerical modeling of acoustic streaming caused by surface acoustic waves in microchannels requires the effect of viscosity to be considered in the equations which complicates the solution.In this paper,it is shown that the major contribution of viscosity and the horizontal component of actuation is concentrated in a narrow region alongside the actuation boundary.Since the inviscid equations are considerably easier to solve,a division into the viscous and inviscid domains would alleviate the computational load significantly.The particles'traces calculated by this approximation are excellently alongside their counterparts from the completely viscous model.It is also shown that the optimum thickness for the viscous strip is about 9-fold the acoustic boundary layer thickness for various flow patterns and amplitudes of actuation.

Scheme of negative acoustic radiation force based on a multiple-layered spherical structure

Acoustic radiation force(ARF), as an important particle manipulation method, has been extensively studied in recent years. With the introduction of the concept of “acoustic tweezers”, negative acoustic radiation has become a research hotspot. In this paper, a scheme of realizing negative ARF based on the multiple-layered spherical structure design is proposed. The specific structure and design idea are presented. Detailed theoretical calculation analysis is carried out.Numerical simulations have been performed to verify the correctness of this prediction. The conjecture that the suppression of backscattering can achieve negative ARF is verified concretely, which greatly expands the application prospect and design ideas of the ARF. This work has laid a theoretical foundation for realizing precise control of the structure.