Practical implementation of Hilbert-Huang Transform algorithm

Hilbert-Huang Transform (HHT) is a newly developed powerful method for nonlinear and non-stationary time series analysis. The empirical mode decomposition is the key part of HHT, while its algorithm was protected by NASA as a US patent, which limits the wide application among the scientific community. Two approaches, mirror periodic and extrema extending methods, have been developed for handling the end effects of empirical mode decomposition. The implementation of the HHT is realized in detail to widen the application. The detailed comparison of the results from two methods with that from Huang et al. (1998, 1999), and the comparison between two methods are presented. Generally, both methods reproduce faithful results as those of Huang et al. For mirror periodic method (MPM), the data are extended once forever. Ideally, it is a way for handling the end effects of the HHT, especially for the signal that has symmetric waveform. The extrema extending method (EEM) behaves as good as MPM, and it is better t     

The Mid-atlantic Ridge (31°S–34°30′S): Temporal and spatial variations of accretionary processes

The ridge located between 31° S and 34°30′S is spreading at a rate of 35 mm yr⁻¹, a transitional velocity between the very slow (≤20 mm yr⁻¹) opening rates of the North Atlantic and Southwest Indian Oceans, and the intermediate rates (60 mm yr⁻¹) of the northern limb of the East Pacific Rise, and the Galapagos and Juan de Fuca Ridges. A synthesis of multi-narrow beam, magnetics and gravity data document that in this area the ridge represents a dynamically evolving system. Here the ridge is partitioned into an ensemble of six distinct segments of variable lengths (12 to 100 km) by two transform faults (first-order discontinuities) and three small offset (< 30 km) discontinuities (second-order discontinuities) that behave non-rigidly creating complex and heterogeneous morphotectonic patterns that are not parallel to flow lines. The offset magnitudes of both the first and second-order discontinuities change in response to differential asymmetric spreading. In addition, along the fossil trace of second-order discontinuities, the lengths of abyssal hills located to either side of a discordant zone are observed to lengthen and shorten creating a saw-toothed pattern. Although the spreading rate remains the same along the length of the ridge studied, the morphology of the spreading segments varies from a deep median valley with characteristics analogous to the rift segments of the North Atlantic to a gently rifted axial bulge that is indistinguishable from the shape and relief of the intermediate rate spreading centers of the East Pacific Rise (i.e., 21°N). Like other carefully surveyed ridge segments at slow and fast rates of accretion, the along-axis profiles of each ridge segment are distinctly convex upwards, and exhibit along-strike changes in relief of 500m to 1500 between the shallowest portion of the segment (approximate center) and the segment ends. Such spatial variations create marked along-axis changes in the morphology and relief of each segment. A relatively low mantle Bouguer anomaly is known to be associated with the ridge segment characterized by a gently rifted axial bulge and is interpreted to indicate the presence of focused mantle upwelling (Kuo and Forsyth, 1988). Moreover, the terrain at the ends of each segment are known to be highly magnetized compared to the centers of each segment (Carbotte et al, 1990). Taken together, these data clearly establish that these profound spatial variations in ridge segment properties between adjoining segments, and along and across each segment, indicate that the upper mantle processes responsible for the formation of this contrasting architecture are not solely related to passive upwelling of the asthenosphere beneath the ridge axis. Rather, there must be differences in the thermal and mechanical structure of the crust and upper mantle between and along the ridge segments to explain these spatial variations in axial topography, crustal structure and magnetization. These results are consistent with the results of investigations from other parts of the ridge and suggest that the emplacement of magma is highly focused along segments and positioned beneath the depth minimum of a given segment. The profound differences between segments indicate that the processes governing the behavior of upwelling mantle are decoupled and the variations in the patterns of axis flanking morphology and rate of accretion indicate that processes controlling upwelling and melt production vary markedly in time as well. At this spreading rate and in this area, the accretionary processes are clearly three-dimensional. In addition, the morphology of a ridge segment is not governed so much by opening rate as by the thermal structure of the mantle which underlies the segment.     

Comparison of methods to model the gravitational gradients from topographic data bases

A number of methods have been developed over the last few decades to model the gravitational gradients using digital elevation data. All methods are based on second-order derivatives of the Newtonian mass integral for the gravitational potential. Foremost are algorithms that divide the topographic masses into prisms or more general polyhedra and sum the corresponding gradient contributions. Other methods are designed for computational speed and make use of the fast Fourier transform (FFT), require a regular rectangular grid of data, and yield gradients on the entire grid, but only at constant altitude. We add to these the ordinary numerical integration (in horizontal coordinates) of the gradient integrals. In total we compare two prism, two FFT and two ordinary numerical integration methods using 1" elevation data in two topographic regimes (rough and moderate terrain). Prism methods depend on the type of finite elements that are generated with the elevation data; in particular, alternative triangulations can yield significant differences in the gradients (up to tens of Eötvös). The FFT methods depend on a series development of the topographic heights, requiring terms up to 14th order in rough terrain; and, one popular method has significant bias errors (e.g. 13 Eötvös in the vertical–vertical gradient) embedded in its practical realization. The straightforward numerical integrations, whether on a rectangular or triangulated grid, yield sub-Eötvös differences in the gradients when compared to the other methods (except near the edges of the integration area) and they are as efficient computationally as the finite element methods.     

HHT方法在地震工程中的应用之初步探讨

