共查询到17条相似文献，搜索用时 500 毫秒
1.
采用新近提出的多次波射线追踪正演算法,结合共轭梯度法求解带约束的阻尼最小二乘最优化反演问题,分析讨论了利用多震相走时资料进行联合反演成像的方法及技术.考虑到不同震相走时的拾取误差不同,反演算法中引入了不同震相种类数据的权系数; 由于同时反演速度模型和反射界面起伏中不同模型参数变化对走时影响程度的不同, Jacobi偏导矩阵元素中引入了不同参数的归一化因子; 另外,为了克服射线密度过大(或过小)区域速度模型的过度(或欠)更新问题,反演算法中引入了等权射线密度的概念.几种数值模拟实例表明(含噪声敏感性试验): 多波走时的联合或同时反演成像技术是一种提高走时成像空间分辨率,进而降低重建模型失真度行之有效的方法. 相似文献
2.
利用地震走时数据,采用联合反演方法获取了江苏地区的一维P波速度模型。与仅采用初至波走时的传统天然地震走时获取方法相比,该方法充分利用了大量存在的续至波参与反演,能有效改进中下地壳的反演能力。针对地震震相目录中常存在震相标识错误的问题,采用的自动判别筛选震相方法能最大限度提高数据走时的精度,可以对不同震相进行有效区分。与其他常用一维速度模型相比,本文反演的模型对Pg、Pn震相走时拟合效果最佳,残差最小。当所用走时数据拥有较高定位精度时,该反演方法能为研究区三维速度结构成像和地震定位提供较可靠的一维速度模型。 相似文献
3.
中国大陆中西部乃至全球造山带普遍具有复杂地壳结构.随着矿产资源勘探和深部探测研究的深入,探测造山带及盆山耦合区下方地壳精细结构正逐渐成为当前面临的巨大挑战.人工源深地震测深方法正越来越清晰地揭示出不同构造域地壳速度结构的基本特征,然而传统的层状结构模型参数化方法难以准确描述复杂地质模型,通常情况下多忽略速度结构的精细间断面且采用层边界平滑处理,难以满足地壳精细结构成像的发展要求.针对上述困难,本文采用最近发展的块状结构建模方案构建三维复杂地壳模型,基于逐段迭代射线追踪正演走时计算方法,推导了走时对三角形界面深度以及网格速度的偏导数,开展了非线性共轭梯度走时反演方法研究.发展了利用直达波和反射波等多震相走时数据对界面深度和网格速度的多参数联合反演方法,并引人不同种类震相数据的权系数和不同类型参数偏导数归一化的方法.数值算例表明,基于块状结构的非线性共轭梯度走时反演方法适用于复杂地壳结构模型,在利用人工源走时数据反演复杂地壳精细结构领域具有良好的应用前景. 相似文献
4.
采用新近研制的分区多步不规则最短路径多震相地震射线追踪正演技术,结合流行的子空间反演算法,提出了一种联合多震相走时资料进行地震三参数 （速度、反射界面和震源位置） 同时反演的方法技术.数值模拟反演实例、以及与双参数 （速度和反射界面或速度和震源位置） 同时反演的对比分析表明：三参数同时反演成像结果大体接近双参数同时反演成像的结果.另外,噪声敏感性试验表明：所提算法对到时数据中可容许的随机误差并不敏感,结果说明多震相走时的联合三参数同时反演成像方法技术不失为一种提高走时成像空间分辨率、进而降低重建模型参数失真度、行之有效的方法技术. 相似文献
5.
球坐标系下多震相走时三参数(速度、震源位置和反射界面)同时反演需要解决两个关键问题:(1)球坐标系下3D速度模型中多次透射、反射(折射)及转换波精确、快速的射线追踪;(2)同时反演时三种不同参数间的强耦合问题.为此,我们将直角坐标系下分区多步不规则最短路径算法推广至球坐标系中,进行区域或者全球尺度的多震相射线追踪.然后将其与适合多参数同时反演的子空间算法相结合,形成一种球坐标系下联合多震相走时三参数同时反演的方法技术.与双参数(速度和反射界面或速度和震源位置)同时反演的数值模拟对比分析显示:三参数与双参数的同时反演结果大体接近,并且它们对到时数据中可容许的随机噪声不太敏感.结果说明本文中的同时反演成像为一种提高成像分辨率,同时反演速度、震源位置和反射界面的有效方法. 相似文献
6.
由多震相走时和波形作三维地震成像 总被引：1，自引：1，他引：0
本项成果包括：提出天然地震走时反演层析成像技术,采用下列方法使得处理结果得以改善：１）利用３２个地震,４０４条射线的Ｐｇ,Ｓｇ,Ｐｍ,Ｓｍ,Ｐｎ,Ｓｎ等震相增加约束条件：２）用已有精度较高的人工地震测深化结果作速度约束;３）用波形反演来修改模型,把波源,介质吸收,散射等物理特征集中反映在记录中,把诸多物理量开发出来互为约束,以修改后的模型再作反演,使解的稳定性大大提高;４）采用最优化过程,选择遗传 相似文献
7.
多震相走时和波形三维地震成象方法及在祁连山地壳结构研究中的应用 总被引：1，自引：0，他引：1
本文的理论方法是以几何射线理论为基础发展起来的、天然地震走时反演技术及天然地震层析成像技术。它存在着震源函数与介质参数的解耦问题。本研究采取了五种方法来改善反演结果。包括,利用,Pg,Pn等震相增大约束条件;用已有精度较高的人工地震测深结果作速度约柬：用波形反演来修改模型,把诸多物理量开发出来互为约束,以修改后的模型再作反演,使解的稳定性大大提高：采用最优化过程,选择遗传算法。可以进行震源定位,走时反演,波形反演：得到任意深度的速度分布及从地表到Moho面的速度剖面。用于在祁连山地区的结果表明,这些层析剖面对认识大地构造、重大深部事件动力学是很有益的。 相似文献
8.
高频假设下的地震射线理论以及相应的地震成像理论表明,在射线稀疏条件下,不可能得到较高分辨率的构造成像;而有限频射线理论更符合实际地震的传播规律,即地震波的走时不仅与中心射线(传统的几何射线)上的速度分布有关,而且与中心射线附近一定范围(称其为第一菲涅耳体)内的速度异常分布有关.鉴于此,本文提出了计算多震相地震波菲涅耳体有限频射线的方法,并定义了走时敏感核函数,同时给出了利用多震相菲涅耳体有限频射线进行速度模型和反射界面同时反演成像的公式.利用多震相走时资料,使用传统射线层析成像方法与有限频射线层析成像方法进行了速度和界面的同时反演成像.结果表明,当射线密度较小时,无论是对速度模型的重建还是对反射界面几何形状的更新,有限频射线层析成像方法均优于传统射线层析成像方法, 而变频有限频射线层析成像则是实际地震层析成像的首选反演算法. 相似文献
9.
WTW(Wave equation traveltime+Waveform inversion)反演是基于波动方程的走时反演（WT反演）和波形反演的联合反演方法.WT反演利用波动方程计算走时和走时关于速度的导数，和传统以射线为基础的走时反演相比，具有不必射线追踪、不必拾取初至、不必高频假设以及初始模型和实际模型差别较大时也能较好收敛等优点，但WT反演与波形反演相比其结果分辨率低.与之互补的是，波形反演的反演结果分辨率高，但是当所给初始模型和实际模型相差太大时，波形反演迭代算法容易陷入局部极小点.可见结合两种方法的WTW反演是一种比较好的联合反演方法.常规WTW迭代算法是首先以WT反演为主反演得到地质模型的整体特征，然后再以波形反演为主反演模型细节，该算法耗时和占用计算机存储空间接近WT反演或波形反演的两倍.为了节省运算耗时和计算机存储空间，往往采取首先单独利用WT反演然后再单独利用波形反演的算法.这样做的缺点是不能紧密结合两种反演方法，使得它们的优缺点在每一次迭代中无法得到互补，从而影响了最终的反演结果.针对以上事实，本文提出一种新的方法实现WTW，使得WTW运算速度和存储空间在任何情况下等同于WT反演或波形反演.模型计算表明新的算法具有更好的收敛性. 相似文献
10.
海洋地球物理勘探中常使用多道地震数据建立浅部地层模型,结合海底地震仪广角地震数据进行射线追踪和走时反演求取地层速度特征.受采集环境和成本限制,部分海底地震仪测线缺少同步采集的多道地震数据,建立的初始模型其浅层结构有较大不确定性,可能会对正演模拟结果造成严重影响.本研究针对海底地震仪广角地震走时模拟中浅部结构的不确定性因素,使用RayInvr软件构建了理论模型和走时数据体,分析讨论了初始模型沉积层厚度与速度变化对正演模拟结果的影响.之后进行实测数据处理,在MCS2019-3测线数据的约束下建立初始模型,对OBS2011-1测线震相重新进行了走时拟合,将得到的模型结果与先前结果进行对比,得到以下结论:沉积层厚度及速度不确定性对走时模拟结果的浅部结构影响较大,深部结构影响较小;使用错误的速度与基底深度组合来拟合沉积层反射震相PsP,会使台站下方的地壳和地幔折射震相Pg和Pn,以及莫霍面反射震相PmP都出现走时小起伏,走时提前对应沉积层速度与厚度偏大,走时延后对应沉积层速度与厚度偏小;PsP震相缺失时,不同的速度-深度组合都可以很好地拟合震相,但会增加结果的不确定性. 相似文献
11.
3-D multi-parameter type traveltime tomography in a spherical coordinate frame: comparison of double and triple class simultaneous inversions 
下载免费PDF全文
下载免费PDF全文 It is now common practice to perform simultaneous traveltime inversion for the velocity field and the reflector geometry in reflection/refraction tomography, or the velocity field and the hypocenter locations in regional earthquake tomography, but seldom are all three classes of model parameters updated simultaneously. This is mainly due to the trade-off between the different types of model parameters and the lack of different seismic phases to constrain the model parameters. Using a spherical-coordinate ray tracing algorithm for first and later(primary reflected) arrival tracing algorithm in combination with a popular linearized inversion solver, it is possible to simultaneously recover the three classes of model parameters in regional or global tomographic studies. In this paper we incorporate the multistage irregular shortest-path ray tracing algorithm(in a spherical coordinate system) with a subspace inversion solver to formulate a simultaneous inversion algorithm for triple model parameters updating using direct and later arrival time information.Comparison tests for two sets of data(noise free and added noise) indicate that the new triple-class parameter inversion algorithm is capable of obtaining nearly the same results as the double-class parameter inversion scheme. Furthermore,the proposed multi-parameter type inversion method is not sensitive to a modest level of picking error in the traveltime data, and also performs well with a relatively large uncertainty in earthquake hypocentral locations. This shows it to be a feasible and promising approach in regional or global tomographic applications. 相似文献
12.
3-D simultaneous inversion for velocity and reflector geometry using multiple classes of arrivals in a spherical coordinate frame 总被引：1，自引：0，他引：1
Guo-jiao Huang Chao-ying Bai Xing-wang Li Stewart Greenhalgh 《Journal of Seismology》2014,18(1):123-135
The two key requirements in conducting 3-D simultaneous traveltime tomography on real data at the regional and global scale with multiple classes of arrival time information are (1) it needs an efficient and accurate arrival tracking algorithm for multiply transmitted, reflected (or refracted) and converted waves in a 3-D variable velocity model with embedded velocity discontinuities (or subsurface interfaces), and (2) a subdimensional inversion solver is required which can easily search for different types of model parameters to balance the trade-off between the different types of model parameter updated in the simultaneous inversion process. For these purposes, we first extend a popular grid/cell-based wavefront expanding ray tracing algorithm (the multistage irregular shortest-path ray tracing method), which previously worked only in Cartesian coordinate at the local scale, to spherical coordinates appropriate to the regional or global scale. We then incorporated a fashionable inversion solver (the subspace method) to formulate a simultaneous inversion algorithm, in which the multiple classes of arrivals (including direct and reflected arrivals from different velocity discontinuities) can be used to simultaneously update both the velocity fields and the reflector geometries. Numerical tests indicate that the new inversion method is both applicable and flexible in terms of computational efficiency and solution accuracy, and is not sensitive to a modest level of noise in the traveltime data. It offers several potential benefits over existing schemes for real data seismic imaging. 相似文献
13.
基于地质统计先验信息的储层物性参数同步反演 总被引：4，自引：1，他引：3
本文提出的储层物性参数同步反演是一种高分辨率的非线性反演方法,该方法综合利用岩石物理和地质统计先验信息,在贝叶斯理论框架下,首先通过变差结构分析得到合理的变差函数,进而利用快速傅里叶滑动平均模拟算法（Fast Fourier TransformMoving Average,FFT-MA）和逐渐变形算法（Gradual Deformation Method,GDM）得到基于地质统计学的储层物性参数先验信息,然后根据统计岩石物理模型建立弹性参数与储层物性参数之间的关系,构建似然函数,最终利用Metropolis算法实现后验概率密度的抽样,得到物性参数反演结果。并将此方法处理了中国陆上探区的一块实际资料,本方法的反演结果具有较高的分辨率,与测井数据吻合度较高;由于可以直接反演储层物性参数,避免了误差的累积,大大减少了不确定性的传递,且计算效率较高。 相似文献
14.
15.
Mrta Kis 《Journal of Applied Geophysics》2002,50(4)
Horizontally layered (1D) earth models are often assumed as a model estimate for the interpretation of geophysical data measured along 2D geological structures. In this process, the individual data sets are usually inverted independently, and it is considered only in a later phase of interpretation that these local (1D) models have common characteristic features. Taking account of these common attributes, instead of the successive independent interpretations, the lateral variations of geometrical and petrophysical parameters can be efficiently determined for the whole 2D structure by applying a series expansion. Using global basis functions, two advantages can be achieved: (i) choosing an appropriate number of basis functions helps us to restrict the complexity of the model; (ii) the integration of all the data sets measured along the profile gives rise to the application of simultaneous or joint inversion methods. This results in a decrease of the number of independent unknowns, a higher stability during the inversion and a more accurate and reliable parameter estimation.In this paper, a joint inversion algorithm is presented using DC geoelectric apparent resistivities and refraction seismic travel times measured along various layouts above a 2D geological model. To describe lateral variations series, expansions are used, and furthermore, to improve the often used approximation of a (locally) 1D forward modelling, the integral mean value of the horizontally changing model parameters (calculated along an appropriately defined interval) is introduced. We call the inversion procedure that combines series expansions and the concept of integral mean Generalised Series Expansion (GSE) inversion. The method was developed and tested for both the simultaneous (integrating data sets of one method or methods on the same physical basis) and the joint inversion (where data sets of methods on different physical bases are joined together), using synthetic and field data sets. It is also demonstrated that the equivalence problem inherent in the independent inversion of DC geoelectric data can efficiently be resolved by the use of the joint GSE inversion method in the cases of conductive and resistive equivalent geological models. 相似文献
16.
跨孔雷达全波形反演是一种使用全波形信息反演两钻孔之间地下信息的层析成像技术.常规的层析成像反演大部分采用射线追踪方法,其中基于初至时的射线追踪方法可以反演出速度剖面（介电常数剖面）,基于最大振幅的层析成像可以反演出衰减剖面（电导率剖面）.常规射线追踪方法有许多不足,究其原因是该方法仅使用了小部分的信号信息.为了进一步提高成像分辨率,本文全面推导了全波形跨孔雷达层析成像反演方法,该方法利用雷达波全幅度相位信息能够反演出地下高分辨率的介电常数和电导率图像.本文通过基于局域网的分布式并行算法,有效地解决了巨量数据正演计算问题.文中首先建立了基于单轴各向异性介质完全匹配层的时间域有限差分二维正演算法,进而通过应用包括时间维度在内的全波场信息与残场逆向传播的全波场信息乘积来计算梯度方向,通过求取以步长为自变量的目标函数的极值确定步长公式,并提出以第一次介电常数反演作为同步反演的初始模型,能够有效提高收敛速度.本文对多组模型进行成像实验,取得了较好的反演效果. 相似文献
17.
为更好地处理与解释复杂海底地形条件下测得的海洋可控源电磁数据,本文提出了一种基于非结构网格剖分的频率域海洋可控源电磁数据三维正则化反演方法.该方法首先对海洋地电模型以非结构四面体单元进行离散,然后基于矢量有限元方法获得海洋可控源电磁响应和灵敏度信息,最后采用共轭梯度法求解高斯-牛顿反演方程计算模型修正量.为提高反演的稳定性,通过在反演过程中采用对数转换方法实现反演模型参数的上下限约束.本文分别测试了单测线水平海底地形反演算例和面积性测量的起伏海底地形反演算例.反演结果表明,本文提出的频率域海洋可控源电磁三维反演能够准确地恢复高阻储油层的位置和电阻率信息,且计算效率较高,可用于实测海洋电磁资料的处理与解释. 相似文献