地震预测中的地电阻率数据处理方法

本文介绍了地电阻率时间序列数据的常用处理方法,即消除年变化法、无量纲法、相对均方差法和差分法共4类8种方法,包括“九五”以来对原方法的改进和新发展的方法,讨论了各方法的原理、数据处理过程、异常分辨能力、异常指标和异常物理机制及其在数据处理中存在的不足.结果显示:① 一般来说,用有效的数据处理方法才能分析、识别出原始数据曲线上的“弱变化”异常;② 文中8种方法的原理简明,异常物理机制清晰或较清晰,各方法的异常识别指标分别为明确、基本明确和定性的,定性的异常在震情研判中仅有参考意义;③ 消除年变化法和无量纲法通常用于识别地震中期、短临异常,而相对均方差法和差分法通常用于识别短临异常;④ 经数据处理得到的异常与原始曲线的“弱变化”异常相协调;⑤ 时间序列数据出现的异常并不等同与地震孕育、发生过程有直接联系的前兆异常,出现数据异常的台站附近不一定会发生显著地震.      

数字化地震前兆台网观测数据质量评价方法

为进一步提升我国地震前兆台网的管理水平,对数字化地震前兆台网观测数据质量进行快速而科学合理的评价,在现有各学科评比办法的基础上,开展了台网数据质量评价方法研究。通过研究获得了1套包含观测数据的完整性、噪声水平、固体潮汐效应、一致性等多方面的数据质量评价指标,在此基础上,采用规范化处理与加权平均等方法建立起台网数据质量评价方法。应用分析显示,评价结果基本能够反映台网的数据质量实际状况。     

Electromagnetic response function for the semi-annual variations estimated by the single-station method using the night-time values of the geomagnetic fields

The geomagnetic night-time values were used to estimate the electromagnetic response function Q1 for half-year period. If the spatial structure of the source field can be described by the approximation, one can estimate the Q1 value using the single-station Z/H method. This technique enables us to carry out regional deep geomagnetic sounding by the method. The data used for analyses are geomagnetic night-time values for about, typically, 26 years from 5 good-quality stations and for several years from 34 stations distributed over the globe. The results indicate that the night-time values yield more reliable response estimates for half-year period compared to the usual estimates obtained from daily means. It implies that the approximation for the night-time fields holds good for the half-year period, but the daily means are not suitable for estimating the response function of the semi-annual variations by using the single-station method. Source field analyses for daily means data and night-time means data have also been carried out in this paper.     

Multiple prediction using the homeomorphic-imaging technique

A new method for predicting different kinds of multiples and peg-leg reflections in unstacked seismic data is discussed. The basis for this method is the fact that kinematic properties of multiples can be represented as a combination of kinematic properties of primary reflections. The prediction is made using a two-step process. In the first step, the values for the angle of emergence and radius of curvature of the wavefront for primary reflections from ‘multiple-generating’ interfaces are obtained. These parameters are estimated directly from unstacked data for every source point using the homeomorphic-imaging technique. The second step consists of prediction of multiples from primary reflections that satisfy a so-called ‘multiple condition’. This condition is the equality of the absolute values of the angles of emergence calculated from the first step. This method is effective even in complex media and information on the subsurface geology is not required. The parameters are estimated directly from the unstacked data and do not require any computational efforts such as in wavefield extrapolation of data.     

Solution of the one-dimensional inverse magnetotelluric problem

Summary Methods of solving the inverse magnetotelluric problem are compared. Basic relations for Newton's method, the least-squares method and Marquardt's method are presented and the convergence properties of these methods are studied. The high effectiveness of Marquardt's method is demonstrated and its application to practical magnetotelluric data is discussed.     

Using recovery data to extend the effective duration of pumping tests

