关于频率电磁测深几个问题的探讨(八)——频率电磁测深与瞬变电磁测深的关系

从探测原理出发,说明人工源频率电磁测深与瞬变电磁测深的关系密切,这是一个问题的两个方面,可以互相转换。两种探测方法实质是一样的,具有共性;但毕竟从频率域和时间域两个侧面反映问题,又具有各自特性。两种方法的原理所依据的麦克斯韦方程,其正演公式可通过傅立叶变换相互转换;从实际探测所发收的信号特征,场区划分与探测深度等主要方面阐明两种探测方法异同,以加深对探测理论与方法的认识,结合实际择优应用。      

基于压缩感知的Curvelet域联合迭代地震数据重建

由于野外采集环境的限制,常常无法采集得到完整规则的野外地震数据,为了后续地震处理、解释工作的顺利进行,地震数据重建工作被广泛的研究.自压缩感知理论的提出,相继出现了基于该理论的多种迭代阈值方法,如CRSI方法（Curvelet Recovery by Sparsity-promoting Inversion method）、Bregman迭代阈值算法（the linearized Bregman method）等.CSRI方法利用地震波形在Curvelet的稀疏特性,通过一种基于最速下降的迭代算法在Curvelet变换域恢复出高信噪比地震数据,该迭代算法稳定,收敛,但其收敛速度慢.Bregman迭代阈值法与CRSI最大区别在于每次迭代时把上一次恢复结果中的阈值前所有能量都保留到本次恢复结果中,从而加快了收敛速度,但随着迭代的进行重构数据中噪声干扰越来越严重,导致最终恢复出的数据信噪比低.综合两种经典方法的优缺点,本文构造了一种新的联合迭代算法框架,在每次迭代中将CRSI和Bregman的恢复量加权并同时加回本次迭代结果中,从而加快了迭代初期的收敛速度,又避免了迭代后期噪声干扰的影响.合成数据和实际数据试算结果表明,我们提出的新方法不仅迭代快速收敛稳定,且能得到高信噪比的重建结果.     

高采样率四分量钻孔应变同震响应分析

In this paper, we adopt the observation data of two 10Hz YRY 4-component borehole strainmeters and one 100Hz FBS-3B broadband seismometer to study the responses of their observations to the April 11, 2012 Mw 8.6 earthquake off the west coast of northern Sumatra, and compare the response characteristics during the coseismic stage after preprocessing of the data. We analyze the spectrum＇s dynamic process during the coseismic stage with the S transformation method, and discuss the spectrum＇s details in different seismic phases. The S transformation results of the data of seismometer and borehole strainmeter at the same station show that the coseismic spectrum evolution of the two kinds of observation are similar to each other and can be cross-referenced in indentifying the wave signals, which proves that the borehole strainmeter measurements at high frequencies are reliable.     

2012年印尼8.6级地震应变地震波的Hilbert-Huang时频分析

本文通过端点效应压制的Hilbert-Huang变换, 对大同及沁源台布置的四分量钻孔应变仪记录的印尼8.6级地震激发的应变地震波形进行时频分析, 结果显示印尼8.6级地震的主震和8.2级余震的应变地震波序列各个震相具有不同的时频特征: ① 地震波到达之前的所谓“环境噪声”部分, 瞬时频率低, 瞬时振幅小; ② P波初至时, 高频成分突然增加, 振幅也随即增强; ③ S波到达时, 频率有所降低而振幅剧烈上升; ④ 面波到达时, 振幅进一步剧烈上升达到整个序列的极大值; ⑤ 尾波部分振幅逐渐降低, 但与噪声部分相比频率依然偏高, 振幅依然偏大。 本文也将应变地震波与地震仪记录的地震波进行对比, 虽然应变地震波与地震波波形和Fourier谱具有极高的相关系数, 但从Hilbert-Huang变换得到的边际谱上看, 应变地震波与地震波有显著的区别, 应变地震波比地震波记录的低频成分相对更多。 通过Hilbert谱, 有助于更好地了解非平稳信号的局部特征, 对于突变信号的地震波, Hilbert-Huang变换是一个较好的时频分析工具。     

