重力固体潮IMF的AM-FM非线性最小二乘拟合

周挚(2005)的重力固体潮研究表明,EMD分解后的重力固体潮IMF(本征模态函数)的信号特征与AM-FM(调频调幅)信号非常接近.AM-FM信号是典犁的非平稳信号.基于AM-FM研究的前沿性,本文借助于函数构造理论和非线性逼近方法,求取重力固体潮IMF信号的数学显式,为进一步的解析工作做好铺垫.首先建立AM-FM数学模型:AM-FM信号y(t)=A(f)*F(t).依据拉格朗日定理、傅立叶级数基本性质和三角函数基本性质,分别用若干正弦函数之和的级数形式来描述AM信号A(t)和FM信号F(t).     

基于EMD的信号瞬时特征的小波分析方法

提出了一种基于EMD(Empirical Mode Decomposition)的信号瞬时特征的小波分析方法。用这种方法提取非平稳信号的瞬时频率和瞬时幅值分三个基本步骤：首先,用EMD把信号分解成IMF(Intrinsic Mode Function)分量;接着,对IMF分量进行小波分析,从小波系数的幅角函数中提取小波脊线;最后,从小波脊线中提取瞬时频率和瞬时幅值。通过对仿真信号的分析,验证了该方法能有效地分析非平稳信号。     

固体潮时频分析新方法

对传统时频分析方法的特点和局限性作了介绍。采用Huang．Hilbert Transformation方法分析了固体潮数据，即利用经验模态分解，将复杂信号分解成有限个本征模态函数之和，求取三维时频空间的Hilbert谱。仿真实验和重力固体潮实例分析表明，该方法适用于调频调幅信号特征的提取。     

基于折线包络的经验模态分解方法

1998年,Huang提出了处理非平稳信号的HHT方法(Hilbert-Huang Transform,简称HILT).该方法包括两个步骤:①任意信号首先经过经验模态分解方法(Empirical Mode Decomposition,简称EMD)被分解为一系列固有模态函数(IntrinsicModeFunction,简称IMF).     

重力固体潮信号独立成分的谱相关解调与地震前兆信息提取

重力固体潮信号是一种多谐波的混合信号,为了提取其中所包含的地震前兆异常信息.本文结合重力固体潮的产生机制,建立了一种重力固体潮正交分解模型.在此基础上,利用独立分量分析算法实现重力固体潮信号的加性分解,然后针对独立分量中的调制关系,利用谱相关方法对其进行乘性解调.从而,完整地提取出了重力固体潮信号中丰富的潮汐谐波信息.进一步,引入理论计算值作为实际测量值的参考背景,在独立分量中凸显出原重力固体潮信号中的异常变化特征.通过对云南地区的实际震例研究表明,重力固体潮独立分量的异常特征与地震事件的时序存在密切的相关性.在地震发生前的1~5个月内,对应于重力固体潮信号长周期谐波系的独立分量在时域波形和循环相关谱方面均有明显的异常变化,而且普遍存在,充分反映了这一异常变化与地震前地壳内部能量的变化有关,很有可能就是重力固体潮信号中隐含的地震前兆信息.     

基于HHT提取昆明、下关重力固体潮的地震前兆信息

在重力固体潮地震前兆分析中引入HHT时频分析新方法.结合HHT的优越性、固体潮的特点和地震的非平稳过程特性,设计重力固体潮地震前兆分析的瞬时频率特征参数;以相应理论计算值作为参照背景,研究固体潮的震前变化特征.昆明、下关的震例分析表明, 的确存在瞬时频率特征参数的震前变化,且具短期、同步正异常特征;瞬时频率特征参数具有明确的物理意义,其震前变化反映了地震非平稳过程对理论重力固体潮的影响.     

基于经验模态分解及独立成分分析的微震信号降噪方法

针对微震信号具有高噪声、突变快、随机性强等特点,基于经验模态分解(EMD)及独立成分分析(ICA)提出一种微震信号降噪方法.首先,对含噪信号进行EMD分解,获得一系列按频率从高到低的内蕴模态函数(IMF),利用原信号与各IMF之间的互相关系数辨识出噪声与信号的分界,将分界之上的高频噪声滤除;其次,为有效去除分界IMF中的模态混叠噪声,基于ICA算法对分界IMF进行盲源分离,提取其中的微震有效信号,并将其与剩余的IMF累加重构,从而得到降噪后的微震信号;最后,利用快速傅里叶变换(FFT)时频谱对比分析降噪前后的信号特征,定性说明本文方法的有效性;引入信噪比和降噪后信号占原信号的能量百分比两个参数,定量说明本文方法能充分保留微震信号的瞬态非平稳特征,降噪效果明显.     