Collection of water-level recovery data is a common practice for pumping tests. The resulting data can provide some of the most useful information from the tests, but are rarely used to their full value. van der Kamp (1989) proposed a general method that uses recovery data to extend the effective duration of pumping. The method is straightforward to implement and applicable for simple or complex hydrogeologic settings. The only assumption invoked is that the response remains linear such that the principle of superposition can be applied. No other assumptions about the properties of the aquifer are required. The method can greatly increase the value of pumping tests by extending the effective duration of the tests for as long as significant residual drawdowns are observed.     

3D simultaneous seismic data reconstruction and noise suppression based on the curvelet transform

Seismic data contain random noise interference and are affected by irregular subsampling. Presently, most of the data reconstruction methods are carried out separately from noise suppression. Moreover, most data reconstruction methods are not ideal for noisy data. In this paper, we choose the multiscale and multidirectional 2D curvelet transform to perform simultaneous data reconstruction and noise suppression of 3D seismic data. We introduce the POCS algorithm, the exponentially decreasing square root threshold, and soft threshold operator to interpolate the data at each time slice. A weighing strategy was introduced to reduce the reconstructed data noise. A 3D simultaneous data reconstruction and noise suppression method based on the curvelet transform was proposed. When compared with data reconstruction followed by denoizing and the Fourier transform, the proposed method is more robust and effective. The proposed method has important implications for data acquisition in complex areas and reconstructing missing traces.     

Accuracy Estimation of the Modern Core Magnetic Field Models Using DMA-Methods for Recognition of the Decreased Geomagnetic Activity in Magnetic Observatory Data

A new approach to identify the signals of the Earth’s main magnetic field (core field) based on the magnetic observatory data processing is suggested. The algorithms implemented in the approach are based on the Discrete mathematical analysis (DMA). The developed method is used to analyze the data from 49 midand low-latitude observatories of the INTERMAGNET network collected in 2011–2015. The results are compared with the classical method for determining the periods of low magnetic activity of external origin which is adopted by the International Association for Geomagnetism and Aeronomy (IAGA). The advantages of the suggested new approach are demonstrated. Based on the data records for the selected time intervals, the time series of the core field components and their secular variations are obtained for each observatory. These data are compared to the values predicted by the most accurate core field models: SIFM, CHAOS-6, and EMM-2015. The accuracy of the models is estimated using a set of statistical parameters: Pearson’s coefficient of linear correlation, Spearman’s and Kendall’s coefficients of rank correlation, the mean and median values over the data sets, and the mean difference between the data obtained by the suggested method from observatory measurements and the model predictions.     

数字化地电阻率观测数据的干扰识别及处理方法

针对辽宁省数字化地电阻率观测数据的特点,找出同步高频变化干扰数据的识别及处理方法;观测资料中缺测数据的干扰识别及处理方法;无物理基础的观测量值识别及处理方法等;通过VB编程快捷、方便地剔除瞬时干扰数据,有利于地震预报人员对数字化观测干扰数据的处理,并对地下介质的动态变化做出正确的分析、判断.     

Curvelet reconstruction of non-uniformly sampled seismic data using the linearized Bregman method

Seismic data reconstruction, as a preconditioning process, is critical to the performance of subsequent data and imaging processing tasks. Often, seismic data are sparsely and non-uniformly sampled due to limitations of economic costs and field conditions. However, most reconstruction processing algorithms are designed for the ideal case of uniformly sampled data. In this paper, we propose the non-equispaced fast discrete curvelet transform-based three-dimensional reconstruction method that can handle and interpolate non-uniformly sampled data effectively along two spatial coordinates. In the procedure, the three-dimensional seismic data sets are organized in a sequence of two-dimensional time slices along the source–receiver domain. By introducing the two-dimensional non-equispaced fast Fourier transform in the conventional fast discrete curvelet transform, we formulate an L1 sparsity regularized problem to invert for the uniformly sampled curvelet coefficients from the non-uniformly sampled data. In order to improve the inversion algorithm efficiency, we employ the linearized Bregman method to solve the L1-norm minimization problem. Once the uniform curvelet coefficients are obtained, uniformly sampled three-dimensional seismic data can be reconstructed via the conventional inverse curvelet transform. The reconstructed results using both synthetic and real data demonstrate that the proposed method can reconstruct not only non-uniformly sampled and aliased data with missing traces, but also the subset of observed data on a non-uniform grid to a specified uniform grid along two spatial coordinates. Also, the results show that the simple linearized Bregman method is superior to the complex spectral projected gradient for L1 norm method in terms of reconstruction accuracy.     

