Seismic wavefield modeling based on time-domain symplectic and Fourier finite-difference method

Seismic wavefield modeling is important for improving seismic data processing and interpretation. Calculations of wavefield propagation are sometimes not stable when forward modeling of seismic wave uses large time steps for long times. Based on the Hamiltonian expression of the acoustic wave equation, we propose a structure-preserving method for seismic wavefield modeling by applying the symplectic finite-difference method on time grids and the Fourier finite-difference method on space grids to solve the acoustic wave equation. The proposed method is called the symplectic Fourier finite-difference （symplectic FFD） method, and offers high computational accuracy and improves the computational stability. Using acoustic approximation, we extend the method to anisotropic media. We discuss the calculations in the symplectic FFD method for seismic wavefield modeling of isotropic and anisotropic media, and use the BP salt model and BP TTI model to test the proposed method. The numerical examples suggest that the proposed method can be used in seismic modeling of strongly variable velocities, offering high computational accuracy and low numerical dispersion. The symplectic FFD method overcomes the residual qSV wave of seismic modeling in anisotropic media and maintains the stability of the wavefield propagation for large time steps.

A stable staggered-grid finite-difference scheme for acoustic modeling beyond conventional stability limit

Staggered-grid finite-difference(SGFD)schemes have been widely used in acoustic wave modeling for geophysical problems.Many improved methods are proposed to enhance the accuracy of numerical modeling.However,these methods are inevitably limited by the maximum Courant-Friedrichs-Lewy(CFL)numbers,making them unstable when modeling with large time sampling intervals or small grid spacings.To solve this problem,we extend a stable SGFD scheme by controlling SGFD dispersion relations and maximizing the maximum CFL numbers.First,to improve modeling stability,we minimize the error between the FD dispersion relation and the exact relation in the given wave-number region,and make the FD dispersion approach a given function outside the given wave-number area,thus breaking the conventional limits of the maximum CFL number.Second,to obtain high modeling accuracy,we use the SGFD scheme based on the Remez algorithm to compute the FD coefficients.In addition,the hybrid absorbing boundary condition is adopted to suppress boundary reflections and we find a suitable weighting coefficient for the proposed scheme.Theoretical derivation and numerical modeling demonstrate that the proposed scheme can maintain high accuracy in the modeling process and the value of the maximum CFL number of the proposed scheme can exceed that of the conventional SGFD scheme when adopting a small maximum effective wavenumber,indicating that the proposed scheme improves stability during the modeling.

An adaptive finite-difference method for seismic traveltime modeling based on 3D eikonal equation

3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.

考虑衬砌截面协调变形约束的既有隧道受盾构下穿施工影响的Fourier能量变分解

目前关于新建盾构开挖诱发既有地铁隧道变形的理论研究,大多将既有隧道简化为不含接头的Euler-Bernoulli梁,而忽视了衬砌截面的剪切变形与管片接头的局部刚度削弱。此外,既有地基模型理论较少考虑衬砌管片与接触土体的协调变形约束。首先,采用Loganathan-Polous公式计算盾构开挖诱发的土体自由位移场,并将既有隧道视为考虑剪切变形的Timoshenko梁,基于沿轴线变化的抗弯与剪切刚度函数来描述管片接头的局部削弱作用;然后,结合能量变分原理与连续弹性体Mindlin解,建立考虑衬砌截面协同变形约束的隧道位移控制方程,将盾构开挖诱发的附加荷载施加于既有隧道,并使用Fourier级数方法求解既有隧道纵向变形;最后,将Fourier能量变分解与3组工程实测数据进行对比验证,取得了较好的一致性。此外,考虑隧道抗弯刚度、剪切刚度、接头刚度削弱作用、土体参数与隧道空间位置等因素的共同影响,使用抗弯系数、剪切系数及接头系数进行敏感性分析。结果表明,考虑截面变形协调约束的理论解更为符合实际,不考虑变形协调约束所得隧道位移结果偏大;衬砌接头局部刚度的削弱作用将造成隧道弯矩沿轴线在接头处发生明显锯齿状突变;抗弯系数与剪切系数的增长均导致隧道结构所受弯矩降低;剪切系数由低到高的增长将导致隧道变形模式发生“弯曲控制-弯剪混合-剪切控制”三阶段变化;隧道变形模式进入“剪切控制”阶段的过程中,抗弯刚度的变化对隧道内力影响逐渐降低;接头系数的增长将造成隧道所受弯矩降低,但隧道抗弯刚度较大时其影响逐渐减弱。

