地震数据实时自适应压缩方法研究

本文探讨了二种地震数据实时自适应压缩方法。自适应采样是我们在总结传统的采样方法的基础上,结合地震信号频率变化特点,对地震信号所采取的一种波形压缩方法。自适应增量编码ADPCM法则在自适应采样压缩的基础上将地震数据又压缩了二倍,进一步提高了压缩效果。两种方法有力配合,使短、中长、长周期地震记录分别压缩了12倍、30倍、10倍以上。压缩恢复波形无论在时域分析(相关系数达0.997以上)还是在频域分析,失真都非常小,基本达到无畸变压缩。上述两种方法算法简单,计算量小,易于地震数据的高速采集、传输及记录。自适应采样法目前已应用于国家地震局地球物理所六室研制的BSO系统中。      

最小二乘反演迭代算法在压制海上变深度缆采集数据虚反射中的应用

海水面的虚反射(鬼波)引起海上拖缆采集数据陷波,导致地震记录频带变窄,而近年发展的变深度缆采集技术,具有多样的陷波特征,通过专门的去虚反射处理方法可获得宽频数据.本文基于已有研究成果,将最小二乘反演迭代压制虚反射算法应用于某海上变深度缆宽频处理.基于频率波数域镜像记录生成方法获得镜像炮集记录,并采用最小二乘解从变深度缆原始和镜像炮集记录中提取上行波.针对镜像炮集记录生成受初始速度模型精度的影响,使得某深度缆接收的上行波和下行波之间的实际延迟时间存在误差,采用最小二乘反演迭代算法最优化计算下行波与上行波之间的平均延迟时间和上行波记录,并采用时空数据窗口滑动克服延迟时间随炮检距和目的层深度变化问题.合成数据及某海上实际变深度缆数据处理测试结果表明,该方法能较好地压制变深度缆由海水面产生的虚反射,能达到拓宽地震记录频带目的.     

地震台阵上信号方位角和慢度的时、频域估计方法比较

地震台阵常用频率-波数分析法来估计地震信号的后方位角和慢度。尽管有关的算法是众所周知的,但不同的实现在某些情况下可能导致不同的结果。如应用在宽频带的澳大利亚WRA台阵的记录时,标准的f-k分析方法往往给出不正确的结果。本文发现错误是由FFT的频谱泄漏效应引起的,如在进行FFT之前对原始数据进行高通滤波则可以有效地避免这样的错误。进一步对在时域中直接估算信号方位角和慢度的一种算法进行了分析,比较发现时域方法的计算速度和精度都不低于频域方法,且在某些低信噪比的情况下,前者可以给出更可靠的结果。     

频率慢度域自适应迭代反演算法压制海上倾斜缆鬼波方法及其应用

海上倾斜缆采集技术具有多样的陷波特征,通过去鬼波处理可获得宽频数据.针对海水面波浪起伏及缆深误差引起的鬼波延迟时间估计误差以及崎岖海底和目的层深度变化使得鬼波和一次反射波的振幅差异系数随偏移距的变化而难以给定一个固定值的问题,本文推导出频率慢度域中鬼波滤波算子以及自适应迭代反演求解上行波算法,该鬼波滤波算子与不同水平慢度对应的鬼波和一次反射波的振幅差异系数以及鬼波延迟时间有关.并基于计算出的理论下行波与实际下行波之间的平方误差最小理论实现自适应反演迭代最优计算该振幅差异系数和鬼波延迟时间.合成的及某海上采集的倾斜缆数据去鬼波处理结果表明,频率慢度域自适应迭代反演算法能较好地去除海上变深度缆鬼波,能达到拓宽地震记录频带目的.     

Crustal structure of northeastern margin of the Tibetan Plateau by receiver function inversion

