共查询到20条相似文献，搜索用时 547 毫秒
1.
面波成像是研究地壳上地幔横波速度结构的一种重要方法．通常,面波相速度或群速度成像假设面波沿大圆路径传播．但是,在地下介质速度结构变化较大时,面波会偏离大圆路径传播,从而导致基于大圆路径假设下的面波成像结果存在一定的误差. 我们采用基于射线追踪的面波成像方法,研究了面波的偏离大圆路径传播对四川西部地区面波相速度成像结果的影响．使用快速行进法(fast marching method)进行面波传播路径的射线追踪,采用子空间反演法(subspace inversion)进行迭代反演,对理论模型合成数据和川西台阵的短周期背景噪声相速度频散数据进行成像分析,并与使用大圆路径传播的成像结果进行对比．对理论模型的测试结果表明,当速度结构变化较大时,基于偏离大圆路径传播的面波成像能够更好地恢复模型异常．对川西台阵的真实数据反演结果显示:在短周期为6 s时,基于偏离大圆路径传播的反演方法较基于大圆路径传播的反演方法所获得的相速度异常的幅度更大些,在四川盆地区域两者的差异接近0.2 km/s;在周期为10 s时,两种反演方法的差异显著减小,基本都在0.1 km/s以内．这主要是因为6 s周期的面波相速度对复杂的上地壳浅层结构更为敏感,从而使得面波的偏离大圆路径传播效应对反演结果的影响更为显著．本文结果表明,当某一周期不同路径的面波相速度测量值变化较大,例如相对于平均相速度的异常超过10％时,则需考虑采用基于偏离大圆路径传播的面波成像方法,否则速度异常较大区域的反演结果可能会造成较大的偏差． 相似文献
2.
A new algorithm for tomographic inversion of traveltimes of reflected and refracted seismic waves is developed. The inversion gives interface configurations and velocity distributions in layers. The important features of the algorithm are: (a) the inclusion of shot time delays in the list of unknown parameters; (b) the regularization is applied in such a way that the most probable model is characterized by the similarity of neighbouring interfaces. As the problem under consideration is non-linear, several iterations are necessary in order to obtain the final model. In the case of a very inexact initial model, a 'layer-by-layer' inversion strategy is recommended as a first inversion step. The inversion program is supplied with a user interface, thanks to which one can: (a) pick interactively and identify seismic traveltimes; (b) build and edit depth/velocity models; and (c) display calculated traveltime curves and compare them with picked traveltimes as well as with the original seismic sections. The efficiency of the inversion software developed is illustrated by a numerical example and a field example in which shallow seismic data are considered. Application to wide-aperture reflection/refraction profiling (WARRP) data is also possible. 相似文献
3.
有别于传统基于梯度信息的反演方法在正则化约束中用总梯度逼近海塞逆矩阵的技术,本文将正则化约束问题的数据拟合项和模型光滑项分开考虑,只利用数据拟合函数的梯度信息对数据拟合项的海塞矩阵进行逼近,通过求解类高斯牛顿下降方向方程得到不依赖前几次迭代正则化因子的更精确下降方向,在求解当前迭代下降方向的过程中,通过保证右端项中两个向量的二范数在同一数量级的原则,实现了正则化因子的自动更新.对理论模型的试算表明这种自适应正则化反演方案可以在拟牛顿反演框架下基本达到OCCAM的算法稳定性,反演结果对初始模型依赖性较小,同时又无需在一次迭代中多次搜索最佳正则化因子.本文还基于此算法讨论了大地电磁各参数对于反演结果的影响,由于本文的反演结果能得到充分的正则化约束,因而在此框架下讨论阻抗和倾子在反演中的作用相对更为客观. 相似文献
4.
深度速度模型的构建仍然是地震成像中的巨大挑战,获得一个精确的深度速度模型和减少深度成像项目周期都是至关重要的。常规层析反演速度建模每次迭代相当于一次线性反演,且需要重新的拾取工作,导致非常耗时,效率低下。本文提出非线性层析反演速度建模技术来建立速度模型。拾取共成像点道集的RMO量,转换到叠前域,作为运动学不变量,通过层析反演迭代进行模型更新。用一个多次的线性反演来逼近一个非线性的物理过程,避免重复的拾取工作,大大提高项目的运转效率。通过实例证明该方法的有效性。 相似文献
5.
6.
T. B. Yanovskaya V. M. Kozhevnikov O. A. Solovei K. R. Akchurin 《Izvestiya Physics of the Solid Earth》2008,44(8):622-630
Dispersion curves of phase velocities of Rayleigh waves are determined by the method of frequency-time analysis in a range of periods of 10–200 s from data of 43 interstation traces in Central Asia. Because the joint use of phase and group velocities significantly decreases the uncertainty in the determination of S wave velocity structures, the same traces were used for calculating group velocities from tomographic reconstructions obtained in [Yanovskaya and Kozhevnikov, 2003, 2006] and determining average velocity structures along these traces. The velocity structures were calculated by the Monte Carlo and linear inversion methods, which gave consistent results. Using velocity values obtained at fixed depths by the 2-D tomography method, lateral variations in velocities at these depths were estimated, which allowed us to construct smoothed vertical velocity structures at some points in the region. The resulting structures were used as initial approximations for constructing local velocity structures solely from previously obtained local dispersion curves of group velocities in the area (32°–56°N, 80°–120°E). Based on these structures, we mapped the lateral distribution of velocity variations at upper mantle depths of 75–400 km and along three vertical profiles. The inferred velocity variations are in good agreement with data on the tectonics of the region. 相似文献
7.
变阻尼约束层析成像及其在VSP资料中的应用(英文) 总被引：1，自引：0，他引：1
初至波走时层析成像已经取得了广泛的应用,然而,由于观测系统的限制,射线在模型中分布不均匀,导致层析结果的分辨能力不足。变阻尼约束方法应用不均匀的先验信息来匹配不均匀的数据分布,可以减小速度模型校正量与射线覆盖程度的相关性。本文将变阻尼约束方法应用于初至波旅行时层析成像中,并将平滑约束方法加入正则化方程组中来避免单独使用变阻尼约束带来的不稳定性,利用阿尔法滤波器对反演中间迭代结果进行平滑和去噪,采用LSQR算法求解线性方程组来提高收敛速度和压制误差传递。本文应用上述层析成像算法对VSP观测系统进行速度反演,分别应用于检测板速度模型数据和实际VSP资料速度反演中,结果表咀变阻尼约束层析成像可以改善射线不均匀覆盖带来的影响,从而提高速度反演结果的质量;VSP资料检波点附近的速度反演结果可靠性高。 相似文献
8.
We develop an inversion procedure using the total variation (TV) regularization method as a stabilizing function to invert surface gravity data to retrieve 3-D density models of geologic structures with sharp boundaries. The developed inversion procedure combines several effective algorithms to solve the TV regularized problem. First, a matrix form of the gradient vector is designed using the Kronecker product to numerically approximate the 3-D TV function. The piecewise polynomial truncated singular value decomposition (PP-TSVD) algorithm is then used to solve the TV regularized inverse problem. To obtain a density model with depth resolution, we use a sensitivity-based depth weighting function. Finally, we apply the Genetic Algorithm (GA) to select the best combination of the PP-TSVD algorithm and the depth weighting function parameters. 3-D simulations conducted with synthetic data show that this approach produces sub-surface images in which the structures are well separated in terms of sharp boundaries, without the need of a priori detailed density model. The method applied to a real dataset from a micro-gravimetry survey of Gotvand Dam, southwestern Iran, clearly delineates subsurface cavities starting from a depth of 40 m within the area of the dam reservoir. 相似文献
9.
The Atmospheric Ultraviolet Radiance Integrated Code (AURIC) is a software package developed by Computational Physics, Inc. (CPI) under the sponsorship of the Air Force Phillips Laboratory/Geophysics Directorate (PL/GP) (currently the Air Force Research Laboratory) for middle and upper atmospheric radiance modeling from the far ultraviolet to the near infrared. The AURIC has been considered as a general model for the radiative transfer simulation of airglow. Based on the theory of MODerate resolution atmospheric TRANsmission (MODTRAN), the AURIC extends calculation to altitudes above 100 km and the wavelength down to 80 nm. A package of AURIC v1.2 was released in 2002, which can be used for single-point simulation from 1947 to 1999. It means that the model is not suitable for atmospheric simulation of large datasets or for atmospheric parameters retrieval from amount of satellite measurements. In this paper, AURIC v1.2 is upgraded to AURIC-2012 based on MATLAB with improvements for modules of the Geomagnetic Parameter (GEOPARM), Atmosphere neutral composition (ATOMS), and Ionospheric electron density (IONOS). The improved AURIC can be used for global automatic airglow simulation and also for automatic retrieval of atmospheric compositions from satellite global observations, such as O/N2 and electron density etc. Besides, the model supplies possibilities for further improvement of airglow radiative mechanism and for substitution of other modules. Based on the AURIC-2012, Limb Column Emission Intensity (L-CEI) and Volume Emission Rate (VER) are calculated. For validation, the results were compared with measurements of the Global Ultraviolet Imager (GUVI) and TIMED Doppler Interferometer (TIDI), respectively. The averaged relative errors of L-CEI and VER at peak altitude are both within 20%. Finally, L-CEI varying with latitude, altitude, solar activity, and geomagnetic activity is simulated, and the distribution characteristics of the simulation and their influencing factors are analyzed subsequently. 相似文献
10.
Electrical imaging provides important subsurface information for the construction of hypervelocity impact models. We here provide an overview and evaluation of the current electrical imaging methods used in impact cratering studies. Although apparent resistivity models are commonly used in the geoelectrical imaging of impact structures, the reliability of these models has not hitherto been determined. In order to assess these imaging approaches in impact cratering, we investigate for the first time the discrepancies between the apparent resistivity and true resistivity models of an impact structure. To this end, we present (1) a new true resistivity model of the Araguainha impact structure in central Brazil by applying L2-norm inversion to previously published data, (2) apparent resistivity model of the impact structure, and (3) models obtained from different stages of the iterative tomographic inversions. Our results show that changes in vertical resistivity gradient are significantly better defined in the true resistivity models than in the apparent resistivity model. On the basis of these results, we outline a new approach that true resistivity models can be effectively assessed by applying both L1- and L2-norm inversion schemes together with the monitoring of intermediate models from iterative inversion. The results of our study highlight the importance of tomographic inversion of resistivity data in impact cratering studies, and they provide a data modeling framework and foundation for cost-effective subsurface imaging of impact structures in the future. 相似文献
11.
地震走时层析成像是反演地层各向异性参数分布的有效方法,但是关于地震各向异性介质走时层析成像的研究并不多,其技术远远没有达到成熟的阶段.在野外数据采集时,地表反射波观测方式相对井间和垂直地震剖面观测方式的成本更低,利用qP反射波走时反演各向异性参数具有更加广泛的实用价值.本文实现的TI介质地震走时层析成像方法结合了TI介质反射波射线追踪算法、走时扰动方程和非线性共轭梯度算法,它可以对任意强度的TI介质模型进行反演,文中尝试利用qP反射波走时重建TI介质模型的参数图像.利用qP反射波对层状介质模型和块状异常体模型进行走时反演,由于qP波相速度对弹性模量参数和Thomsen参数的偏微分不同,所以可以分别反演弹性模量参数和Thomsen参数.数值模拟结果表明:利用qP反射波可以反演出TI介质模型的弹性模量参数与Thomsen参数,不同模型的走时迭代反演达到了较好的收敛效果,与各向同性介质走时反演结果相比较,各向异性介质走时反演结果具有较好的识别能力. 相似文献
12.
本文提出了一种初至纵波(P波)与瑞雷面波的交叉梯度联合反演策略.通过对初至P波进行全波形反演可以获得近地表P波速度结构;通过对仅含瑞雷面波信息的地震数据转换到频率-波数域进行加窗振幅波形反演(Windowed-Amplitude Waveform Inversion,w-AWI)可获得近地表横波(S波)速度结构.在二者反演的目标函数中均加入P波速度和S波速度的交叉梯度作为正则化约束项,使得在反演过程中P波速度和S波速度相互制约,相互约束,从而实现对地震初至P波与瑞雷面波的联合反演.数值模拟结果表明交叉梯度联合反演可以提高S波速度反演分辨率,而P波速度反演结果并没有得到提高.实际资料的反演结果表明,交叉梯度联合反演能够获得更加可信的近地表速度结构. 相似文献
13.
Raffaella De Matteis Annalisa Romeo Giuseppe Pasquale Giovanni Iannaccone Aldo Zollo 《Studia Geophysica et Geodaetica》2010,54(3):367-387
A new dataset of first P-wave arrival times is used to derive the 3D tomographic model of the Campania-Lucania region in the
southern Apennines (Italy). We address the issue related to the non-uniqueness of the tomographic inversion solution through
massive numerical experimentation based on the global exploration of the model parameter space starting from a large variety
of physically plausible initial models. The average of all the realizations is adopted as the best-fit solution and the uncertainty
of the model parameters is studied using a statistical approach based on a Monte Carlo-type analysis. How the uncertainty
in the initial model, earthquake locations, and data influences the inversion result is studied by considering separately
the individual effects. Checkerboard tests are performed to estimate the resolving power of the dataset. Re-located seismicity
in a reliable new 3D tomographic model allows us to correlate the earthquake distribution with the main seismogenic structures
present in the area. 相似文献
14.
在波数域中用重力反演莫霍面深度时通常假定壳幔密度差为一常数,但这只是一种近似的密度模型,本文采用了密度随深度呈指数变化的变密度模型来反演莫霍面深度,给出了利用指数密度模型在波数域中计算重力异常的正演公式及界面深度的反演公式.利用指数密度模型及重力资料反演了青藏高原莫霍面的深度,分析了莫霍面的特征.结果表明,青藏高原莫霍面呈现出边缘浅、中部深的特点,边缘变化快、梯度大,中间变化梯度趋缓.中心地带的羌塘地体莫霍面深度达74 km,向四周慢慢变浅至67 km左右,边缘地区突然变浅至50km左右.通过常密度模型、变密度模型及地震反演得到的莫霍面的比较,证实变密度模型更适合于莫霍面结构的反演. 相似文献
15.
16.
17.
利用走时成像中的正则化反演方法对天津—北京的高分辨地震折射剖面获得的初至波资料进行成像.分析了不同的初始模型和参数变化对反演结果的影响,并采用检测板方法对最终模型的分辨进行评估,得到了天津—北京剖面的上地壳精细速度结构,给出了该地区的隐伏断层.结果显示,最终模型分辨达到5 km×1 km,部分地区更加细致;研究区域坳陷区与隆起区相间,大坳陷内部有小坳陷;结晶基底的深度变化约2～8 km,其起伏形态符合地质构造情况;断裂活动和第三系以来的壳内运动对速度结构影响较大. 相似文献
18.
无线电波透视法是常用的工作面地质构造探测方法之一,目前普遍使用的SIRT方法层析分辨率不高。本文采用约束正则化方法,推导Tikhonov正则化和全变差正则化的最小化问题表达式,讨论影响层析结果的主要因素,对典型理论模型进行了层析成像实验。结果表明:正则化方法具有比SIRT方法更好的分辨率;射线条数越多、噪声水平越低,层析分辨率越高;Tikhonov正则化在正则参数增大时层析结果更光滑,减小时则更贴近异常,全变差正则化与其相反。最后对实际坑透数据进行层析,识别出的异常构造基本吻合已知疑似构造位置,从而说明正则化方法在无线电波透视应用中的可行性。 相似文献
19.
针对地表剧烈起伏,速度纵、横向变化大的复杂区,层析静校正较以往的折射波静校正方法有明显的优势,但是受初至时间拾取精度、炮检距的选择、近地表模型约束等问题的影响,层析反演的精度还不能满足低幅度构造预测的需求,为此提出微测井约束分步层析的静校正方法,即将浅层速度模型与最终近地表模型分步进行层析反演,并应用微测井信息约束浅层速度模型层析反演,有效地提高了近地表速度模型反演的精度,通过在古峰庄地区的应用较好地解决了复杂地表条件下的静校正问题,低幅度构造预测精度得到提高. 相似文献
20.
《Journal of Atmospheric and Solar》2002,64(3):367-376
Observations of winds in mesospheric airglow layers have been made at Mt. John (44°S,170°E), New Zealand for some years. We present a modelling study of airglow emissions which shows that the properties of wind detection based on airglow emission means that high-frequency gravity waves are effectively filtered from the wind spectrum observed. This filtering means that any waves with periods of the order of hours should be detectable in the record (as they will not be hidden in the noise of the higher-frequency waves ubiquitous at these heights). One example of such a wave is shown. As part of the analysis, we show that because the airglow layers differ in width, some waves might be observed in only one airglow layer, even when present in both. 相似文献