共查询到18条相似文献，搜索用时 234 毫秒
1.
初至波走时反演难以得到复杂的地下构造,而早至波波形反演在处理远偏移距地震记录时易产生周期跳跃现象,只能利用近偏移距数据进行反演.但近偏移距的早至波所携带的信息较少,难以实现对地下介质的高精度速度建模.为解决上述问题,并充分利用远偏移距的早至波信息,本文提出基于初至波波场重构的早至波反演方法.该方法将早至波波动方程作为惩罚项,加入到传统的全波形反演的目标函数中.新的反演问题转化为一个选择优化问题,可以用变量映射法对该问题进行求解:首先给定初始速度,在早至波波形空间寻找合适的早至波解;在准确重构早至波场之后,利用早至波场与近地表速度的近线性关系,实现对近地表速度场的高精度建模.该方法不包含观测数据的匹配过程,可以更好的处理远偏移距信息并避免周期跳跃.模型测试验证了本文方法的准确性和有效性. 相似文献
2.
利用天然地震记录直达P波后续20s内的波形信息,研究了地壳、土地慢速度结构反演的方法、唯一性、精度及应注意的问题,并将求解有条件极值的惩罚数法和求解无条件极值的单纯形法引入到波形反演.数值计算证明,波形反演方法对记录的误差有压制作用,在深源、远震、各向同性水平层状近似合理的前提下,可用一个地震事件在射线平面内的二分量记录反演台站区地壳、上地幔P、S波速度结构.在新疆天山地区转换波流动台站观测中,选择出较好的记录进行波形反演,给出了该转换被测深剖面上6个台站下方的地壳、上地慢P、S波速度结构.通过塔里木盆地内一测点上波形反演和地震勘探结果的对比证明,波形反演方法具有满意的精度. 相似文献
3.
关于区域性的地震面波层析反演方法的讨论 总被引：4，自引：0，他引：4
讨论了在地球表面局部区域进行面波层析反演时所用的各种反演方法，指出各种方法用于区域性层析反演时的优缺点，并就由面波资料对区域性的地壳和上地幔结构进行反演的方法──分块反演方法、本征函数展开法、球面雷当变换方法、Ｔａｒａｎｔｏｌａ概率方法和波形反演的方法进行了讨论。指出了各种方法的优缺点及适用范围，并对一些方法进行了数值模拟，给出了计算实例。关键词 相似文献
4.
高精度的成像需要准确的速度场信息,波形反演被认为是目前具有最高分辨率的速度反演方法之一.计算效率是目前全波形反演需要考虑的一个主要问题.为了很好地解决计算效率的问题,本文引入了一种高效的无串扰同时源反演方法,并详细介绍了其原理与计算流程.同时,基于此多震源同时反演方法,本文拓展出实波场反演、虚波场反演及复波场反演的三种反演策略,进而分析复波场的实部和虚部对全波形反演的影响.实验表明,相比于复波场全波形反演,无论是实波场反演还是虚波场反演,反演分辨率有所降低;相比于实波场反演,虚波场反演对初始模型的依赖性较小,目标函数的非线性较弱;最后,通过使用组合反演策略,即初始阶段采用虚波场反演,中后期阶段采用复波场反演,不仅可以降低反演的非线性,而且能够保证高精度建模. 相似文献
5.
为了提高表层速度反演精度,本文提出了一种新的波形反演方法.该方法只利用初至波波形信息以减少波形反演对初始模型的依赖性,降低反演多解性与稳定性.由于只利用初至波波形信息,所以该方法利用高斯束计算格林函数和正演波场,以减少正演计算量.为了避免庞大核函数的存储,该方法基于Born波路径,利用矩阵分解算法实现方向与步长的累加计算.将此基于Born波路径的初至波波形反演方法应用于理论模型实验,并与声波方程全波形反演和初至波射线走时层析方法相对比,发现该方法的反演效果略低于全波形反演方法,但明显优于传统初至波射线走时层析方法,而计算效率却与射线走时层析相当.同时,相对于全波形反演,本文方法对初始模型的依赖性也有所降低. 相似文献
6.
地面资料全波形反演采用低波数回折波敏感核和高波数偏移等时线恢复地下模型的长波长和短波长分量.当回折波的穿透深度有限时,很容易陷入局部极值.波动方程反射波走时反演采用反射波的透射敏感核更新速度模型.当浅层速度存在较大误差时,仍无法得到满意的反演结果.VSP资料中直达波是一种透射波,穿透深度大,可用于井旁背景速度场建模.本文发展了针对VSP观测方式的波动方程走时联合反演(直达波和反射波)方法.采用全局优化参数反演法实现VSP上下行波和纵横波分离,构建了基于下行直达波和上行反射波走时残差的混合目标函数,推导了相应的伴随状态方程和梯度公式,给出了背景速度和反射系数分步反演流程.理论和实际VSP资料应用表明：直达波和反射波走时联合反演可以得到运动学特征精确的速度模型,为后续全波形反演提供可靠的初始模型;联合反演比单独直达波/反射波走时反演的精度高,比直达波和反射波走时级联反演的耗时少. 相似文献
7.
通过对波场的时间二阶积分运算以增强地震数据中的低频成分,提出了一种可有效减小对初始速度模型依赖性的地震数据全波形反演方法—时间二阶积分波场的全波形反演方法.根据散射理论中的散射波场传播方程,推导出时间二阶积分散射波场的传播方程,再利用一阶Born近似对时间二阶积分散射波场传播方程进行线性化.在时间二阶积分散射波场传播方程的基础上,利用散射波场反演地下散射源分布,再利用波场模拟的方法构建地下入射波场,然后根据时间二阶积分散射波场线性传播方程中散射波场与入射波场、速度扰动间的线性关系,应用类似偏移成像的公式得到速度扰动的估计,以此建立时间二阶积分波场的全波形迭代反演方法.最后把时间二阶积分波场的全波形反演结果作为常规全波形反演的初始模型可有效地减小地震波场全波形反演对初始模型的依赖性.应用于Marmousi模型的全频带合成数据和缺失4Hz以下频谱成分的缺低频合成数据验证所提出的全波形反演方法的正确性和有效性,数值试验显示缺失4Hz以下频谱成分数据的反演结果与全频带数据的反演结果没有明显差异. 相似文献
8.
采用弹性波全波形反演方法精确重建深部金属矿多参数模型,建模过程采用基于地震照明的反演策略.首先给出基于照明理论的观测系统可视性定义,利用可视性分析构建新的目标函数,对反演目标可视性较高的炮检对接收到的地震记录在波场匹配时占有更高的权重,确保了参与反演计算中的地震数据的有效性;其次将给定观测系统对地下介质的弹性波场照明强度作为优化因子,根据地震波在波阻抗界面处的能量分配特点,自适应补偿波场能量分布和优化速度梯度,以提高弹性波全波形反演过程的稳定性和反演结果的精度.理论模型和金属矿模型反演试验结果表明,基于可视性分析和能量补偿的反演策略可以使弹性波全波形反演更快地收敛到目标函数的全局极小值,获得适用于金属矿高分辨率地震偏移成像的多参数模型. 相似文献
9.
低频成分缺失和地下速度强烈变化会导致严重的周期跳现象,是地震数据全波形反演的难题.通过对地震数据加时间阻尼和时间积分降主频处理,提出了一种可有效去除周期跳现象的多主频波场时间阻尼全波形反演方法.由浅到深的速度不准确会造成波形走时失配和走时失配的累积.浅部速度的准确反演可有效地减小深部波形走时失配与周期跳现象.对地震数据施加时间阻尼得到时间阻尼数据,利用不同阻尼值的时间阻尼地震数据实现由浅到深的全波形反演.低主频波场的周期跳现象相对高主频波场的要弱.对地震波场进行不同阶的时间积分以得到不同主频的波场,把低主频波场的全波形反演结果作为高主频波场全波形反演的初始模型.应用缺失4 Hz以下频谱成分的二维盐丘模型合成数据验证所提出的全波形反演方法的正确性和有效性,数值试验结果显示多主频波场的时间阻尼全波形反演方法对缺失低频成分地震数据和地下速度强烈变化具有很好的适应性. 相似文献
10.
利用CDSN以及境外的数字地震台网观测的宽频带体波波形资料,采用体波波形反演方法,对中国北部及蒙古地区的上地幔平均P波速度结构以及部分地区的横向不均匀性特征进行了研究．结果表明,该地区上地幔盖层的P波速度较低(约7.8～8.0 km/s),平均盖层厚度约60 km,在410和665 km附近存在速度跳跃分别为0.29和0.55 km/s左右的速度间断面．准噶尔盆地上地幔顶部P波速度约7.7 km/s．上地幔盖层具有较高的速度梯度(平均速度梯度>0.005 5/s)和较大的厚度(90～100 km),在140 km深处P波速度可达8.2 km/s左右．贝加尔湖附近上地幔盖层的平均P波速度介于8.0～8.05 km/s之间,上地幔盖层厚度约30 km． 相似文献
11.
IntroductionHistorically, surface waves provided much of our understanding of the lateral heterogeneity in the crust and the upper mantle, especially in oceans or desert areas. Traditional method utilizes the great circle theorem for surface waves. It assumed that surface wave propagates along the great circle connecting epicenter with station. Firstly, traditional surface wave analysis was carried out by using multi-frequency-filter dispersion analysis to acquire the group or phase velocity d… 相似文献
12.
We develop a two‐dimensional full waveform inversion approach for the simultaneous determination of S‐wave velocity and density models from SH ‐ and Love‐wave data. We illustrate the advantages of the SH/Love full waveform inversion with a simple synthetic example and demonstrate the method's applicability to a near‐surface dataset, recorded in the village ?achtice in Northwestern Slovakia. Goal of the survey was to map remains of historical building foundations in a highly heterogeneous subsurface. The seismic survey comprises two parallel SH‐profiles with maximum offsets of 24 m and covers a frequency range from 5 Hz to 80 Hz with high signal‐to‐noise ratio well suited for full waveform inversion. Using the Wiechert–Herglotz method, we determined a one‐dimensional gradient velocity model as a starting model for full waveform inversion. The two‐dimensional waveform inversion approach uses the global correlation norm as objective function in combination with a sequential inversion of low‐pass filtered field data. This mitigates the non‐linearity of the multi‐parameter inverse problem. Test computations show that the influence of visco‐elastic effects on the waveform inversion result is rather small. Further tests using a mono‐parameter shear modulus inversion reveal that the inversion of the density model has no significant impact on the final data fit. The final full waveform inversion S‐wave velocity and density models show a prominent low‐velocity weathering layer. Below this layer, the subsurface is highly heterogeneous. Minimum anomaly sizes correspond to approximately half of the dominant Love‐wavelength. The results demonstrate the ability of two‐dimensional SH waveform inversion to image shallow small‐scale soil structure. However, they do not show any evidence of foundation walls. 相似文献
13.
波动方程深度偏移是解决复杂地质体成像的关键技术,基于波动方程的速度建模为其提供更为精确的速度模型.频率域波形反演是目前研究最为广泛的波动方程速度建模方法之一,它推动了波形反演在勘探尺度下的应用.本文通过对频率域波形反演的实现,分析对比了其有效执行过程中与频率相关的影响因素.介绍了时间域的多尺度反演方法在频率域的一种实现方式,对比分析了输入数据的频点带宽和应用的子波频带范围不同时对反演结果的影响.本文通过设计的山地地质模型对频率域波形反演进行了测试和对比,得到的结论为频率域波形反演的有效计算提供了依据和参考. 相似文献
14.
Crustal velocity structure of central Gansu Province from regional seismic waveform inversion using firework algorithm 
下载免费PDF全文
下载免费PDF全文 The firework algorithm (FWA) is a novel swarm intelligence-based method recently proposed for the optimization of multi-parameter, nonlinear functions. Numerical waveform inversion experiments using a synthetic model show that the FWA performs well in both solution quality and efficiency. We apply the FWA in this study to crustal velocity structure inversion using regional seismic waveform data of central Gansu on the northeastern margin of the Qinghai-Tibet plateau. Seismograms recorded from the moment magnitude (M W) 5.4 Minxian earthquake enable obtaining an average crustal velocity model for this region. We initially carried out a series of FWA robustness tests in regional waveform inversion at the same earthquake and station positions across the study region, inverting two velocity structure models, with and without a low-velocity crustal layer; the accuracy of our average inversion results and their standard deviations reveal the advantages of the FWA for the inversion of regional seismic waveforms. We applied the FWA across our study area using three component waveform data recorded by nine broadband permanent seismic stations with epicentral distances ranging between 146 and 437 km. These inversion results show that the average thickness of the crust in this region is 46.75 km, while thicknesses of the sedimentary layer, and the upper, middle, and lower crust are 3.15, 15.69, 13.08, and 14.83 km, respectively. Results also show that the P-wave velocities of these layers and the upper mantle are 4.47, 6.07, 6.12, 6.87, and 8.18 km/s, respectively. 相似文献
15.
Guido Baeten Jan Willem de Maag René‐Edouard Plessix Rini Klaassen Tahira Qureshi Maren Kleemeyer Fons ten Kroode Zhang Rujie 《Geophysical Prospecting》2013,61(4):701-711
Velocity model building and impedance inversion generally suffer from a lack of intermediate wavenumber content in seismic data. Intermediate wavenumbers may be retrieved directly from seismic data sets if enough low frequencies are recorded. Over the past years, improvements in acquisition have allowed us to obtain seismic data with a broader frequency spectrum. To illustrate the benefits of broadband acquisition, notably the recording of low frequencies, we discuss the inversion of land seismic data acquired in Inner Mongolia, China. This data set contains frequencies from 1.5–80 Hz. We show that the velocity estimate based on an acoustic full‐waveform inversion approach is superior to one obtained from reflection traveltime inversion because after full‐waveform inversion the background velocity conforms to geology. We also illustrate the added value of low frequencies in an impedance estimate. 相似文献
16.
Elastic waves, such as Rayleigh and mode‐converted waves, together with amplitude versus offset variations, serve as noise in full waveform inversion using the acoustic approximation. Heavy preprocessing must be applied to remove elastic effects to invert land or marine data using the acoustic inversion method in the time or frequency domains. Full waveform inversion using the elastic wave equation should be one alternative; however, multi‐parameter inversion is expensive and sensitive to the starting velocity model. We implement full acoustic waveform inversion of synthetic land and marine data in the Laplace domain with minimum preprocessing (i.e., muting) to remove elastic effects. The damping in the Laplace transform can be thought of as an automatic time windowing. Numerical examples show that Laplace‐domain acoustic inversion can yield correct smooth velocity models even with the noise originating from elastic waves. This offers the opportunity to develop an accurate smooth starting model for subsequent inversion in the frequency domain. 相似文献
17.
Liangguo Dong Zhongyi Fan Hongzhi Wang Benxin Chi Yuzhu Liu 《Geophysical Prospecting》2018,66(8):1503-1520
Reflection full waveform inversion can update subsurface velocity structure of the deeper part, but tends to get stuck in the local minima associated with the waveform misfit function. These local minima cause cycle skipping if the initial background velocity model is far from the true model. Since conventional reflection full waveform inversion using two‐way wave equation in time domain is computationally expensive and consumes a large amount of memory, we implement a correlation‐based reflection waveform inversion using one‐way wave equations to retrieve the background velocity. In this method, one‐way wave equations are used for the seismic wave forward modelling, migration/de‐migration and the gradient computation of objective function in frequency domain. Compared with the method using two‐way wave equation, the proposed method benefits from the lower computational cost of one‐way wave equations without significant accuracy reduction in the cases without steep dips. It also largely reduces the memory requirement by an order of magnitude than implementation using two‐way wave equation both for two‐ and three‐dimensional situations. Through numerical analysis, we also find that one‐way wave equations can better construct the low wavenumber reflection wavepath without producing high‐amplitude short‐wavelength components near the image points in the reflection full waveform inversion gradient. Synthetic test and real data application show that the proposed method efficiently updates the background velocity model. 相似文献
18.
Offset‐variable density improves acoustic full‐waveform inversion: a shallow marine case study 下载免费PDF全文
下载免费PDF全文 Akela Silverton Michael Warner Joanna Morgan Adrian Umpleby 《Geophysical Prospecting》2016,64(5):1201-1214
We have previously applied three‐dimensional acoustic, anisotropic, full‐waveform inversion to a shallow‐water, wide‐angle, ocean‐bottom‐cable dataset to obtain a high‐resolution velocity model. This velocity model produced an improved match between synthetic and field data, better flattening of common‐image gathers, a closer fit to well logs, and an improvement in the pre‐stack depth‐migrated image. Nevertheless, close examination reveals that there is a systematic mismatch between the observed and predicted data from this full‐waveform inversion model, with the predicted data being consistently delayed in time. We demonstrate that this mismatch cannot be produced by systematic errors in the starting model, by errors in the assumed source wavelet, by incomplete convergence, or by the use of an insufficiently fine finite‐difference mesh. Throughout these tests, the mismatch is remarkably robust with the significant exception that we do not see an analogous mismatch when inverting synthetic acoustic data. We suspect therefore that the mismatch arises because of inadequacies in the physics that are used during inversion. For ocean‐bottom‐cable data in shallow water at low frequency, apparent observed arrival times, in wide‐angle turning‐ray data, result from the characteristics of the detailed interference pattern between primary refractions, surface ghosts, and a large suite of wide‐angle multiple reflected and/or multiple refracted arrivals. In these circumstances, the dynamics of individual arrivals can strongly influence the apparent arrival times of the resultant compound waveforms. In acoustic full‐waveform inversion, we do not normally know the density of the seabed, and we do not properly account for finite shear velocity, finite attenuation, and fine‐scale anisotropy variation, all of which can influence the relative amplitudes of different interfering arrivals, which in their turn influence the apparent kinematics. Here, we demonstrate that the introduction of a non‐physical offset‐variable water density during acoustic full‐waveform inversion of this ocean‐bottom‐cable field dataset can compensate efficiently and heuristically for these inaccuracies. This approach improves the travel‐time match and consequently increases both the accuracy and resolution of the final velocity model that is obtained using purely acoustic full‐waveform inversion at minimal additional cost. 相似文献