共查询到18条相似文献，搜索用时 906 毫秒
1.
混杂的全局最优化算法是一种模拟退火法和下坡的单纯形法相结合的方法，该法用来反演远震体波求解震源参数，在时间域对震源时间函数加约束，反演双力偶参数而不是反演矩张互通主考虑多重震源导致非线性和多重模型问题，在这些问题中，目标函数包括许多局部极小值。 相似文献
2.
震源时间函数与震源破裂过程 总被引：2，自引：0，他引：2
简要介绍了现代数字地震学的两个重要概念。一个是震源时间函数，另一个是地震震源的时空破裂过程。首先，从地震断层位移表示定量出发，介绍震源时间函数和震源破裂过程的基本概念。然后，分别介绍两种从远场地震记录中提取震源时间函数和获取有限断层面上时空破裂过程图像的反演方法。 相似文献
3.
在用矩张量反演方法确定地震的震源机制和震源的破裂过程时，由于测定震源深度的误差，很难保证计算格林函数时所采用的震源深度(理论震源深度)和实际震源深度完全相等．理论震源深度和实际震源深度的差异(震源深度误差)对矩张量的反演或多或少会造成影响．借助于合成地震图较系统地讨论了震源深度误差对3种基本类型(正断层、逆断层和走滑断层)的震源反演造成的影响．对于正断层和逆断层，震源深度误差主要影响各向同性成分和补偿线性向量偶极成分．在理论震源深度大于实际震源深度时:对于正断层，震源深度误差导致出现虚假的正的各向同性分量和负的补偿线性向量偶极分量；随着差异的增大其绝对值逐渐增大；对于逆断层，震源深度误差导致出现虚假的负的各向同性分量和正的补偿线性向量偶极分量，其绝对值随着震源深度误差的增加而增加．在理论震源深度小于实际震源深度时，结果恰好相反:对于正断层，震源深度误差产生了虚假的负的各向同性分量和正的补偿线性向量偶极分量；对于逆断层，震源深度误差产生了虚假的正的各向同性分量和负的补偿线性向量偶极分量．同样，随着深度误差的增大，它们的绝对值也逐渐增大．对于走滑断层，震源深度误差对矩张量反演的影响则不同于正断层和逆断层两种情况，受影响最大的是震源时间函数的形状．无论理论震源深度小于还是大于实际震源深度，震源时间函数的形状都发生不同程度的拖尾现象．数字试验表明，当震源深度误差小于20 km时对于地震的总体机制的反演没有明显影响．另外，当理论震源深度大于实际震源深度时，震源深度误差对矩张量反演的影响相对较小． 相似文献
4.
震源深度是核试验以及塌陷等浅源地震研究中的关键参数,可以为事件成因分析提供关键信息.然而朝鲜核试验区域地形起伏较大,地形效应可能对震源深度反演的结果造成影响.本文基于理论地震图进行测试,研究了地形起伏对震源深度反演的影响.发现震源深度小于2km时,不考虑地形影响,反演得到的震源深度会系统偏浅0.2km左右.然后利用MDJ2速度结构模型,我们反演了2017年9月3日朝鲜M6.3事件的震源参数,结果显示震源深度约为0.8km.进一步基于带地形的格林函数重新反演了该事件的震源深度,发现在1km处波形拟合结果较好.不同速度模型测试结果显示该事件的震源深度反演误差约为1km.案例研究表明,基于层状均匀速度模型,利用区域地震波形资料反演的震源深度可以为浅源事件成因分析提供关键约束. 相似文献
5.
近震源破裂过程反演研究--Ⅱ.9.21中国台湾集集地震破裂过程的近场反演 总被引：1，自引：0，他引：1
根据非线性规划研究的最新成果所设计的一种全新的震源破裂过程的反演方法, 用近场地震波观测资料反演了1999年9月21日发生在中国台湾省集集Mw7.6地震震源破裂过程.为了使反演中设置的断层模型与集集地震实际破裂面尽可能一致, 以尽可能减小由于断层模型设置的不确定性对震源破裂过程反演结果的影响, 设立的断层模型为与集集地震造成的主要地表破裂尽可能拟合的弯曲面模型.反演结果显示: (1) 集集地震震源的破裂大体持续了32 s, 其中主要破裂发生在第6～27 s间, 破裂主要集中发生在断层北段向东拐弯处.(2) 震源破裂以逆冲为主, 平均滑动角为64.5°, 与USGS, Harvard及CWB(台湾中央气象局)的结果相当.标量地震矩为7.76×1020牛顿米, 稍大于USGS和Harvard反演的标量地震矩.(3) 集集地震震源的破裂存在清晰的成核过程, 成核过程经历6 s后, 地震矩释放明显加速.起始破裂从断层南段开始, 10 s后破裂主要集中在断层北段发生.最后将反演结果与震后GPS观测结果进行了对照分析, 并对反演结果的科学意义进行了讨论. 相似文献
6.
由于深部金属矿埋深和自身的复杂性,利用重、磁、电方法和一般的地震方法很难有效地对其进行高精度定位.全波形反演通过最小化模拟数据与观测数据的差异使深部金属矿的高精度探测成为可能,但全波形反演是一个局部优化过程,需要准确的低频数据作为起始,而这在一般的地震数据采集中难以做到.本文先在频域中使用伴随状态震源函数反演方法,通过震源附近的直达波能精确地反演出震源函数的形态.然后利用得到的高精度震源子波结合褶积与反褶积思想及相似性现象重构含有低频成分的自激自收数据.将该数据应用到全波形反演中,有效缓减了反演过程中出现的周波跳跃现象,并提高了模型反演的正确性.Marmousi模型和金属矿模型的数值模拟实验证明了新方法改善了在没有低频数据时的全波形反演结果,并有较好抗噪性. 相似文献
7.
近期发展的震源反演方法并没有充分利用区域地震图的宽频带性质,原因是用平行层模型计算的格林函数和真实地壳的复杂传播现象不同。这种反演通常去掉了短周期成分,并且只反演记录中选定的部分。在本文中,我们介绍一种估算震源参数的方法,它可以较好地应用整个宽频带记录,而且只采用目前仅有的不完善的格林函数即可。这个方法通过分别拟合格林函数的各个部分来使得记录图与合成图主要的地壳到达波到时差别不大。确定震源参数的原 相似文献
8.
根据非线性规划研究的最新成果所设计的一种全新的震源破裂过程的反演方法, 用近场地震波观测资料反演了1999年9月21日发生在中国台湾省集集Mw7.6地震震源破裂过程. 为了使反演中设置的断层模型与集集地震实际破裂面尽可能一致, 以尽可能减小由于断层模型设置的不确定性对震源破裂过程反演结果的影响, 设立的断层模型为与集集地震造成的主要地表破裂尽可能拟合的弯曲面模型. 反演结果显示: (1) 集集地震震源的破裂大体持续了32 s, 其中主要破裂发生在第6~27 s间, 破裂主要集中发生在断层北段向东拐弯处. (2) 震源破裂以逆冲为主, 平均滑动角为64.5°, 与USGS, Harvard及CWB(台湾中央气象局)的结果相当. 标量地震矩为7.76×1020牛顿米, 稍大于USGS和Harvard反演的标量地震矩. (3) 集集地震震源的破裂存在清晰的成核过程, 成核过程经历6 s后, 地震矩释放明显加速. 起始破裂从断层南段开始, 10 s后破裂主要集中在断层北段发生. 最后将反演结果与震后GPS观测结果进行了对照分析, 并对反演结果的科学意义进行了讨论. 相似文献
9.
用近震源波形资料拟合反演地震的震源破裂过程,所包含的一些不确定因素将对反演结果的精度及可靠性产生影响,文中的数值实验分析了所假定的反演断层模型参数的某些不确定性对反演结果的影响程度,并对观测波形的截取长度对反演精度的影响进行了讨论．结果表明：(1)近震源地震波形资料能较好地分辨断层浅部的破裂过程.然而对断层深部的位错分布的约束和反演能力较差．联合使用近、远场地震波资料进行反演,能反演出一个更为完全的整个断层破裂过程的图像.(2)用近震源地震波资料反演时,反演结果对所假定的反演断层的走向和倾角非常敏感.断层走向偏离真实值2°或倾角偏离真实值5°都会导致一个虚假的反演结果.(3)反演中所使用的介质速度结构模型的不确定性,也会对反演结果产生影响. 相似文献
10.
本文提出解调包络方法来重构地震记录中缺失的低频信号,同时该方法能够降低全波形反演的非线性程度;提出伴随状态震源函数反演方法来得到精确的震源函数,并推导了梯度计算公式;解调包络方法结合低通滤波技术,实现了从低频到高频的多尺度反演策略,有效缓解了全波形反演的周波跳跃问题.数值算例证明了解调包络、伴随状态震源函数反演方法和低通滤波多尺度反演策略的可行性及优越性.震源函数反演精度测试结果表明:即使观测记录在缺失低频信息的情况下,也能反演得到精确的震源函数.缺失低频测试和抗噪能力测试结果表明:即使地震数据中缺失9Hz以下的低频信号或者信噪比极低的情况下,利用反演得到的精确震源函数进行解调包络多尺度全波形反演,同样可以得到高精度的全波形反演结果.与Hilbert包络全波形反演对比结果表明:解调包络在重构低频和降低伴随震源主频方面具有一定优势. 相似文献
11.
基于L2范数的常规全波形反演目标函数是一个强非线性泛函,在反演过程中容易陷入局部极小值.本文提出归一化能量谱目标函数来缓解全波形反演过程中的强非线性问题,同时能够有效地缓解噪声和震源子波不准等因素的影响.能量谱目标函数是通过匹配观测数据与模拟数据随频率分布的能量信息来实现最小二乘反演的,其忽略了地震数据波形与相位变化的细节特征,这在反演的过程中能够有效缓解波形匹配错位等问题.数值测试结果表明,基于归一化能量谱目标函数在构建初始速度模型、抗噪性和缓解震源子波依赖等方面都优于归一化全波形反演目标函数.金属矿模型测试结果表明,即使地震数据缺失低频分量,基于归一化能量谱目标函数的全波形反演方法在像金属矿这样的强散射介质反演问题上同样具有一定的优势. 相似文献
12.
用平面波延拓方程进行地震数据的叠前速度反演 总被引：1，自引：1，他引：1
本文讨论地震勘探数据的叠前速度反演方法及其在海洋地震勘探数据上的反演试验.反演主要的计算步骤是:1.采用Fourier-Hankel变换把球面波分解为平面谐波;2.用平面谐波的延拓方程将上行波与下行波同时向下延拓,并计算每一层底部的反射系数和下一层的波阻抗;3.用最小二乘法从波阻抗中确定该层的声波速度.重复第2步与第3步,直到某一预定深度时结束.通过反演试验,对地震振幅比例的改变,子波变形,以及第1层速度和密度的误差对反演方法的稳定性及其精度的影响进行了分析.还通过实际海洋地震勘探数据的反演试验,对这一方法在地震勘探中的应用前景作了论述. 相似文献
13.
由于地震资料是带限的,而井资料约束反演过程中的地震反演频带的展宽也是有一定限度的,使得反演地震波阻抗也是带限的。如何进一步拓宽反演地震波阻抗的频带宽度,提高地震资料的分辨率是文中研究的内容。利用井约束下的地震反演波阻抗资料、地震资料及井资料,借鉴地质模拟技术,研究如何利用地层的空间相对稳定性特征,通过相关函数的分析,估计反演更高分辨率的地震波阻抗资料。考虑到地质体空间的相关性及分辨率,相关半径取相关函数从最大点到第1个极值点的范围。在目标函数的计算过程中,既考虑当前计算点的权系数信息,又考虑井外推时的权系数信息,这样使计算结果更加稳定可靠。通过综合模拟计算,能够识别较小砂体 相似文献
14.
15.
B. URSIN 《Geophysical Prospecting》1986,34(8):1213-1218
The one-dimensional seismic inverse problem consists of recovering the acoustic impedance (or reflectivity function) as a function of traveltime from the reflection response of a horizontally layered medium excited by a plane-wave impulsive source. Most seismic sources behave like point sources, and the data must be corrected for geometrical spreading before the inversion procedure is applied. This correction is usually not exact because the geometrical spreading is different for primary and multiple reflections. An improved algorithm is proposed which takes the geometrical spreading from a point source into account. The zero-offset reflection response from a stack of homogeneous layers of variable thickness is used to compute the thickness, velocity and density of each layer. This is possible because the geometrical spreading contains additional information about the velocities. 相似文献
16.
17.
JAN HELGESEN 《Geophysical Prospecting》1991,39(3):313-336
Three methods for least-squares inversion of receiver array-filtered seismic data are investigated: (1) point receiver inversion where array effects are neglected; (2) preprocessing of the data with an inverse array filter, followed by point receiver inversion; (3) array inversion, where the array effects are included in the forward modelling. The methods are tested on synthetic data generated using the acoustic wave equation and a horizontally stratified earth model. It is assumed that the group length and the group interval are identical. For arrays that are shorter than the minimum wavelength of the emitted wavefield, and when the data are appropriately muted, point receiver inversion (first method) gives satisfactory results. For longer arrays, array inversion (third method) should be used. The failure of the inverse array filter (second method) is due to aliasing problems in the data. 相似文献
18.
应用地震资料约束下大地电磁资料反演古潜山或基底内部结构的新方法. 首先利用地震及井资料,反演得到古潜山或基底之上各层的厚度及电阻率,建立了古潜山或基底以上的地电模型;并利用大地电磁一维反演资料,建立了古潜山或基底及其以下的地电模型. 在初始模型建立的基础上,应用高效快速的二维大地电磁正演算法——模式匹配方法,进行正演计算. 利用共扼梯度迭代方法,进行正反演迭代计算. 通过对胜利油田南北618剖面反演, 基底以下在地震资料没有显示的构造信息,在反演结果中得到了较好地反映. 相似文献