基于深度学习的重力异常与重力梯度异常联合反演   总被引：1，自引：0，他引：1  

张志厚  廖晓龙  曹云勇  侯振隆  范祥泰  徐正宣  路润琪  冯涛  姚禹  石泽玉 《地球物理学报》2021,64(4):1435-1452

高效高精度的反演算法在重力大数据时代背景下显得尤为重要,受深度学习卓越的非线性映射能力的启发,本文提出了一种基于深度学习的重力异常及重力梯度异常的联合反演方法.文中首先提出了一种基于网格点几何格架的重力异常及重力梯度异常的空间域快速正演算法,这为本文深度学习反演算法的实现奠定了基础;其次对大量的不同密度模型进行正演计算获得样本数据集;然后设计了一种端到端的深度学习网络结构(GraInvNet),再利用样本数据对该网络结构进行训练;最后进行反演预测.组合模型试验表明,多维度数据联合反演相比单一分量反演其结果更“聚焦”,且与模型边界高度吻合,并且对于复杂模型的姿态与物性预测具有极为显著的优势,以及对于含噪声数据的反演,其质量也不会降低;Vinton岩丘实测重力数据也验证了文中方法的有效性;从而证明了深度学习在重力数据的高效高精度反演方面具有的巨大潜力.  相似文献   

Determining the focal mechanisms and depths of relatively small earthquakes using a few stations by full‐waveform modelling     

Hussam Busfar  M. Nafi Toksöz 《Geophysical Prospecting》2013,61(4):712-724

Determining the focal mechanism of earthquakes helps us to better define faults and understand the stress regime. This technique can be helpful in the oil and gas industry where it can be applied to microseismic events. The objective of this paper is to find double couple focal mechanisms, excluding scalar seismic moments, and the depths of small earthquakes using data from relatively few local stations. This objective is met by generating three‐component synthetic seismograms to match the observed normalized velocity seismograms. We first calculate Green's functions given an initial estimate of the earthquake's hypocentre, the locations of the seismic recording stations and a 1D velocity model of the region for a series of depths. Then, we calculate the moment tensor for different combinations of strikes, dips and rakes for each depth. These moment tensors are combined with the Green's functions and then convolved with a source time function to produce synthetic seismograms. We use a grid search to find the synthetic seismogram with the largest objective function that best fits all three components of the observed velocity seismogram. These parameters define the focal mechanism solution of an earthquake. We tested the method using three earthquakes in Southern California with moment magnitudes of 5.0, 5.1 and 4.4 using the frequency range 0.1–2.0 Hz. The source mechanisms of the events were determined independently using data from a multitude of stations. Our results obtained, from as few as three stations, generally match those obtained by the Southern California Earthquake Data Center. The main advantage of this method is that we use relatively high‐frequency full‐waveforms, including those from short‐period instruments, which makes it possible to find the focal mechanism and depth of earthquakes using as few as three stations when the velocity structure is known.  相似文献   

P‐wave velocity distribution in basalt flows of the Enni Formation in the Faroe Islands from refraction seismic analysis     

Uni K. Petersen  R. James Brown  Morten S. Andersen 《Geophysical Prospecting》2013,61(1):168-186

The main objective of this work is to establish the applicability of shallow surface‐seismic traveltime tomography in basalt‐covered areas. A densely sampled ～1300‐m long surface seismic profile, acquired as part of the SeiFaBa project in 2003 ( Japsen et al. 2006 ) at Glyvursnes in the Faroe Islands, served as the basis to evaluate the performance of the tomographic method in basalt‐covered areas. The profile is centred at a ～700‐m deep well. V_P, V_S and density logs, a zero‐offset VSP, downhole‐geophone recordings and geological mapping in the area provided good means of control. The inversion was performed with facilities of the Wide Angle Reflection/Refraction Profiling program package ( Ditmar et al. 1999 ). We tested many inversion sequences while varying the inversion parameters. Modelled traveltimes were verified by full‐waveform modelling. Typically an inversion sequence consists in several iterations that proceed until a satisfactory solution is reached. However, in the present case with high velocity contrasts in the subsurface we obtained the best result with two iterations: first obtaining a smooth starting model with small traveltime residuals by inverting with a high smoothing constraint and then inverting with the lowest possible smoothing constraint to allow the inversion to have the full benefit of the traveltime residuals. The tomogram gives usable velocity information for the near‐surface geology in the area but fails to reproduce the expected velocity distribution of the layered basalt flows. Based on the analysis of the tomogram and geological mapping in the area, a model was defined that correctly models first arrivals from both surface seismic data and downhole‐geophone data.  相似文献   

