EMD-Wavelet-ICA耦合模型及其在GPS坐标序列降噪中的应用

在分析比较经验模态分解(EMD)、小波变换(Wavelet)和独立分量分析(ICA)优缺点的基础上,提出一种新的EMD-Wavelet-ICA耦合模型。该模型充分利用了EMD的自适应性,对原始信号进行分解获得不同频率的模态函数(IMF),采用标准化模量的累计均值对IMF进行尺度划分;进而分别采用Wavelet和ICA对高频和低频IMF进行降噪,将降噪后的IMF进行多尺度重构,获得降噪后的信号;采用信噪比、标准差、偏差和相关系数等指标对降噪效果进行评价。仿真数据和GPS坐标序列的处理结果表明:与EMD模型和EMD-ICA模型相比,新模型的标准差、偏差均有不同程度的减小;信噪比和相关系数有一定程度的增大,可以获得更好的降噪效果。     

Impact of sub-pixel heterogeneity on modelled brightness temperature for an agricultural region

Knowledge of sub-pixel heterogeneity, particularly at the passive microwave scale, can improve the brightness temperature (and ultimately the soil moisture) estimation. However, the impact of surface heterogeneity (in terms of soil moisture, soil temperature and vegetation water content) on brightness temperature in an agricultural setting is relatively unknown. The Soil Moisture Active Passive Validation Experiment 2012 (SMAPVEX12) provided an opportunity to evaluate sub-pixel heterogeneity at the scale of a Soil Moisture Ocean Salinity (SMOS) or the Soil Moisture Active Passive (SMAP) radiometer footprint using field measured data. The first objective of this study was to determine if accounting for surface heterogeneity reduced the error between estimated brightness temperature (T_b) and T_b measured by SMOS. It was found that when accounting for variation in surface soil moisture, temperature and vegetation water content within the pixel footprint, the error between the modelled T_b and the measured T_b was less than if a homogeneous pixel were modelled. The correlation between the surface parameters and the error associated with not accounting for surface heterogeneity were investigated. It was found that there was low to moderate correlation between the error and the coefficient of variance associated with the measured soil moisture, soil temperature and vegetation volumetric water content during the field campaign. However, it was found that the correlations changed depending on the stage of vegetation growth and the amount of time following a precipitation event. At the start of the field campaign (following a precipitation event), there was strong correlation between the error and all three surface parameters (r ≥ 0.75). Following a precipitation event close to the middle of the field campaign (during which there was rapid growth in vegetation), there was strong correlation between the error and the variability in vegetation water content (r = 0.89), moderate correlation with soil moisture (r = 0.61) and low correlation with soil temperature (r = 0.26).     

多尺度细节增强的遥感影像超分辨率重建

鉴于现有超分辨率重建方法难以突显重建影像细节信息的问题,提出多尺度细节增强的遥感影像超分辨率重建模型框架。首先,通过最小二乘滤波方法将序列影像分解成包含大尺度边缘的平滑信息和包含中小型尺度的细节信息;其次,利用插值方法得到相应的高分辨率细节信息和平滑信息,构造纹理细节增强函数,提升中小型细节的增强幅度;最终,融合细节信息和平滑信息,得到初始的超分辨率重建结果,并利用局部优化模型进一步改善重建影像质量。选取同时相和多时相遥感影像作为试验数据。试验结果表明,本文重建结果与插值方法、TV方法和MAP方法相比,在客观评价指标上均有显著提高,明显改善了重建影像的纹理细节。论文提出的多尺度细节增强的超分辨率重建方法,可以使重建影像提供更多高频细节信息,具有较好鲁棒性和普适性。     

结合像元分解和STARFM模型的遥感数据融合

