The W transform and its improved methods for time-frequency analysis of seismic data

The conventional linear time-frequency analysis method cannot achieve high resolution and energy focusing in the time and frequency dimensions at the same time,especially in the low frequency region.In order to improve the resolution of the linear time-frequency analysis method in the low-frequency region,we have proposed a W transform method,in which the instantaneous frequency is introduced as a parameter into the linear transformation,and the analysis time window is constructed which matches the instantaneous frequency of the seismic data.In this paper,the W transform method is compared with the Wigner-Ville distribution(WVD),a typical nonlinear time-frequency analysis method.The WVD method that shows the energy distribution in the time-frequency domain clearly indicates the gravitational center of time and the gravitational center of frequency of a wavelet,while the time-frequency spectrum of the W transform also has a clear gravitational center of energy focusing,because the instantaneous frequency corresponding to any time position is introduced as the transformation parameter.Therefore,the W transform can be benchmarked directly by the WVD method.We summarize the development of the W transform and three improved methods in recent years,and elaborate on the evolution of the standard W transform,the chirp-modulated W transform,the fractional-order W transform,and the linear canonical W transform.Through three application examples of W transform in fluvial sand body identification and reservoir prediction,it is verified that W transform can improve the resolution and energy focusing of time-frequency spectra.

Novel Doppler Frequency Extraction Method Based on Time-Frequency Analysis and Morphological Operation

A novel method of Doppler frequency extraction is proposed for Doppler radar scoring systems. The idea is that the time-frequency map can show how the Doppler frequency varies along the time-line, so the Doppler frequency extraction becomes curve detection in the image-view. A set of morphological operations are used to implement curve detection. And a map fusion scheme is presented to eliminate the influence of strong direct current （DC） component of echo signal during curve detection. The radar real-life data are used to illustrate the performance of the new approach. Experimental results show that the proposed method can overcome the shortcomings of piecewise-processing-based FFT method and can improve the measuring precision of miss distance.

A Recursive Method of Time-Frequency Analysis for the Signal Processing of Flutter Test with Progression Variable Speed

Focused on the non-statlonarity and real-time analysis of signal in flutter test with progression variable speed （FTPVS）, a new method of recursive time-frequency analysis is presented. The time-varying system is tracked on-line by building a time-varying parameter model, and then the relevant parameter spectrum can be obtained. The feasibility and advantages of the method are examined by digital simulation. The results of FTPVS at low-speed wind-tunnel promise the engineering application perspective of the method.

NEW SAR IMAGE INTERPRETATION METHOD OF AIRCRAFT BASED ON JOINT TIME-FREQUENCY ANALYSIS

With the continuous improvement of Synthetic Aperture Radar(SAR) resolution, interpreting the small targets like aircraft in SAR images becomes possible and turn out to be a hot spot in SAR application research. However, due to the complexity of SAR imaging mechanism, interpreting targets in SAR images is a tough problem. This paper presents a new aircraft interpretation method based on the joint time-frequency analysis and multi-dimensional contrasting of basic structures. Moreover, SAR data acquisition experiment is designed for interpreting the aircraft. Analyzing the experiment data with our method, the result shows that the proposed method largely makes use of the SAR data information. The reasonable results can provide some auxiliary support for the SAR images manual interpretation.

Application of Wavelet Packet De-noising in Time-Frequency Analysis of the Local Wave Method

The local wave method is a very good time-frequency method for nonstationaryvibration signal analysis. But the interfering noise has a big influence on the accuracy oftime-frequency analysis. The wavelet packet de-noising method can eliminate the interference ofnoise and improve the signal-noise-ratio. This paper uses the local wave method to decompose thede-noising signal and perform a time-frequency analysis. We can get better characteristics. Finally,an example of wavelet packet de-noising and a local wave time-frequency spectrum application ofdiesel engine surface vibration signal is put forward.

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.

Generalized nth-Order Perturbation Method Based on Loop Subdivision Surface Boundary Element Method for Three-Dimensional Broadband Structural Acoustic Uncertainty Analysis

