地震资料分频处理

在地震资料信噪比较低、叠加次数较低的情况下,常规处理效果不好,可以采用分频处理。根据频谱分析及能量扫描,把地震记录分成几个频带,分别进行自动剩余静校正、互相关静校正、地表一致性振幅补偿、反Q滤波、消除噪音、自适应倾斜面元等处理。然后再将各频带的地震反射波能量加起来做最终修饰,获得最终叠加剖面。这样可以充分利用各个频带的反射波能量,提高处理质量。     

九江-瑞昌5.7级地震余震的S波分裂研究

利用在九江-瑞昌MS5.7地震震中附近架设的丁家山台(DJS)、狮子洞台(SZD)和武蛟台(WUJ)的地震波形资料,采用S波分裂系统分析方法,对余震进行了S波分裂分析。结果表明,震中距较小、台站附近断层分布复杂的丁家山台(DJS)的慢波时间延迟相对较大,快波偏振方向不太集中。台站附近断裂分布单一的武蛟台(WUJ)快波偏振优势方向与断裂走向角度相差约35°,与区域主压应力方向也不一致。台站附近无断层通过的狮子洞台(SZD)快波偏振优势方向接近区域主压应力方向。慢波时间延迟大小与震源深度没有明显的规律性关系。     

Inversion of multicomponent seismic time shifts for reservoir pressure and length: a feasibility study

Pressure drops associated with reservoir production generate excess stress and strain that cause travel‐time shifts of reflected waves. Here, we invert time shifts of P‐, S‐, and PS‐waves measured between baseline and monitor surveys for pressure reduction and reservoir length. The inversion results can be used to estimate compaction‐induced stress and strain changes around the reservoir. We implement a hybrid inversion algorithm that incorporates elements of gradient, global/genetic, and nearest neighbour methods and permits exploration of the parameter space while simultaneously following local misfit gradients. Our synthetic examples indicate that optimal estimates of reservoir pressure from P‐wave data can be obtained using the reflections from the reservoir top. For S‐waves, time shifts from the top of the reservoir can be accurately inverted for pressure if the noise level is low. However, if noise contamination is significant, it is preferable to use S‐wave data (or combined shifts of all three modes) from reflectors beneath the reservoir. Joint wave type inversions demonstrate improvements over any single pure mode. Reservoir length can be estimated using the time shifts of any mode from the reservoir top or deeper reflectors. We also evaluate the differences between the actual strain field and those corresponding to the best‐case inversion results obtained using P‐ and S‐wave data. Another series of tests addresses the inversion of the time shifts for the pressure drops in two‐compartment reservoirs, as well as for the associated strain field. Numerical testing shows that a potentially serious source of error in the inversion is a distortion in the strain‐sensitivity coefficients, which govern the magnitude of stiffness changes. This feasibility study suggests which wave types and reflector locations may provide the most accurate estimates of reservoir parameters from compaction‐induced time shifts.     

亚暴期间等离子体对流速度及热层纬向风速度变化的统计学分析

本文利用2001年至2005年间CHAMP卫星及DMSP(F13,F15)卫星观测数据,对亚暴发生前后,亚暴初始位置所在磁地方时(Magnetic Local Time, MLT)东侧2 h(+2 h,MLT)至西侧4 h(-4 h,MLT)范围内等离子体对流速度(Vy)及热层纬向风速度(Uy)的变化进行了统计学分析.研究发现在亚暴发生后的1.5 h内,所有MLT区间的Vy均明显增大,且峰值位置的地磁纬度向赤道侧移动,1.5 h后,Vy减小,峰值的纬度向极区移动,表明亚暴的发生能显著增强等离子体对流速度;在亚暴发生位置的西侧0~2 h内,Vy增幅最大,这表明亚暴对热层的影响主要在西侧,影响最大的区域是西侧0~2 h MLT区间;Uy在亚暴发生后3 h持续增大,其对亚暴的响应相较于Vy有1.5 h的延迟.     

Application of a cascading filter implemented using morphological filtering and time–frequency peak filtering for seismic signal enhancement

Inspired by the idea of the iterative time–frequency peak filtering, which applies time–frequency peak filtering several times to improve the ability of random noise reduction, this article proposes a new cascading filter implemented using mathematic morphological filtering and the time–frequency peak filtering, which we call here morphological time–frequency peak filtering for convenience. This new method will be used mainly for seismic signal enhancement and random noise reduction in which the advantages of the morphological algorithm in processing nonlinear signals and the time–frequency peak filtering in processing nonstationary signals are utilized. Structurally, the scheme of the proposed method adopts mathematic morphological operation to first preprocess the signal and then applies the time–frequency peak filtering method to ultimately extract the valid signal. Through experiments on synthetic seismic signals and field seismic data, this paper demonstrates that the morphological time–frequency peak filtering method is superior to the time–frequency peak filtering method and its iterative form in terms of valid signal enhancement and random noise reduction.     