应用多种来源重力异常编制中国海陆及邻区空间重力异常图及重力场解读

本文通过分析陆地实测空间重力异常数据、海洋船载测量空间重力异常数据、卫星测高重力异常,布格重力异常数据、EGM2008地球重力模型数据等多种来源数据的性质和精度,并对相关数据进行对比,研究了编制1:500万中国海陆空间重力异常图的数据使用方案和技术方法.在地形较为平坦、实测数据分布均匀的陆区,使用实测数据,在地形复杂,实测数据稀少以及没有实测数据的陆区或岛屿,利用布格重力异常反推空间异常的方法合成平均空间重力数据,西藏地区的数据对比实验证明合成平均空间重力异常数据是一种有效的数据补充.利用三观测列方差分解法在南海地区对船载测量空间重力数据和美国SS系列及丹麦DNSC08GRA卫星重力数据进行了方差分解计算,结果表明不同来源的卫星测高重力数据具有很大的一致性,数据精度较以往有了很大的提高.海区空间重力数据使用原则是在船载重力测量数据校准下,全面使用卫星测高重力数据进行编图.海陆过渡区的异常处理应以EGM2008地球重力模型重力场为基准参考场,实现海陆异常平缓过渡,无缝连接.对中国海陆空间重力异常场进行了小波变换处理,对空间重力异常场进行了解读,勾画出三横四竖的一级重力梯级带及其所围限的8个一级重力异常区,并划分了二级重力异常区和梯级带,为块体构造学体系中大地构造格架的建立提供了地球物理证据.     

频率域偶层位曲面位场处理和转换方法研究

在空间域偶层位法的基础上,研究了完整的频率域偶层位曲面位场处理和转换方法.该法可应用于平面或曲面、规则网或非规则网的位场数据处理和转换.通过对偶层面z坐标和计算面z坐标平移不同的量来加速正演快速收敛和保证反演稳定、快速收敛;提出了适合于不规则网曲面处理和转换的核心算法——单点快速Fourier变换;提出了频率域不规则网曲面处理和转换方法技术.通过以上技术措施解决了大数据量特别是曲面不规则网的位场处理和转换问题,模型试算以及实际资料处理验证了该方法的应用效果.     

Direct analytic signal (DAS) method in the interpretation of magnetic data

We presented a new method for interpreting 2D magnetic data, called direct analytic signal (DAS) method, which directly used the analytic signal of magnetic anomaly to compute the depth and the structural index of the source. The DAS method needs only the computation of the first order derivatives of magnetic anomaly, so that the inversion results are more stable than the results obtained by the other existing analytic signal methods. The DAS method is tested on synthetic magnetic data with and without noise, and the DAS method can successfully obtain the depth and the structural index of the source. We also applied the DAS method to interpret a real magnetic data over a shallow geological source whose source parameters are known from closely drilling information, and the inversion results are in accord with the true values.     

利用零偏移VSP资料估计介质品质因子方法研究

利用峰值频率移动法估算零偏VSP资料的品质因子Q．该方法用Ricker子波和匹配地震子波分别逼近零相位和混合相位的震源子波,得到了峰值频率移动法估计Q值的公式．进而针对常规方法估计的地震子波峰值频率精度不高的问题,提出了估计地震子波峰值频率的特征结构法．通过合成零偏VSP资料的仿真试验,验证了峰值频率移动法估计Q值的正确性．仿真结果表明,与快速Fourier变换和Burg最大熵方法相比较,特征结构法得到的峰值频率和Q值精度高一些．仿真结果也表明,用峰值频率移动法估计Q值时需要选取恰当的子波参数,否则影响Q值估计的精度．     

