一种ARMA模型地震子波提取方法研究

在地震子波因果、混合相位的假设下,本文突破突破滑动平均(MA)型假设的地震子波提取技术,提出一种阶数吝啬的自回归滑动平均(ARMA)模型对地震子波进行参数化准确建模的方法.针对地震子波ARMA模型定阶困难,超定阶容易造成计算量大、运算速度慢,欠定阶不能满足精确子波的要求,本文采用基于自相关函数的奇异值分解(SVD)法确定AR模型阶数,同时将信息量准则法与高阶累积量法相结合,提出了一种新的MA模型定阶法.数值仿真和实际地震数据处理结果表明,本文所用方法可以有效地压制加性高斯色噪声,信息量准则法可有效提高MA定阶的准确率,在保证子波精度的同时尽可能降低模型阶数,实现运算高效率.     

用对称映射ARMA模型的零极点研究子波相位对反射系数序列反演的影响

为研究地震子波相位对反射系数序列反演的影响,在自回归滑动平均(ARMA)模型描述子波的基础上,提出采用z域对称映射ARMA模型零极点的方法构造了一系列相同振幅谱、不同相位谱的地震子波,并结合谱除法对人工合成地震记录进行反射系数序列反演.理论分析表明,子波相位估计不准时反射系数序列反演结果中残留一个纯相位滤波器,该纯相位滤波器的相位谱为真实子波和构造子波的相位谱之差.采用丰度和变分作为评价方法,在反演结果中确定出真实的或准确的反射系数序列.仿真实验和实际数据处理结果也验证了子波相位对反射系数序列反演的影响规律和评价方法的有效性,为进一步提高反射系数序列反演结果精度指明了研究方向.     

基于ARMA模型拟合分析方法的地磁监测数据研 究——以长春地磁台观测数据为例

ARMA模型也叫做自回归滑动平均模型,它是研究时间序列的重要方法,其中AR为自回归模型,MA为移动平均模型。长春地磁数据F分量的ARMA模型拟合分析是基于MATLAB平台,运用时间序列分析的研究方法,对2007—2018年长春地磁数据F分量进行ARMA模型拟合。按照不同ARMA模型相异的数据特性确定模型类别,判断其属于ARMA模型、AR模型或是MA模型。运用AIC准则法,找到使得AIC值最小的参数并确定为模型参数,最后对残差序列运用Q值检验法,根据检验结果判断模型的拟合优度。本文中的ARMA模型拟合方法是针对F分量的预处理分数据建立数学模型,运用一阶差分方法去除年变趋势,最终得出F分量变化幅度的近似拟合回归方程。通过ARMA模型拟合方程可初步预测F值增减幅度的变化趋势及极值范围,这将为地磁异常乃至其它地震异常定量分析和预测提供一种新思路,也将为地磁异常的判定核实及地震分析预报起到推动作用。     

6级左右历史地震震源参数的模拟估计

本文借鉴直接拟合烈度数据点和枚举震源参数的做法,设计了一种利用烈度资料估计6级左右历史地震震源参数的方法.该方法对震源参数所有可能的组合进行枚举,采用地震波场模拟计算转换的理论烈度值,利用模型选择方法评估各可能的震源参数组合模型与历史破坏记录推断的地震烈度数据点的拟合程度,对震源参数做出估计.该方法充分考虑到历史资料相对稀少对震源参数估计的影响,以多种震源参数估计结果和相应权重值来定量化表示估计的不确定性.通过对给定震中位置、震源深度和滑动角的Bootstrap数值恢复检测与2006年美国Parkfield 6.0级地震实例的测试,表明该方法得出的震源参数估计结果具有统计一致性和一定的无偏性.将该方法应用于1882年河北深县6级地震的震源参数估计,结果显示东西向旧城北断层或何庄断层及北东东走向的深西断层为深县地震的发震构造的可能性较大.     

参数化地震子波估计模型定阶方法研究综述

参数化模型地震子波提取中模型阶数的确定对子波提取的精度至关重要,笔者针对因果地震子波MA、ARMA模型,总结了目前比较成熟的四类定阶方法：相关分析法、信息量准则法、线性代数法及基于优化算法的模型定阶方法,分析比较发现这些算法都存在不同程度的缺陷,只能在各自的适用范围内得到较好的定阶结果,针对算法在子波建模应用中的不足提出了改进建议;同时分析比较了基于高阶累积量的非因果系统模型定阶的各种方法,并对非因果地震子波模型定阶进行了展望.     