基于全波形采样的激电多信息提取方法研究与应用

激电法是地质勘查中的一种重要方法,时、频域测量各具优势,但传统的时、频域测量实现方法相对独立.根据频域和时域可相互转换的理论,本文提出一种基于全波形采样的时频激电多参数提取的地球物理数据处理方法.采取时间域激电法的观测模式,接收机高精度同步记录整周期电压-电流全部采样点的波形数据,通过该文提出的处理方法,仅一次供电和测量即可提取多种时、频域激电参数:时域激电的视电阻率、视极化率和频域激电的多频视电阻率、视相位、视频散率、去耦后的视相位等参数.在云南保山某典型矿区开展方法的试验应用,提取了地质体大量的相关激电信息,总结了目标体的激电参数响应规律.应用结果表明,该方法便捷高效,抗干扰能力强,相较于传统的时域激电测量,多种参数组合可提升对目标地质体的认识,增强激电法的应用效果.     

基于背景噪声互相关技术的雅浦俯冲带海底地震仪时间校正

由于电磁波无法在水中传播,导致海底地震仪(OBS)与全球定位系统(GPS)无法进行实时对钟,因此会产生一定的时钟漂移.而受船时、海况等影响,被动源OBS在回收后往往电量亏损,致使设备回收时并没有完成对钟,无法直接利用对钟前后的时钟漂移量进行线性校正.本文选取中国科学院海洋研究所2016—2017年间在雅浦俯冲带进行的长周期被动源OBS实验所获得的数据,利用背景噪声互相关方法对OBS进行了时间校正.结果发现,各OBS在经过一年左右时间的布放后,时钟漂移量在1～3 s之间,且基本呈线性趋势,将该结果平均分配到每天从而完成线性校正.另外,探讨了时钟的非线性漂移、误差源、参考台站的选取等因素对时钟准确度的影响,为OBS的数据处理进行了有益的探索.     

利用剪切波分裂研究青藏高原东北缘及东南缘上地幔各向异性

利用从IRIS上下载的青藏高原东北缘238个台站以及中国地震科学探测台阵喜马拉雅一期350个台站记录到的远震波形数据,通过采用剪切波分裂方法,获取各个台站下方各向异性分裂参数——快波偏振方向（φ）和分裂时差（δt）,从而得到青藏高原东北缘、东南缘上地幔的各向异性特征。研究结果表明,祁连块体、阿拉善块体和鄂尔多斯块体北部,快波方向为NNW-SSE,明显不同于GPS测量得到的近NE-SW的地表位移场方向,延迟时间平均~0.85 s;羌塘块体以及松潘-甘孜褶皱带的西部,快波方向呈现沿顺时针方向旋转的趋势,并且与GPS测量得到的地表位移场方向一致,延迟时间平均为~1.24 s;松潘-甘孜褶皱带东部、秦岭造山带与鄂尔多斯块体南部的交界处,快波方向呈现无序分布,与GPS方向表现出不一致的分布,延迟时间平均~1.08 s;川滇块体北部,快波方向近似N-S方向,与GPS测量得到的地表位移场方向相同,平均延迟时间为~0.925 s;位于北纬26°以南的川滇块体南部,快波方向近E-W方向分布,明显不同于GPS地表位移场方向,平均延迟时间~1.065 s。综合分析推测,羌塘块体、松潘-甘孜褶皱带的西部以及川滇块体北部,地表形变和深部之间的变形是相互耦合的;祁连块体、阿拉善块体和鄂尔多斯块体北部,松潘-甘孜褶皱带东部、秦岭造山带与鄂尔多斯块体南部的交界处以及川滇块体南部,壳幔之间可能存在解耦现象。     

Performance of a flux splitting scheme when solving the single-phase pressure equation discretized by MPFA

In this article we present a series of tests to study how well suited the TPFA coefficient matrix is as a preconditioner for the MPFA discrete system of equations in an iterative solver, using a flux splitting method. These tests have been conducted for single-phase flow for a wide range of anisotropy, heterogeneity, and grid skewness (mainly parallelogram grids). We use the K-orthogonal part of the MPFA transmissibilities for a parallelogram grid to govern the TPFA transmissibilities. The convergence of the flux splitting method is for each test case measured by the spectral radius of the iteration matrix.