In this paper,a generalized nth-order perturbation method based on the isogeometric boundary element method is proposed for the uncertainty analysis of broadband structural acoustic scattering problems.The Burton-Miller method is employed to solve the problem of non-unique solutions that may be encountered in the external acoustic field,and the nth-order discretization formulation of the boundary integral equation is derived.In addition,the computation of loop subdivision surfaces and the subdivision rules are introduced.In order to confirm the effectiveness of the algorithm,the computed results are contrasted and analyzed with the results under Monte Carlo simulations(MCs)through several numerical examples.

Dynamic Coupling Analysis of Semisubmersible Platform Float-over Method for Docking Case

In this paper,the multi-body coupled dynamic characteristics of a semisubmersible platform and an HYSY 229 barge were investigated.First,coupled hydrodynamic analysis of the HYSY 229 barge and the semisubmersible platform was performed.Relevant hydrodynamic parameters were obtained using the retardation function method of three-dimensional frequency-domain potential flow theory.The results of the hydrodynamic analysis were highly consistent with the test findings,verifying the accuracy of the multifloating hydrodynamic coupling analysis,and key hydrodynamic parameters were solved for different water depths and the coupling effect.According to the obtained results,the hydrodynamic influence was the largest in shallow waters when the coupling effect was considered.Furthermore,the coupled motion equation combined with viscous damping,fender system,and mooring system was established,and the hydrodynamics,floating body motion,and dynamic response of the fender system were analyzed.Motion analysis revealed good agreement among the surge,sway,and yaw motions of the two floating bodies.However,when the wave period reached 10 s,the motion of the two floating bodies showed severe shock,and a relative motion was also observed.Therefore,excessive constraints should be added between the two floating bodies during construction to ensure construction safety.The numerical analysis and model test results of the semisubmersible platform and HYSY 229 barge at a water depth of 42 m and sea conditions of 0°,45°,and 90° were in good agreement,and the error was less than 5%.The maximum movement of the HYSY 229 barge reached 2.61 m in the sway direction,whereas that of the semisubmersible platform was 2.11 m.During construction,excessive constraints should be added between the two floating bodies to limit their relative movement and ensure construction safety.

Simplified quantitative analysis method and its application in the insitu synthesized copper-based azide chips

Copper-based azide(Cu(N_(3))2 or CuN_(3),CA)chips synthesized by in-situ azide reaction and utilized in miniaturized explosive systems has become a hot research topic in recent years.However,the advantages of in-situ synthesis method,including small size and low dosage,bring about difficulties in quantitative analysis and differences in ignition capabilities of CA chips.The aim of present work is to develop a simplified quantitative analysis method for accurate and safe analysis of components in CA chips to evaluate and investigate the corresponding ignition ability.In this work,Cu(N_(3))2 and CuN_(3)components in CA chips were separated through dissolution and distillation by utilizing the difference in solubility and corresponding content was obtained by measuring N_(3)-concentration through spectrophotometry.The spectrophotometry method was optimized by studying influencing factors and the recovery rate of different separation methods was studied,ensuring the accuracy and reproducibility of test results.The optimized method is linear in range from 1.0-25.0 mg/L,with a correlation coefficient R^(2)=0.9998,which meets the requirements of CA chips with a milligram-level content test.Compared with the existing ICP method,component analysis results of CA chips obtained by spectrophotometry are closer to real component content in samples and have satisfactory accuracy.Moreover,as its application in miniaturized explosive systems,the ignition ability of CA chips with different component contents for direct ink writing CL-20 and the corresponding mechanism was studied.This study provided a basis and idea for the design and performance evaluation of CA chips in miniaturized explosive systems.

Reliability analysis of carbon fiber rod-reinforced umbilical cable under tension using an improved sampling method

The umbilical cable is a vital component of subsea production systems that provide power,chemical agents,control signals et al.,and its requirement for reliability is exceedingly high.However,as the umbilical cable is a composite structure comprising multiple functional units,the reliability analysis of such cables involves numerous parameters that can impact calculation efficiency.In this paper,the reliability analysis of a new kind of umbilical cable with carbon fiber rod under tension is analyzed.The global dynamic analytical model is first established to determine the maximum tension load,then the local analytical model of umbilical cable including each unit are constructed by finite element method(FEM).Based on the mechanical analytical model,the reliability of umbilical cable under tension load is studied using response surface method(RSM)and Monte Carlo method.During the calculation process,a new tangent plane sampling method to calculate the response surface function(RSF)is proposed in this paper,which could make sampling points faster come close to the RSF curve,and it is proved that the calculation efficiency increases about 33%comparing with traditional method.