基于Fourier变换红外光谱分析咖啡豆成分

采用Fourier变换红外光谱对不同产区、不同海拔的5组浅中烘焙的咖啡豆样品进行分析,根据官能团的红外光谱特征峰分析样品主要成分;为进一步分析5组样品的成分差异,对原始光谱进行二阶导数处理,结合聚类分析理论建立平均偏离分析方法,计算出5组样品间的平均偏离大小并进行分析.结果表明:5组样品的红外光谱特征峰峰形相似,即主要成分相同;根据红外光谱得到平均偏离大小与产区海拔差值呈正相关.研究结果为分析咖啡豆产区和海拔提供鉴别依据,并对红外光谱研究提供一定参考.

分数阶Fourier变换的测不准原理

该文参考Fourier变换的性质研究了离散分数阶Fourier变换的测不准原理以及连续分数阶Fourier变换在Lebesgue测度下的测不准原理,使得分数阶Fourier变换的测不准原理性质更一般化.

周期边界条件下四阶特征值问题的一种有效的Fourier谱逼近

文章提出了周期边界条件下四阶特征值问题的一种有效的Fourier谱逼近方法.首先,根据周期边界条件引入了适当的Sobolev空间和相应的逼近空间,建立了原问题的一种弱形式及其离散格式,并推导了等价的算子形式.其次,定义了正交投影算子,并证明了其逼近性质,结合紧算子的谱理论证明了逼近特征值的误差估计.另外,构造了逼近空间中的一组基函数,推导了离散格式基于张量积的矩阵形式.最后,文章给出了一些数值算例,数值结果表明其算法是有效的和谱精度的.

Fourier正则化同时反演热传导方程的初值和源项

研究边值问题中从已知的给定数据中同时识别初始条件和未知源的问题。用Fourier变换求出问题的精确解,发现该反问题是严重不适定的。为了解决所提出的问题,本文中用Fourier正则化方法构造出正则化解近似解,并在边界条件的假设下,得到先验参数选取下的精确解与其正则化近似解之间的对数型误差估计。

Super-Fast Approximation Algorithms Using Classical Fourier Tools

In the author’s recent publications, a parametric system biorthogonal to the corresponding segment of the exponential Fourier system was unusually effective. On its basis, it was discovered that knowledge of a finite number of Fourier coefficients of function f from an infinite-dimensional set of elementary functions allows f to be accurately restored (the phenomenon of over-convergence). Below, parametric biorthogonal systems are constructed for classical trigonometric Fourier series, and the corresponding phenomena of over-convergence are discovered. The decisive role here was played by representing the space L2 as an orthogonal sum of two corresponding subspaces. As a result, fast parallel algorithms for reconstructing a function from its truncated trigonometric Fourier series are proposed. The presented numerical experiments confirm the high efficiency of these convergence accelerations for smooth functions. In conclusion, the main results of the work are summarized, and some prospects for the development and generalization of the proposed approaches are discussed.

Enhanced Fourier Transform Using Wavelet Packet Decomposition

Many domains, including communication, signal processing, and image processing, use the Fourier Transform as a mathematical tool for signal analysis. Although it can analyze signals with steady and transitory properties, it has limits. The Wavelet Packet Decomposition (WPD) is a novel technique that we suggest in this study as a way to improve the Fourier Transform and get beyond these drawbacks. In this experiment, we specifically considered the utilization of Daubechies level 4 for the wavelet transformation. The choice of Daubechies level 4 was motivated by several reasons. Daubechies wavelets are known for their compact support, orthogonality, and good time-frequency localization. By choosing Daubechies level 4, we aimed to strike a balance between preserving important transient information and avoiding excessive noise or oversmoothing in the transformed signal. Then we compared the outcomes of our suggested approach to the conventional Fourier Transform using a non-stationary signal. The findings demonstrated that the suggested method offered a more accurate representation of non-stationary and transient signals in the frequency domain. Our method precisely showed a 12% reduction in MSE and a 3% rise in PSNR for the standard Fourier transform, as well as a 35% decrease in MSE and an 8% increase in PSNR for voice signals when compared to the traditional wavelet packet decomposition method.

基于Fourier拟合的光伏跟踪系统设计