基于测井、VSP和地面地震数据最佳拟合的子波估计(英文)

本文提出了基于测井、VSP和地震数据拟合的子波估计方法,从输入输出都包含随机噪声的统计模型出发,采用相关性拟合技术来提取子波。拟合度和误差分析为整个过程提供了定量的质量控制手段,可以评估数据拟合和子波估计的可靠性。实际数据试算表明,该方法在含有噪声的实际数据中稳定而有效,在地震频带内的子波估计和数据拟合是可靠的。该方法无需对子波相位和振幅谱进行任何假设,其主要优点在于确定相位的能力。     

基于保角变换的双谱域相位谱估计方法及其在地震子波估计中的应用(英文)

地震子波估计是地震资料处理与解释中的重要环节,它的准确与否直接关系到反褶积及反演等结果的好坏。高阶谱(双谱和三谱)地震子波估计方法是一类重要的、新兴的子波估计方法,然而基于高阶谱的地震子波估计往往因为高阶相位谱卷绕的原因,导致子波相位谱求解产生偏差,进而影响了混合相位子波估计的效果。针对这一问题,本文在双谱域提出了一种基于保角变换的相位谱求解方法。通过缩小傅里叶相位谱的取值范围,有效避免了双谱相位发生卷绕的情况,从而消除了原相位谱估计中双谱相位卷绕的影响。该方法与最小二乘法相位谱估计相结合,构成了基于保角变换的最小二乘地震子波相位谱估计方法,并与最小二乘地震子波振幅谱估计方法一起,应用到了地震资料混合相位子波估计中。理论模型和实际资料验证了该方法的有效性。同时本文将双谱域地震子波相位谱估计中保角变换的思想推广到三谱域地震子波相位谱估计中。     

基于矩阵Toeplitz稀疏分解的相对波阻抗反演方法

针对利用地震道进行相对波阻抗反演中遇到的横向连续性难以保持、初始子波容错度差以及随机噪声干扰影响反演结果等问题,提出了一种基于矩阵Toeplitz稀疏分解的相对波阻抗反演方法.该方法将地震数据剖面的Toeplitz稀疏分解问题分解为两个子反演问题,其一以Toeplitz子波矩阵元素为待反演的参数,用Fused Lasso方法求解,可保证子波具有紧支集且是光滑的;其二以稀疏反射系数矩阵元素为待反演参数,用基于回溯的快速萎缩阈值迭代算法求解,大大降低了目标函数中参数选择的难度.通过交替迭代求解上述两个子反演问题可将地震数据剖面因式分解为一个Toeplitz子波矩阵和一个稀疏反射系数矩阵;然后由反射系数矩阵递推反演可以得到高分辨率的相对波阻抗剖面;利用测井资料加入低频分量后,也可得到高分辨率的绝对波阻抗剖面.Marmousi2模型生成的合成记录算例和实际地震资料算例均表明:本文方法可以从带限地震数据中有效地反演相对波阻抗,反演结果分辨率高并且能够很好地保持地震数据的横向连续性;即使在初始估计子波存在误差和地震数据被随机噪声污染的情况下也能取得较好的效果.     

基于ARMA模型非因果空间预测滤波(英文)

常规频域预测滤波方法是建立在自回归(autoregressive,AR)模型基础上的,这导致滤波过程中前后假设的不一致,即首先利用源噪声的假设计算误差剖面,却又将其作为可加噪声而从原始剖面中减去来得到有效信号。本文通过建立自回归-滑动平均(autoregres sive/moving-average,ARMA)模型,首先求解非因果预测误差滤波算子,然后利用自反褶积形式投影滤波过程估计可加噪声,进而达到去除随机噪声目的。此过程有效避免了基于AR模型产生的不一致性。在此基础上,将一维ARMA模型扩展到二维空间域,实现了基于二维ARMA模型频域非因果空间预测滤波在三维地震资料随机噪声衰减中的应用。模型试验与实际资料处理表明该方法在很好保留反射信息同时,压制随机噪声更加彻底,明显优于常规频域预测去噪方法。     

研究各向异性介质参数的合成地震图方法

基于EDA各向异性理论, 利用合成地震图方法求取各向异性介质参数. 该方法利用求取克利斯多夫方程特征向量和特征值的解析表达式, 避免了矩阵运算过程中奇异解的产生, 也减少了大量的运算时间, 具有一定的优越性. 同时, 运用凝聚函数法对合成记录与原始记录进行检验, 结果显示在S波的主频范围内二者具有很好的线性关系. 表明用该方法求得的各向异性介质参数, 可以客观反映和描述研究区域的介质各向异性特性.      