AHermitian C^(2) Differential Reproducing Kernel Interpolation Meshless Method for the 3D Microstructure-Dependent Static Flexural Analysis of Simply Supported and Functionally Graded Microplates

This work develops a Hermitian C^(2) differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredependent static flexural behavior of a functionally graded(FG)microplate subjected to mechanical loads and placed under full simple supports.In the formulation,we select the transverse stress and displacement components and their first-and second-order derivatives as primary variables.Then,we set up the differential reproducing conditions(DRCs)to obtain the shape functions of the Hermitian C^(2) differential reproducing kernel(DRK)interpolant’s derivatives without using direct differentiation.The interpolant’s shape function is combined with a primitive function that possesses Kronecker delta properties and an enrichment function that constituents DRCs.As a result,the primary variables and their first-and second-order derivatives satisfy the nodal interpolation properties.Subsequently,incorporating ourHermitianC^(2)DRKinterpolant intothe strong formof the3DCCST,we develop a DRKIM method to analyze the FG microplate’s 3D microstructure-dependent static flexural behavior.The Hermitian C^(2) DRKIM method is confirmed to be accurate and fast in its convergence rate by comparing the solutions it produces with the relevant 3D solutions available in the literature.Finally,the impact of essential factors on the transverse stresses,in-plane stresses,displacements,and couple stresses that are induced in the loaded microplate is examined.These factors include the length-to-thickness ratio,the material length-scale parameter,and the inhomogeneity index,which appear to be significant.

Superconvergence of Direct Discontinuous Galerkin Methods:Eigen-structure Analysis Based on Fourier Approach

This paper investigates superconvergence properties of the direct discontinuous Galerkin(DDG)method with interface corrections and the symmetric DDG method for diffusion equations.We apply the Fourier analysis technique to symbolically compute eigenvalues and eigenvectors of the amplification matrices for both DDG methods with different coefficient settings in the numerical fluxes.Based on the eigen-structure analysis,we carry out error estimates of the DDG solutions,which can be decomposed into three parts:(i)dissipation errors of the physically relevant eigenvalue,which grow linearly with the time and are of order 2k for P^(k)(k=2,3)approximations;(ii)projection error from a special projection of the exact solution,which is decreasing over the time and is related to the eigenvector corresponding to the physically relevant eigenvalue;(iii)dissipative errors of non-physically relevant eigenvalues,which decay exponentially with respect to the spatial mesh sizeΔx.We observe that the errors are sensitive to the choice of the numerical flux coefficient for even degree P^(2)approximations,but are not for odd degree P^(3)approximations.Numerical experiments are provided to verify the theoretical results.

A novel box-counting method for quantitative fractal analysis of threedimensional pore characteristics in sandstone

Fractal theory offers a powerful tool for the precise description and quantification of the complex pore structures in reservoir rocks,crucial for understanding the storage and migration characteristics of media within these rocks.Faced with the challenge of calculating the three-dimensional fractal dimensions of rock porosity,this study proposes an innovative computational process that directly calculates the three-dimensional fractal dimensions from a geometric perspective.By employing a composite denoising approach that integrates Fourier transform(FT)and wavelet transform(WT),coupled with multimodal pore extraction techniques such as threshold segmentation,top-hat transformation,and membrane enhancement,we successfully crafted accurate digital rock models.The improved box-counting method was then applied to analyze the voxel data of these digital rocks,accurately calculating the fractal dimensions of the rock pore distribution.Further numerical simulations of permeability experiments were conducted to explore the physical correlations between the rock pore fractal dimensions,porosity,and absolute permeability.The results reveal that rocks with higher fractal dimensions exhibit more complex pore connectivity pathways and a wider,more uneven pore distribution,suggesting that the ideal rock samples should possess lower fractal dimensions and higher effective porosity rates to achieve optimal fluid transmission properties.The methodology and conclusions of this study provide new tools and insights for the quantitative analysis of complex pores in rocks and contribute to the exploration of the fractal transport properties of media within rocks.

Numerical Analysis of Bacterial Meningitis Stochastic Delayed Epidemic Model through Computational Methods

