共查询到20条相似文献，搜索用时 125 毫秒
1.
《应用地球物理》2021,(1)
常规的地震随机反演方法一般用变差函数表征地下地层的空间结构信息,但是变差函数很难有效表征实际复杂地层的非均质特性。本文利用多参数的思想充分发掘已知地震和测井资料中蕴含的地下地层信息,利用多个统计特征参数更全面地描述地下复杂储层的空间结构特征,并与地震随机反演相结合,提出了基于统计特征参数先验信息的叠前地震随机反演方法。该方法基于随机介质理论从已知的地震和测井资料中求取多个统计特征参数,选取符合实际地质特征的椭圆自相关函数,构建满足地层空间结构特征的先验信息模型,并将多参数约束融入似然函数中构建目标函数,用非常快速量子退火算法进行优化更新,实现目标函数的求解得到最终的反演结果。由模型试算可知,相比于传统的先验模型构建方法,本文基于多个参数构建的先验模型包含更丰富的地层信息,能更好地实现对地下复杂储层的描述。实际资料分析表明,基于统计特征参数先验信息的随机反演方法能够有效地进行地下复杂储层的地球物理预测,且具有较高的分辨率。 相似文献
2.
基于遗传算法和统计学的电阻率测深二维反演研究 总被引：2，自引：0，他引：2
本文介绍了一种以遗传算法和统计学为基础的反演方法,该方法利用一维反演结果的统计学特征确定地层电性结构并设计初始模型,使用遗传算法修改模型,经过二维有限单元法正演达到拟合误差最小,从而完成电阻率测深数据的二维反演.这个方法的优点在于不需要解大型线性方程组,降低计算量,先判断地层电性结构,反演后给出用地电断面表示的直观的解释结果.理论模型试算和模型实验证明该方法效果较好. 相似文献
3.
三维重力反演是地质工作者了解地球深部构造,认知地下结构的重要手段.按照反演单元划分,三维重力反演有离散多面体(Discrete)反演和网格节点(Voxels)反演两种方式.离散多面体反演由于易于吸收先验地质信息得到的理论场能够很好地拟合观测场,因此,在实际重力反演中更受欢迎.目前离散多面体重力反演中初始模型的建立方法繁杂不一,实际应用受到很大的限制.本文本着充分挖掘利用先验信息和重力观测数据得到丰富可靠的反演结果这一原则,以离散多面体反演技术为基础,改进建模过程.在初始模型的建立中,吸收贝叶斯算法优势,采用隐马尔科夫链改善朴素贝叶斯方法的分类效果,通过最大似然函数算法求解,再采取模型降阶技术,固定所建模型中几何体的形态或密度,达到在几何体形态(x,y,z)、密度(σ)和重力值(g)五个参数中降低维数目的,从而减小高维不确定性和正演的计算量,由此反演计算的地质体密度和分布范围相对更准确,更利于重现重力模型结构.通过单位球体和任意形态几何体模拟实验,以及安徽省泥河矿区三维重力反演实践,得到非常接近实际的密度或重力值,大幅提高了三维重力反演的精度和效率,说明该方法是有效、实用的. 相似文献
4.
基于地质统计先验信息的储层物性参数同步反演 总被引：4，自引：1，他引：3
本文提出的储层物性参数同步反演是一种高分辨率的非线性反演方法,该方法综合利用岩石物理和地质统计先验信息,在贝叶斯理论框架下,首先通过变差结构分析得到合理的变差函数,进而利用快速傅里叶滑动平均模拟算法（Fast Fourier TransformMoving Average,FFT-MA）和逐渐变形算法（Gradual Deformation Method,GDM）得到基于地质统计学的储层物性参数先验信息,然后根据统计岩石物理模型建立弹性参数与储层物性参数之间的关系,构建似然函数,最终利用Metropolis算法实现后验概率密度的抽样,得到物性参数反演结果。并将此方法处理了中国陆上探区的一块实际资料,本方法的反演结果具有较高的分辨率,与测井数据吻合度较高;由于可以直接反演储层物性参数,避免了误差的累积,大大减少了不确定性的传递,且计算效率较高。 相似文献
5.
本文提出的储层物性参数同步反演是一种高分辨率的非线性反演方法,该方法综合利用岩石物理和地质统计先验信息,在贝叶斯理论框架下,首先通过变差结构分析得到合理的变差函数,进而利用快速傅里叶滑动平均模拟算法(Fast Fourier TransformMoving Average,FFT-MA)和逐渐变形算法(Gradual Deformation Method,GDM)得到基于地质统计学的储层物性参数先验信息,然后根据统计岩石物理模型建立弹性参数与储层物性参数之间的关系,构建似然函数,最终利用Metropolis算法实现后验概率密度的抽样,得到物性参数反演结果。并将此方法处理了中国陆上探区的一块实际资料,本方法的反演结果具有较高的分辨率,与测井数据吻合度较高;由于可以直接反演储层物性参数,避免了误差的累积,大大减少了不确定性的传递,且计算效率较高。 相似文献
6.
针对重力与地震联合反演存在的问题,结合已有的研究成果,本文研究实现了速度和密度随机分布共网格单元模型的建模技术,以适应密度和速度剧烈变化的复杂模型及联合反演的计算要求.重力正演利用了该网格的二度半体模型,并进一步改进了地震走时的二维射线追踪计算方法,以适用于速度随机分布的网格介质.结合改进的模拟退火算法,实现了这种共网格条件下的重力与地震资料的同步联合反演.模型试验证明了重力与地震联合反演可以准确确定复杂物性界面的密度和速度结构,适用于物性界面不完全一致和物性变化剧烈的复杂模型,并且联合反演结果要优于单独的重力反演.带先验信息约束下的实际资料的联合反演,进一步证明了该方法的适用性和效果,可提高反演精度并减少多解性. 相似文献
7.
利用能够整合测井信息与井间地震信息的地质统计学随机模拟方法,结合传统的地质统计学反演思路,得到了一种能够同时整合测井、井间地震与地面地震三种先验信息的地质统计学反演与储层建模方法.由于井间射线信息、测井信息与地面地震数据在随机反演与建模过程当中都得到了尊重,因此与传统地质统计学反演仅利用了测井与地面地震数据相比,本文的地质统计学反演与建模方法更充分地利用了先验信息,有效提高了反演的精度,降低了随机建模中的多解性.基于理论数据的测试证实了上述观点. 相似文献
8.
地磁测深研究的周期范围通常为10~5~10~7s,缺少反映浅部电性结构的短周期信息,而C-响应受浅部电阻率影响明显,因此本文提出在反演中增加浅部(约200 km)电阻率约束以提高深部反演的稳定性和可靠性.在磁层环状电流满足P_1~0假设的条件下,球坐标系中一维导电薄球层状地球的C-响应和电导率分布关系由边界条件通过递推的方法计算得到.反演采用有限内存拟牛顿(L-BFGS)法;浅部电阻率约束通过将目标函数对模型参数的梯度设为零来实现;通过置信区间分析评价约束反演结果的可靠性.合成数据的无约束反演虽然最终的拟合效果很好,但浅部电阻率受初始模型影响,差异较大;采用浅部约束后,反演结果对初始模型依赖性明显减小,同时还能显著提高200~600 km范围内反演结果的准确性.对全球近地轨道卫星观测的C-响应数据约束反演后结果与前人一致,表现为地幔电导率整体上随着深度的增加而增加.参数置信区间分析表明,由于约束反演加入了浅部信息,电阻率的变化范围更加紧致,说明反演结果更加可靠.因此,有必要通过其他地球物理方法,如长周期大地电磁测深等获得浅部电阻率分布,作为先验信息参加反演,进行浅部约束的C-响应反演,获得更可靠的一维全深度电性结构,为地磁测深数据解释奠定基础. 相似文献
9.
在已有研究成果的基础上,为了适应物性参数剧烈变化的复杂模型并满足联合反演的要求,开发了速度和电阻率随机分布共网格单元模型的建模技术.基于这种统一的物性随机分布的网格介质模型,利用有限元方法和改进的射线追踪法分别正演计算大地电磁场和地震走时,结合改进的模拟退火算法,研究实现了电阻率和速度随机分布条件下的大地电磁与地震资料的同步联合反演.对物性界面不完全一致和物性变化剧烈的带地形复杂模型的试验,表明了该方法在精细反演复杂电阻率和速度结构方面的效果,克服了以往研究局限于简单模型的不足.对地震资料品质差的地区开展的实际资料联合反演,表明了方法的适用性,先验信息约束下的联合反演提高了反演精度. 相似文献
10.
大地电磁资料反演方法较多,但构建二维反演初始模型常规方法主要有两种:以均匀半空间为初始模型、直接以一维反演结果为初始模型.本文利用典型实测剖面资料对这两种方法进行反复试验对比,结果表明:在地形及构造复杂地区,由于实测资料存在严重的静态效应,利用常规方法构建反演初始模型进行二维反演,难以达到理想的反演效果.基于此,本文提出了一种构建二维反演初始模型与反演新的思路方法,其核心为用于反演与构建初始模型数据为同一基础数据,但对基础数据分别进行处理,得到两个处理及校正程度不同的数据,分别用于初始模型构建与二维反演输入数据,这样得到的反演结果既能较好地保留浅部至深部有用异常信息,又能够最大程度上消除实测资料中的静态效应.本文进一步对实测剖面反演结果与已知地质信息进行对比,验证了本方法的有效性和可靠性,认为此方法为地形及地质条件复杂地区大地电磁资料处理反演提供了新的思路,具有一定的参考意义和借鉴作用. 相似文献
11.
在频率域弹性波有限元正演方程的基础上,依据匹配函数（也就是观测数据和正演数据残差的二次范数）最小的准则,用矩阵压缩存储与LU分解技术来存储和求解频率域正演方程中的大型稀疏复系数矩阵、用可调阻尼因子的Levenberg Marquard方法求解反演方程组,直接求取地下介质的弹性波速度,导出了频率域弹性波有限元最小二乘反演算法. 为了利用地下地质体的分布规律,减少反演所求的未知数个数,本文又提出了规则地质块体建模方法引入到反演中来. 经数值模型验证,在噪声干扰很大（噪声达到50髎）或初始模型与真实模型相差很大的情况下,反演也能取得很满意的效果,证明本方法具有很好的抗噪性与“强壮性”. 相似文献
12.
GPS AND LEVELING CONSTRAINED CO-SEISMIC SOURCE AND SLIP DISTRIBUTION OF THE LUSHAN MS7.0 EARTHQUAKE ON 20 APRIL 2013 
下载免费PDF全文
下载免费PDF全文 Co-seismic deformation results calculated by the observations of GPS continuous sites and precise leveling are used as three-dimensional constraints for the deformation field of Lushan MS7.0 earthquake on 20 April 2013.The inversion of seismic source parameters are calculated based on the consideration of different value-taking schemes and fault models. Seven data type selection schemes, three fault models and two data coverage selection schemes are designed in order to discuss the effect of data selection and fault model selection on the inversion results.
The results show that the fault model using blind reverse fault for the inversion is superior to the model using the fault that ruptures from its upper boundary to the earth'surface, which may indicate that the Lushan earthquake fault is most likely a blind reverse fault; there are no obvious differences in the inversion results between the blind listric reverse fault models and the rectangle blind reverse fault models.
The best inversion result of Lushan earthquake seismic moment is MW6.7.And it also shows that the distribution of dislocations on the fault plane is concentrated in the range of 30km×30km, the northern flank of the seismogenic fault of Lushun earthquake is of dextral faulting and the southern flank of the fault is of sinistral faulting, the sinistral component is larger that dextral component, showing a wedge deformation mode. 相似文献
13.
Real time 2D inversion for an induction logging instrument may be achieved using a fast forward modeling and special inversion strategy. The fast forward modeling employs a low-frequency approximation of an induction response known as Doll's geometric factor. Modeling geometric factors is much faster than modeling the electromagnetic field in the frequency domain. To transform real data into the Doll's limit, multi-frequency skin-effect correction is applied. The correction technique involves an asymptotic theory of the integral equation for a 2D boundary value problem. The inversion is based on separating the parameter space into subspaces of lower dimension. Initially, adaptive overlapping windows split logging data into manageable portions. Each window consists of three subwindows: the predictor, corrector and upgrader. Further separation of parameters is introduced by Doll's approximation: the low-frequency response is linear with respect to formation conductivity. This allows us to split inversion for conductivity and geometric parameters. The next level of splitting inversion is achieved by independently determining parameters of the near borehole zone and remote formation areas. This is done by utilizing different subsets of sensors. The inversion does not require initial guess: layers are introduced dynamically, if necessary. The resolution is improved in sequential iterations by adding finer details to the previously obtained models. The final selection of parameters satisfies a variety of a priori constraints formulated as target resistivity distributions. The technique for imposing constraints is based on the analysis of data mapping into the model space. Interpretation of synthetic and real data confirms the viability of the method. 相似文献
14.
高频电磁波测井响应对地层纵向边界相当敏感 ,在其资料反演中必须重视地层纵向边界位置的提取和反演精度 .首先应用特征识别技术由测井曲线提取地层参数 ,在此基础上建立反演初始地层模型 ,通过分别构造反演地层电参数和纵向边界位置的目标函数 ,采用改进的阻尼型高斯 -牛顿算法 ,实现了对地层电参数和纵向边界位置的交替循环反演 .对含1 7个水平分层的典型地层模型的模拟测井资料进行反演计算 ,测井响应的反演精度达到2 .0 % (平均相对误差 ) ,电参数的误差为 4.4% . 相似文献
15.
在界面两侧地层的弹性参数弱反差的假设难以成立的情况下,本文提出用平均入射角道集进行PP波与PS波的联合反演.首先,在PP波与PS波AVA(amplitude versus angle,振幅随入射角变化)道集的基础上,分别选择小入射角范围与大入射角范围的AVA道集进行局部加权叠加,以获得由两个角度组成的平均入射角道集,并作为后续反演的输入数据.然后,再通过最小二乘原理建立了PP波与PS波联合反演目标函数,推导了模型修改量的向量公式,建立了平均入射角道集联合反演的流程.模型数据与实际数据的测试结果表明:在信噪比较低、地层弹性参数反差较大、层厚较薄的情况下,该反演方法的精度在很大程度上超过了基于近似反射系数的反演方法,为复杂油气藏勘探提供了新的思路. 相似文献
16.
17.
本文介绍了用GPS资料反演断层运动中存在的一些问题。GPS资料区域构造信息的提取是反演中非常重要的环节,结合一定区域反演时选择符合断层运动的位错模型是必要的,断层运动参数的选取问题,最后是反演的算法问题。 相似文献
18.
跨孔雷达全波形反演是一种使用全波形信息反演两钻孔之间地下信息的层析成像技术.常规的层析成像反演大部分采用射线追踪方法,其中基于初至时的射线追踪方法可以反演出速度剖面（介电常数剖面）,基于最大振幅的层析成像可以反演出衰减剖面（电导率剖面）.常规射线追踪方法有许多不足,究其原因是该方法仅使用了小部分的信号信息.为了进一步提高成像分辨率,本文全面推导了全波形跨孔雷达层析成像反演方法,该方法利用雷达波全幅度相位信息能够反演出地下高分辨率的介电常数和电导率图像.本文通过基于局域网的分布式并行算法,有效地解决了巨量数据正演计算问题.文中首先建立了基于单轴各向异性介质完全匹配层的时间域有限差分二维正演算法,进而通过应用包括时间维度在内的全波场信息与残场逆向传播的全波场信息乘积来计算梯度方向,通过求取以步长为自变量的目标函数的极值确定步长公式,并提出以第一次介电常数反演作为同步反演的初始模型,能够有效提高收敛速度.本文对多组模型进行成像实验,取得了较好的反演效果. 相似文献
19.
中国大陆中西部乃至全球造山带普遍具有复杂地壳结构.随着矿产资源勘探和深部探测研究的深入,探测造山带及盆山耦合区下方地壳精细结构正逐渐成为当前面临的巨大挑战.人工源深地震测深方法正越来越清晰地揭示出不同构造域地壳速度结构的基本特征,然而传统的层状结构模型参数化方法难以准确描述复杂地质模型,通常情况下多忽略速度结构的精细间断面且采用层边界平滑处理,难以满足地壳精细结构成像的发展要求.针对上述困难,本文采用最近发展的块状结构建模方案构建三维复杂地壳模型,基于逐段迭代射线追踪正演走时计算方法,推导了走时对三角形界面深度以及网格速度的偏导数,开展了非线性共轭梯度走时反演方法研究.发展了利用直达波和反射波等多震相走时数据对界面深度和网格速度的多参数联合反演方法,并引人不同种类震相数据的权系数和不同类型参数偏导数归一化的方法.数值算例表明,基于块状结构的非线性共轭梯度走时反演方法适用于复杂地壳结构模型,在利用人工源走时数据反演复杂地壳精细结构领域具有良好的应用前景. 相似文献
20.
从各向异性弹性波的有限元正演方程出发 ,导出了反问题中时间域雅可比矩阵求解的计算公式 .它具有与时间域有限元正演方程相同的表达形式 ,故可通过有限元正演计算来获得雅可比矩阵 .研究了有限元正演算法的效率和精度、吸收边界条件等方面的问题 ,以提高反演系统的效率和精度 .在此基础上 ,实现了叠前全波场各向异性弹性参数反演 .计算表明 ,在初始模型偏离真实模型较大的情况下 ,层状模型和横向不均匀模型的反演结果均能准确地收敛到真实模型上 . 相似文献