傅里叶变换在深州井资料处理中的应用

深州井水位日常观测中,经常会因为仪器供电不稳、雷击等原因出现一些幅度很大的脉冲、突跳及数据缺测,严重影响了观测资料的内在质量。本文采用数字信号处理技术[1]对深州井水位资料进行分析处理,首先使用分段线性插值法对缺测数据进行插值处理,然后用快速傅里叶变换对水位数据进行频谱分析,绘制振幅谱,分析水位数据中的不同周期成分,选取恰当的截止频率进行滤波处理,去除了干扰,获得了井水位的变化动态,提高了资料质量。最后应用滤波后的水位数据做一些震例分析。     

基于基本多项式的GPS精密星历插值方法研究

采用基本多项式插值方法对GPS精密星历进行了插值,该方法利用Householder变换来避免对病态的法方程系数矩阵直接求逆而导致的插值计算错误,并对高阶基本多项式插值过程中出现的龙格现象进行了讨论,给出了多项式拟合区间长度、有效区间长度和基本插值多项式阶数的选取准则,最后对插值的结果进行比较分析。结果表明,对于15分钟采样率的精密星历,基本多项式插值精度达到毫米级,能够满足精密定位的需要。     

资源三号卫星不同空间分辨率图像斑状植被检测比较研究

黄河三角洲滨海区域大量分布着由柽柳、芦苇、翅碱蓬、二色补血草、白茅等等组成的斑状植被,多呈棋盘状规则分布,胜利油田地震勘探被认为是产生这一现象的原因之一。深入研究其分布格局并探讨其形成机理,对于区域湿地生态系统植被恢复至关重要。为此,首先就需要对斑状植被进行有效地检测和制图。本文评价了我国首颗民用高分辨率光学传输型立体测图卫星（资源三号卫星）不同空间分辨率图像检测黄河三角洲斑状植被的能力,结果表明：在利用霍夫曼圆检测方法下,前视、后视、正视全色图像和正视多光谱图像都能很好地检测斑状植被,空间分辨率高的图像有更高的检测精度。     

北斗非GEO卫星信号码并行捕获

2012年12月,北斗二代的接口控制文件正式公布。北斗B1I频点的非GEO卫星上调制20bit的NH码,因此每1ms信号可能存在比特跳变。借鉴传统GPS并行码相位捕获算法,本文研究一种基于补零FFT的并行码相位算法。该算法能够实现1-20ms相干积分,利用非相干累加提高信噪比,并采用一种对信号强度和采样率不敏感的判决变量。仿真表明,通过相干积分和非相干累加,该算法使用400ms数据可以捕获24dBHz的弱信号。     

三维密度界面的正反演研究和应用

重力位场的界面反演是位场处理解释中的重要问题.本文将基于快速傅里叶变换的频率域界面反演方法Parker-Oldenburg公式推广到物性可随深度变化的三维情况,得出了密度可以横向、纵向任意变化的重力界面正反演公式.该方法在计算时可以合理地选取地面下某一深度作为基准面以减小界面起伏,使迭代易于收敛.理论模型试验表明该方法反演精度高,收敛速度快,在密度界面反演中具有广泛的实用价值.最后利用该方法反演华北地区莫霍面的深度,反演结果得到了地震测深数据的验证.     

多次波压制技术在南海北部陆缘深水区的应用

南海北部洋陆过渡带水深变化剧烈,海底地形崎岖,构造复杂多变,地震勘探中多次波十分发育,常规处理难以取得很好的成像效果.由于方法的局限性,任何一种多次波压制方法都很难彻底地压制全部多次波.本文通过对国家基金委共享航次在南海北部陆缘采集的二维反射地震资料的分析,总结了不同多次波压制技术的优缺点和适用范围,结合滤波类方法和波动理论类方法,综合应用了表面相关多次波衰减法(SRME,Surface Related Multiple Elimination)、τ-p域预测反褶积、高分辨率抛物线Radon变换以及叠后压制残余多次波等方法去除多次波,有效突出了一次反射波,证明了方法的有效性和实用性,同时也形成了针对南海深水资料多次波压制的一种有效的处理流程和方法组合.