共查询到20条相似文献，搜索用时 31 毫秒
1.
为了有效解决目前大地电磁和地震走时资料单方法反演结果一致性不好的问题,同时克服基于岩石不同物性参数间关系耦合约束联合反演的局限性,本文研究了基于交叉梯度耦合约束的大地电磁与地震走时资料的三维联合反演算法.以较为成熟的天然地震走时资料三维正反演和大地电磁三维正反演算法为基础,实现了具有共同的反演网格,以交叉梯度结构耦合约束,并能同时获得电阻率和速度模型的三维联合反演算法.分别利用单棱柱体模型和双棱柱体模型合成数据进行了联合反演试算.结果表明:无论是单棱柱体模型还是双棱柱体模型,联合反演结果比单独反演对异常体的空间形态都有更好的恢复,其中单棱柱体模型反演的异常体电阻率更接近于真实电阻率,双棱柱体模型的联合反演结果不仅消除了围岩的部分电阻率假异常,而且增强了对异常体深部速度结构特征的恢复程度.联合反演还能同时改善电阻率和速度反向变化异常体的单独反演结果,进一步证明交叉梯度耦合不依赖于岩石物性关系,而强调地下结构的相似性,具有更普遍的适用性. 相似文献
2.
岩石物性约束是地球物理联合反演重要的实现方式.在以往联合反演研究中,这种约束通常以点对点的映射方式出现,要求经验性的岩石物性关联是相对精确的,无形中提高了岩石物性约束应用于联合反演的条件.基于此,本文结合全局寻优的模拟退火算法,提出了宽范围岩石物性约束技术.该技术将岩石物性约束与模拟退火模型扰动、模型接受准则相结合,实现了空间对空间的映射,将电阻率和速度耦合在一起.文章以MT和地震正则化同步联合反演为例,检验了不同的先验信息条件下的宽范围岩石物性约束技术的适用性.模型试验表明:新技术容易实现,容错性高,可降低"不精确"岩石物性关联信息引入带来的风险,提高岩石物性约束应用的灵活性. 相似文献
3.
The seismic reflection method provides high-resolution data that are especially useful for discovering mineral deposits under deep cover. A hindrance to the wider adoption of the seismic reflection method in mineral exploration is that the data are often interpreted differently and independently of other geophysical data unless common earth models are used to link the methods during geological interpretation. Model-based inversion of post-stack seismic data allows rock units with common petrophysical properties to be identified and permits increased bandwidth to enhance the spatial resolution of the acoustic-impedance model. However, as seismic reflection data are naturally bandlimited, any inversion scheme depends upon an initial model, and must deal with non-unique solutions for the inversion. Both issues can be largely overcome by using constraints and integrating prior information. We exploit the abilities of fuzzy c-means clustering to constrain and to include prior information in the inversion. The use of a clustering constraint for petrophysical values pushes the inversion process to select models that are primarily composed of several discrete rock units and the fuzzy c-means algorithm allows some properties to overlap by varying degrees. Imposing the fuzzy clustering techniques in the inversion process allows solutions that are similar to the natural geologic patterns that often have a few rock units represented by distinct combinations of petrophysical characteristics. Our tests on synthetic models, with clear and distinct boundaries, show that our methodology effectively recovers the true model. Accurate model recovery can be obtained even when the data are highly contaminated by random noise, where the initial model is homogeneous, or there is minimal prior petrophysical information available. We demonstrate the abilities of fuzzy c-means clustering to constrain and to include prior information in the acoustic-impedance inversion of a challenging magnetotelluric/seismic data set from the Carlin Gold District, USA. Using fuzzy c-means guided inversion of magnetotelluric data to create a starting model for acoustic-impedance proved important in obtaining the best result. Our inversion results correlate with borehole data and provided a better basis for geological interpretation than the seismic reflection images alone. Low values of the acoustic impedance in the basement rocks were shown to be prospective by geochemical analysis of rock cores, as would be predicted for later gold mineralization in weak, decalcified rocks. 相似文献
4.
储层弹性与物性参数可直接应用于储层岩性预测和流体识别,是储层综合评价和油气藏精细描述的基本要素之一.现有的储层弹性与物性参数地震同步反演方法大都基于Gassmann方程,使用地震叠前数据,通过随机优化方法反演储层弹性与物性参数;或基于Wyllie方程,使用地震叠后数据,通过确定性优化方法反演储层弹性与物性参数.本文提出一种基于Gassmann方程、通过确定性优化方法开展储层弹性和物性参数地震叠前反演的方法,该方法利用Gassmann方程建立储层物性参数与叠前地震观测数据之间的联系,在贝叶斯反演框架下以储层弹性与物性参数的联合后验概率为目标函数,通过将目标函数的梯度用泰勒公式展开得到储层弹性与物性参数联合的方程组,其中储层弹性参数对物性参数的梯度用差分形式表示,最后通过共轭梯度算法迭代求解得到储层弹性与物性参数的最优解.理论试算与实际资料反演结果证明了方法的可行性. 相似文献
5.
重力反演方法是研究地壳结构和物性界面起伏的有效地球物理手段之一.本文收集了南北地震带南段67个已有的固定台站接收函数反演的Moho面深度结果,并使用基于EGM2008重力异常模型计算的布格重力异常,验证了本文提出的重震联合密度界面反演方法的有效性.利用接收函数对台站下方Moho面深度估计作为先验约束,定义了一类评价函数,通过对重力反演算法中尺度因子,平移因子和稳定性因子的最优选择,最小化重力反演结果与接收函数模型之间的差异.结果表明,本文提出的方法,可以有效地同化不同地球物理方法获得的反演模型,且通过重震联合反演可以改进由于对空间分布不均匀的接收函数结果插值可能而引起的误差.本文还通过引入Crust1.0的Moho面深度为初值,同时考虑地壳密度的横向不均匀分布,通过模型之间的联合反演有效改善了地球物理反演模型间的不一致性问题.本文反演得到的最优化Moho面深度模型与已知67个台站位置接收函数模型之间的标准差约1.9km,小于Crust1.0与接收函数结果模型之间标准差为3.73km的统计结果.本文研究结果对于同化重震反演结果、精化地壳密度界面模型,都具有十分重要的参考意义. 相似文献
6.
长波长假设条件下，各向同性背景地层中发育一组平行排列的垂直裂缝可等效为具有水平对称轴的横向各向同性（HTI）介质.基于不同观测方位的岩石地震响应特征变化，宽方位地震数据不仅可实现裂缝岩石弹性参数与各向异性参数的预测，同时也蕴含着丰富的孔隙度等储层物性参数信息.本文结合实际地震资料提出了贝叶斯框架下岩石物理驱动的储层裂缝参数与物性参数概率地震联合反演方法，首先基于AVAZ反演裂缝岩石的弹性参数与各向异性参数，并在此基础上通过统计岩石物理模型表征孔隙度、裂缝密度等各向异性介质储层参数与裂缝岩石参数的相互关联，并采用马尔科夫链蒙特卡洛（MCMC）抽样方法进行大量样本的随机模拟，使用期望最大化（EM）算法估计后验条件概率分布，最终寻找最大后验条件概率对应的孔隙度、裂缝密度等HTI裂缝介质储层参数即为反演结果.测井及实际地震数据处理表明，该方法能够稳定合理地从方位地震资料中获取裂缝岩石弹性参数与各向异性参数，并提供了一种较为可靠的孔隙度、裂缝密度等裂缝介质储层参数概率地震反演方法. 相似文献
7.
In geophysical inverse problems, the posterior model can be analytically assessed only in case of linear forward operators, Gaussian, Gaussian mixture, or generalized Gaussian prior models, continuous model properties, and Gaussian-distributed noise contaminating the observed data. For this reason, one of the major challenges of seismic inversion is to derive reliable uncertainty appraisals in cases of complex prior models, non-linear forward operators and mixed discrete-continuous model parameters. We present two amplitude versus angle inversion strategies for the joint estimation of elastic properties and litho-fluid facies from pre-stack seismic data in case of non-parametric mixture prior distributions and non-linear forward modellings. The first strategy is a two-dimensional target-oriented inversion that inverts the amplitude versus angle responses of the target reflections by adopting the single-interface full Zoeppritz equations. The second is an interval-oriented approach that inverts the pre-stack seismic responses along a given time interval using a one-dimensional convolutional forward modelling still based on the Zoeppritz equations. In both approaches, the model vector includes the facies sequence and the elastic properties of P-wave velocity, S-wave velocity and density. The distribution of the elastic properties at each common-mid-point location (for the target-oriented approach) or at each time-sample position (for the time-interval approach) is assumed to be multimodal with as many modes as the number of litho-fluid facies considered. In this context, an analytical expression of the posterior model is no more available. For this reason, we adopt a Markov chain Monte Carlo algorithm to numerically evaluate the posterior uncertainties. With the aim of speeding up the convergence of the probabilistic sampling, we adopt a specific recipe that includes multiple chains, a parallel tempering strategy, a delayed rejection updating scheme and hybridizes the standard Metropolis–Hasting algorithm with the more advanced differential evolution Markov chain method. For the lack of available field seismic data, we validate the two implemented algorithms by inverting synthetic seismic data derived on the basis of realistic subsurface models and actual well log data. The two approaches are also benchmarked against two analytical inversion approaches that assume Gaussian-mixture-distributed elastic parameters. The final predictions and the convergence analysis of the two implemented methods proved that our approaches retrieve reliable estimations and accurate uncertainties quantifications with a reasonable computational effort. 相似文献
8.
Estimating elastic parameters from prestack seismic data remains a subject of interest for the exploration and development of hydrocarbon reservoirs. In geophysical inverse problems, data and models are in general non‐linearly related. Linearized inversion methods often have the disadvantage of strong dependence on the initial model. When the initial model is far from the global minimum, inversion iteration is likely to converge to the local minimum. This problem can be avoided by using global optimization methods. In this paper, we implemented and tested a prestack seismic inversion scheme based on a quantum‐behaved particle swarm optimization (QPSO) algorithm aided by an edge‐preserving smoothing ( EPS) operator. We applied the algorithm to estimate elastic parameters from prestack seismic data. Its performance on both synthetic data and real seismic data indicates that QPSO optimization with the EPS operator yields an accurate solution. 相似文献
9.
基于地质统计先验信息的储层物性参数同步反演 总被引：4，自引：1，他引：3
本文提出的储层物性参数同步反演是一种高分辨率的非线性反演方法,该方法综合利用岩石物理和地质统计先验信息,在贝叶斯理论框架下,首先通过变差结构分析得到合理的变差函数,进而利用快速傅里叶滑动平均模拟算法（Fast Fourier TransformMoving Average,FFT-MA）和逐渐变形算法（Gradual Deformation Method,GDM）得到基于地质统计学的储层物性参数先验信息,然后根据统计岩石物理模型建立弹性参数与储层物性参数之间的关系,构建似然函数,最终利用Metropolis算法实现后验概率密度的抽样,得到物性参数反演结果。并将此方法处理了中国陆上探区的一块实际资料,本方法的反演结果具有较高的分辨率,与测井数据吻合度较高;由于可以直接反演储层物性参数,避免了误差的累积,大大减少了不确定性的传递,且计算效率较高。 相似文献
10.
为了降低单一地球物理方法反演的多解性及受噪声的影响程度,本文围绕重力、磁法和大地电磁法开展了三维联合反演的研究.重、磁采用基于对数障碍法的正则化反演算法,大地电磁使用limited-memory BroydenFletcher-Goldfarb-Shanno(L-BFGS)反演算法,引入交叉梯度函数实现了三种物性结构的相互耦合,最终开发出一套重磁电三维联合反演算法,并实现MPI并行加速计算.通过理论模型算例验证了算法的准确性,结果表明:不论是单棱柱体模型还是组合棱柱体模型,联合反演结果相较单独反演对于异常体的空间形态刻画以及物性数值恢复具有较好的提升;单棱柱体模型算例使得异常体的物性参数(密度、磁化率和电阻率)更加接近于真实的物性参数;组合棱柱体模型的联合反演结果不仅仅消除了围岩物性参数的假异常,而且还增强了异常体边界结构的恢复程度. 相似文献
11.
Full waveform inversion aims to use all information provided by seismic data to deliver high-resolution models of subsurface parameters. However, multiparameter full waveform inversion suffers from an inherent trade-off between parameters and from ill-posedness due to the highly non-linear nature of full waveform inversion. Also, the models recovered using elastic full waveform inversion are subject to local minima if the initial models are far from the optimal solution. In addition, an objective function purely based on the misfit between recorded and modelled data may honour the seismic data, but disregard the geological context. Hence, the inverted models may be geologically inconsistent, and not represent feasible lithological units. We propose that all the aforementioned difficulties can be alleviated by explicitly incorporating petrophysical information into the inversion through a penalty function based on multiple probability density functions, where each probability density function represents a different lithology with distinct properties. We treat lithological units as clusters and use unsupervised K-means clustering to separate the petrophysical information into different units of distinct lithologies that are not easily distinguishable. Through several synthetic examples, we demonstrate that the proposed framework leads full waveform inversion to elastic models that are superior to models obtained either without incorporating petrophysical information, or with a probabilistic penalty function based on a single probability density function. 相似文献
12.
We describe a linear Bayesian inversion method to estimate the relevant petrophysical properties of the media forming a reflecting interface from the observations of amplitude variation with incidence angle. Three main steps characterize the proposed approach:
– information from borehole logs are statistically analysed to estimate the empirical models that describe the functional relationship between petrophysical (e.g. porosity, saturation, pressure or depth) and seismic variable(P and S velocities and density);
– the pure-mode (PP) reflection coefficient is parameterized in terms of the relevant petrophysical variables and is linearized in order to implement the linear inversion;
– the sought petrophysical parameters are estimated from the seismic reflected amplitudes by applying the linearized inversion where a priori information, data and model errors and solutions are described by probability density functions.
We test the method on synthetic and real data relative to reflections from a shale/gas-sand interface where the amplitude versus angle response, besides the lithological contrast, is mainly controlled by the saturation and porosity of the sand layer. The outcomes of the linearized inversion are almost identical to those obtained by a previously developed non-linear inversion method demonstrating the applicability of the linear inversion. It turns out that the gas-sand saturation in the range 0%–95% is a poorly resolved parameter while the porosity is the best resolved parameter. The issues of robustness and resolution of the inversion are discussed either through singular value decomposition analysis or the observation of the a posteriori probability density functions.
The linear inversion algorithm, compared with the previously developed non-linear method, reduces significantly the computation time allowing for more extensive applications. 相似文献
– information from borehole logs are statistically analysed to estimate the empirical models that describe the functional relationship between petrophysical (e.g. porosity, saturation, pressure or depth) and seismic variable(P and S velocities and density);
– the pure-mode (PP) reflection coefficient is parameterized in terms of the relevant petrophysical variables and is linearized in order to implement the linear inversion;
– the sought petrophysical parameters are estimated from the seismic reflected amplitudes by applying the linearized inversion where a priori information, data and model errors and solutions are described by probability density functions.
We test the method on synthetic and real data relative to reflections from a shale/gas-sand interface where the amplitude versus angle response, besides the lithological contrast, is mainly controlled by the saturation and porosity of the sand layer. The outcomes of the linearized inversion are almost identical to those obtained by a previously developed non-linear inversion method demonstrating the applicability of the linear inversion. It turns out that the gas-sand saturation in the range 0%–95% is a poorly resolved parameter while the porosity is the best resolved parameter. The issues of robustness and resolution of the inversion are discussed either through singular value decomposition analysis or the observation of the a posteriori probability density functions.
The linear inversion algorithm, compared with the previously developed non-linear method, reduces significantly the computation time allowing for more extensive applications. 相似文献
13.
We explore the link between basin modelling and seismic inversion by applying different rock physics models. This study uses the E‐Dragon II data in the Gulf of Mexico. To investigate the impact of different rock physics models on the link between basin modelling and seismic inversion, we first model relationships between seismic velocities and both (1) porosity and (2) effective stress for well‐log data using published rock physics models. Then, we build 1D basin models to predict seismic velocities derived from basin modelling with different rock physics models, in a comparison with average sonic velocities measured in the wells. Finally, we examine how basin modelling outputs can be used to aid seismic inversion by providing constraints for the background low‐frequency model. For this, we run different scenarios of inverting near angle partial stack seismic data into elastic impedances to test the impact of the background model on the quality of the inversion results. The results of the study suggest that the link between basin modelling and seismic technology is a two‐way interaction in terms of potential applications, and the key to refine it is establishing a rock physics models that properly describes changes in seismic signatures reflecting changes in rock properties. 相似文献
14.
重力异常对地壳横向密度变化敏感,而无约束重力反演得到的密度模型其垂向分辨能力往往不理想.为了改善反演结果的垂向分辨率,本文参考已有先验分层模型,基于贝叶斯原理,提出了一种重震联合反演的新策略,可实现多种参考模型和复杂加权参数条件下的最大后验概率估计.理论模型测试结果表明,对于深度加权、多参考模型约束等多种问题,本文提出的新方法都可以稳健地获得最优化的模型参数.本文同时以中国地震科学台阵在龙门山地区及周边的一维接收函数分层模型和地震层析成像结果为参考,通过此方法对该区的重力异常进行反演,获得了该区的高精度三维密度结构,其水平分辨率优于10 km,垂直分辨率优于5 km.结合四条通过汶川和芦山地震震中的剖面进行分析后发现,反演得到的密度结构模型在过强震震源区位置横向变形显著,其揭示的分层地壳结构和变形模式与地表已知断裂构造具有相关性.本文提出的重震联合反演新策略,可为研究潜在强震风险源区的地壳结构和物性特征提供有效的科技方法支撑. 相似文献
15.
Mattia Aleardi 《Geophysical Prospecting》2019,67(2):331-348
We use Legendre polynomials to reparameterize geophysical inversions solved through a particle swarm optimization. The subsurface model is expanded into series of Legendre polynomials that are used as basis functions. In this framework, the unknown parameters become the series of expansion coefficients associated with each polynomial. The aim of this peculiar parameterization is threefold: efficiently decreasing the number of unknowns, inherently imposing a 1D spatial correlation to the recovered subsurface model and searching for maximally decoupled parameters. The proposed approach is applied to two highly non-linear geophysical optimization problems: seismic-petrophysical inversion and 1D elastic full-waveform inversion. In this work, with the aim to maintain the discussion at a simple level, we limit the attention to synthetic seismic data. This strategy allows us to draw general conclusions about the suitability of this peculiar parameterization for solving geophysical problems. The results demonstrate that the proposed approach ensures fast convergence rates together with accurate and stable final model predictions. In particular, the proposed parameterization reveals to be effective in reducing the ill conditioning of the optimization problem and in circumventing the so-called curse-of-dimensionality issue. We also demonstrate that the implemented algorithm greatly outperforms the outcomes of the more standard approach to global inversion in which each subsurface parameter is considered as an independent unknown. 相似文献
16.
Fresnel zone inversion for lateral heterogeneities in the earth 总被引：2，自引：0，他引：2
Kiyoshi Yomogida 《Pure and Applied Geophysics》1992,138(3):391-406
We propose a different kind of seismic inversion from travel-time or waveform inversion for lateral heterogeneities in the earth: Fresnel zone inversion. Amplitude and phase delay of data in several frequency ranges are inverted for model space around ray paths with a width corresponding to the considered frequency so that primary effect of finiteness of wavelength be included. For vertically heterogeneous media, Fréchet derivatives for inversion are obtained very efficiently using the paraxial ray approximation, with nearly similar amounts of computation compared to travel-time inversion. As an example, Fréchet derivatives are computed for a teleseismic observation system for a three-dimensional structure in the lithosphere beneath an array of seismic stations. Even if the used frequency is around 2 Hz, the width of Fréchet derivatives cannot be neglected, particularly near the bottom of the lithosphere. Sensitivity of model parameters to observations is, moreover, different in our approach from conventional travel-time inversion: it is zero along ray paths but large slightly away from them. Some model calculations show that travel-time inversion, particularly with models divided into very fine meshes or blocks, might give misleading results. An example of inversion for a simple Camembert model, in the event that travel-time inversion gives no reliable results, shows how this technique works with much smaller data sets and computation than waveform inversions. 相似文献
17.
Characterization of a reservoir model requires determination of its petrophysical parameters, such as porosity and saturation.
We propose a new method to determine these parameters directly from seismic data. The method consists of the computation and
inversion of seismic waveforms. A high frequency method is presented to model wave propagation through an attenuative and
dispersive poroelastic medium. The high frequency approximation makes it possible to efficiently compute sensitivity functions.
This enables the inversion of seismic waveforms for porosity and saturation. The waveform inversion algorithm is applied to
two laboratory crosswell datasets of a water saturated sand. The starting models were obtained using travel time tomography.
The first dataset is inverted for porosity. The misfit reduction for this dataset is approximately 50%. The second dataset
was obtained after injection of a nonaqueous-phase liquid (NAPL), possibly with some air, which made the medium more heterogeneous.
This dataset was inverted for NAPL and air saturation using the porosity model obtained from the first inversion. The misfit
reduction of the second experiment was 70%. Regions of high NAPL and high air saturation were found at the same location.
These areas correlate well with the position of one of the injection points as well as regions of higher NAPL concentrations
found after excavation of the sand. It is therefore possible to directly invert waveforms for pore fluid saturation by taking
into account the attenuation and dispersion caused by the poroelasticity. 相似文献
18.
宽范围物性约束技术容易实现、具有一定容错性,目前已在大地电磁测深（MT）和地震、MT和重力联合反演中实现,但该技术是结合模拟退火算法实现的.差分进化算法（DE）是一种全局优化算法,但该算法在地球物理联合反演领域应用较少.基于此,本文以双种群设置方案为框架改进了DE算法,并提出了基于改进DE算法的宽范围物性约束技术.MT和重力联合反演的模型试验表明：与传统的DE算法相比,改进的DE算法收敛速度更快,寻优能力更强;基于改进DE算法的宽范围物性约束技术可以促进不同岩石物性参数在一定"范围"内实现耦合,既可以利用岩石物性关联的导向作用,又可以发挥优化算法的寻优能力,进而降低地球物理联合反演对先验信息的要求;此外,该技术的实现也验证了宽范围物性约束思想在联合反演领域中的适用性,具有进一步推广至其他优化算法中的潜质. 相似文献
19.
对于陆相沉积环境下的复杂隐蔽岩性储层,由于观测信息不准确,如信息重叠、信息缺失和噪音污染,以及岩石物理关系模糊等原因,储层横向预测存在不惟一性、不稳定性和不确定性.基于线性假定的常规储层横向预测技术已不适用于复杂隐蔽岩性储层的勘探.本文采用一种非线性储层岩性物性褶积模型,建立波阻抗与孔隙度/泥质含量的函数关系;通过多级结构分解和双向边沿子波检测来刻画复杂岩石物理关系;通过Caianiello褶积神经网络实现确定性反演、统计反演和非线性理论三者有机结合;最后联合应用基于逆算子的反演方法和基于正算子的重建算法实现了综合地质、测井和地震波阻抗信息进行高分辨率储层物性参数反演.非线性储层物性参数反演采用多井约束机制和分频反演方式,在陆相和近海油气勘探资料的实际应用中,取得了明显应用效果. 相似文献
20.
Seismic envelope inversion: reduction of local minima and noise resistance 总被引：3，自引：0，他引：3 
下载免费PDF全文
下载免费PDF全文 Waveform inversion met severe challenge in retrieving long‐wavelength background structure. We have proposed to use envelope inversion to recover the large‐scale component of the model. Using the large‐scale background recovered by envelope inversion as new starting model, we can get much better result than the conventional full waveform inversion. By comparing the configurations of the misfit functional between the envelope inversion and the conventional waveform inversion, we show that envelope inversion can greatly reduce the local minimum problem. The combination of envelope inversion and waveform inversion can deliver more faithful and accurate final result with almost no extra computation cost compared to the conventional full waveform inversion. We also tested the noise resistance ability of envelope inversion to Gaussian noise and seismic interference noise. The results showed that envelope inversion is insensitive to Gaussian noise and, to a certain extent, insensitive to seismic interference noise. This indicates the robustness of this method and its potential use for noisy data. 相似文献