[目的]为提高光伏支架的发电效率,以基于Fourier拟合算法的视日运动轨迹光伏跟踪器为研究设计对象。[方法]借助MATLAB分析了传统视日运动轨迹算法的跟踪精度情况,更进一步地在此基础上提出了精度更优的8阶傅里叶拟合视日运动轨迹跟踪算法,并结合GNSS定位技术,得出太阳实时运行轨迹。系统设计以STM32单片机作为微处理器,采用6轴加速度传感器,通过重力矢量在其轴上的投影确定光伏跟踪支架的实时倾斜度,以此二者建立跟踪支架三维运行模型,通过分析太阳光线在支架上的余弦效应即可得到电机需要运行的动作角。同时为解决晨昏时刻光伏组件间的阴影遮挡问题,系统增加逆跟踪算法。并考虑到大规模机组的应用,提出主从电机同步控制策略。[结果]傅-2里叶拟合算法拥有更高的跟踪精度,精度可达到10数量级,高于传统算法一个数量级。同时采用逆跟踪技术可节约24.3%光伏阵列用地面积,大幅提高了土地利用率。[结论]本研究采用更为精准的视日运动轨迹跟踪模型,较好地解决了太阳辐射利用的余弦效应,提升了电站的发电效率,实现了安全高效绿色能源体系的建设,促进了我国“双碳”目标的推进与达成。

A Deepfake Detection Algorithm Based on Fourier Transform of Biological Signal

Deepfake-generated fake faces,commonly utilized in identity-related activities such as political propaganda,celebrity impersonations,evidence forgery,and familiar fraud,pose new societal threats.Although current deepfake generators strive for high realism in visual effects,they do not replicate biometric signals indicative of cardiac activity.Addressing this gap,many researchers have developed detection methods focusing on biometric characteristics.These methods utilize classification networks to analyze both temporal and spectral domain features of the remote photoplethysmography(rPPG)signal,resulting in high detection accuracy.However,in the spectral analysis,existing approaches often only consider the power spectral density and neglect the amplitude spectrum—both crucial for assessing cardiac activity.We introduce a novel method that extracts rPPG signals from multiple regions of interest through remote photoplethysmography and processes them using Fast Fourier Transform(FFT).The resultant time-frequency domain signal samples are organized into matrices to create Matrix Visualization Heatmaps(MVHM),which are then utilized to train an image classification network.Additionally,we explored various combinations of time-frequency domain representations of rPPG signals and the impact of attention mechanisms.Our experimental results show that our algorithm achieves a remarkable detection accuracy of 99.22%in identifying fake videos,significantly outperforming mainstream algorithms and demonstrating the effectiveness of Fourier Transform and attention mechanisms in detecting fake faces.

Mid-infrared Optical Frequency Comb-based Fourier Transform Spectrometer for Broadband Molecular Spectroscopy

Optical frequency combbased Fourier transform spectroscopy has the features of broad spectral bandwidth,high sensitivity,andmultiplexed trace gas detection,which has valuable application potential in the fields of precision spectroscopy and trace gas detection.Here,we report the development of a mid-infrared Fourier transform spectrometer based on an optical frequency comb combined with a Herriott-type multipass cell.Using this instrument,the broadband absorption spectra of several important molecules,including methane,acetylene,water molecules and nitrous oxide,are measured by near real-time data acquisition in the 2800-3500 cm^(-1)spectral region.The achieved minimum detectable absorption of the instrument is 4.4×10^(-8)cm^(-1)·Hz^(-1/2)per spectral element.Broadband spectra of H_(2)0 are fited using the Voigt profile multispectral fitting technique and the consistency of the concentration inversion is 1%.Our system also enables precise spectroscopic measurements,and it allows the determination of the spectral line positions and upper state constants of N_(2)O in the(0002)-(1000)band,with results in good agreement with those reported by Toth[Appl.Opt.30,5289(1991)].

An improved non-uniform fast Fourier transform method for radio imaging of coronal mass ejections

Radioheliographs can obtain solar images at high temporal and spatial resolution,with a high dynamic range.These are among the most important instruments for studying solar radio bursts,understanding solar eruption events,and conducting space weather forecasting.This study aims to explore the effective use of radioheliographs for solar observations,specifically for imaging coronal mass ejections(CME),to track their evolution and provide space weather warnings.We have developed an imaging simulation program based on the principle of aperture synthesis imaging,covering the entire data processing flow from antenna configuration to dirty map generation.For grid processing,we propose an improved non-uniform fast Fourier transform(NUFFT)method to provide superior image quality.Using simulated imaging of radio coronal mass ejections,we provide practical recommendations for the performance of radioheliographs.This study provides important support for the validation and calibration of radioheliograph data processing,and is expected to profoundly enhance our understanding of solar activities.