Hilbert-Huang变换在提取阵列声波信号动力特性中的应用

阵列声波信号是典型的非线性、非平稳信号,其动力特性的量化提取对于进行地层结构构造分析提供了必要的基础资料.而Hilbert-Huang变换(HHT)是一种处理非线性、非平稳信号的新方法.它通过经验模态分解(EMD)将信号分解为有限个固有模态函数(IMF),并对每个固有模态函数进行Hilbert变换得到Hilbert谱.本文将这种方法应用于阵列声波信号动力特性的提取,有效地获得了信号能量的时频分布,瞬时能量、Hilbert能量、最大振幅对应的时频分布等动力特性,显示了HHT的优势以及对于进一步实现地层结构构造分析的重要意义.     

过渡内蕴模态函数对经验模态分解去噪结果的影响研究及改进算法

经验模态分解(Empirical Mode Decomposition, EMD)是一种具有较大应用潜力的去噪算法.目前,该算法存在的一个较大问题是过渡内蕴模态函数(Intrinsic Mode Function, IMF)中混叠噪声不能有效处理.过渡内蕴模态函数中混叠噪声不易剔除,限制了该算法的应用.本文针对此问题,通过研究过渡IMF的特点,首次提出一种有效去除过渡IMF中混叠噪声的方法.该方法首先对原信号进行一次EMD处理,得到包含过渡IMF的初步去噪结果,并将其与合适的余弦信号结合,改变其包络分布,然后对其结果再次进行EMD处理,仿真实验表明该方法在保留有效信号的同时,可以有效的去除过渡IMF中混叠的噪声,并将该方法用于实际地震资料随机噪声压制,处理效果令人满意.     

基于分布估计算法对重力固体潮信号独立主成分分析及谐波提取

固体潮信号是地球自转在月球和太阳作用下产生的混合复杂信号.本文提出一个重力固体潮信号的正交分解模型,通过将重力固体潮信号谐波分量,分解在两个正交的方向上,以提取独立的谐波分量.为了实现这一分解,本文采用独立成分分析算法进行验证分析.独立成分分析主要用于提取信号中的独立成分,但一般的独立成分分析算法具有对初始值敏感,收敛较早熟等问题.本文利用分布估计算法的学习概率模型全局寻优,克服独立成分分析算法缺点.通过实验仿真分析,改进算法可将重力固体潮信号成功分离出以下独立的谐波状态:半日波信号、日波信号、长周期波信号.并且,这些分离出的独立成分与本文提出的重力固体潮信号正交分解模型一致,分别对应于模型中的正交谐波分量.     

Contribution to the knowledge of the magmatic evolution by the study of the variation of the coefficient of viscosity

The determination of the coefficient of viscosity of eruptive products gives useful elements to the knowledge of possible variations of composition and physical conditions of the magmas.     

Pictorial Series of the Manifestations of the Dynamics of the Earth

Studia Geophysica et Geodaetica -     

Evolution of the Source of the Active Region of the Sun

Geomagnetism and Aeronomy - The magnetic evolution of the active region of the Sun at the level of the photosphere has been studied. Magnetic synoptic maps of the 2007 Kitt Peak Observatory were...     

Magnetic stratigraphy of the Ordovician in the lower reach of the Kotuy River: the age of the Bysy-Yuryakh stratum and the rate of geomagnetic reversals on the eve of the superchron