Seismic wavelet estimation via a system identification method

On the assumption that the wavelet is causal and nonminimum phase, an autoregressive moving average (ARMA) model is introduced to fit the seismic trace. Seismic wavelet extraction is converted to parameters estimation of the ARMA model. Singular value decomposition (SVD) of an appropriate matrix formed by autocorrelation is exploited to determine the autoregressive (AR) order, and the cumulant-based SVD-TLS (total least squares) approach is proposed to obtain the AR parameters. The author proposes a new moving average (MA) model order determination method via combining the information theoretic criteria method and higher-order cumulant method. The cumulant approach is used to achieve the MA parameters. Theoretical analysis and numerical simulations demonstrate the feasibility of the wavelet extraction approach.     

Time-frequency decomposition of seismic signals via quantum swarm evolutionary matching pursuit

Matching pursuit belongs to the category of spectral decomposition approaches that use a pre-defined discrete wavelet dictionary in order to decompose a signal adaptively. Although disengaged from windowing issues, matching point demands high computational costs as extraction of all local structure of signal requires a large size dictionary. Thus in order to find the best match wavelet, it is required to search the whole space. To reduce the computational cost of greedy matching pursuit, two artificial intelligence methods, (1) quantum inspired evolutionary algorithm and (2) particle swarm optimization, are introduced for two successive steps: (a) initial estimation and (b) optimization of wavelet parameters. We call this algorithm quantum swarm evolutionary matching pursuit. Quantum swarm evolutionary matching pursuit starts with a small colony of population at which each individual, is potentially a transformed form of a time-frequency atom. To attain maximum pursuit of the potential candidate wavelets with the residual, the colony members are adjusted in an evolutionary way. In addition, the quantum computing concepts such as quantum bit, quantum gate, and superposition of states are introduced into the method. The algorithm parameters such as social and cognitive learning factors, population size and global migration period are optimized using seismic signals. In applying matching pursuit to geophysical data, typically complex trace attributes are used for initial estimation of wavelet parameters, however, in this study it was shown that using complex trace attributes are sensitive to noisy data and would have lower rate of convergence. The algorithm performance over noisy signals, using non-orthogonal dictionaries are investigated and compared with other methods such as orthogonal matching pursuit. The results illustrate that quantum swarm evolutionary matching pursuit has the least sensitivity to noise and higher rate of convergence. Finally, the algorithm is applied to both modelled seismograms and real data for detection of low frequency anomalies to validate the findings.     

SEISMIC PHASE UNWRAPPING: METHODS,RESULTS, PROBLEMS1

Six known methods of seismic phase unwrapping (or phase restoration) are compared. All the methods tested unwrap the phase satisfactorily if the initial function is a simple theoretical wavelet. None of the methods restore the phase of a synthetic trace exactly. An initial validity test of the phase-unwrapping method is that the sum of the restored wavelet phase spectrum and the restored pulse-trace phase spectrum (assuming the convolutional model for the seismic trace) must be equal to the restored phase spectrum of the synthetic trace. Results show that none of the tested methods satisfy this test. Quantitative estimation of the phase-unwrapping accuracy by correlation analysis of the phase deconvolution results separated these methods, according to their efficiency, into three groups. The first group consists of methods using a priori wavelet information. These methods make the wavelet phase estimation more effective than the minimum-phase approach, if the wavelet is non-minimum-phase. The second group consists of methods using the phase increment Δø(Δω) between two adjacent frequencies. These methods help to decrease the time shift of the initial synthetic trace relative to the model of the medium. At the same time they degrade the trace correlation with the medium model. The third group consists of methods using an integration of the phase derivative. These methods do not lead to any improvement of the initial seismic trace. The main problem in the phase unwrapping of a seismic trace is the random character of the pulse trace. For this reason methods based on an analysis of the value of Δø(Δω) only, or using an adaptive approach (i.e. as Δω decreases) are not effective. In addition, methods based on integration of the phase derivative are unreliable, due to errors in numerical integration and differentiation.     

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

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

Estimation of the seismic wavelet through homomorphic deconvolution and well log data: application on well-to-seismic tie procedure