利用余弦变换的重力梯度Parker正演方法及其应用     

刘繁明  张迎发荆心 李艳 《测绘学报》2013,42(2):177-190

重力测量数据存在地形数据产生的高频分量的影响,高精度地形数据正演重力梯度也能较好地反映重力局部高频特征。为获得高精度重力梯度数据,实现基准梯度数据库精确快速构建,研究了利用数字高程模型正演重力梯度的频率域快速计算方法,推导出基于余弦变换的Parker正演重力梯度理论公式。数值实验结果表明,余弦变换频率域正演方法平均绝对误差可达到0.5E左右精度要求,与傅里叶变换正演方法相比误差可减小3dB左右,与棱柱法等空间域正演方法相比,该方法计算规模小,速度优势明显。  相似文献   

Blocky regularization schemes for Full‐Waveform Inversion★     

Antoine Guitton 《Geophysical Prospecting》2012,60(5):870-884

With ill‐posed inverse problems such as Full‐Waveform Inversion, regularization schemes are needed to constrain the solution. Whereas many regularization schemes end up smoothing the model, an undesirable effect with FWI where high‐resolution maps are sought, blocky regularization does not: it identifies and preserves strong velocity contrasts leading to step‐like functions. These models might be needed for imaging with wave‐equation based techniques such as Reverse Time Migration or for reservoir characterization. Enforcing blockiness in the model space amounts to enforcing a sparse representation of discontinuities in the model. Sparseness can be obtained using the ?¹ norm or Cauchy function which are related to long‐tailed probability density functions. Detecting these discontinuities with vertical and horizontal gradient operators helps constraining the model in both directions. Blocky regularization can also help recovering higher wavenumbers that the data used for inversion would allow, thus helping controlling the cost of FWI. While the Cauchy function yields blockier models, both ?¹ and Cauchy attenuate illumination and inversion artifacts.  相似文献   

基于Born敏感核函数的速度、密度双参数全波形反演          下载免费PDF全文

杨积忠  刘玉柱  董良国 《地球物理学报》2016,59(3):1082-1094

速度、密度之间的相互耦合使得密度在多参数全波形反演中较难获得.本文将截断高斯-牛顿法用于声介质速度、密度双参数全波形反演,通过考虑近似Hessian矩阵中反映速度、密度相互作用的非主对角块元素,有效解决了多参数全波形反演中速度、密度之间的耦合问题,在不采用反演策略的情况下,仍能够获得精度较高的速度、密度反演结果.常规的截断牛顿类全波形反演通常利用一阶伴随状态法求取目标函数对模型参数的梯度,利用二阶伴随状态法或有限差分法求解Hessian-向量乘,在每一步内循环迭代过程中需要额外求解两次正演问题,计算量较大.本文基于Born近似,将梯度计算中的核函数-向量乘表示为具有明确物理意义的向量-标量乘的累加运算,同时将Hessian-向量乘转化为两次核函数-向量乘,无需额外求解正演问题,有效降低了计算量.数值实验证明了本文提出的方法的有效性.  相似文献   

Modelling of rutting of two flexible pavements with the shakedown theory and the finite element method   总被引：1，自引：0，他引：1  

Cyrille Chazallon  Georg Koval  Pierre Hornych  Fatima Allou  Saida Mouhoubi 《Computers and Geotechnics》2009

This paper presents a finite element program, for the modelling of rutting of flexible pavements. In its present version, the program incorporates a permanent deformation model for unbound granular materials based on the concept of the shakedown theory developed by Zarka for metallic structures under cyclic loadings and has been used to estimate the permanent deformations of unbound granular materials (UGM) subjected to traffic loading. The calculation is performed in two steps: the first step consists in modelling the resilient behaviour of the pavement in 3D, using non-linear elastic models, to determine the stress field in the pavement. Then stress paths are derived and used to calculate the permanent deformations and the displacements, using a Drucker–Prager yield surface. An application to the prediction of the permanent deformations of experimental pavements with an unbound granular base, tested on the LCPC pavement testing facility is presented.  相似文献   

Pore structure characteristics of tight sandstones in the northern Songliao Basin,China     

《Marine and Petroleum Geology》2017