Using seismic data of about one year recorded by 18 broadband stations of ASCENT project,we obtained 2547 receiver functions in the northeastern Tibetan Plateau.The Moho depths under 14 stations were calculated by applying the H-κ domain search algorithm.The Moho depths under the stations with lower signal-noise ratio(SNR) were estimated by the time delay of the PS conversion.Results show that the Moho depth varies in a range of ~40–60 km.The Moho near the Haiyuan fault is vague,and its depth is larger than those on its two sides.In the Qinling-Qilian Block,the Moho becomes shallower gradually from west to east.To the east of 105°E,the average depth of the Moho is 45 km,whereas the west is 50 km or even deeper.Combining our results with surface wave research,we suggest a boundary between the Qinling and the Qilian Mountains at around 105°E.S wave velocities beneath 15 stations have been obtained through a linear inversion by using Crust2.0 as an initial model,and the crustal thickness that was derived by H-κ domain search algorithm was also taken into account.The results are very similar to the results of previous active source studies.The resulting figure indicates that low velocity layers developed in the middle and lower crust beneath the transition zone of the Tibet Block and western Qinling,which may be related to regional faults and deep earth dynamics.The velocity of the middle and lower crust increases from the Songpan Block to the northeastern margin of Tibetan Plateau.Based on the velocity of the crust,the distribution of the low velocity zone and the composition of the curst(Poisson's ratio),we infer that the crust thickening results from the crust shortening along the direction of compression.     

Seismic data compression based on integer wavelet transform

Introduction With the development of the seismological observation technique and deep-going of seismicdata application fields, especially the digitization of data in earthquake station networks, theimprovement of the precision, the data quantity increases as geometric order, which bringdifficulty to saving and transfering these data. To keep all information, seismic data, like medicalimages, should be compressed without error in many applications. In generally, traditionalcompression meth…     

Seismic data compression based on integer wavelet transform

Due to the particularity of the seismic data, they must be treated by lossless compression algorithm in some cases. In the paper, based on the integer wavelet transform, the lossless compression algorithm is studied. Comparing with the traditional algorithm, it can better improve the compression rate. CDF (2, n) biorthogonal wavelet family can lead to better compression ratio than other CDF family, SWE and CRF, which is owe to its capability in canceling data redundancies and focusing data characteristics. CDF (2, n) family is suitable as the wavelet function of the lossless compression seismic data. Contribution No.04FE1019, Institute of Geophysics, China Earthquake Administration.     

应力张量反演的遗传算法及其在青藏高原东北缘的应用

为准确而快速地求解区域构造应力张量,构建了基于震源机制解反演应力张量的遗传算法策略,详细介绍了计算原理,分析了以断层面上剪应力方向和滑动方向的偏差为约束反演应力张量而忽略剪应力大小偏差对反演结果的影响,表明仅考虑剪应力方向和滑动方向的偏差就可以得到正确结果.利用不同应力状态下人工合成的包含不同噪声水平的震源机制数据对该方法进行检验,并与网格搜索法得到的结果比较.表明本文方法所得结果的拟合差均小于网格搜索法结果的拟合差,并且二者相差较小,体现了本文方法的稳健性.将该方法应用于青藏高原东北缘(103-106°E,34.5-37.5°N)应力张量的估计,结果显示,该区域主要受控于青藏高原近东西向的挤压,从而导致阿拉善块体以及华南块体方向的拉张,并且向华南块体方向的拉张作用强于向阿拉善块体的拉张作用.     

A numerical comparison of time and frequency‐domain marine electromagnetic methods for hydrocarbon exploration in shallow water

In shallow water the frequency domain controlled source electromagnetic method is subject to airwave saturation that strongly limits the sensitivity to resistive hydrocarbon targets at depth. It has been suggested that time‐domain CSEM may offer an improved sensitivity and resolution of these deep targets in the presence of the airwave. In order to examine and test these claims, this work presents a side‐by‐side investigation of both methods with a main focus on practical considerations, and how these effect the resolution of a hydrocarbon reservoir. Synthetic noisy data for both time‐domain and frequency domain methods are simulated using a realistic frequency dependent noise model and frequency dependent scaling for representative source waveforms. The synthetic data studied here include the frequency domain response from a compact broadband waveform, the time‐domain step‐response from a low‐frequency square wave and the time‐domain impulse response obtained from pseudo‐random binary sequences. These data are used in a systematic resolution study of each method as a function of water‐depth, relative noise and stacking length. The results indicate that the broadband frequency domain data have the best resolution for a given stacking time, whereas the time‐domain data require prohibitively longer stacking times to achieve similar resolution.     

承德地震台VP型宽频带垂直摆运行效能分析

与VS型垂直摆各项参数进行对比,分析VP型宽频带垂直摆具有的优势。采用VS、VP型垂直摆记录资料,对比数据连续率、完整率及固体潮记录能力、同震变化等;分析VP型垂直摆精度、潮汐因子绝对误差、相位误差等,以偏离正常年变背景为基准,判断其数据记录是否异常,并分析年变特征,拟合观测数据分析异常特征,通过R值评分,进一步分析承德地震台VP型宽频带垂直摆映震能力。分析发现,VP型宽频带垂直摆记录固体潮和同震曲线的能力较好,精度较高,潮汐因子误差及绝对误差较小,产出数据更为可靠,异常信息更为明显,映震能力更高。     

VP宽频带倾斜仪现场自动校准装置设计

为了对VP宽频带倾斜仪进行现场校准,设计了一款基于单片机AT89C51的新型现场自动校准装置。经验证该装置在+8V激励电压情况下与高精度激光干涉仪结合使用,其校准相对误差为0.36%,满足小于1%的校准相对误差要求。     

24道宽频带程控数字记录器硬件的研制

24道宽频带程控数字记录器为数字化测振仪器．该仪器的放大增益、频带、记录方式用单片机控制,可以广泛用于浅层地震信息数字化数据的采集记录．     

使用双线性变换构造实时地震波形仿真滤波器的方法研究

介绍了直接使用双线性变换把模拟系统的传递函数转换为离散系统的传递函数,构造出由宽带数字地震记录,仿真出传统模拟地震仪器记录波形的递归滤波器的方法,从而使在线实时波形仿真成为可能。并通过比较利用实际宽带数字地震记录,分别采用频域仿真和本文算法构造的仿真滤波器仿真得到的波形,验证了本文算法的正确性。     

Well tie for broadband seismic data

The seismic industry is increasingly acquiring broadband data in order to reap the benefits of extra low‐ and high‐frequency contents. At the low end, as the sharp low‐cut decay gets closer to zero frequency, it becomes harder for a well tie to estimate the low‐frequency response correctly. The fundamental difficulty is that well logs are too short to allow accurate estimation of the long‐period content of the data. Three distinctive techniques, namely parametric constant phase, frequency‐domain least squares with multi‐tapering, and Bayesian time domain with broadband priors, are introduced in this paper to provide a robust solution to the wavelet estimation problem for broadband seismic data. Each of these techniques has a different mathematical foundation that would enable one to explore a wide range of solutions that could be used on a case‐by‐case basis depending on the problem at hand. A case study from the North West Shelf Australia is used to analyse the performance of the proposed techniques. Cross‐validation is proposed as a robust quality control measure for evaluating well‐tie applications. It is observed that when the seismic data are carefully processed, then the constant phase approach would likely offer a good solution. The frequency‐domain method does not assume a constant phase. This flexibility makes it prone to over‐fitting when the phase is approximately constant. Broadband priors for the time‐domain least‐squares method are found to perform well in defining low‐frequency side lobes to the wavelet.     

GM3磁通门磁力仪数据实时压缩方法的改进

研究了GM3磁通门磁力仪数据实时压缩方法,将国际通用的测震数据SEED格式中采用的Steim—1压缩算法引入到地磁数据的实时处理系统中,通过改进压缩算法,提高了压缩率,节省了仪器数据存储空间,减少了台站数据传输时间。     

探地雷达小波域三维波动方程偏移

阐述了矩阵多分辨分析理论中的标准形式与非标准形式,并以Hilbert算子为例,说明了算子多分辨表示的压缩效果,为小波域偏移算法奠定了理论基础.从三维雷达波动方程出发,利用爆炸反射原理和浮动坐标变换,推导出三维探地雷达波动方程差分格式,并通过方程分裂算法及小波多分辨算法,在小波域求解波场外推矩阵,进而得到探地雷达小波域三维波动方程偏移算法,在此基础上,开发了探地雷达小波域偏移处理程序,并把该程序应用于三个球体空洞的3-D正演结果及实际的雷达数据中,通过对比偏移处理前后的雷达资料,得知该三维偏移算法能使3-D正演剖面中的反射波归位、绕射波收敛,极大地提高了雷达剖面的分辨率,有利于探地雷达资料的地质解释.     

DPCM-based vibration sensor data compression and its effect on structural system identification

Due to the large scale and complexity of civil infrastructures, structural health monitoring typically requires a substantial number of sensors, which consequently generate huge volumes of sensor data. Innovative sensor data compression techniques are highly desired to facilitate efficient data storage and remote retrieval of sensor data. This paper presents a vibration sensor data compression algorithm based on the Differential Pulse Code Modulation (DPCM) method and the consideration of effects of signal distortion due to lossy data compression on structural system identification. The DPCM system concerned consists of two primary components: linear predictor and quantizer. For the DPCM system considered in this study, the Least Square method is used to derive the linear predictor coefficients and Jayant quantizer is used for scalar quantization. A 5-DOF model structure is used as the prototype structure in numerical study. Numerical simulation was carried out to study the performance of the proposed DPCM-based data compression algorithm as well as its effect on the accuracy of structural identification including modal parameters and second order structural parameters such as stiffness and damping coefficients. It is found that the DPCM-based sensor data compression method is capable of reducing the raw sensor data size to a significant extent while having a minor effect on the modal parameters as well as second order structural parameters identified from reconstructed sensor data.     

逆散射级数法预测层间多次波的算法改进及其策略

波动方程预测层间多次波有两类算法.一是由Berkhout和Verschuur（1998）提出的基于CFP（共聚焦点）延拓的算法,另一是由Weglein（1997）等人提出的基于逆散射级数法的算法.ISS（逆散射级数）算法具有不依赖速度模型的优点,但是其计算成本很高.本文改进了Weglein提出的1-D预测公式,提高了计算效率.从理论方面,改进算法的计算速度提高了大约12倍.实际结果表明其计算速度更快,提高近80倍左右,主要原因在于后者的空间复杂度小.此外,由于时间-空间域比伪深度-波数域的层间多次波预测有更多的优点,本文推导了1.5-D时间-空间域ISS层间多次波预测算法,该算法有计算速度快、预测噪音小、适应观测系统能力强和不依赖地震子波等优点.在实际应用中,计算成本高是ISS预测层间多次波算法的主要障碍.因此,本文还讨论了预测层间多次波的策略.结合本文提出的改进算法和应用策略,能够在实际地震资料处理中较好地应用ISS预测层间多次波.     

Efficient seismic response data storage and transmission using ARX model‐based sensor data compression algorithm

This paper presents a linear predictor (LP)‐based lossless sensor data compression algorithm for efficient transmission, storage and retrieval of seismic data. Auto‐Regressive with eXogenous input (ARX) model is selected as the model structure of LP. Since earthquake ground motion is typically measured at the base of monitored structures, the ARX model parameters are calculated in a system identification framework using sensor network data and measured input signals. In this way, sensor data compression takes advantage of structural system information to maximize the sensor data compression performance. Numerical simulation results show that several factors including LP order, measurement noise, input and limited sensor number affect the performance of the proposed lossless sensor data compression algorithm concerned. Generally, the lossless data compression algorithm is capable of reducing the size of raw sensor data while causing no information loss in the sensor data. Copyright © 2005 John Wiley & Sons, Ltd.     

Best Practice in Earthquake Location Using Broadband Three-component Seismic Waveform Data

—?We present an earthquake location algorithm, the Broadband Waveform Regional Earthquake Location Program (BW_RELP), which utilizes phase onset times and wave azimuths recorded by three-component broadband seismic stations and an adaptive migrating grid search algorithm to find the global minimum in an arbitrary normed misfit parameter. The performance of BW_RELP is demonstrated using regional (300–800?km distant) broadband recordings to locate events in the 1995 Ridgecrest, California earthquake sequence. The purpose of this study is to introduce the BW_RELP algorithm in detail and to expand on the previous paper by Deger et?al. (BSSA, 88, 1353–1362, 1998), using one Berkeley Digital Seismic Network (BDSN) station (YBH) and two USNSN stations (ELK and MNV) which span 300–800?km in distance and 55 degrees in azimuth, to further investigate the capability of a sparse broadband network of three-component stations at monitoring a region located outside of the network, as will be the case in the monitoring of the Comprehensive Test-Ban-Treaty (CTBT) for low magnitude seismic events. We assess the capability of this sparse three-station broadband network and we compare locations estimated from phase onset time and wave azimuth measurements to a ground-truth catalog of high-quality earthquake locations derived from data recorded by the Southern California Seismic Network (SCSN). The results indicate that in the regional distance range it is possible, when an appropriate calibration event is available, to obtain absolute event locations to within 18?km as is prescribed by the CTBT.