Based on theWorld Health Organization(WHO),Meningitis is a severe infection of the meninges,the membranes covering the brain and spinal cord.It is a devastating disease and remains a significant public health challenge.This study investigates a bacterial meningitis model through deterministic and stochastic versions.Four-compartment population dynamics explain the concept,particularly the susceptible population,carrier,infected,and recovered.The model predicts the nonnegative equilibrium points and reproduction number,i.e.,the Meningitis-Free Equilibrium(MFE),and Meningitis-Existing Equilibrium(MEE).For the stochastic version of the existing deterministicmodel,the twomethodologies studied are transition probabilities and non-parametric perturbations.Also,positivity,boundedness,extinction,and disease persistence are studiedrigorouslywiththe helpofwell-known theorems.Standard and nonstandard techniques such as EulerMaruyama,stochastic Euler,stochastic Runge Kutta,and stochastic nonstandard finite difference in the sense of delay have been presented for computational analysis of the stochastic model.Unfortunately,standard methods fail to restore the biological properties of the model,so the stochastic nonstandard finite difference approximation is offered as an efficient,low-cost,and independent of time step size.In addition,the convergence,local,and global stability around the equilibria of the nonstandard computational method is studied by assuming the perturbation effect is zero.The simulations and comparison of the methods are presented to support the theoretical results and for the best visualization of results.

A vector sum analysis method for stability evolution of expansive soil slope considering shear zone damage softening

Slope stability analysis is a classical mechanical problem in geotechnical engineering and engineering geology.It is of great significance to study the stability evolution of expansive soil slopes for engineering construction in expansive soil areas.Most of the existing studies evaluate the slope stability by analyzing the limit equilibrium state of the slope,and the analysis method for the stability evolution considering the damage softening of the shear zone is lacking.In this study,the large deformation shear mechanical behavior of expansive soil was investigated by ring shear test.The damage softening characteristic of expansive soil in the shear zone was analyzed,and a shear damage model reflecting the damage softening behavior of expansive soil was derived based on the damage theory.Finally,by skillfully combining the vector sum method and the shear damage model,an analysis method for the stability evolution of the expansive soil slope considering the shear zone damage softening was proposed.The results show that the shear zone subjected to large displacement shear deformation exhibits an obvious damage softening phenomenon.The damage variable equation based on the logistic function can be well used to describe the shear damage characteristics of expansive soil,and the proposed shear damage model is in good agreement with the ring shear test results.The vector sum method considering the damage softening behavior of the shear zone can be well applied to analyze the stability evolution characteristics of the expansive soil slope.The stability factor of the expansive soil slope decreases with the increase of shear displacement,showing an obvious progressive failure behavior.

Application of Choquet Integral-Importance-Performance Analysis and TOPSIS Methods in Approaching the Preference Factors of Calligraphy and Seal Engraving Imagery

Classical Chinese characters,presented through calligraphy,seal engraving,or painting,can exhibit different aesthetics and essences of Chinese characters,making them the most important asset of the Chinese people.Calligraphy and seal engraving,as two closely related systems in traditional Chinese art,have developed through the ages.Due to changes in lifestyle and advancements in modern technology,their original functions of daily writing and verification have gradually diminished.Instead,they have increasingly played a significant role in commercial art.This study utilizes the Evaluation Grid Method(EGM)and the Analytic Hierarchy Process(AHP)to research the key preference factors in the application of calligraphy and seal engraving imagery.Different from the traditional 5-point equal interval semantic questionnaire,this study employs a non-equal interval semantic questionnaire with a golden ratio scale,distinguishing the importance ratio of adjacent semantic meanings and highlighting the weighted emphasis on visual aesthetics.Additionally,the study uses Importance-Performance Analysis(IPA)and Technique for Order Preference by Similarity to Ideal Solution(TOPSIS)to obtain the key preference sequence of calligraphy and seal engraving culture.Plus,the Choquet integral comprehensive evaluation is used as a reference for IPA comparison.It is hoped that this study can provide cultural imagery references and research methods,injecting further creativity into industrial design.

Time-frequency analysis of Li solid-phase diffusion in spherical active particles under typical discharge modes

