共查询到20条相似文献，搜索用时 718 毫秒
1.
2.
地震层析成像反演中3种广义解的误差分析与评价 总被引：3，自引：3，他引：0
在地震波速度和地震波 值层析成像中,对于非线性反演问题解的可靠性研究非常重要．为了讨论ＡＲＴ（ＡｌｇｅｂｒａｉｃＲｅｃｏｎｓｔｒｕｃｔｉｏｎ Ｔｅｃｈｎｉｑｕｅ）、ＳＶＤ（Ｓｉｎｇｕｌａｒ ＶａｌｕｅＤｅｃｏｍｐｏｓｉｔｉｏｎ）、ＬＳＱＷ（ＬｅａｓｔＳｕａｒｅ）３种方法的算法稳定性以及反演效果,本研究将这３种算法用于求解非常病态问题和一般稳定问题,根据数值计算结果,分别给出了在非常病态问题和一般稳定问题中３种算法的相对误差及其算法稳定性评价．进而针对３种算法在同一地质模型上的速度成像结果进行了比较,为地震层析成像反演效果分析提供了定量参数． 相似文献
3.
在地震波速度和地震波 值层析成像中,对于非线性反演问题解的可靠性研究非常重要．为了讨论ＡＲＴ（ＡｌｇｅｂｒａｉｃＲｅｃｏｎｓｔｒｕｃｔｉｏｎ Ｔｅｃｈｎｉｑｕｅ）、ＳＶＤ（Ｓｉｎｇｕｌａｒ ＶａｌｕｅＤｅｃｏｍｐｏｓｉｔｉｏｎ）、ＬＳＱＷ（ＬｅａｓｔＳｕａｒｅ）３种方法的算法稳定性以及反演效果,本研究将这３种算法用于求解非常病态问题和一般稳定问题,根据数值计算结果,分别给出了在非常病态问题和一般稳定问题中３种算法的相对误差及其算法稳定性评价．进而针对３种算法在同一地质模型上的速度成像结果进行了比较,为地震层析成像反演效果分析提供了定量参数． 相似文献
4.
Andrzej Kijko 《Pure and Applied Geophysics》1977,115(4):999-1009
Summary An algorithm for the optimum distribution of any number of regional seismic stations is formulated. The information needed from a given region is the relative distribution of seismic intensity and the travel times of seismic waves used in earthquake location procedure. To define the optimum distribution of a regional network, some generalization of D-planning (minimization of the ellipsoid volume of earthquake location errors) was applied. The D-optimum criterion is generalized for a case where equations describing the process are nonlinear and when the most probable values of the sought parameters (origin time and earthquake focus coordinates) are not known.As an example of the algorithm application, the optimum distribution of seismic stations in the central district of the Lublin Coal Basin in Poland is given. 相似文献
5.
通过引入带有二阶正则算子的正则化项,建立了一种双参数混合正则化方法.为确定最佳正则化参数,这里主要应用L-曲线法、偏差原理和广义交叉校验准则的优化组合来确定.首先对理论模型进行了数值模拟,通过与截断奇异值分解法、共轭梯度法及标准Tikhonov正则化法的结果比较,表明该方法不仅精度高,而且对于数据的随机扰动具有稳定性.... 相似文献
6.
主要介绍各种地震定位方法,概述各种地震定位方法的基本原理,重点介绍Geiger的经典方法以及在此基础上建立的各种线性方法:单一地震事件定位法与多个地震事件定位法以及对每一种地震定位方法的应用情况,尤其是国内的应用情况做了总结;同时指出各种方法的特点,并进行相应的比较. 相似文献
7.
福建台网负责监测中国台湾地区地震。对于中深源地震使用何种定位方法能获得较好的地震参数,这直接影响到地震定位精度。利用JOPENS系统中交互分析软件MSDP提供的定位方法,对同一地震进行两次定位,即不使用和使用接入的台湾台站,将福建台网得出的两次结果与中国台湾公布的地震参数进行对比,分析定位精度,进而找出适用于台湾地区中深源地震的定位方法,以便进一步判断在地震速报中使用这些台站进行辅助定位的可行性,并给出相关的操作方法及建议。 相似文献
8.
烃源岩的定量地震刻画对于勘探开发区块的优选、盆地油气资源量的估算都具有重要意义.陆相沉积环境下的浅湖或半深湖相的烃源岩横向变化快,其空间展布需要依靠钻井约束下的反射地震进行刻画,但是其地震弹性特征与岩性和有机质含量的映射关系呈现高度非线性化,因而很难利用传统基于地震岩石物理模型驱动的烃源岩地震预测方法进行有效刻画.本文以低勘探区的东海盆地长江坳陷为例,提出了一种在数据驱动的机器学习框架下,综合利用地质约束、钻井录井、测井、地球化学和叠前地震数据进行烃源岩的定量地震刻画的工作流程.其核心思想是利用随机森林集成学习算法对小样本数据表现优异的特征,以井位处的测井弹性数据(纵波速度和密度)、岩性、地球化学标定的总有机碳含量(TOC)为样本标签数据,在地质导向约束下通过随机森林算法生成学习网络,并将该网络与叠前地震反演结果相结合,采取先预测泥岩再预测总有机碳含量的“两步走”策略,完成对烃源岩空间分布及其非均质性的定量地震刻画,并对预测结果的不确定性进行评价.测试结果显示,随机森林算法相较于其他的机器学习算法能够更准确的识别陆相沉积地层的泥岩,并比传统的利用阻抗转化方法获得更可靠的总有机碳含量预测结果. 相似文献
9.
10.
根据大震速报和快速预警实际需要,首先对Inglada线性单事件定位方法进行了适当的改进, 使其在仅有P波到时数据的情况下也能快速定位,且求解过程仅需简单迭代而不用奇异值分解; 其次, 尝试将改进后的方法从单层均匀模型引入到分层均匀模型中的近源台网定位情形,并通过单层均匀和分层均匀两种不同模型的理论实验讨论了该方法的可行性和适用范围; 最后整合了质量高且分布较好的距离2008年汶川MS8.0地震震中最近的强震、微震, 以及川西流动台阵等观测记录资料,对汶川MS8.0主震初始破裂点的时空参数进行了多种模型的定位实验. 结果表明, 改进后的线性单事件定位方法简单、快捷、易用,可广泛应用于近源地震定位,尤其是用于无法得到S波到时的中强以上直至巨大地震的速报、地震现场流动台网的快速定位以及地震的快速预警等. 相似文献
11.
基于谐波入射下的波动理论频域精确解,导出线性时域精确解,并在Matlab环境中编制相应的计算程序;选取8个简化的中硬场地剖面,用LSSRLI-1、精确解和SHAKE2000三种方法计算各场地在不同输入条件下的地震反应。结果表明:程序计算所得地表反应谱和土体剪应变分布与SHAKE2000结果一致;LSSRLI-1方法得到的地表反应谱与前二者结果一致;LSSRLI-1方法在某些情况下得到的土体剪应变分布与另外二者结果存在较大偏差,该偏差对地表反应谱有着不可忽略甚至非常显著的影响。 相似文献
12.
Cauchy priori distribution-based Bayesian AVO reflectivity inversion may lead to sparse estimates that are sensitive to large reflectivities. For the inversion, the computation of the covariance matrix and regularized terms requires prior estimation of model parameters, which makes the iterative inversion weakly nonlinear. At the same time, the relations among the model parameters are assumed linear. Furthermore, the reflectivities, the results of the inversion, or the elastic parameters with cumulative error recovered by integrating reflectivities are not well suited for detecting hydrocarbons and fuids. In contrast, in Bayesian linear AVO inversion, the elastic parameters can be directly extracted from prestack seismic data without linear assumptions for the model parameters. Considering the advantages of the abovementioned methods, the Bayesian AVO reflectivity inversion process is modified and Cauchy distribution is explored as a prior probability distribution and the time-variant covariance is also considered. Finally, we propose a new method for the weakly nonlinear AVO waveform inversion. Furthermore, the linear assumptions are abandoned and elastic parameters, such as P-wave velocity, S-wave velocity, and density, can be directly recovered from seismic data especially for interfaces with large reflectivities. Numerical analysis demonstrates that all the elastic parameters can be estimated from prestack seismic data even when the signal-to-noise ratio of the seismic data is low. 相似文献
13.
以日本KiK-net强震观测台网中硬场地井下记录为样本,对传统等效线性化方法LSSRLI-1（频域）、SHAKE2000,时域非线性方法DEEPSOIL和频域一致等效线性化方法SOILQUAKE等几种计算程序在硬土场地地震反应分析中的可靠性进行对比检验。检验工况包括KiK-net井下台网中地表峰值加速度不小于0.05g的水平硬场地的总计344台次的加速度记录,涉及5个台站,土层厚度6~50m,地表峰值加速度范围0.050~0.805g。结果表明:认为硬土场地以往计算方法能够体现土层放大的认识是片面的,在中强地震动情况下,现有流行方法计算出的地表响应会偏小,强地震动下会严重偏小,会给工程抗震设计提供偏于危险的输入;硬场地中烈度8度以下（地表PGA在0.19g以下）,SHAKE2000和DEEPSOIL、LSSRLI-1（频域）计算结果与实际记录差距可以接受,但烈度8度以上,计算出地表响应均较实际记录明显偏小,且随地震动强度增加差距急剧增大;硬场地中,无论何种烈度,SOILQUAKE16计算的地表加速度响应与实测相当,可体现出土层放大作用。 相似文献
14.
研究了六层层状介质模型瑞利波基阶和高模式波相速度对横波速度、深度的敏感性,结果表明：基阶波较高模式波对7 m以内浅部地层的横波速度更敏感,敏感性频带在10～25 Hz范围内,峰值频带集中在18 Hz左右;高模式波较基阶波对深部地层的横波速度更敏感,敏感性频带宽,峰值分散.基阶波对浅层的敏感性和高模式波穿透深度更深的特点为近地表岩土层二维横波速度结构的联合反演提供了前提条件.利用阻尼最小二乘SVD（Singular Value Decomposition）算法联合基阶与高模式波对理论模型和实例数据进行横波速度反演,反演结果表明联合反演增强了反演的稳定性,提高了反演的精度. 相似文献
15.
— A P-wave tomographic method for 3-D complex media (3-D distribution of elastic parameters and curved interfaces) with orthorhombic symmetry is presented in this paper. The technique uses an iterative linear approach to the nonlinear travel-time inversion problem. The hypothesis of orthorhombic anisotropy and 3-D inhomogeneity increases the set of parameters describing the model dramatically compared to the isotropic case. Assuming a Factorized Anisotropic Inhomogeneous (FAI) medium and weak anisotropy, we solve the forward problem by a perturbation approach. We use a finite element approach in which the FAI medium is divided into a set of elements with polynomial elastic parameter distributions. Inside each element, analytical expressions for rays and travel times, valid to first-order, are given for P waves in orthorhombic inhomogeneous media. More complex media can be modeled by introducing interfaces separating FAI media with different elastic properties. Simple formulae are given for the Fréchet derivatives of the travel time with respect to the elastic parameters and the interface parameters. In the weak anisotropy hypothesis the P-wave travel times are sensitive only to a subset of the orthorhombic parameters: the six P-wave elastic parameters and the three Euler angles defining the orientation of the mirror planes of symmetry. The P-wave travel times are inverted by minimizing in terms of least-squares the misfit between the observed and calculated travel times. The solution is approached using a Singular Value Decomposition (SVD). The stability of the inversion is ensured by making use of suitable a priori information and/or by applying regularization. The technique is applied to two synthetic data sets, simulating simple Vertical Seismic Profile (VSP) experiments. The examples demonstrate the necessity of good 3-D ray coverage when considering complex anisotropic symmetry. 相似文献
16.
从瞬变电磁响应中提取IP信息的研究(英文) 总被引：1，自引：0，他引：1
在对浙江省某铜矿进行瞬变电磁勘查时,由于感应激发极化效应的影响,造成瞬变电磁晚期测道的数据发生倒转。采用常规的瞬变电磁数据处理方法难以处理,造成晚期测道数据不可用,影响TEM的探测深度及精度。针对这个情况,采用Cole.Cole模型分析了均匀半空问模型中直流电阻率、充电率、时间常数以及频率相关系数对磁性源瞬变电磁响应的影响规律。利用奇异值分解法(The Singular Value Decomposition,简称SVD)对实测TEM数据进行分析及反演,并且从瞬变电磁响应数据中分离出勘探区瞬变电磁测量数据Cole—Cole模型的各参数,将分离出来的参数用于探测结果的辅助解释。经过分析得出充电率和频率相关系数对瞬变电磁响应的影响较大,而直流电阻率和时间常数对其影响较小的结论。通过常规方法和奇异值分解法对实测数据分别进行处理,发现后者视电阻率断面图的异常更加突出,从而较准确的圈定了铜矿体地范围,此结论得到_『钻孔的验证,与实际地质情况吻合较好。 相似文献
17.
The current analysis of earthquakes is typically based on linear mathematical models that may fail to describe and forecast particular behaviors, because in many cases the data complexity may induce a highly non linear behavior. In this paper the implementation of an alternative method for seismic time series analysis is presented. The RPs (Recurrence Plots) enables recognition and treatment of measured accelerations. An RP obtained from seismic data allows a more efficient interpretation of the ground motions and this explanation contributes to characterize materials and responses. The nonlinear attributes from RPs analysis can be used as filters to reveal patterns or be combined to predict a seismic property. Automated seismic data characterization, based on nonlinear seismic attributes, could rewrite the rules of earthquake phenomena interpretation. The objective of this work is to establish a new methodology for practical application of nonlinear dynamics in seismic pattern/attributes recognition, an evolving and challenging engineering field. 相似文献
18.
This paper introduces the eigenspace structural identification technique for tall buildings subjected to ambient excitations that are stationary and where only the response time histories are measured. Based on the forward innovation model of the Kalman filter sequence, the actual response can be constructed as a function of the measured response time history with contamination of either displacement or velocity. The response time history is decomposed into subspace matrices using QR decomposition and Quotient Singular Value Decomposition (QSVD) techniques. These are then substituted into the least-square formulation to obtain the solution which is non-unique. Similarity transformation is applied to arrive at the desired solution employing the fact that eigenvalues of self-similar systems are identical. The advantages of this eigenspace technique are that it is non-iterative, initial estimates of the parameters to the identified are not required, well-established numerical algorithm of the decomposition techniques employed are available, and the method can handle MDOF systems efficiently. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
19.
20.
Multi-source seismic technology is an efficient seismic acquisition method that requires a group of blended seismic data to be separated into single-source seismic data for subsequent processing. The separation of blended seismic data is a linear inverse problem. According to the relationship between the shooting number and the simultaneous source number of the acquisition system, this separation of blended seismic data is divided into an easily determined or overdetermined linear inverse problem and an underdetermined linear inverse problem that is difficult to solve. For the latter, this paper presents an optimization method that imposes the sparsity constraint on wavefields to construct the object function of inversion, and the problem is solved by using the iterative thresholding method. For the most extremely underdetermined separation problem with single-shooting and multiple sources, this paper presents a method of pseudo-deblending with random noise filtering. In this method, approximate common shot gathers are received through the pseudo-deblending process, and the random noises that appear when the approximate common shot gathers are sorted into common receiver gathers are eliminated through filtering methods. The separation methods proposed in this paper are applied to three types of numerical simulation data, including pure data without noise, data with random noise, and data with linear regular noise to obtain satisfactory results. The noise suppression effects of these methods are sufficient, particularly with single-shooting blended seismic data, which verifies the effectiveness of the proposed methods. 相似文献