数学形态滤波与大地电磁噪声压制

数学形态滤波是一种新型的非线性滤波方法.介绍了形态学滤波的基本原理,针对大地电磁信号表现出的非线性、非平稳性和非最小相位特性,综合结构元素特征我们提出一种基于数学形态滤波的大地电磁噪声压制方法.为了有效抑制目标信号中的噪声干扰及修正统计偏倚现象,通过选用合理的结构元素及形态开-闭、闭-开组合,将正、负结构元素级联构造组合广义形态滤波器对实测大地电磁信号进行噪声压制.实验结果表明,该方法是切实可行的,有效地剔除了大尺度干扰及基线漂移,较好地还原了大地电磁信号的原始特征,修正了标准形态算子所产生的统计偏倚现象,去噪精度高.该方法计算速度快,具有潜在优势,为矿集区海量大地电磁信号与强干扰的分离提供了一条新的途径,应用前景广阔.     

漳州地震台阵背景噪声的特性研究

对2007年10~12月在福建省漳州地区进行的勘址观测数据进行了分析,采用了均方根、功率谱、相干函数等方法对漳州地区的地震背景噪声水平进行了研究。研究结果表明,漳州台阵内环半径参考值为510m,外环半径参考值为1700m。滨临海边的台阵勘址测点噪声水平偏高,并在1.5、5.0Hz处发现明显噪声源。1.5Hz噪声源主要来自水库地区的大风引起的水波,5.0Hz噪声源来自台阵场地东北方向的发电厂。     

Interference effects of aircraft on earth's electromagnetic response at very low frequency and low frequency

Over the last few decades, very low frequency electromagnetics has been widely and successfully applied in mineral exploration and groundwater exploration. Many radio transmitters with strong signal‐to‐noise ratios are scattered in the very low frequency band and low frequency band. Based on experiences gained from ground measurements with the radio‐magnetotelluric technique operating in the frequency interval 1–250 kHz, broadband magnetometers have been used to cover both very low frequency (3–30 kHz) and low frequency (30–300 kHz) bands to increase the resolution of the near‐surface structure. The metallic aircraft as a conductive body will distort the magnetic signal to some extent, and thus it is important to investigate aircraft interference on the electromagnetic signal. We studied noise caused by rotation of an aircraft and the aircraft itself as a metallic conductive body with three methods: 3D wave polarization, determination of transmitter direction and full tipper estimation. Both very low frequency and low frequency bands were investigated. The results show that the magnetic field is independent of the aircraft at low frequencies in the very low frequency band and part of the low frequency band (below 100 kHz). At high frequencies (above 100 kHz), the signals are more greatly influenced by the aircraft, and the wave polarization directions are more scattered, as observed when the aircraft turned. Some aircraft generated noise mixed with radio transmitter signals, detected as ‘dummy’ signals by the 3D wave polarization method. The estimated scalar magnetic transfer functions are dependent on the aircraft flight directions at high frequencies, because of aircraft interference. The aircraft eigenresponse in the transfer functions (tippers) between vertical and horizontal magnetic field components was compensated for in the real part of the estimated tippers, but some unknown effect was still observed in the imaginary parts.     

Eliminating 3D pre‐stack migration artefacts by 6D filtering

State‐of‐the‐art 3D seismic acquisition geometries have poor sampling along at least one dimension. This results in coherent migration noise that always contaminates pre‐stack migrated data, including high‐fold surveys, if prior‐to‐migration interpolation was not applied. We present a method for effective noise suppression in migrated gathers, competing with data interpolation before pre‐stack migration. The proposed technique is based on a dip decomposition of common‐offset volumes and a semblance‐type measure computation via offset for all constant‐dip gathers. Thus the processing engages six dimensions: offset, inline, crossline, depth, inline dip, and crossline dip. To reduce computational costs, we apply a two‐pass (4D in each pass) noise suppression: inline processing and then crossline processing (or vice versa). Synthetic and real‐data examples verify that the technique preserves signal amplitudes, including amplitude‐versus‐offset dependence, and that faults are not smeared.     

Non‐convex compressed sensing with frequency mask for seismic data reconstruction and denoising

Compressed Sensing has recently proved itself as a successful tool to help address the challenges of acquisition and processing seismic data sets. Compressed sensing shows that the information contained in sparse signals can be recovered accurately from a small number of linear measurements using a sparsity‐promoting regularization. This paper investigates two aspects of compressed sensing in seismic exploration: (i) using a general non‐convex regularizer instead of the conventional one‐norm minimization for sparsity promotion and (ii) using a frequency mask to additionally subsample the acquired traces in the frequency‐space () domain. The proposed non‐convex regularizer has better sparse recovery performance compared with one‐norm minimization and the additional frequency mask allows us to incorporate a priori information about the events contained in the wavefields into the reconstruction. For example, (i) seismic data are band‐limited; therefore one can use only a partial set of frequency coefficients in the range of reflections band, where the signal‐to‐noise ratio is high and spatial aliasing is low, to reconstruct the original wavefield, and (ii) low‐frequency characteristics of the coherent ground rolls allow direct elimination of them during reconstruction by disregarding the corresponding frequency coefficients (usually bellow 10 Hz) via a frequency mask. The results of this paper show that some challenges of reconstruction and denoising in seismic exploration can be addressed under a unified formulation. It is illustrated numerically that the compressed sensing performance for seismic data interpolation is improved significantly when an additional coherent subsampling is performed in the domain compared with the domain case. Numerical experiments from both simulated and real field data are included to illustrate the effectiveness of the presented method.     

