共查询到20条相似文献，搜索用时 155 毫秒
1.
作为全局非线性优化的新方法之一的遗传算法,近年来已从生物工程流行到大地电磁测深资料解释中.然而,大地电磁反演问题具有不适定性,解的非唯一性.通过结合求解不适定问题的Tikhonov正则化方法,本文采用实数编码遗传算法求解大地电磁二维反演问题.此算法在构建目标函数时引入正则化的思想,利用遗传算法求解最优化问题.常规的基于局部线性化的最优化反演方法易使解陷入局部极小值,而且严重的依赖初始模型的选择.与传统线性化的迭代反演方法相比,实数编码遗传算法能够克服传统方法的不足且能获得更好的反演结果.通过对大地电磁测深理论模型进行计算,结果表明:该算法具有收敛速度快、解的精度高和避免出现早熟等优点,可用于大地电磁资料解释. 相似文献
2.
3.
4.
5.
常规非线性反演方法虽然对初始模型的依赖大为减弱，但局部收敛现象和计算速度慢仍然是瓶颈. 本文提出了一种新的反演方法——量子路径积分算法(Quantum Path Integral Algorithm,简称QPIA). 该方法引入量子力学的横向场、传播子等概念，并充分利用量子隧穿效应，大大提高反演的效率，具体是通过对反演目标函数的构建，并以Feynman的传播子来构成模型的接收概率来实现. 在对一维大地电磁模型和实际数据进行试验后，表明该方法比常规反演方法更能够精确、稳定和快速地逼近真实模型. 相似文献
6.
7.
8.
大地电磁与地震联合反演研究现状与展望 总被引：2，自引：0，他引：2
由于大地电磁探测和地震勘探各自存在优势和局限,因此联合反演比单独一种地球物理资料反演更优越.本文剖析了联合反演的两类分类方法,指出大地电磁与地震方法具备联合反演的基础-理论上相似,物性上有联系,并重点分析国内外大地电磁与地震数据联合反演的研究现状与存在问题,总结层状均匀和非均匀两种介质结构的电震联合反演算法.对于如何耦合不同类型的物性参数这一联合反演的核心问题,指出目前广泛使用的岩石物理法和结构法两类方法的适用范围.最后提出大地电磁与地震联合反演的发展方向:层状均匀介质结构要优化线性迭代技术;非均匀的复杂地质结构,需加紧非线性反演理论向非线性反演发展,其中结构法反演前景广阔,要加强实际资料的运用;此外,电震联合反演也需与岩石物性等先验地质信息相结合. 相似文献
9.
本文从大地电磁扬的色散关系出发,阐述了进行大地电磁资料一维联合反演的物理意义。论证了对于以色散关系相联系的一对大地电磁响应函数,特别是阻抗实部和虚部视电阻率进行联合反演的必要性和可能的应用前景。对比研究了大地电磁阻抗各单参量及其联合反演的效果。结果表明,在加快反演迭代速度,限定反演解的非唯一性、保证反演迭代的稳定进行等方面,联合反演比单参量反演具有明显的优越性;在大多数情况下,阻抗实部视电阻率和阻抗虚部视电阻率联合反演的效果要比通常所采用的常规视电阻率和阻抗相位的联合反演效果好。实际大地电磁资料的试验表明,各种不同平均阻抗的利用和综合反演解释将有助于获得更为可靠的结果。 相似文献
10.
适用于近地表(2000m以内)勘探的频率域电磁法主要有音频大地电磁法(audio-frequency magnetotellurics,AMT),无线电大地电磁法(radio-magnetotellurics,RMT),可控源音频大地电磁法(controlled source audiofrequency magnetotellurics,CSAMT),广域电磁法(Wide Field Electromagnetic Method,WFEM).本文拟从最新的数据采集技术、数据处理技术、正反演算法、实例等四个方面,论述适用于浅部勘探的AMT,RMT,CSAMT和WFEM方法的国内外最新进展,总结目前AMT,RMT,CSAMT和WFEM方法遇到的困难,对潜在的发展方向提出建议.综述表明:(1)张量测量、多站阵列、多站叠加可提高AMT、RMT和CSAMT数据的质量.利用近区数据WFEM法可获得良好的效果.国产与国外仪器在质量方面的差距正在逐步缩小.(2)数学形态滤波技术、Hilbert-Huang变换等可有效分离出有用的数据,局部畸变仍然是亟待解决的难题,需要更为深入的研究.(3)矢量有限元与非结构网格的出现大幅度提高了有限元处理复杂电磁问题模拟的精度与应用范围,成为目前电磁正演的首选工具.完全非线性反演算法仍然局限于1D、2D问题,共轭梯度法和高斯牛顿算法等为解决3D问题的发展趋势.地质约束的引入和多数据联合反演可以减小反演的非唯一性.各向异性的反演为目前反演研究的热点之一.(4)野外数据解释的正确性严重依赖于对地下结构先期的维性判别,在2D特性不明显、3D特性明显时,需要采用3D进行反演解释. 相似文献
11.
本文将大地电磁反演中的RRI方法应用于线源频率测深模型资料的反演中.当背景电导率变化很小时,可用变化前的电场来近似代替变化后的电场,从而在反演方程的推导过程中将有源电磁场中的源项消掉,得到和大地电磁场相同的反演方程,使有源电磁波的反演也可以应用RRI方法.反演过程中所需要的模型资料通过有限元方法得到,该资料不需做近场校正,直接用适合于有源电磁场的RRI方法反演,避免了近场校正带来的误差.数值模型结果证明该方法是可行的.最后,用RRI方法讨论了当源和目标区间存在低阻异常体时只对目标区反演的可行性,对实际工作的解释有一定的指导意义. 相似文献
12.
在已有研究成果的基础上,为了适应物性参数剧烈变化的复杂模型并满足联合反演的要求,开发了速度和电阻率随机分布共网格单元模型的建模技术.基于这种统一的物性随机分布的网格介质模型,利用有限元方法和改进的射线追踪法分别正演计算大地电磁场和地震走时,结合改进的模拟退火算法,研究实现了电阻率和速度随机分布条件下的大地电磁与地震资料的同步联合反演.对物性界面不完全一致和物性变化剧烈的带地形复杂模型的试验,表明了该方法在精细反演复杂电阻率和速度结构方面的效果,克服了以往研究局限于简单模型的不足.对地震资料品质差的地区开展的实际资料联合反演,表明了方法的适用性,先验信息约束下的联合反演提高了反演精度. 相似文献
13.
实现大地电磁法快速三维反演的关键在于找到快速计算灵敏度矩阵的方案. 本文在对大地电磁三维张量阻抗表达式进行深入分析的基础上，获得了三维快速松弛反演算法的灵敏度函数表达式，解决了三维快速松弛反演的核心问题；为了减小反演解的多解性，针对三维问题定义了最小构造函数，实现了求最小构造的三维快速松弛反演算法. 对二维棱柱体、三维棱柱体理论模型的大地电磁合成数据进行了反演试算，反演结果与理论模型相吻合. 对日本Kayabe地区实测资料进行了XY模式反演，反演得到的地电模型较好地反映了地热资源的分布. 反演试算结果表明大地电磁三维快速松弛反演具有计算速度快、结果稳定可靠等特点；使三维快速松弛反演在普通微机上得以实现. 相似文献
14.
针对常规大地电磁（Magnetotelluric,MT）反演方法对电阻率异常体边界不太敏感的问题,本文尝试基于贝叶斯理论开展二维大地电磁电阻率尖锐边界反演研究.在反演中,模型参数由边界位置及内部电阻率组成,通过贝叶斯理论将模型参数与数据相联系,采用Markov Chain Monte Carlo（MCMC）的Metropolis-Hastings（MH）方法对后验概率密度函数（Posteriori Probability Density,PDD）进行采样.采样过程中无罚值函数约束,完全以数据自身所包含的信息对模型进行约束,同时与有限约束进行比较,并考虑不同起始采样点对结果的影响.以接受率为参考,用模型算例说明MH方法中建议分布函数选择的重要性.当模型参数间相关性较弱时,使用边缘概率分布对采样结果进行分析.该方法能给出模型参数的分布范围,并给出该模型参数范围对应的数据范围.通过与已知模型的对比及数据拟合情况分析检验了该反演方法的有效性.该方法有助于提高大地电磁尖锐边界反演的分辨能力. 相似文献
15.
Two-dimensional joint inversion of radiomagnetotelluric and direct current resistivity data 总被引：2，自引：0，他引：2
An algorithm for the two-dimensional (2D) joint inversion of radiomagnetotelluric and direct current resistivity data was developed. This algorithm can be used for the 2D inversion of apparent resistivity data sets collected by multi-electrode direct current resistivity systems for various classical electrode arrays (Wenner, Schlumberger, dipole-diplole, pole-dipole) and radiomagnetotelluric measurements jointly. We use a finite difference technique to solve the Helmoltz and Poisson equations for radiomagnetotelluric and direct current resistivity methods respectively. A regularized inversion with a smoothness constrained stabilizer was employed to invert both data sets. The radiomagnetotelluric method is not particularly sensitive when attempting to resolve near-surface resistivity blocks because it uses a limited range of frequencies. On the other hand, the direct current resistivity method can resolve these near-surface blocks with relatively greater accuracy. Initially, individual and joint inversions of synthetic radiomagnetotelluric and direct current resistivity data were compared and we demonstrated that the joint inversion result based on this synthetic data simulates the real model more accurately than the inversion results of each individual method. The developed 2D joint inversion algorithm was also applied on a field data set observed across an active fault located close to the city of Kerpen in Germany. The location and depth of this fault were successfully determined by the 2D joint inversion of the radiomagnetotelluric and direct current resistivity data. This inversion result from the field data further validated the synthetic data inversion results. 相似文献
16.
17.
常规三维大地电磁反演的正则项为L2范数,它以电阻率空间分布函数处处光滑为模型期望,弱化了算法对电性突变界面的分辨能力.本文实现了正则项为L1范数的三维大地电磁反演算法,让模型空间梯度向量更有机会取得稀疏解,在充分正则的迭代下能够有效突出模型真实电性界面.为避免L1范数零点不可导带来的求解困难,使用迭代重加权最小二乘法把原问题转换为一系列L2正则子问题迭代求解.每个子问题的极小方法使用改进型拟牛顿法,其下降方向既能保证正则项海塞矩阵的精确性,又能允许反演过程随迭代灵活更新正则因子.使用比值法或分段衰减法自适应更新正则因子以避免迭代早期陷入奇异解,从而提升反演收敛的稳定性并降低初始模型依赖度.合成的无噪数据反演表明L1正则算法的模型恢复效果优于L2正则;不同噪声水平的合成数据反演表明本文的算法具有稳健性;实测数据反演对比表明在合理的正则因子调整策略下,L1正则反演结果的模型分辨率优于L2正则.另外,不同初始模型的反演测试还表明,正则因子选取不合理时L1正则可能造成方块状假异常. 相似文献
18.
Yanling Shi Zhanxiang He Zuzhi Hu Qiang Wei Dechun Li Cuixian Meng Liansheng Ji Sheng Zhang 《Studia Geophysica et Geodaetica》2018,62(1):167-185
Mapping deep geological hydrocarbon targets is of significant importance in basin exploration. In areas lacking reliable seismic data, magnetotelluric (MT) and gravity explorations are helpful to delineate the distribution of potential deep geological hydrocarbon targets. Here we investigate the effectiveness of the integrated 3D MT and gravity explorations for mapping the potential deep hydrocarbon source rocks. The result based on the data from the W Basin (part of the Ordes Basin) of China demonstrates that the method is efficient and economical for basin exploration. The method is particularly useful in target areas which are of great interest for oil and gas exploration but lack high quality seismic data. In our method, we first use the high-precision 3D small-bin MT data acquisition to improve the data accuracy. Then we perform datum static correction method and apply 3D inversion to obtain the3D resistivity distribution. We also develop a layered resistivity model based on resistivity logging to assist the interpretation of the inverted 3D resistivity data so as to derive an initial 3D geological model. Starting from the initial model, we use 2D gravity data to update the model via 2D inversion line by line, and then pass the updated model for the next round of the 3D MT inversion. The integrated inversion is implemented iteratively so the model converges to satisfy the need of final geological analysis. The application to the W Basin shows that we could successfully delineate the geological distribution of the potential deep hydrocarbon source rocks within the basin and map the thickness of the upper Paleozoic. 相似文献
19.
Nuree Han Myung Jin Nam Hee Joon Kim Tae Jong Lee Yoonho Song Jung Hee Suh 《Journal of Applied Geophysics》2009,68(4):533-545
Magnetotelluric (MT) surveys were conducted in Pohang, Korea, for low-temperature geothermal exploration in 2002 and 2003. Pohang is located in the southeastern part of the Korean Peninsula and close to the East Sea. In the interpretation of MT data from a coastal environment, sea effects must be correctly included because seawater is a strong conductor. We first constructed a five-layered earth model with a realistic coastline and bathymetry to investigate sea effects on MT data measured in Pohang. This model clearly shows that the Pohang data are significantly influenced by sea water at frequencies blow 1 Hz at the whole measurement sites. Next, we utilized a three-dimensional inversion algorithm based on the Gauss–Newton approach to produce a reliable resistivity model. Seawater is excluded from the inversion domain to fix the resistivity, while included in the modeling domain to simulate sea effects on MT responses. Blocks for the sub-seafloor are included in unknown parameters since they are sufficiently close to the survey area to affect MT responses in Pohang. Static shifts are also considered in inversion for more accurate interpretation. The rms data misfit is smoothly reduced from 11.2 to 1.87 after 7 iterations. The resulting resistivity model shows a pattern of low–high–low resistivity with depth. The model is compatible with resistivity logs obtained from four boreholes in the survey area, and can explain major geological features in Pohang. 相似文献
20.
Si-Yu Li Shu-Ming Wang Peng-Fei Wang Xiao-Lu Su Xin-Song Zhang Zhi-Hui Dong 《Acta Geophysica》2018,66(4):607-621
The grey wolf optimizer (GWO) is a novel bionics algorithm inspired by the social rank and prey-seeking behaviors of grey wolves. The GWO algorithm is easy to implement because of its basic concept, simple formula, and small number of parameters. This paper develops a GWO algorithm with a nonlinear convergence factor and an adaptive location updating strategy and applies this improved grey wolf optimizer (improved grey wolf optimizer, IGWO) algorithm to geophysical inversion problems using magnetotelluric (MT), DC resistivity and induced polarization (IP) methods. Numerical tests in MATLAB 2010b for the forward modeling data and the observed data show that the IGWO algorithm can find the global minimum and rarely sinks to the local minima. For further study, inverted results using the IGWO are contrasted with particle swarm optimization (PSO) and the simulated annealing (SA) algorithm. The outcomes of the comparison reveal that the IGWO and PSO similarly perform better in counterpoising exploration and exploitation with a given number of iterations than the SA. 相似文献