Understanding the pore structure characteristics of tight gas sandstones is the primary purpose of reservoir evaluation and efforts to characterize tight gas transport and storage mechanisms and their controls. Due to the various pore types and multi-scale pore sizes in tight reservoirs, it is essential to combine several techniques to characterize pore structure. Scanning electron microscopy (SEM), nitrogen gas adsorption (N₂GA), mercury intrusion porosimetry (MIP) and nuclear magnetic resonance (NMR) were conducted on tight sandstones from the Lower Cretaceous Shahezi Formation in the northern Songliao Basin to investigate pore structure characteristics systematically (e.g., type and size distribution of pores) and to establish how significant porosity and permeability are for different pore types. The studied tight sandstones are composed of intergranular pores, dissolution pores and intercrystalline pores. The integration of N₂GA and NMR can be used as an efficient method to uncover full pore size distribution (PSD) of tight sandstones, with pore sizes ranging from 2 nm to dozens of microns. The full PSDs indicate that the pore sizes of tight sandstones are primarily distributed within 1.0 μm. With an increase in porosity and permeability, pores with larger sizes contribute more to porosity. Intercrystalline pores and intergranular/dissolution pores can be clearly distinguished on the basis of mercury intrusion and surface fractal. The relative contribution of intercrystalline pores to porosity ranges from 58.43% to 91.74% with an average of 79.74%. The intercrystalline pores are the primary contributor to pore space, whereas intergranular/dissolution pores make a considerably greater contribution to permeability. A specific quantity of intergranular/dissolution pores is the key to producing high porosity and permeability in tight sandstone reservoirs. The new two permeability estimation models show an applicable estimation of permeability with R² values of 0.955 and 0.962 for models using D_max (pore diameter corresponding to displacement pressure) and D_f (pore diameter at inflection point), respectively. These results indicate that both D_max and D_f are key factors in determining permeability.  相似文献   

基于全波形反演的探地雷达数据逆时偏移成像   总被引：1，自引：1，他引：0       下载免费PDF全文

雷林林  刘四新  傅磊  吴俊军 《地球物理学报》2015,58(9):3346-3355

逆时偏移成像(RTM)常用来处理复杂速度模型,包括陡倾角及横向速度变化剧烈的模型.与常规偏移成像方法(如Kirchhoff偏移)相比,逆时偏移成像能提供更好的偏移成像结果,近些年逆时偏移成像越来越广泛地应用到勘探地震中,它逐渐成为石油地震勘探中的一种行业标准.电磁波和弹性波在动力学和运动学上存在相似性,故本文开发了基于麦克斯韦方程组的电磁波逆时偏移成像算法,并将其应用到探地雷达数据处理中.时间域有限差分(FDTD)用于模拟电磁波正向和逆向传播过程,互相关成像条件用于获得最终偏移结果.逆时偏移成像算法中,偏移成像结果受初始模型影响较大,而其中决定电磁波传播速度的介电常数的影响尤为重要.本文基于时间域全波形反演(FWI)算法反演获得了更为精确的地下介电常数模型,并将其反演结果作为逆时偏移成像的初始介电常数模型.为了验证此算法的有效性,首先构建了一个复杂地质结构模型,合成了共偏移距及共炮点探地雷达数据,分别应用常规Kirchhoff偏移算法及逆时偏移成像算法进行偏移处理,成像结果显示由逆时偏移成像算法得到的偏移结果与实际模型具有较高的一致性;此外本文在室内沙槽中进行了相关的物理模拟实验,采集了共偏移距及共炮点探地雷达数据,分别应用Kirchhoff和叠前逆时偏移成像算法进行处理,结果表明叠前逆时偏移成像在实际应用中能获得更好的成像效果.  相似文献   

不依赖子波、基于包络的FWI初始模型建立方法研究   总被引：3，自引：3，他引：0       下载免费PDF全文

敖瑞德  董良国  迟本鑫 《地球物理学报》2015,58(6):1998-2010

地震全波形反演(FWI)从理论走向实际面临着诸多难题,其中之一就是需要一个较高精度的初始模型,另一个难题就是需要一个较为精确的震源子波,初始模型和震源子波的准确程度严重影响着全波形反演的最终结果.为此,本文提出了不依赖子波、基于包络的FWI初始模型建立的方法,建立了相应的目标函数,推导出了反演的梯度,给出了伴随震源的表达式,理论上分析了不依赖子波FWI的可行性.在数值试验中,讨论了参考道的选取方式,通过分析归一化目标函数收敛速率,认为近偏移距参考道优于远偏移距参考道,在地震数据含干扰噪音时,平均道作为参考道要优于最小偏移距参考道.通过包络、包络对数、包络平方三种目标函数反演结果的比较,发现包络对数目标函数对深层的反演效果最好.通过不同子波的试验进一步验证了本方法的正确性.  相似文献