高空间、时间分辨率遥感数据在监测地表快速变化方面具有重要的作用。然而,对于特定传感器获取的遥感影像在空间分辨率和时间分辨率上存在不可调和的矛盾,遥感数据时空融合技术是解决这一矛盾的有效方法。本文利用像元分解降尺方法(Downscaling mixed pixel)和STARFM模型(Spatial and Temporal Adaptive Reflectance Fusion Model)相结合的CDSTARFM算法(Combination of Downscaling Mixed Pixel Algorithm and Spatial and Temporal Adaptive Reflectance Fusion Model)进行遥感数据融合。首先,利用像元分解降尺度方法对参与融合的MODIS数据进行分解降尺度处理;其次,利用分解降尺度的MODIS数据替代STARFM模型中直接重采样的MODIS数据进行数据融合;最后以Landsat 8和MODIS遥感影像数据对该方法进行了实验。结果表明:(1)CDSTARFM算法比STARFM和像元分解降尺度算法具有更高的融合精度;(2)CDSTARFM能够在较小的窗口下获得更高的融合精度,在相同的窗口下其融合精度也高于STARFM;(3)CDSTARFM融合的影像更接近真实影像,消除了像元分解降尺度影像中的"图斑"和STARFM模型融合影像中的"MODIS像元边界"。     

Seismic attenuation estimation using a complete ensemble empirical mode decomposition-based method

The frequency attenuation gradient method can provide important information for hydrocarbon detection. In this paper, a method using Complete Ensemble Empirical Mode Decomposition (CEEMD), Hilbert transform and the least-squares curve-fitting is proposed for seismic attenuation estimation as an effective frequency attenuation gradient estimation approach. We first use CEEMD to obtain the different Intrinsic Mode Functions (IMFs), which have a narrow band and can enhance the physical meaning of instantaneous attributes trace by trace. The time-frequency spectrum, which is computed using a Hilbert transform of each IMF, is represented as a spectrum with a single-peak that has narrow side lobes, which is conducive to frequency attenuation gradient estimation. Second, for each time sample, the frequency-amplitude spectrum of each IMF trace is extracted from the time-frequency spectrum to conduct the attenuation gradient computation. Then, the logarithm operation is performed for each IMF trace. Due to the very narrow bands of some IMFs in some seismic traces, a variable frequency window is adopted along the IMF trace according to the local data characteristics. Finally, the attenuation gradient for each IMF in a seismic trace can be computed using least-squares fitting. A different IMF reflects a seismic trace with a different spatiotemporal scale and can highlight different geologic and stratigraphic information. The correlation weighted average operation is used to highlight some useful details in seismic trace and obtains the attenuation gradient for each seismic trace. Field data examples demonstrate our method and its effectiveness. The proposed method can stably estimate the frequency attenuation gradient.     

Application of multi-seismic attributes analysis in the study of distributary channels

During the exploration of stratigraphic reservoirs, the key to locating these reservoirs is to identify the sandstone distribution. Seismic data can be used to recognize large-scale distributary-channel sedimentary bodies; however, depiction of sedimentary bodies in small distributary channels using traditional interpretation data from seismic profiles is extremely difficult. In the Upper Cretaceous first sandstone bed of the third member of the Nenjiang Formation (called the Nen 3 member for simplicity) in the XB area of the southern Songliao basin, distributary channels are the main sedimentary facies. The distributary channels migrate frequently; therefore, it is difficult to precisely depict the morphology of the channel and the distribution of sandstone. In this study, we investigated the deposition of distributary channels using equal-scale formation plane seismic attributes such as post-stack amplitude, instantaneous amplitude and seismic waveform classification, analyzed well logging data from target formations, and mapped the distribution of distributary channels. This study shows that seismic sedimentological study and stacking multiple attribute analysis can identify the distribution of distributary channels precisely and effectively. The attributes of equal-scale strata slices are evidently superior to those of time and horizon slices; Attribute extraction and selection of equal-scale formation slices is an essential step. The comprehensive selection of seismic attributes that show a good correlation to a single-well can be used to clearly depict channel bodies. The overlap of the 40-Hz single-frequency energy and the RMS amplitude depicts the sedimentary characteristics and shapes of sandstone bodies in main channels as well as small-scale distributary channels more clearly than a single attribute. These attributes show that the first sandstone bed of the Nen 3 member was from distributary channels with complex shapes. The distributary channels cut across each other vertically and laterally. It indicates that, during the deposition of the Nen 3 member in XB area, the direction of the main channel is from north to south, and the provenance could be from north and east. The method of this paper may provide some helpful suggestions to the geologist using seismic attributes to undertake research on sedimentary environment at other places in the world.     

