The rapidly advancing Convolutional Neural Networks(CNNs)have brought about a paradigm shift in various computer vision tasks,while also garnering increasing interest and application in sensor-based Human Activity Rec...The rapidly advancing Convolutional Neural Networks(CNNs)have brought about a paradigm shift in various computer vision tasks,while also garnering increasing interest and application in sensor-based Human Activity Recognition(HAR)efforts.However,the significant computational demands and memory requirements hinder the practical deployment of deep networks in resource-constrained systems.This paper introduces a novel network pruning method based on the energy spectral density of data in the frequency domain,which reduces the model’s depth and accelerates activity inference.Unlike traditional pruning methods that focus on the spatial domain and the importance of filters,this method converts sensor data,such as HAR data,to the frequency domain for analysis.It emphasizes the low-frequency components by calculating their energy spectral density values.Subsequently,filters that meet the predefined thresholds are retained,and redundant filters are removed,leading to a significant reduction in model size without compromising performance or incurring additional computational costs.Notably,the proposed algorithm’s effectiveness is empirically validated on a standard five-layer CNNs backbone architecture.The computational feasibility and data sensitivity of the proposed scheme are thoroughly examined.Impressively,the classification accuracy on three benchmark HAR datasets UCI-HAR,WISDM,and PAMAP2 reaches 96.20%,98.40%,and 92.38%,respectively.Concurrently,our strategy achieves a reduction in Floating Point Operations(FLOPs)by 90.73%,93.70%,and 90.74%,respectively,along with a corresponding decrease in memory consumption by 90.53%,93.43%,and 90.05%.展开更多
This paper presents a new technique for measuring the bunch length of a high-energy electron beam at a bunch-by-bunch rate in storage rings.This technique uses the time–frequency-domain joint analysis of the bunch si...This paper presents a new technique for measuring the bunch length of a high-energy electron beam at a bunch-by-bunch rate in storage rings.This technique uses the time–frequency-domain joint analysis of the bunch signal to obtain bunch-by-bunch and turn-by-turn longitudinal parameters,such as bunch length and synchronous phase.The bunch signal is obtained using a button electrode with a bandwidth of several gigahertz.The data acquisition device was a high-speed digital oscilloscope with a sampling rate of more than 10 GS/s,and the single-shot sampling data buffer covered thousands of turns.The bunch-length and synchronous phase information were extracted via offline calculations using Python scripts.The calibration coefficient of the system was determined using a commercial streak camera.Moreover,this technique was tested on two different storage rings and successfully captured various longitudinal transient processes during the harmonic cavity debugging process at the Shanghai Synchrotron Radiation Facility(SSRF),and longitudinal instabilities were observed during the single-bunch accumulation process at Hefei Light Source(HLS).For Gaussian-distribution bunches,the uncertainty of the bunch phase obtained using this technique was better than 0.2 ps,and the bunch-length uncertainty was better than 1 ps.The dynamic range exceeded 10 ms.This technology is a powerful and versatile beam diagnostic tool that can be conveniently deployed in high-energy electron storage rings.展开更多
The traditional calculation method of frequency-domain Green function mainly utilizes series or asymptotic expansion to carry out numerical approximation, however, this method requires very careful zoning, thus the co...The traditional calculation method of frequency-domain Green function mainly utilizes series or asymptotic expansion to carry out numerical approximation, however, this method requires very careful zoning, thus the computing process is complex with many cycles, which has greatly affected the computing efficiency. To improve the computing efficiency, this paper introduces Gaussian integral to the numerical calculation of the frequency-domain Green function and its partial derivatives. It then compares the calculation result with that in existing references. The comparison results demonstrate that, on the basis of its sufficient accuracy, the method has greatly simplified the computing process, reduced the zoning and improved the computing efficiency.展开更多
This compendium review focuses on the spatial distribution of sensitivity to localized absorption changes in optically diffuse media,particularly for measurements relevant to near-infrared spectroscopy.The three tempo...This compendium review focuses on the spatial distribution of sensitivity to localized absorption changes in optically diffuse media,particularly for measurements relevant to near-infrared spectroscopy.The three temporal domains,continuous wave,frequency domain,and time domain,each obtain different optical data types whose changes may be related to effective homogeneous changes in the absorption coefficient.Sensitivity is the relationship between a localized perturbation and the recovered effective homogeneous absorption change.Therefore,spatial sensitivity maps representing the perturbation location can be generated for the numerous optical data types in the three temporal domains.The review first presents a history of the past 30 years of work investigating this sensitivity in optically diffuse media.These works are experimental and theoretical,presenting one-,two-,and three-dimensional sensitivity maps for different Near-Infrared Spectroscopy methods,domains,and data types.Following this history,we present a compendium of sensitivity maps organized by temporal domain and then data type.This compendium provides a valuable tool to compare the spatial sensitivity of various measurement methods and parameters in one document.Methods for one to generate these maps are provided in Appendix A,including the code.This historical review and comprehensive sensitivity map compendium provides a single source researchers may use to visualize,investigate,compare,and generate sensitivity to localized absorption change maps.展开更多
The complexity of an elastic wavefield increases the nonlinearity of inversion, To some extent, multiscale inversion decreases the nonlinearity of inversion and prevents it from falling into local extremes. A multisca...The complexity of an elastic wavefield increases the nonlinearity of inversion, To some extent, multiscale inversion decreases the nonlinearity of inversion and prevents it from falling into local extremes. A multiscale strategy based on the simultaneous use of frequency groups and layer stripping method based on damped wave field improves the stability of inversion. A dual-level parallel algorithm is then used to decrease the computational cost and improve practicability. The seismic wave modeling of a single frequency and inversion in a frequency group are computed in parallel by multiple nodes based on multifrontal massively parallel sparse direct solver and MPI. Numerical tests using an overthrust model show that the proposed inversion algorithm can effectively improve the stability and accuracy of inversion by selecting the appropriate inversion frequency and damping factor in low- frequency seismic data.展开更多
Frequency-domain airborne electromagnetics is a proven geophysical exploration method.Presently,the interpretation is mainly based on resistivity-depth imaging and onedimensional layered inversion;nevertheless,it is d...Frequency-domain airborne electromagnetics is a proven geophysical exploration method.Presently,the interpretation is mainly based on resistivity-depth imaging and onedimensional layered inversion;nevertheless,it is difficult to obtain satisfactory results for two- or three-dimensional complex earth structures using 1D methods.3D forward modeling and inversion can be used but are hampered by computational limitations because of the large number of data.Thus,we developed a 2.5D frequency-domain airborne electromagnetic forward modeling and inversion algorithm.To eliminate the source singularities in the numerical simulations,we split the fields into primary and secondary fields.The primary fields are calculated using homogeneous or layered models with analytical solutions,and the secondary(scattered) fields are solved by the finite-element method.The linear system of equations is solved by using the large-scale sparse matrix parallel direct solver,which greatly improves the computational efficiency.The inversion algorithm was based on damping leastsquares and singular value decomposition and combined the pseudo forward modeling and reciprocity principle to compute the Jacobian matrix.Synthetic and field data were used to test the effectiveness of the proposed method.展开更多
As a high quality seismic imaging method, full waveform inversion (FWI) can accurately reconstruct the physical parameter model for the subsurface medium. However, application of the FWI in seismic data processing i...As a high quality seismic imaging method, full waveform inversion (FWI) can accurately reconstruct the physical parameter model for the subsurface medium. However, application of the FWI in seismic data processing is computationally expensive, especially for the three-dimension complex medium inversion. Introducing blended source technology into the frequency-domain FWI can greatly reduce the computational burden and improve the efficiency of the inversion. However, this method has two issues: first, crosstalk noise is caused by interference between the sources involved in the encoding, resulting in an inversion result with some artifacts; second, it is more sensitive to ambient noise compared to conventional FWI, therefore noisy data results in a poor inversion. This paper introduces a frequency-group encoding method to suppress crosstalk noise, and presents a frequency- domain auto-adapting FWI based on source-encoding technology. The conventional FWI method and source-encoding based FWI method are combined using an auto-adapting mechanism. This improvement can both guarantee the quality of the inversion result and maximize the inversion efficiency.展开更多
For real-time jamming signal generation in deceiving inverse synthetic aperture radar(ISAR),the target characteristics modulation is always processed in the expensive field programmable gate array(FPGA).Due to the...For real-time jamming signal generation in deceiving inverse synthetic aperture radar(ISAR),the target characteristics modulation is always processed in the expensive field programmable gate array(FPGA).Due to the large computational complexity of the traditional modulating operation,the size and structure of simulated false-target are limited.With regard to the principle of dechirping in range compression of linear frequency modulated(LFM) radar,a novel algorithm named "inverse dechirping" is proposed for target characteristics modulation.This algorithm only needs one complex multiplier in the FPGA to generate the jamming signal when the radar signal is intercepted,which can be obtained by multiplication of radar signal samplings and the equivalent dechirped target echo in the time domain.As the complex synthesis of dechirped target echo can be realized by cheap digital signal processor(DSP) within the interpulse time,the overall cost of the jamming equipment will be reduced and the false-target size will not be limited by the scale of FPGA.Numerical simulations are performed to verify the correctness and effectiveness of the proposed algorithm.展开更多
The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the tim...The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.展开更多
AIM To establish whether frequency-domain optical coherence tomography(FD-OCT) is safe and effective in the evaluation and treatment of angiographicallyintermediate coronary lesions(ICL) METHODS Sixty-four patients wi...AIM To establish whether frequency-domain optical coherence tomography(FD-OCT) is safe and effective in the evaluation and treatment of angiographicallyintermediate coronary lesions(ICL) METHODS Sixty-four patients with 2-dimensional quantitativecoronary angiography(2D-QCA) demonstrating ICL were included. OCT imaging was performed. According to predetermined OCT criteria, patients were assigned to either of 2 groups: OCT-guided percutaneous coronary intervention(PCI) or OCT-guided optimal medical therapy(OMT). The primary efficacy endpoint was to demonstrate the superiority and higher accuracy of FD-OCT compared to 2D-QCA in evaluating stenosis severity in patients with ICL. The primary safety endpoint was the incidence of 30-d major adverse cardiac events(MACE). Secondary endpoints included MACE at 12 mo and other clinical events.RESULTS Analysis of the primary efficacy endpoint demonstrates that 2D-QCA overestimates the stenosis severity of ICL in both the OCT-guided PCI and OMT groups, proving FD-OCT to be superior to and more precise than 2D-QCA in treating this subset of lesions. The primary safety endpoint was fully met with the incidence of 30-d MACE being nil in both the OCT-guided PCI and OCTguided OMT groups. Incidences of secondary endpoints were found to be low in both arms, the only exception being the relatively high incidence of recurrent episodes of angina which was, however, very similar in the 2 groups.CONCLUSION FD-OCT is safe and effective in the evaluation and treatment of ICL. Larger studies are needed to firmly establish the efficacy and safety of FD-OCT in treating ICL across all coronary artery disease population subgroups.展开更多
Earth medium is not completely elastic, with its viscosity resulting in attenuation and dispersion of seismic waves. Most viscoelastic numerical simulations are based on the finite-difference and finite-element method...Earth medium is not completely elastic, with its viscosity resulting in attenuation and dispersion of seismic waves. Most viscoelastic numerical simulations are based on the finite-difference and finite-element methods. Targeted at viscoelastic numerical modeling for multilayered media, the constant-Q acoustic wave equation is transformed into the corresponding wave integral representation with its Green's function accounting for viscoelastic coefficients. An efficient alternative for full-waveform solution to the integral equation is proposed in this article by extending conventional frequency-domain boundary element methods to viscoelastic media. The viscoelastic boundary element method enjoys a distinct characteristic of the explicit use of boundary continuity conditions of displacement and traction, leading to a semi-analytical solution with sufficient accuracy for simulating the viscoelastic effect across irregular interfaces. Numerical experiments to study the viscoelastic absorption of different Q values demonstrate the accuracy and applicability of the method.展开更多
基金supported by National Natural Science Foundation of China(Nos.61902158 and 62202210).
文摘The rapidly advancing Convolutional Neural Networks(CNNs)have brought about a paradigm shift in various computer vision tasks,while also garnering increasing interest and application in sensor-based Human Activity Recognition(HAR)efforts.However,the significant computational demands and memory requirements hinder the practical deployment of deep networks in resource-constrained systems.This paper introduces a novel network pruning method based on the energy spectral density of data in the frequency domain,which reduces the model’s depth and accelerates activity inference.Unlike traditional pruning methods that focus on the spatial domain and the importance of filters,this method converts sensor data,such as HAR data,to the frequency domain for analysis.It emphasizes the low-frequency components by calculating their energy spectral density values.Subsequently,filters that meet the predefined thresholds are retained,and redundant filters are removed,leading to a significant reduction in model size without compromising performance or incurring additional computational costs.Notably,the proposed algorithm’s effectiveness is empirically validated on a standard five-layer CNNs backbone architecture.The computational feasibility and data sensitivity of the proposed scheme are thoroughly examined.Impressively,the classification accuracy on three benchmark HAR datasets UCI-HAR,WISDM,and PAMAP2 reaches 96.20%,98.40%,and 92.38%,respectively.Concurrently,our strategy achieves a reduction in Floating Point Operations(FLOPs)by 90.73%,93.70%,and 90.74%,respectively,along with a corresponding decrease in memory consumption by 90.53%,93.43%,and 90.05%.
基金supported by the National Key R&D Program(No.2022YFA1602201)。
文摘This paper presents a new technique for measuring the bunch length of a high-energy electron beam at a bunch-by-bunch rate in storage rings.This technique uses the time–frequency-domain joint analysis of the bunch signal to obtain bunch-by-bunch and turn-by-turn longitudinal parameters,such as bunch length and synchronous phase.The bunch signal is obtained using a button electrode with a bandwidth of several gigahertz.The data acquisition device was a high-speed digital oscilloscope with a sampling rate of more than 10 GS/s,and the single-shot sampling data buffer covered thousands of turns.The bunch-length and synchronous phase information were extracted via offline calculations using Python scripts.The calibration coefficient of the system was determined using a commercial streak camera.Moreover,this technique was tested on two different storage rings and successfully captured various longitudinal transient processes during the harmonic cavity debugging process at the Shanghai Synchrotron Radiation Facility(SSRF),and longitudinal instabilities were observed during the single-bunch accumulation process at Hefei Light Source(HLS).For Gaussian-distribution bunches,the uncertainty of the bunch phase obtained using this technique was better than 0.2 ps,and the bunch-length uncertainty was better than 1 ps.The dynamic range exceeded 10 ms.This technology is a powerful and versatile beam diagnostic tool that can be conveniently deployed in high-energy electron storage rings.
基金Supported by the National Natural Science Foundation of China under Grant No.50779007the National Science Foundation for Young Scientists of China under Grant No.50809018+2 种基金the Specialized Research Fund for the Doctoral Program of Higher Education of China under Grant No.20070217074the Defence Advance Research Program of Science and Technology of Ship Industry under Grant No.07J1.1.6Harbin Engineering University Foundation under Grant No.HEUFT07069
文摘The traditional calculation method of frequency-domain Green function mainly utilizes series or asymptotic expansion to carry out numerical approximation, however, this method requires very careful zoning, thus the computing process is complex with many cycles, which has greatly affected the computing efficiency. To improve the computing efficiency, this paper introduces Gaussian integral to the numerical calculation of the frequency-domain Green function and its partial derivatives. It then compares the calculation result with that in existing references. The comparison results demonstrate that, on the basis of its sufficient accuracy, the method has greatly simplified the computing process, reduced the zoning and improved the computing efficiency.
文摘This compendium review focuses on the spatial distribution of sensitivity to localized absorption changes in optically diffuse media,particularly for measurements relevant to near-infrared spectroscopy.The three temporal domains,continuous wave,frequency domain,and time domain,each obtain different optical data types whose changes may be related to effective homogeneous changes in the absorption coefficient.Sensitivity is the relationship between a localized perturbation and the recovered effective homogeneous absorption change.Therefore,spatial sensitivity maps representing the perturbation location can be generated for the numerous optical data types in the three temporal domains.The review first presents a history of the past 30 years of work investigating this sensitivity in optically diffuse media.These works are experimental and theoretical,presenting one-,two-,and three-dimensional sensitivity maps for different Near-Infrared Spectroscopy methods,domains,and data types.Following this history,we present a compendium of sensitivity maps organized by temporal domain and then data type.This compendium provides a valuable tool to compare the spatial sensitivity of various measurement methods and parameters in one document.Methods for one to generate these maps are provided in Appendix A,including the code.This historical review and comprehensive sensitivity map compendium provides a single source researchers may use to visualize,investigate,compare,and generate sensitivity to localized absorption change maps.
基金supported by the Natural Science Foundation of China(No.41374122)
文摘The complexity of an elastic wavefield increases the nonlinearity of inversion, To some extent, multiscale inversion decreases the nonlinearity of inversion and prevents it from falling into local extremes. A multiscale strategy based on the simultaneous use of frequency groups and layer stripping method based on damped wave field improves the stability of inversion. A dual-level parallel algorithm is then used to decrease the computational cost and improve practicability. The seismic wave modeling of a single frequency and inversion in a frequency group are computed in parallel by multiple nodes based on multifrontal massively parallel sparse direct solver and MPI. Numerical tests using an overthrust model show that the proposed inversion algorithm can effectively improve the stability and accuracy of inversion by selecting the appropriate inversion frequency and damping factor in low- frequency seismic data.
基金supported by the Doctoral Fund Project of the Ministry of Education(No.20130061110060 class tutors)the National Natural Science Foundation of China(No.41504083)National Basic Research Program of China(973Program)(No.2013CB429805)
文摘Frequency-domain airborne electromagnetics is a proven geophysical exploration method.Presently,the interpretation is mainly based on resistivity-depth imaging and onedimensional layered inversion;nevertheless,it is difficult to obtain satisfactory results for two- or three-dimensional complex earth structures using 1D methods.3D forward modeling and inversion can be used but are hampered by computational limitations because of the large number of data.Thus,we developed a 2.5D frequency-domain airborne electromagnetic forward modeling and inversion algorithm.To eliminate the source singularities in the numerical simulations,we split the fields into primary and secondary fields.The primary fields are calculated using homogeneous or layered models with analytical solutions,and the secondary(scattered) fields are solved by the finite-element method.The linear system of equations is solved by using the large-scale sparse matrix parallel direct solver,which greatly improves the computational efficiency.The inversion algorithm was based on damping leastsquares and singular value decomposition and combined the pseudo forward modeling and reciprocity principle to compute the Jacobian matrix.Synthetic and field data were used to test the effectiveness of the proposed method.
基金financially supported by the National Natural Science Foundation of China(No.41074075/D0409)the National Science and Technology Major Project(No.2011ZX05025-001-04)
文摘As a high quality seismic imaging method, full waveform inversion (FWI) can accurately reconstruct the physical parameter model for the subsurface medium. However, application of the FWI in seismic data processing is computationally expensive, especially for the three-dimension complex medium inversion. Introducing blended source technology into the frequency-domain FWI can greatly reduce the computational burden and improve the efficiency of the inversion. However, this method has two issues: first, crosstalk noise is caused by interference between the sources involved in the encoding, resulting in an inversion result with some artifacts; second, it is more sensitive to ambient noise compared to conventional FWI, therefore noisy data results in a poor inversion. This paper introduces a frequency-group encoding method to suppress crosstalk noise, and presents a frequency- domain auto-adapting FWI based on source-encoding technology. The conventional FWI method and source-encoding based FWI method are combined using an auto-adapting mechanism. This improvement can both guarantee the quality of the inversion result and maximize the inversion efficiency.
基金supported by the National Natural Science Foundation of China(6127144261401481)
文摘For real-time jamming signal generation in deceiving inverse synthetic aperture radar(ISAR),the target characteristics modulation is always processed in the expensive field programmable gate array(FPGA).Due to the large computational complexity of the traditional modulating operation,the size and structure of simulated false-target are limited.With regard to the principle of dechirping in range compression of linear frequency modulated(LFM) radar,a novel algorithm named "inverse dechirping" is proposed for target characteristics modulation.This algorithm only needs one complex multiplier in the FPGA to generate the jamming signal when the radar signal is intercepted,which can be obtained by multiplication of radar signal samplings and the equivalent dechirped target echo in the time domain.As the complex synthesis of dechirped target echo can be realized by cheap digital signal processor(DSP) within the interpulse time,the overall cost of the jamming equipment will be reduced and the false-target size will not be limited by the scale of FPGA.Numerical simulations are performed to verify the correctness and effectiveness of the proposed algorithm.
基金supported by the Joint Fund of Seismological Science(Grant No.U1839206)the National R&D Program on Monitoring,Early Warning and Prevention of Major Natural Disaster(Grant No.2017YFC1500301)+2 种基金supported by IGGCAS Research Start-up Funds(Grant No.E0515402)National Natural Science Foundation of China(Grant No.E1115401)supported by National Natural Science Foundation of China(Grant No.11971258).
文摘The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.
基金Supported by the National Natural Science Foundation of China,No.81570363the National Key Research and Development Program of China,No.2016YFA0201304the Priority Academic Program Development of Jiangsu Higher Education Institutions,No.KYZZ15_0263
文摘AIM To establish whether frequency-domain optical coherence tomography(FD-OCT) is safe and effective in the evaluation and treatment of angiographicallyintermediate coronary lesions(ICL) METHODS Sixty-four patients with 2-dimensional quantitativecoronary angiography(2D-QCA) demonstrating ICL were included. OCT imaging was performed. According to predetermined OCT criteria, patients were assigned to either of 2 groups: OCT-guided percutaneous coronary intervention(PCI) or OCT-guided optimal medical therapy(OMT). The primary efficacy endpoint was to demonstrate the superiority and higher accuracy of FD-OCT compared to 2D-QCA in evaluating stenosis severity in patients with ICL. The primary safety endpoint was the incidence of 30-d major adverse cardiac events(MACE). Secondary endpoints included MACE at 12 mo and other clinical events.RESULTS Analysis of the primary efficacy endpoint demonstrates that 2D-QCA overestimates the stenosis severity of ICL in both the OCT-guided PCI and OMT groups, proving FD-OCT to be superior to and more precise than 2D-QCA in treating this subset of lesions. The primary safety endpoint was fully met with the incidence of 30-d MACE being nil in both the OCT-guided PCI and OCTguided OMT groups. Incidences of secondary endpoints were found to be low in both arms, the only exception being the relatively high incidence of recurrent episodes of angina which was, however, very similar in the 2 groups.CONCLUSION FD-OCT is safe and effective in the evaluation and treatment of ICL. Larger studies are needed to firmly establish the efficacy and safety of FD-OCT in treating ICL across all coronary artery disease population subgroups.
基金supported by the National Natural Science Foundation of China (No. 41130418)the Strategic Leading Science and Technology Programme (Class B) of the Chinese Academy of Sciences (No. XDB10010400)
文摘Earth medium is not completely elastic, with its viscosity resulting in attenuation and dispersion of seismic waves. Most viscoelastic numerical simulations are based on the finite-difference and finite-element methods. Targeted at viscoelastic numerical modeling for multilayered media, the constant-Q acoustic wave equation is transformed into the corresponding wave integral representation with its Green's function accounting for viscoelastic coefficients. An efficient alternative for full-waveform solution to the integral equation is proposed in this article by extending conventional frequency-domain boundary element methods to viscoelastic media. The viscoelastic boundary element method enjoys a distinct characteristic of the explicit use of boundary continuity conditions of displacement and traction, leading to a semi-analytical solution with sufficient accuracy for simulating the viscoelastic effect across irregular interfaces. Numerical experiments to study the viscoelastic absorption of different Q values demonstrate the accuracy and applicability of the method.