共查询到19条相似文献，搜索用时 200 毫秒
1.
利用混合法构建三维数字岩心(英文) 总被引：9，自引：0，他引：9
在弹性波有三维数字岩心描述了岩石的微观孔隙结构。X射线CT扫描是获取三维数字岩心最准确和直接的方法,但实验成本高。本文结合沉积过程模拟和模拟退火算法,提出了重建三维数字岩心的混合法,基于岩石二维图像重建三维数字岩心。利用岩石颗粒沉积算法构建初始数字岩心,作为模拟退火算法的初始状态。运用模拟退火算法调整岩石颗粒和孔隙的相对位置,使三维数字岩心与岩心二维图像具有相似的自相关函数,从而建立三维数字岩心。与传统模拟退火算法相比,该方法计算时间明显减小。运用局部孔隙度理论定量比较了重建数字岩心和岩心X射线微CT图像,两种数字岩心具有相似的均质性和孔隙连通性。利用有限元方法和格子玻尔兹曼方法分别模拟了重建三维数字岩心的地层因素和渗透率,数值模拟结果与实验结果相符。相比传统模拟退火算法,混合法重建数字岩心的传导特性更接近真实岩心的传导特性。 相似文献
2.
数字岩心微观孔隙结构十分复杂,有限元模拟物性参数与弹性参数之间关系是非线性的,直接反演其物性参数准确度低、稳定性差.本文发展了一种数字岩石物理逆建模方法,实现了基于数字岩心的储层参数有效预测.从数字岩心基函数的构建出发,基于有限元方法,计算了一系列具有等间距物性参数值(孔隙度、泥质含量和含水饱和度)的数字岩心弹性参数(体积模量、剪切模量和密度),通过插值算法建立了数字岩心弹性参数三维数据集,从而实现了弹性模量的有限元数值解的快速构建;然后搜索弹性参数的单值等值面,通过等值面的空间交会得到交点,完成储层参数预测.测试结果表明:基于数字岩心逆建模理论的储层参数预测结果与实际模型一致,具有可行性,并且可以通过增加插值点数目提高预测的准确性;孔隙度和泥质含量预测结果稳定性很好,而含水饱和度对噪声的加入较为敏感. 相似文献
3.
地层纵横向非均质性强,工区间数据分布存在差异.这导致基于已有工区数据构建的机器学习储层参数预测模型,推广到新工区会存在较大预测误差.常规地质方法是在岩心与测井响应特征分析基础上建模,利用测井资料计算储层参数,流程复杂.该方法需要岩心校准模型,同样难以快速推广到新的工区.考虑地层纵横向非均质性,本文设计了一种深度Transformer迁移学习网络,通过已有工区的测井与岩心资料构建预测模型,实现未取心新工区储层参数快速准确预测.首先利用无监督学习算法-孤立森林剔除测井数据中存在的异常噪声数据.然后设计Transformer特征提取网络,提高网络特征提取能力,以此深入挖掘测井数据与储层参数的内在联系.最后设计深度迁移学习网络,构建网络损失函数,利用随机梯度下降算法优化网络参数,实现储层参数准确预测.本方案应用于四川南部地区五峰组—龙马溪组页岩储层参数孔隙度、总有机碳含量和总含气量预测.实验结果与工区校正后计算结果、主流机器学习模型预测结果对比,本方案结果与岩心数据具有更高的一致性.应用结果表明：本文方案具有实用性、有效性和可推广性. 相似文献
4.
《地球物理学进展》2017,(3)
以川西地区的须家河组致密储层岩石为研究对象,利用微CT技术结合Avizo软件先进的数学算法构建了三维数字岩心模型,可以表征砂砾岩储层岩石的孔隙结构特征,并将数字岩心和有限元软件Comsol结合,实现了基于数字岩心的水驱气模拟过程的可视化.并在此基础上开展了水驱气模拟,研究微观孔隙结构特征对岩心中气水两相流的影响.研究结果表明:致密砂岩岩心的孔喉分布状态主要呈连片状和孤立状,其中连片状孔隙在空间上连通性好,主要与残余粒间孔或粒间溶蚀孔有关,而孤立状孔隙在空间上多呈孤立分布,主要与粒内溶蚀孔有关;致密砂岩样品等效孔径主要分布范围在0.5μm以下,储层物性差的样品孔隙结构要比储层物性好的样品复杂,且前者的孤立孔隙多且小孔隙占比高,连通孔隙较少,其对渗透率贡献性少;在水驱气的过程中,岩心的微观孔隙结构将改变驱替前缘形状以及造成气水两相流中舌进现象;随着岩心孔隙度和渗透率增大,水驱气的驱替效率增大,残余气饱和度降低. 相似文献
5.
6.
利用水平与竖向谱比（HVSR）方法反演场地速度结构是国际上迅速发展的研究领域.HVSR反演计算实质是一个土层场地模型空间搜索的全局优化问题,当模型搜索空间的复杂程度增大时,目前常用的搜索算法收敛速度慢,计算效率较低.本文实现了一种结合遗传和模拟退火方法优点的混合全局优化HVSR反演算法,通过理论模型和竖向台阵实测数据的检验,表明该算法能获得很好的反演效果,较好地解决了蒙特卡罗方法收敛速度慢,遗传算法收敛早熟和模拟退火算法搜索效率低的问题.本文在此基础上讨论了单台加速度S波记录用于场地速度结构HVSR反演的适用性,为基于单个地震台的地震观测记录反演浅层速度结构提供了一种高效且较为准确的反演方法. 相似文献
7.
8.
最优化测井解释方法能充分利用各种测井资料及地质信息,可以有效地评价复杂岩性油气藏.优化算法的选择是最优化测井解释方法的关键,影响着测井解释结果的准确性.细菌觅食算法(BFA)是新兴的一种智能优化算法,具有较强的全局搜索能力,但在寻优后期收敛速度较慢.复合形算法(CM)局部搜索能力极强,将其与BFA算法相结合构成BFA-CM混合算法,既提高了搜索精度又提高了搜索效率.利用BFA-CM最优化测井解释方法对苏里格致密砂岩储层实际资料进行了处理,计算结果与岩心及薄片分析资料吻合得很好. 相似文献
9.
实验室声发射三维定位软件 总被引：5，自引：2，他引：5
简要介绍了地震动力学实验室新一代声发射定位软件的主要功能、定位算法的实现方案以及定位效果的评测结果。新软件配有多种通用数据库接口,具有与设备无关的处理能力。它除了可以进行三维定位之外,还可以进行面定位与线定位。定位时,可以给定样品的波速,也可以把波速当未知数求解。定位算法的主要改进有两个方面:1)以稳健的走时绝对偏差最小作为定位目标函数;2)方程组求解优化过程采用了模拟退火法直接搜索目标函数的全局最小值。新软件显著提高了定位的精度,同时提供了一个方便快捷的数据处理平台 相似文献
10.
《地球物理学进展》2017,(5)
在数字岩石物理中数字岩心是弹性参数模拟的基础,现有的数字岩心建模方法虽然较多,但一般不能满足弹性参数模拟对骨架所提出的要求.过程法以岩石粒径分布曲线为约束条件构建数字岩心,是一种比较灵活的数字岩心建模方法.本文根据前人研究结果在过程法的压实过程中用一面状像素层表示颗粒与颗粒间的接触边界,使相邻的骨架颗粒能够相互区分,提高了数字岩心骨架结构的精细程度.在弹性参数模拟中,当变化数字岩心中的颗粒、颗粒接触边界和胶结物的弹性模量时,数字岩心的弹性模量发生相应的变化.数字岩心弹性模量的变化量与组成物质的体积含量及其弹性模量赋值有关.这说明能够利用数字岩石物理研究各种地质和环境因素对岩石弹性性质的影响.进一步分析认为,目前的建模方法对数字岩心的定义还不够充分,利用数字岩石物理解决实际地质问题尚有难度. 相似文献
11.
《Advances in water resources》1998,21(6):499-508
Global optimization methods such as simulated annealing, genetic algorithms and tabu search are being increasingly used to solve groundwater remediation design and parameter identification problems. While these methods enjoy some unique advantages over traditional gradient based methods, they typically require thousands to tens of thousands of forward simulation runs before reaching optimal or near-optimal solutions. Thus, one severe limitation associated with these global optimization methods is very long computation time. To mitigate this limitation, this paper presents a new approach for obtaining, repeatedly and efficiently, the solutions of a linear forward simulation model subject to successive perturbations. The proposed approach takes advantage of the fact that successive forward simulation runs, as required by a global optimization procedure, usually involve only slight changes in the coefficient matrices of the resultant linear equations. As a result, the new solution to a system of linear equations perturbed by the changes in aquifer properties and/or sinks/sources can be obtained as the sum of a non-perturbed base solution and the solution to the perturbed portion of the linear equations. The computational efficiency of the proposed approach arises from the fact that the perturbed solution can be derived directly without solving the linear equations again. A two-dimensional test problem with 20 by 30 nodes demonstrates that the proposed approach is much more efficient than repeatedly running the simulation model, by more than 15 times after a fixed number of model evaluations. The ratio of speedup increases with the number of model evaluations and also the size of simulation model. The main limitation of the proposed approach is the large amount of computer memory required to store the inverse matrix. Effective ways for limiting the storage requirement are briefly discussed. 相似文献
12.
非线性地球物理反演方法:回顾与展望 总被引：19，自引：13，他引：19
自20世纪90年代以来，非线性地球物理反演方法已走向成熟，这些方法包括线性化迭代法，仿真淬火法，遗传算法及联合反演方法等，线性化迭代反演基于开放物理系统状态发生相变的原理，要进一步改善模型参数化的技术和迭代过程的自组织；仿真淬火法与遗传算法基于自然过程的指数率或生物演化的优生率，可以相互结合以提高解估计的分辨率与置信度；联合反演要结合岩石物理性质的统计规律，才能取得兼容地质与综合方法的应用效果，地球动力学中的反问题不仅涉及偏微分方程系数项的求解，而且涉及初始条件或初始边界的求解，对地球动力学作用过程研究特别重要。 相似文献
13.
It is well known that seismic inversion based on local model optimization methods, such as iterative use of linear optimization, may fail when prior information is sparse. Where the seismic events corresponding to reflectors of interest remain to be identified, a global optimization technique is required. We investigate the use of a global, stochastic optimization method, that of simulated annealing, to solve the seismic trace inversion problem, in which the two-way traveltimes and reflection coefficients are to be determined. The simulated annealing method is based on an analogy between the model-algorithm system and a statistical mechanical system. We exploit this analogy to produce improved annealing schedules. It is shown that even in cases of virtually no prior information about two-way traveltimes and reflection coefficients, the method is capable of producing reliable results. 相似文献
14.
The problem of equivalence in direct current (DC) resistivity and electromagnetic methods for a thin resistive and conducting layer is well‐known. Attempts have been made in the past to resolve this problem through joint inversion. However, equivalence still remains an unresolved problem. In the present study, an effort is made to reduce non‐uniqueness due to equivalence using global optimization and joint inversion by successive refinement of the model space. A number of solutions derived for DC resistivity data using very fast simulated annealing global inversion that fits the observations equally well, follow the equivalence principle and show a definite trend. For a thin conductive layer, the quotient between resistivity and thickness is constant, while for a resistive one, the product between these magnitudes is constant. Three approaches to obtain very fast simulated annealing solutions are tested. In the first one, layer resistivities and thicknesses are optimized in a linear domain. In the second, layer resistivities are optimized in the logarithmic domain and thicknesses in the linear domain. Lastly, both layer resistivities and thicknesses are optimized in the logarithmic domain. Only model data from the mean models, corresponding to very fast simulated annealing solutions obtained for approach three, always fit the observations. The mean model defined by multiple very fast simulated annealing solutions shows extremely large uncertainty (almost 100%) in the final solution after inversion of individual DC resistivity or electromagnetic (EM) data sets. Uncertainty associated with the intermediate resistive and conducting layers after global optimization and joint inversion is still large. In order to reduce the large uncertainty associated with the intermediate layer, global optimization is performed over several iterations by reducing and redefining the search limits of model parameters according to the uncertainty in the solution. The new minimum and maximum limits are obtained from the uncertainty in the previous iteration. Though the misfit error reduces in the solution after successive refinement of the model space in individual inversion, it is observed that the mean model drifts away from the actual model. However, successive refinement of the model space using global optimization and joint inversion reduces uncertainty to a very low level in 4–5 iterations. This approach works very well in resolving the problem of equivalence for resistive as well as for conducting layers. The efficacy of the approach has been demonstrated using DC resistivity and EM data, however, it can be applied to any geophysical data to solve the inherent ambiguities in the interpretations. 相似文献
15.
在数值模拟中,隐式有限差分具有较高的精度和稳定性.然而,传统隐式有限差分算法大多由于需要求解大型矩阵方程而存在计算效率偏低的局限性.本文针对一阶速度-应力弹性波方程,构建了一种优化隐式交错网格有限差分格式,然后将改进格式由时间-空间域转换为时间-波数域,利用二范数原理建立目标函数,再利用模拟退火法求取优化系数.通过对均匀模型以及复杂介质模型进行一阶速度-应力弹性波方程数值模拟所得单炮记录、波场快照分析表明:这种优化隐式交错网格差分算法与传统的几种显式和隐式交错网格有限差分算法相比不但降低了计算量,而且能有效的压制网格频散,使弹性波数值模拟的精度得到有效的提高. 相似文献
16.
We derive a governing second-order acoustic wave equation in the time domain with a perfectly matched layer absorbing boundary condition for general inhomogeneous media. Besides, a new scheme to solve the perfectly matched layer equation for absorbing reflections from the model boundaries based on the rapid expansion method is proposed. The suggested scheme can be easily applied to a wide class of wave equations and numerical methods for seismic modelling. The absorbing boundary condition method is formulated based on the split perfectly matched layer method and we employ the rapid expansion method to solve the derived new perfectly matched layer equation. The use of the rapid expansion method allows us to extrapolate wavefields with a time step larger than the ones commonly used by traditional finite-difference schemes in a stable way and free of dispersion noise. Furthermore, in order to demonstrate the efficiency and applicability of the proposed perfectly matched layer scheme, numerical modelling examples are also presented. The numerical results obtained with the put forward perfectly matched layer scheme are compared with results from traditional attenuation absorbing boundary conditions and enlarged models as well. The analysis of the numerical results indicates that the proposed perfectly matched layer scheme is significantly effective and more efficient in absorbing spurious reflections from the model boundaries. 相似文献
17.
18.
2D多尺度混合优化地球物理反演方法及其应用(英文) 总被引：1，自引：0，他引：1
局部优化和全局优化方法广泛应用到地球物理反演,但是两者各有其优缺点。将两类方法结合起来可以取长补短。将退火遗传算法（SAGA）和单纯形算法相结合,得到了一种高效、健全的2D非线性混合地震走时反演方法。首先,利用SAGA进行大范围的全局搜索,然后由单纯形方法进行快速局部搜索。为了降低层析成像的多解性,我们采用了多尺度逐次逼近的技巧。把速度场划分为不同的空间尺度,定义网格节点上的速度作为待反演参数,采用双三次样条函数模型参数化,正问题采用有限差分走时计算方法,反问题采用多尺度混合反演方法。一个低速度异常体的数值模拟试验和抗走时扰动试验表明该方法是有效和健全的。我们将该方法应用到青藏高原东北缘阿尼玛卿rlet,Meyer,Marr,缝合带东段上部地壳速度结构研究中。数字模型试验和实际资料的应用表明了方法的有效性和健全性。 相似文献
19.
Use of VFSA for resolution, sensitivity and uncertainty analysis in 1D DC resistivity and IP inversion 总被引：1，自引：0，他引：1
We present results from the resolution and sensitivity analysis of 1D DC resistivity and IP sounding data using a non-linear inversion. The inversion scheme uses a theoretically correct Metropolis–Gibbs' sampling technique and an approximate method using numerous models sampled by a global optimization algorithm called very fast simulated annealing (VFSA). VFSA has recently been found to be computationally efficient in several geophysical parameter estimation problems. Unlike conventional simulated annealing (SA), in VFSA the perturbations are generated from the model parameters according to a Cauchy-like distribution whose shape changes with each iteration. This results in an algorithm that converges much faster than a standard SA. In the course of finding the optimal solution, VFSA samples several models from the search space. All these models can be used to obtain estimates of uncertainty in the derived solution. This method makes no assumptions about the shape of an a posteriori probability density function in the model space. Here, we carry out a VFSA-based sensitivity analysis with several synthetic and field sounding data sets for resistivity and IP. The resolution capability of the VFSA algorithm as seen from the sensitivity analysis is satisfactory. The interpretation of VES and IP sounding data by VFSA, incorporating resolution, sensitivity and uncertainty of layer parameters, would generally be more useful than the conventional best-fit techniques. 相似文献