Until recently, the existing data prevented the geophysicists from accurately dating the Bysy-Yuryakh stratum, which outcrops in the middle reach of the Kotuy River, constraining the time of its formation to a wide interval from the end of the Late Cambrian to the beginning of the Silurian. The obtained paleomagnetic data unambiguously correlate the Bysy-Yuryakh stratum to the Nyaian regional stage and constrain its formation, at least a considerable part of it, by the Tremadocian. This result perfectly agrees with the data on the Bysy-Yuryakh conodonts studied in this work and yields a spectacular example of the successful application of paleomagnetic studies in solving important tasks of stratigraphy and, correspondingly, petroleum geology. Within the Bysy-Yuryakh stratum, we revealed a large normal-polarity interval corresponding to the long (>1 Ma) period when the geomagnetic reversals were absent. This result, in combination with the data for the Tremadocian and Middle–Upper Cambrian sequences of the other regions, indicates that (1) the rate of occurrence of the geomagnetic reversals on the eve of the Ordovician Moyero superchron of reversed polarity was at most one reversal per Ma; (2) the superchron does not switch on instantaneously but is preceded by a certain gradual change in the operation conditions of the dynamo mechanism which, inter alia, manifests itself by the reduction of the frequency of geomagnetic reversals with the approach of the superchron. This finding supports the views according to which a process preparing the establishment of the superchrons takes place at the core–mantle boundary.     

On the determination of the non-linearity of the gravimeter scale

Summary The results of two different methods of the determination of the non-linearity of the gravimeter scale are compared; the economy is discussed as well.     

Influence of the topography of the surface of the Earth on vertical sounding of the temperature profile

Atmospheric temperature and humidity fields as well as information on other meteorological parameters are nowadays retrieved from radiance measurements recorded by operational meteorological satellites. Up to now, the inversion procedures used only take into account crude information on the topography of the Earth’s surface. However, the applied radiative transfer codes have to consider the Earth’s surface as the lower boundary of the atmospheric model and, therefore, need a more precise mean elevation and a classification of the roughness of the Earth’s surface. The influence of the topography of the Earth surface on retrieved temperature profiles is studied by using a physico-statistical inversion method. An objective analysis is made of the more precise mean elevation and derivation of roughness parameters using a new high-resolution digital elevation model (DEM) with a resolution of 500 m×500 m. By means of a geomorphological process and a newly developed topography rejection test, areas with a high surface roughness are localized and singled out. The influence of topography on the retrieved temperature profiles is illustrated by case studies. Changes are found predominantly in areas with a high variation of topography. Using the new high-resolution DEM and the topography rejection test, the geographical position of the calculated temperature profiles tends to be shifted towards areas with a small vertical variation of topography. The mean elevation determined by the new elevation model better characterizes the area observed. Hence, the temperature profiles can be calculated down to lower atmospheric levels. Furthermore, a guess profile better describing the atmospheric situation is selected by the more precise elevation. In addition, the temperature profiles obtained near the coast are improved considerably by the more precise determination of the surface property ‘sea’ and ‘land,’ respectively. Integration of an independent physical information such as topography leads, on average, to a slight improvement of the results of the physico-statistical inversion procedure. In some cases, however, significant improvements have been achieved regarding the desired accuracy of temperature profiles of the order of 1 K. In future, the spatial resolution of new high-resolution sounding instrumentation on the next generation of operational meteorological satellites will be increased. To exploit the resolving power of this new instrumentation, the different variation of the topography of the Earth surface, especially in regions with a high variation of topography, can be taken into account more precisely by using a high-resolution DEM.     

The effect of the pollution of the ground layer of the atmosphere on the albedo

Summary A simple model is proposed suitable for studying the effect of the ground layer of the atmosphere, polluted by aerosol, on the albedo. This model is founded on solving the equation of transfer of radiative energy. The numerical results are discussed, particular attention being paid to the analysis of errors due to neglecting the multiple reflection of solar radiation on the aerosol particles. A method which would also include the multiple reflection is proposed, and the conditions under which the presence of the aerosol is responsible for an increase or decrease of the solar radiation balance on the Earth's surface, are analysed.     

On the solution of the reverse problem of geodesics with the help of the chord

Резюме Дается вывод формул для вычислений длины хорды, длины и азимутов нормаляных сечений и длины геодезической линии эллипсоида. Применяются тригонометрические функции данных величин и постоянные принятого ргфгргнц-эллипсоида, что является целесообразным при использовании вычислителQjных машин, особенно автоматов. Рабочие формулы рекомендуются в форме (8–11).

---

---

Dedicated to Professor František Fiala on His 85th Birthday

---

Address: Veveři 95, Brno.     

Evolutionary Changes in the Configuration of the Magnetic Field of the Solar Corona in the Epoch of the Minimum

Geomagnetism and Aeronomy - The evolution of the magnetic configuration in the solar corona is studied. The curvature of the K-corona helmets obtained from an analytical approximation of the...