High thresholds for avoidance of sonar by free-ranging long-finned pilot whales (Globicephala melas)

The potential effects of exposing marine mammals to military sonar is a current concern. Dose–response relationships are useful for predicting potential environmental impacts of specific operations. To reveal behavioral response thresholds of exposure to sonar, we conducted 18 exposure/control approaches to 6 long-finned pilot whales. Source level and proximity of sonar transmitting one of two frequency bands (1–2 kHz and 6–7 kHz) were increased during exposure sessions. The 2-dimensional movement tracks were analyzed using a changepoint method to identify the avoidance response thresholds which were used to estimate dose–response relationships. No support for an effect of sonar frequency or previous exposures on the probability of response was found. Estimated response thresholds at which 50% of population show avoidance (SPL_max = 170 dB re 1 μPa, SEL_cum = 173 dB re 1 μPa² s) were higher than previously found for other cetaceans. The US Navy currently uses a generic dose–response relationship to predict the responses of cetaceans to naval active sonar, which has been found to underestimate behavioural impacts on killer whales and beaked whales. The navy curve appears to match more closely our results with long-finned pilot whales, though it might underestimate the probability of avoidance for pilot-whales at long distances from sonar sources.     

基于背景噪声研究云南地区面波速度非均匀性和方位各向异性

基于中国地震科学探测台阵项目一期在南北地震带南段架设的300多个地震台站,利用2011年10月至2012年9月的连续观测记录,采用基于背景噪声互相关函数的面波层析成像技术,研究了青藏高原东南缘的云南地区面波群速度和方位各向异性分布.结果显示,地壳的面波快波方向呈现近南北向,整体表现出围绕东喜马拉雅构造结顺时针旋转的趋势,和地表GPS速度场以及S波分裂的快波方位较为一致.小江断裂东西两侧的快波方位有一定差异.对反映深度大概在下地壳和上地幔顶部的长周期面波,快波方向从近南北向逐渐向北西向过渡,在菱形块体附近,26°N以南,快波方向和红河断裂的走向趋于一致;其他区域相比上地壳的快波方向也有较大变化,这种结果较为支持青藏高原东南缘的云南地区壳幔变形的解耦.在滇西南,澜沧江向东弧形展布区域,中下地壳快波方位呈现局部的圆周旋转趋势,结合该区地震分布规律和应力主轴方向,推测这种现象和块体挤压及旋转具有一定相关性.     

噪声对磁场向下延拓迭代法的计算误差影响分析

磁场向下延拓可以增强数据的分辨率.在现有的向下延拓方法中,迭代法具有较好的计算效果和较大的向下延拓高度.首先介绍了迭代法的基本原理,分析迭代法的特点;然后利用迭代法对不同噪声水平的模型数据和实测数据进行向下延拓计算,研究噪声对迭代法的计算误差影响.通过延拓结果和实际数据比较分析,客观评价了迭代法的计算精度、可靠性及适应性,对推广应用该方法具有实质性意义.     

自适应噪声消除器在地电场观测数据处理中的应用

在地电场和电阻率同场地观测的台站, 地电场观测必然会受到地电阻率观测人工供电的干扰, 尤其是当地电场观测的采样率提高到1次/s以上时, 这种干扰更是无可避免, 而且不易消除。 该文阐述了应用自适应噪声消除器对受种干扰的地电场观测数据进行处理, 消除观测数据中的干扰噪声, 从而获得地电场真实信号波形。 对实际观测数据的处理结果表明, 采用此方法可以很好的消除秒采样地电场观测数据中地电阻率供电产生的干扰。 这种方法还可用于消除其他观测方法中存在的有规律的干扰。     

基于谐波建模和自相关的磁共振信号消噪与提取方法研究

磁共振探测技术（Magnetic Resonance Sounding,MRS）以其无损、定量、直接等优势,被广泛应用于地下水调查、水文环境评价以及灾害水源预警等领域.在实际应用中,强工频谐波和随机噪声等严重影响MRS信号的质量,导致后续水文参数解释不准确.针对这一问题,提出谐波建模和自相关相结合的方法进行消噪以及信号特征参数提取.首先构建工频谐波模型,针对建模算法严重依赖工频基频精度的问题,采用自适应扫描方式搜索方案,大幅提高搜索准确度和速度;其次推导了MRS信号自相关表达式,提出了自相关参数提取的非线性拟合方法.仿真数据结果表明,建模消噪方法有效消除了工频谐波,信噪比平均提升了17.03 dB;自相关处理后,信噪比进一步提升了16.10 dB,初始振幅和弛豫时间参数提取结果的准确度比处理前分别提高了3.8倍和2.8倍.通过不同信噪比和弛豫时间的重复实验,得到当噪声水平小于200 nV和弛豫时间大于200 ms时,自相关参数提取具有较高的稳定性.最后,通过野外实测数据处理实验,进一步验证了联合消噪和参数提取方法的有效性.