非负矩阵分解和Curvelet在遥感图像融合中的应用

非负矩阵分解是一种提取图像原始信息局部特征的新方法,第二代Curvelet变换是一种效果较好的多尺度变换分析方法。结合两者特征提出一种基于NMF和Curvelet的遥感图像的融合方法,首先对已配准的多光谱图像和全色图像进行Curevelet分解,得到各层系数（Coarse、Detail和Fine尺度层）。然后对Coarse尺度层（低频系数）进行NMF分解,提取出包含特征基的低频系数;对Detail和Fine尺度层（高频系数）采用方差为测度参数进行邻域融合。最后进行Curevelet逆变换得到融合图像。实验结果表明,该方法的融合图像能较好地保留光谱信息,并在空间细节信息上得到改善,优于小波方法、Curvelet等方法。     

希尔伯特-黄变换技术及在边界层湍流研究中的应用

大气边界层中的湍流运动可以视为由许多不同时空尺度的湍涡叠加形成的,谱分析方法可以将这些不同尺度的涡旋分离开来。在介绍希尔伯特-黄变换（HHT）的基础上,选取2011年4月16日00—01时（北京时）、05—06时、12—13时不同时段的大气湍流观测数据,采用希尔伯特-黄变换技术分析了内蒙古科尔沁沙地地区的谱分布特征,得到的各内在模函数对应的中心瞬时频率可以表现出不同尺度上的湍流波动,其中,包含周期为3—5和10—13 min的阵风和重力内波,并在希尔伯特边际谱中体现。三维希尔伯特谱分布表明湍流能量多集中在频率小于0.5 Hz的低频段;大气湍流脉动较强的时段,其内在模函数和希尔伯特谱也对应着较大的能量。对比传统的快速傅里叶变换（FFT）方法,希尔伯特-黄变换可以更清晰地表现出湍流的频谱分布。     

基于二进小波的数字形变资料趋势异常分析探讨

利用二进正交离散小波，通过Mallat快速算法对形变倾斜资料进行了处理，并将计算结果与该地区的地震活动进行了对比分析。综合分析认为，小波分解的本质是窗口加权滑动平均的结果，能有效地对各种不同频率成分的形变序列按不同的频段分离和重构，可以起到明显放大低频段趋势成份的作用，为有效识别形变趋势异常提供一种可行方法。     

超道相干体技术在复杂断块断裂解释中的应用

针对复杂断块的断裂地震特征,笔者利用超道构建技术,并与基于小波变换的相干技术集成,再进行相干计算。超道构建技术具有不改变地层的原有信息而保留了地层的倾角特征、降低属性分析的平均效应的优势,试算结果证明,该方法不仅计算速度快,而且抑制噪声能力强。通过与传统相干算法实际应用对比,该算法有效地突出了地层的高连续性特征,提高了识别断层的分辨率,地质特征信息丰富,有利于断层解释,尤其适用于低信噪比地震资料地区。     

Dipping-interface Mapping Using Mode-separated Rayleigh Waves

Multichannel analysis of surface waves (MASW) method is a non-invasive geophysical technique that uses the dispersive characteristic of Rayleigh waves to estimate a vertical shear (S)-wave velocity profile. A pseudo-2D S-wave velocity section is constructed by aligning 1D S-wave velocity profiles at the midpoint of each receiver spread that are contoured using a spatial interpolation scheme. The horizontal resolution of the section is therefore most influenced by the receiver spread length and the source interval. Based on the assumption that a dipping-layer model can be regarded as stepped flat layers, high-resolution linear Radon transform (LRT) has been proposed to image Rayleigh-wave dispersive energy and separate modes of Rayleigh waves from a multichannel record. With the mode-separation technique, therefore, a dispersion curve that possesses satisfactory accuracy can be calculated using a pair of consecutive traces within a mode-separated shot gather. In this study, using synthetic models containing a dipping layer with a slope of 5, 10, 15, 20, or 30 degrees and a real-world example, we assess the ability of using high-resolution LRT to image and separate fundamental-mode Rayleigh waves from raw surface-wave data and accuracy of dispersion curves generated by a pair of consecutive traces within a mode-separated shot gather. Results of synthetic and real-world examples demonstrate that a dipping interface with a slope smaller than 15 degrees can be successfully mapped by separated fundamental waves using high-resolution LRT.     