Li transient concentration distribution in spherical active material particles can affect the maximum power density and the safe operating regime of the electric vehicles(EVs). On one hand, the quasiexact/exact solution obtained in the time/frequency domain is time-consuming and just as a reference value for approximate solutions;on the other hand, calculation errors and application range of approximate solutions not only rely on approximate algorithms but also on discharge modes. For the purpose to track the transient dynamics for Li solid-phase diffusion in spherical active particles with a tolerable error range and for a wide applicable range, it is necessary to choose optimal approximate algorithms in terms of discharge modes and the nature of active material particles. In this study, approximation methods,such as diffusion length method, polynomial profile approximation method, Padé approximation method,pseudo steady state method, eigenfunction-based Galerkin collocation method, and separation of variables method for solving Li solid-phase diffusion in spherical active particles are compared from calculation fundamentals to algorithm implementation. Furthermore, these approximate solutions are quantitatively compared to the quasi-exact/exact solution in the time/frequency domain under typical discharge modes, i.e., start-up, slow-down, and speed-up. The results obtained from the viewpoint of time-frequency analysis offer a theoretical foundation on how to track Li transient concentration profile in spherical active particles with a high precision and for a wide application range. In turn, optimal solutions of Li solid diffusion equations for spherical active particles can improve the reliability in predicting safe operating regime and estimating maximum power for automotive batteries.

Parameterized time-frequency analysis to separate multi-radar signals

Multi-radar signal separation is a critical process in modern reconnaissance systems. However, the complicated battlefield is typically confronted with increasing electronic equipment and complex radar waveforms. The intercepted signal is difficult to separate with conventional parameters because of severe overlapping in both time and frequency domains. On the contrary, time-frequency analysis maps the 1D signal into a 2D time-frequency plane, which provides a better insight into the signal than traditional methods. Particularly, the parameterized time-frequency analysis (PTFA) shows great potential in processing such non stationary signals. Five procedures for the PTFA are proposed to separate the overlapped multi-radar signal, including initiation, instantaneous frequency estimation with PTFA, signal demodulation, signal separation with adaptive filter and signal recovery. The proposed method is verified with both simulated and real signals, which shows good performance in the application on multi-radar signal separation.

Time-Frequency Analysis of EEG Signals Evoked by Voluntary, Stimulated and Imaginary Motions

In order to investigate the characteristics of sensorimotor cortex during motor execution(ME), voluntary, stimulated and imaginary finger flexions were performed by ten volunteer subjects. Electroencephalogram(EEG) data were recorded according to the modified 10-20 International EEG System. The patterns were compared by the analysis of the motion-evoked EEG signals focusing on the contralateral(C3) and ipsilateral(C4) channels for hemispheric differences. The EEG energy distributions at alpha(8—13 Hz), beta(14—30 Hz) and gamma(30—50 Hz) bands were computed by wavelet transform(WT) and compared by the analysis of variance(ANOVA). The timefrequency(TF) analysis indicated that there existed a contralateral dominance of alpha post-movement event-related synchronization(ERS) pattern during the voluntary task, and that the energy of alpha band increased in the ipsilateral area during the stimulated(median nerve of wrist) task. Besides, the contralateral alpha and beta event-related desynchronization(ERD) patterns were observed in both stimulated and imaginary tasks. Another significant difference was found in the mean power values of gamma band(p<0.01)between the imaginary and other tasks. The results show that significant hemispheric differences such as alpha and beta band EEG energy distributions and TF changing phenomena(ERS/ERD) were found between C3 and C4 areas during all of the three patterns. The largest energy distribution was always at the alpha band for each task.

The combined time-frequency analysis of the acoustic signals backscattered from ultrasonic contrast agents in the evaluation of the blood perfusion

Objective: To analyze the non-periodic, unstable and even chaotic echoes scattered from microbubbles which are extremely sensitive and may easily collapse, fragment or shrink when ultrasound contrast agents are exposed to ultrasound (US) irradiation. Methods: The combined time-frequency analysis was applied to the original signals instead of the traditional Fourier spectral analysis technique. Results: The results obtained from simulation as well as experiment showed that the subharmonic, 2nd harmonic and ultra harmonic of the microbubbles occurred during the oscillation and varied with time. The dependence on the incident ultrasonic amplitude and microbubble parameters were established. Conclusion: The transient echoes backscattered from the ultrasound agent in the evaluation of the blood perfusion can be analyzed thoroughly by the technique of combined-frequency analysis and the time detail of the frequency contents can be revealed.