详细介绍了一种新的非线性、非平稳信号的处理方法-HHT，讨论了其有效性和优越性，简单介绍了其在地震工程中的应用前景，并利用HHT处理了地震工程中常用El Centro地震记录，得到了该记录的振幅-频率-时间分布特性，并将得到的边际谱与傅立叶谱做了对比。     

Hilbert-Huang变换在密频结构阻尼识别中的应用

Hilbert—Huang变换是一种新的数据处理方法，由经验模分解(Empirical Mode Decomposition)技术及Hilbert变换两部分组成。本文研究此方法对于密频结构阻尼识别的应用。首先对于两自由度系统模型，说明该方法用于阻尼识别的步骤。进而研究存在频率密集现象的高层建筑的阻尼识别问题。上述结果与理论值及由半功率带宽法的识别值进行了比较，对比显示Hilbert．Huang方法较传统方法具有良好的识别密频结构阻尼的性能，适用于大型结构的系统识别。     

SH Wave Propagation in Viscoelastic Media

When comparing solutions for the propagation of SH waves in plane parallel layered elastic and viscoelastic (anelastic) media, one of the first things that becomes apparent is that in the elastic case the location of the saddle points required to obtain a high frequency approximation are located on the real p axis. This is true of the branch points also. In a viscoelastic medium this is not typical. The saddle point corresponding to an arrival lies in the first quadrant of the complex p-plane as do the branch points. Additionally, in the elastic case the saddle point and branch points lie on a straight line drawn through the origin (the positive real axis in the complex p-plane), while in the viscoelastic case this is generally not the case and the saddle point and branch points lie in such a manner as to indicate the degree of their complex values.In this paper simple SH reflected and transmitted particle displacement arrivals due to a point torque source at the surface in a viscoelastic medium composed of a layer over a half space will be considered. The path of steepest descent defining the saddle point in the first quadrant will be parameterized in terms of a real variable and the high frequency solutions and intermediate analytic results obtained will be used to formulate more specific constraints and observations regarding saddle point location relative to branch point locations in the complex p-plane.As saddle point determination for an arrival is, in general, the solution of a non-linear equation in two unknowns (the real and imaginary parts of the complex saddle point p ₀), which must be solved numerically, the use of analytical methods for investigating this problem type is somewhat limited.Numerical experimentation using well documented solution methods, such as Newton's method, was undertaken and some observations were made. Although fairly basic, they did provide for the design of algorithms for the computation of synthetic traces that displayed more efficient convergence and accuracy than those previously employed. This was the primary motivation for this work and the results from the SH problem may be used with minimal modifications to address the more complicated subject of coupled P-SV wave propagation in viscoelastic media.Another reason for revisiting a problem that has received some attention in the literature was to approach it in a fairly comprehensive manner so that a number of specific observations may be made regarding the location of the saddle point in the complex p-plane and to incorporate these into computer software. These have been found to result in more efficient algorithms for the SH wave propagation and a significant enhancement of the comparable software in the P-SV problem.     

Analysis of wave groups by wave envelope-phase and the Hilbert Huang transform methods

The wave groups are studied by both conventional wave analysis methods and by the non-stationary Hilbert Huang Transform (HHT) method. Full-scale wave records containing abnormal waves are used. Instantaneous quantities, such as envelope, phase and frequency, are adopted to study the wave grouping. A refined definition of wave group is proposed considering that the wave process is simultaneously amplitude and frequency modulated. The validation of the proposed definition is conducted by analysis of numerical simulation data. Group parameters are proposed based on the time-frequency distribution of energy. An attempt is made to find the relationship between the characteristics of abnormal waves and the group characteristics.     

希尔伯特-黄变换的端点延拓

希尔伯特-黄变换(HHT)是近年来发展起来的一种新的时间序列信号分析方法.该文在对HHT深入研究与充分肯定的基础上,发展了信号的镜像闭合延拓和包络的极值延拓两种方法.通过几个典型的例子检验了两种方法,并与Huang等(1998,1999)进行了比较,得到了令人满意的结果.镜像闭合延拓法根据信号端点的分布特性,把镜子放在具有对称性的极值位置,通过镜像法把镜内信号映射成一个周期性的信号,不存在端点,从根本上避免了经验模态分解和希尔伯特变换的端点问题.极值延拓法简单易行,具有与镜像闭合法相当的效果,在处理非对称波形信号时更显其优越性.     

HHT-EEMD用于IGS站高程时间序列分析

利用经验模态分解（EMD）和整体经验模态分解（EEMD）方法,将BJFS站2000~2015年高程时间序列进行分解,发现其不仅存在1 a、0.5 a、0.25 a、2个月、1个月以及长周期等周期项,还存在以往方法很难探测出来的近似2 a周期信号。与EMD分解结果对比,整体经验模态分解可以明显减弱模式混叠现象。对各分量进行Hilbert 变换（HHT）,得到时间-频率-能量的Hilbert 频谱图。结果表明,年周期和2 a周期变化是高程运动的主要贡献项。利用小波变换方法对比验证EEMD的分解结果表明,与小波分析相比,EEMD重构信号与高程序列差异的RMS更小,证明了HHT-EEMD方法在数据资料分析过程中的有效性。对环境负载及GRACE负载造成的测站位移进行功率谱分析得出,环境负载确实会造成IGS站高程时间序列的1 a、0.5 a以及季节性运动,GRACE负载还验证了2 a信号的存在。