无限弹性介质中隧道内爆炸动力响应的解析解

基于圆柱形薄壳结构理论,建立爆炸荷载作用下隧道结构的动力平衡方程.采用积分变换法求得Fourier变换域内爆炸荷栽作用下隧道位移场的积分表达式.结合具体算例,采用离散Fourier逆变换的数值计算方法,给出了爆心处隧道位移时程曲线、最大正应力时程曲线及某一时刻隧道位移幅值和最大正应力在x方向的变化曲线.其结果与数值计算结果相一致,验证了方法的可行性.     

Analysis of earthquake ground motions using an improved Hilbert–Huang transform

Some limitations of the Hilbert–Huang transform (HHT) for nonlinear and nonstationary signal processing are remarked. As an enhancement to the HHT, a time varying vector autoregressive moving average (VARMA) model based method is proposed to calculate the instantaneous frequencies of the intrinsic mode functions (IMFs) obtained from the empirical mode decomposition (EMD) of a signal. By representing the IMFs as time varying VARMA model and using the Kalman filter to estimate the time varying model parameters, the instantaneous frequencies are calculated according to the time varying parameters, then the instantaneous frequencies and the envelopes derived from the cubic spline interpolation of the maxima of IMFs are used to yield the Hilbert spectrum. The analysis of the length of day dataset and the ground motion record El Centro (1940, N–S) shows that the proposed method offers advantages in frequency resolution, and produces more physically meaningful and readable Hilbert spectrum than the original HHT method, short-time Fourier transform (STFT) and wavelet transform (WT). The analysis of the seismic response of a building during the 1994 Northridge earthquake shows that the proposed method is a powerful tool for structural damage detection, which is expected as the promising area for future research.     

Effect of wellbore storage and finite thickness skin on flow to a partially penetrating well in a phreatic aquifer

An analytical model is presented for the analysis of constant flux tests conducted in a phreatic aquifer having a partially penetrating well with a finite thickness skin. The solution is derived in the Laplace transform domain for the drawdown in the pumping well, skin and formation regions. The time-domain solution in terms of the aquifer drawdown is then obtained from the numerical inversion of the Laplace transform and presented as dimensionless drawdown–time curves. The derived solution is used to investigate the effects of the hydraulic conductivity contrast between the skin and formation, in addition to wellbore storage, skin thickness, delayed yield, partial penetration and distance to the observation well. The results of the developed solution were compared with those from an existing solution for the case of an infinitesimally thin skin. The latter solution can never approximate that for the developed finite skin. Dimensionless drawdown–time curves were compared with the other published results for a confined aquifer. Positive skin effects are reflected in the early time and disappear in the intermediate and late time aquifer responses. But in the case of negative skin this is reversed and the negative skin also tends to disguise the wellbore storage effect. A thick negative skin lowers the overall drawdown in the aquifer and leads to more persistent delayed drainage. Partial penetration increases the drawdown in the case of a positive skin; however its effect is masked by the negative skin. The influence of a negative skin is pronounced over a broad range of radial distances. At distant observation points the influence of a positive skin is too small to be reflected in early and intermediate time pumping test data and consequently the type curve takes its asymptotic form.     

小波变换方法在辽宁跨断层形变资料分析中的应用

运用二进小波变换方法,分析了辽宁省跨断层观测资料不同尺度的小波变换的细节信号变换特征。研究结果表明,不同尺度小波变换分解的细节信号显示出不同的时变特征,可以识别和提取原始测值中不易分辨的震前异常。当尺度为3时的细节信号的年周期特征被打破,或者当尺度为3时的细节信号变化超出2倍均方差时,测点周围250km左右范围内发生中强地震的可能性很大。     

Radon变换在非规则观测系统中波场分离的应用

Radon变换是数据处理中广泛应用的一种方法技术.本文介绍了Radon变换在规则及其非规则观测系统中波场分离中应用,讨论了Radon变换的公式、实现方法和有关计算参数的选择.通过理论模型的试算,阐明了Radon变换在缺失道、道距变化情况下依然可以高精度的分离波场,并对线性、抛物Radon变换进行了比较,对隧道超前预报反射波资料进行了处理,给出了隧道超前预报处理流程中提取反射波的算例.