Wavelet estimation and well-tie procedures are important tasks in seismic processing and interpretation. Deconvolutional statistical methods to estimate the proper wavelet, in general, are based on the assumptions of the classical convolutional model, which implies a random process reflectivity and a minimum-phase wavelet. The homomorphic deconvolution, however, does not take these premises into account. In this work, we propose an approach to estimate the seismic wavelet using the advantages of the homomorphic deconvolution and the deterministic estimation of the wavelet, which uses both seismic and well log data. The feasibility of this approach is verified on well-to-seismic tie from a real data set from Viking Graben Field, North Sea, Norway. The results show that the wavelet estimated through this methodology produced a higher quality well tie when compared to methods of estimation of the wavelet that consider the classical assumptions of the convolutional model.     

Evaluating Lower Computational Burden Approaches for Calibration of Large Environmental Models

Realistic environmental models used for decision making typically require a highly parameterized approach. Calibration of such models is computationally intensive because widely used parameter estimation approaches require individual forward runs for each parameter adjusted. These runs construct a parameter-to-observation sensitivity, or Jacobian, matrix used to develop candidate parameter upgrades. Parameter estimation algorithms are also commonly adversely affected by numerical noise in the calculated sensitivities within the Jacobian matrix, which can result in unnecessary parameter estimation iterations and less model-to-measurement fit. Ideally, approaches to reduce the computational burden of parameter estimation will also increase the signal-to-noise ratio related to observations influential to the parameter estimation even as the number of forward runs decrease. In this work a simultaneous increments, an iterative ensemble smoother (IES), and a randomized Jacobian approach were compared to a traditional approach that uses a full Jacobian matrix. All approaches were applied to the same model developed for decision making in the Mississippi Alluvial Plain, USA. Both the IES and randomized Jacobian approach achieved a desirable fit and similar parameter fields in many fewer forward runs than the traditional approach; in both cases the fit was obtained in fewer runs than the number of adjustable parameters. The simultaneous increments approach did not perform as well as the other methods due to inability to overcome suboptimal dropping of parameter sensitivities. This work indicates that use of highly efficient algorithms can greatly speed parameter estimation, which in turn increases calibration vetting and utility of realistic models used for decision making.     

新型随机地震动模型

在研究结构的随机地震反应时，要用大量的符合场地条件的地震记录作为输入数据。但强震历史记录却不是每个地区都有的，因此根据符合场地条件的现有地震记录建立随机地震动模型具有重要意义。本文利用中国抗震规范2001版修正选取的样本波作为目标波，考虑了幅值和频率的双重非平稳性，建立了新型随机地震动模型——改进的时变ARMA模型随机地震动模型。通过使用残差的卡方检验法，对多种非平稳ARMA模型生成的模拟波进行检验；同时又比较丁模拟波与目标波的功率谱密度图和反应谱图。结果证明：此法能够更精确地反映不同场地条件地震动的频谱和幅值的真实内容，从而建立符合目标场地条件的更为有效的模拟地震动，为相关研究与工程设计架起一座桥梁。     

子波相位不准对反演结果的影响(英文)

本文重点讨论在振幅谱估计准确的情况下,采用不同相位谱子波作为实际估计子波进行线性最小二乘反演,并对结果进行分析。除子波相位外,所有其它影响反演结果的因素均忽略。稀疏反射系数模型(块状波阻抗模型)反演结果表明:(1)使用不同相位谱子波进行反演,其反演结果合成的记录与原始记录都非常匹配,但反演的反射系数和声波阻抗结果与真实模型有差异;(2)反演结果的可靠程度主要与不同相位子波z变换的根的分布有关,当估计子波与真实子波Z变换的根的分布仅在单位圆附近有差异时,反演的反射系数和声波阻抗与真实模型很接近;(3)尽管反演前后地震记录都匹配了,并且评价反演结果好坏的柯西准则或改进柯西准则(反演参数没有进行自适应处理)已经达到了最优(最小),但反演结果与真实模型仍存在较大差异。最后,针对子波相位估计不准可能导致反演效果较差这个问题,我们提出采用求L1范数、丰度、变分、柯西准则(反演参数进行了自适应处理)或/和改进柯西准则(反演参数进行了自适应处理)的最优值或次优值作为评价准则的一种解决办法,理论上得到了好的效果。     

用遗传算法实现地震信号反褶积

遗传算法作为寻优手段具有全局优化和很好的稳定性．本文将遗传算法用于地震信号反褶积处理，与已往方法相比它具有更好的分辨率和稳定性我们采用Bernoulli－Gaussian模型和ARMA模型分别描述地震反射系数序列和地震子波，用最大似然和最小预测误差准则分别构造用于估计反射系数序列和地震子波的目标函数，用遗传算法优化目标函数，以实现地震信号反褶积．