频谱分解技术在泥页岩储层含气性预测中的应用——以柴达木盆地鱼卡凹陷为例

针对柴达木盆地北缘鱼卡凹陷中侏罗统大煤沟组泥页岩,利用地震频谱曲线在不同流体介质储层中的形态、斜率、主频与中心频率的差值等属性参数特征,识别陆相泥页岩层系的地震频谱属性特征,预测含气性。研究结果表明,中侏罗统大煤沟组上部能量系数较低,平均为0.6719,下部能量系数较高,平均为0.7521,结合现场含气量和全烃特征认为,可以以能量系数0.7为界线,作为识别下部和上部含气性的标志,下部明显比上部含气性好。平面上进一步圈定了2个有利目标区,得到柴页1井勘探验证。因此,频谱分解技术可应用于中国陆相页岩气的勘查开发。     

微波消解-电感耦合等离子体质谱法测定海泡石中的稀土元素

海泡石是一种纤维状含水的富镁硅酸盐黏土矿,其中的稀土元素含量在1×10~(-7)～1×10~(-5)之间,目前还没有建立海泡石中稀土元素的国家标准分析方法。测定岩石中的稀土元素主要是采用电感耦合等离子体质谱法(ICP-MS),样品前处理一般采用封闭溶矿和碱熔,但这两种处理方法耗时较长,效率不高。本文通过比较硝酸-氢氟酸-过氧化氢、硝酸-氢氟酸、硝酸-过氧化氢三种样品前处理方法,确定使用硝酸-氢氟酸溶矿,然后进行微波消解同时赶去氢氟酸,避免氢氟酸与稀土元素生成难溶的氟化物,再采用ICP-MS法测定15种稀土元素的含量。由于海泡石中的镁含量较高,为降低基体效应,以~(103)Rh和~(185)Re作内标补偿基体效应和校正灵敏度漂移,各元素测定值的准确性显著提高,回收率为91. 2%～110. 9%,检出限为0. 002～0. 011μg/L,精密度≤2. 79%。本方法与封闭酸溶ICP-MS法的分析结果吻合较好,且用酸量少(7 mL),溶矿效率高(1 h),检出限更低。     

2000-2014年西藏雅鲁藏布江流域积雪时空变化分析及对气候的响应研究

利用2000-2014年MODIS逐日无云积雪产品对雅鲁藏布江流域积雪特征的空间分布及变化、积雪随高程变化的规律进行了分析,并采用被动微波数据SMMR （1979-1987年）和SSM/I （1988-2008年）以及中国地面降水和气温0.5°×0.5°日值格点数据集,研究了雅鲁藏布江流域关键积雪参数对气候要素的响应等。结果表明：流域下游积雪日较大且变化剧烈;流域整体上呈显著减少的趋势;积雪日随高程的上升而增加;流域内降水呈不显著的增加趋势,而气温呈显著的增加趋势,最高气温对积雪变化影响最大;气温对积雪终日的影响明显高于积雪初日;在积雪消融期降水的增多促进了积雪的消融。     

微波辐射计在对流天气预报中的应用

选取2016—2018年每年4—9月份RPG-HATPRO型42通道微波辐射计观测的不稳定指数参数(K、SI、CAPE、LI)及水汽参数(IWV、LWP),研究得出各参数触发雷雨大风、短时强降水的阈值条件为K>37℃、SI<-1℃、IWV>60 kg/m~2、LWP>400 g/m~2,而LI、CAPE无法对3种天气类型进行区分。利用费舍判别分析方法,将不稳定指数参数及水汽参数作为预报因子,建立预报方程并进行检验,结果表明:二级判别方程预测对流天气的准确率为76%,可以作为预报对流天气的辅助工具;多级判别方程不能很好地区分3种天气类型,但将其作为修正后的二级判别方程使用,能提高对流天气的测中概率。