对称幂平方L-函数Fourier系数的高次均值估计

设f(z)为全模群SL(2,Z)上权为偶数k的Hecke特征形式,L(s,sym^(2)f)为其对应的对称幂平方L-函数,λsym^(2)f(n)为对称幂平方L-函数L(s,sym^(2)f)的Fourier展开式中第n个正规化Fourier系数.本文借助Newton和Thorne近期关于函数Symjf的重要工作以及自守L-函数的亚凸界和均值估计结果,研究了λsym^(2)f(n)在四个整数平方和序列上的三次均值估计和四次均值估计,得到了均值估计的渐近公式,改进了之前的结果.

Performance of Continuous Wavelet Transform over Fourier Transform in Features Resolutions

This study presents a comparative analysis of two image enhancement techniques, Continuous Wavelet Transform (CWT) and Fast Fourier Transform (FFT), in the context of improving the clarity of high-quality 3D seismic data obtained from the Tano Basin in West Africa, Ghana. The research focuses on a comparative analysis of image clarity in seismic attribute analysis to facilitate the identification of reservoir features within the subsurface structures. The findings of the study indicate that CWT has a significant advantage over FFT in terms of image quality and identifying subsurface structures. The results demonstrate the superior performance of CWT in providing a better representation, making it more effective for seismic attribute analysis. The study highlights the importance of choosing the appropriate image enhancement technique based on the specific application needs and the broader context of the study. While CWT provides high-quality images and superior performance in identifying subsurface structures, the selection between these methods should be made judiciously, taking into account the objectives of the study and the characteristics of the signals being analyzed. The research provides valuable insights into the decision-making process for selecting image enhancement techniques in seismic data analysis, helping researchers and practitioners make informed choices that cater to the unique requirements of their studies. Ultimately, this study contributes to the advancement of the field of subsurface imaging and geological feature identification.

Design and implementation of code acquisition using sparse Fourier transform

Code acquisition is the kernel operation for signal synchronization in the spread-spectrum receiver.To reduce the computational complexity and latency of code acquisition,this paper proposes an efficient scheme employing sparse Fourier transform(SFT)and the relevant hardware architecture for field programmable gate array(FPGA)and application-specific integrated circuit(ASIC)implementation.Efforts are made at both the algorithmic level and the implementation level to enable merged searching of code phase and Doppler frequency without incurring massive hardware expenditure.Compared with the existing code acquisition approaches,it is shown from theoretical analysis and experimental results that the proposed design can shorten processing latency and reduce hardware complexity without degrading the acquisition probability.

Modeling a Periodic Signal Using Fourier Series

This paper covers the concept of Fourier series and its application for a periodic signal. A periodic signal is a signal that repeats its pattern over time at regular intervals. The idea inspiring is to approximate a regular periodic signal, under Dirichlet conditions, via a linear superposition of trigonometric functions, thus Fourier polynomials are constructed. The Dirichlet conditions, are a set of mathematical conditions, providing a foundational framework for the validity of the Fourier series representation. By understanding and applying these conditions, we can accurately represent and process periodic signals, leading to advancements in various areas of signal processing. The resulting Fourier approximation allows complex periodic signals to be expressed as a sum of simpler sinusoidal functions, making it easier to analyze and manipulate such signals.

A Physical Security Technology Based upon Doubly Multiple Parameters Weighted Fractional Fourier Transform

Enhancing the security of the wireless communication is necessary to guarantee the reliable of the data transmission, due to the broadcast nature of wireless channels. In this paper, we provide a novel technology referred to as doubly multiple parameters weighted fractional Fourier transform(DMWFRFT), which can strengthen the physical layer security of wireless communication. This paper introduces the concept of DM-WFRFT based on multiple parameters WFRFT(MP-WFRFT), and then presents its four properties. Based on these properties, the parameters decryption probability is analyzed in terms of the number of parameters. The number of parameters for DM-WFRFT is more than that of the MP-WFRFT,which indicates that the proposed scheme can further strengthen the the physical layer security. Lastly, some numerical simulations are carried out to illustrate that the efficiency of proposed DM-WFRFT is related to preventing eavesdropping, and the effect of parameters variety on the system performance is associated with the bit error ratio(BER).

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.