Depth Estimation from Gravity Data Using the Maximum Entropy Method (MEM) and the Multi Taper Method (MTM)

The depth determination from the gravity data in frequency domain is carried out using the classical fast Fourier transform (FFT) method utilizing scaling properties of ensemble of anomalous source. The problem of calculating power spectrum from the FFT is well described in the literature. Here, the application of other high-resolution methods of power spectrum calculation, such as maximum entropy method (MEM) and multi-taper method (MTM) are explored to estimate depth to anomalous sources. At the outset, the FFT, the MEM and the MTM are tested on synthetic gravity data, generated for different types of synthetic models and then all these methods are applied to the field gravity data of the Bengal basin. The MTM with scaling is found to be superior for providing the detailed subsurface information rather than the MEM and the FFT methods in the case of synthetic as well as field examples.     

希尔伯特-黄变换在地震资料处理中的应用

希尔伯特-黄变换是最新发展起来的处理非线性非平稳信号的时频分析方法.本文介绍了该方法的基本原理和几个关键技术,以及其优势所在,并将它应用于含随机噪音的地震信号和实际地震资料的分解.分解后的信号能够重构原始信号,重构后在信噪比和分辨率方面有了一定的提高,对地震资料的去噪有重要意义.     

四分位法在地磁日变化观测数据分析中的应用

四分位法可用于描述包括偏态分布、分布不明和分布末端无确定值等任何分布资料的特征,并且可以去除数据中大数的影响,较均值、标准差等更为稳定.针对地磁日变化数据的复杂性,不同时段数据分布特征的差异性和不确定性,在此引进四分位法加以分析.通过分析,提取地磁日变化中来自地磁场扰动、仪器记录错误、人为干扰以及可能与地震的孕育有关因...     

在反演过程中利用不同范数准则压制噪声对反演结果的影响

本文对误差数据在各种范数准则下的反演结果进行了详细的分析,在此基础上提出了在反演过程中利用不同范数准则压制噪声的方法,即在反演过程中先用一种范数准则建立目标泛函并求取梯度、再利用另一种范数准则计算修正步长,并在建立目标泛函和求取梯度时引入了非整数次范数准则.该方法有效地避免了在压制噪声过程中的人为干扰,提高了数据的使用率.从合成数据的反演结果可以看出,该方法对噪声具有很好的压制能力.     

环工作面三维直流电阻率法研究及应用

本文在矿井直流电法理论基础上,推导环工作面三维直流电法测量深度公式,得到环工作面三维直流电法记录点坐标位置.基于以上理论,建立了环工作面三维直流电法异常体模型,分析了环工作面三维直流电法二极装置、三极装置、四极装置异常体响应特征,得出环工作面三维直流电法各装置敏感度低的主要原因是缺少必要的数据采集点.为了提高环工作面三维直流电法分辨率,最大限度增加工作面内侧的采集数据点.基于以上原因本文提出了对角偶极装置采集模式及多装置的数据叠加反演处理方法,以增加工作面内采集数据点,提高三维勘探的精度.本文通过某矿8802工作面底板富水性探查实例,说明该方法的有效性.     

大地电磁测深资料的二次函数逼近非线性反演

将二次函数逼近非线性优化首次应用于大地电磁测深反演问题,该反演方法利用二次函数有唯一最小值的特点进行逼近大地电磁反演模型,从而避免了常规的迭代反演过程中陷入局部极小问题,实现了对目标函数求全局极小,较好地解决了非唯一性问题;同时该方法不用求灵敏度矩阵,且对初始模型无任何要求。通过理论模型检验、井旁MT点反演结果与测井曲线的对比及MT测线的反演电阻率深度剖面与地震测线的时间剖面对比均表明,本文方法取得较好的应用效果。