共查询到20条相似文献，搜索用时 46 毫秒
1.
2.
Zhdanov提出的拟线性(QL)近似方法有效解决了积分方程计算速度慢、占用内存多的缺陷,但因为三维频谱激电(SIP)换源频繁,每次换源都需要重新计算剖分单元的并矢格林函数和一次场,所以用该方法模拟三维SIP的速度仍然缓慢.本文根据一次场及格林函数的空间对称性,提出了一种适用于多源电磁法的快速QL近似正演方法.在此基础上,实现了Cole-Cole参数范围约束的三维SIP共轭梯度反演方法.理论试算结果表明:快速QL近似方法适用于三维SIP正演模拟,计算速度较QL近似方法有了很大提高;Cole-Cole参数范围约束的三维SIP共轭梯度反演方法对零频电阻率、极化率等参数的反演效果良好,将异常区域剖分成2000块时迭代一次只需约0.2s的时间,占用内存大约为140 MB. 相似文献
3.
为推进大地电磁三维反演的实用化，本文实现了基于L-BFGS算法的带地形大地电磁三维反演.首先推导了大地电磁法三维反演的Tikhonov正则化目标函数以及Hessian矩阵逆矩阵近似表达式和计算方法，然后设计了一种既能保证空气电阻率固定不变又能保证模型平滑约束的协方差矩阵统一表达式，解决带地形反演问题.在反演算法中采用正则化因子冷却法以及基于Wolf条件的步长搜索策略，提升了反演的稳定性.利用开发的算法对多个带地形地电模型（山峰地形下的单个异常模型、峰-谷地形下的棋盘模型）的合成数据进行了三维反演，并与已有大地电磁三维反演程序（ModEM）进行对比，验证了本文开发的三维反演算法的正确性和可靠性.最后，利用该算法反演了华南某山区大地电磁实测数据，得到该区三维电性结构，揭示了研究区以高阻介质为基底，中间以低阻不整合面和相对低阻介质连续分布，浅部覆盖高阻介质的电性结构特征，进一步验证了本文算法的实用性. 相似文献
4.
对地震波偏移成像而言,由于观测采集系统不规则、不完备采样以及地下复杂构造对波场传播的影响,在成像点处通常存在照明不均匀的现象,对偏移成像结果进行照明补偿是提高地质体成像精度的有效方法.本文给出了一种基于逆时偏移实现的高效照明补偿成像方法.在高频近似下,本文首先定义了衡量地下照明均匀程度的雅克比矩阵,并以此为基础给出照明补偿算子.通过波场的边界积分表达式,进一步将照明补偿算子利用外推波场进行表达,避免了照明补偿算子中高频格林函数的显式计算,以此构造出高效的照明补偿成像算法.合成数据和实际数据的数值实验结果验证了本文提出方法的正确性. 相似文献
5.
层析反演是速度建模中最重要的方法之一,结合偏移成像在成像域进行走时层析速度反演是当前比较成熟有效且广泛应用的技术.本文从高斯束偏移成像条件出发,在波动方程的一阶Born近似和Rytov近似下,推导了成像域走时扰动与速度扰动的线性关系,建立了成像域走时层析方程及其显式表达的层析核函数.该核函数的本质是有限频层析核函数,利用该核函数替换常规射线层析核函数可以明显提高层析反演精度.该核函数的计算关键是背景波场格林函数的计算,本文利用高斯束传播算子计算格林函数进而得到走时层析核函数,实现方式灵活高效且计算精度较高.基于高斯束传播算子的偏移成像与层析成像相结合进行深度域建模迭代,体现了速度建模与偏移成像一体化的思想.数值计算及实际数据应用证明了基于高斯束传播算子的成像域走时层析方法的有效性. 相似文献
6.
时间域航空电磁快速成像研究 总被引：4，自引：0，他引：4
基于电磁法理论,编写了时间域中心回线航空电磁正演算法程序,并通过查表法求取了视电阻率,时间域航空电磁的快速成像采用了经本文改进的快速一阶与二阶近似成像算法,通过四种典型地电模型对两种成像方法进行了对比,分析了两种算法在航空物探应用中的优缺点.试算结果表明,本文的两种快速成像算法均得到了比较理想的拟合效果,二阶近似成像算法对地层电阻率变化灵敏性尤为突出. 相似文献
7.
提出了针对轴对称二维井间模型的一维、二维联合反演方法.该方法基于纵向成层背景地层的Green函数，以一维反演结果作为联合反演的迭代初始值，同时反演一维背景地层和二维异常剖面的电导率.采用递推技术计算Green函数的系数，可以很方便地同时得到Green函数对空间坐标的积分和对背景层电导率的微分，大大加快了计算雅可比矩阵〖WTHX〗M〖WTBZ〗的速度并使计算方便、准确.反演实例说明了该方法的有效性.反演实例还显示，背景层电导率较二维成像剖面的电导率收敛速度要快，测量数据的精度对二维剖面电导率成像质量的影响大于对背景层电导率分布的影响. 相似文献
8.
三维反演解释是电磁法勘探发展的重要趋势,而如何提高三维反演的可靠性、稳定性和计算效率是算法开发者们目前的研究重点.本文实现了一种频率域可控源电磁(CSEM)三维反演算法.其中正演基于拟态有限体积法离散化,利用直接矩阵分解技术来求解大型线性系统方程,不仅准确、稳定,而且特别有利于含有大量发射场源位置的CSEM勘探情况;对目标函数的最优化采用高斯牛顿法(GN),具有近似二次的收敛性;使用预条件共轭梯度法(PCG)求解每次GN迭代所得到的法方程,避免了显式求解和存储灵敏度矩阵,减小了计算量.以上这些方法的结合应用,使得本文的三维反演算法准确、稳定且高效.通过陆地和海洋CSEM勘探场景中的典型理论模型的反演测试,验证了本文算法的有效性. 相似文献
9.
利用积分方程方法以及半空间并矢格林函数的快速算法对埋地目标体矢量电磁散射进行正演计算。首先,利用半空间电磁并矢格林函数建立起埋地目标体的体积分方程。然后通过将空间偏导转移至格林函数谱域积分的积分号之外,并采用离散复镜像方法来近似余下的零阶索末菲积分,进一步得到并矢格林函数各个分量的闭合形式。由于避免了对索末菲积分的繁琐数值计算,使得生成反应矩阵和计算散射场时由半空间并矢格林函数计算带来的瓶颈问题得到较好的克服,因而极大地提高了埋地目标体电磁响应正演计算效率,同时也能保证足够的精度。 相似文献
10.
利用加速差分进化算法反演非均匀介质电磁成像 总被引：1，自引：1，他引：0
本文研究了差分进化算法在地球物理反演中的几种应用.利用双频电磁波电导率成像原理建立成像方程后,根据其严重病态性质,将Tikhonov正则化方法与差分进化算法结合,反演其成像方程.为加速差分进化算法的收敛速度,提出了将种群熵的自适应差分进化(ARDE)算法以及粒子群差分进化混合(PSODE)算法分别与Tikhonov正则化方法结合.在大型反演计算中,这两种方法可以在不影响反演效果的前提下,不同程度地提高收敛速度,降低时间成本.适宜于在正则化参数选取困难情况时的地球物理反演问题的求解. 相似文献
11.
实现大地电磁法快速三维反演的关键在于找到快速计算灵敏度矩阵的方案. 本文在对大地电磁三维张量阻抗表达式进行深入分析的基础上，获得了三维快速松弛反演算法的灵敏度函数表达式，解决了三维快速松弛反演的核心问题；为了减小反演解的多解性，针对三维问题定义了最小构造函数，实现了求最小构造的三维快速松弛反演算法. 对二维棱柱体、三维棱柱体理论模型的大地电磁合成数据进行了反演试算，反演结果与理论模型相吻合. 对日本Kayabe地区实测资料进行了XY模式反演，反演得到的地电模型较好地反映了地热资源的分布. 反演试算结果表明大地电磁三维快速松弛反演具有计算速度快、结果稳定可靠等特点；使三维快速松弛反演在普通微机上得以实现. 相似文献
12.
13.
地层宏观电性各向异性会对可控源电磁响应产生重要影响.由于海底地层电性结构常表现为电导率各向异性，若仅对海洋可控源电磁（MCSEM）数据进行常规各向同性反演，有可能无法获得准确的反演解释结果，从而削弱MCSEM技术的可靠性.本文实现了电导率垂直各向异性（VTI）条件下频率域海洋可控源电磁数据三维反演算法.其中，三维正演采用基于二次场控制方程的交错网格有限体积法，并利用直接矩阵分解技术来求解离散所得的大型线性方程组，有利于快速计算多场源的响应.反演采用具有近似二次收敛性的高斯牛顿算法对目标函数进行最优化.最后，对具有VTI电性各向异性特征的盐丘构造模型的MCSEM合成数据分别进行了电导率各向同性和垂直各向异性三维反演，结果表明：各向同性三维反演算法无法对受VTI介质影响的MCSEM数据进行正确的反演解释，而垂直各向异性三维反演能够获得更为可靠的地下电阻率结构和异常体分布，展现出对海底电性各向异性结构更为优良的反演解释能力. 相似文献
14.
过套管电阻率测井是通过测量套管壁电势实现测量地层的视电阻率,基于传输线方程理论,针对层状地层,给出了套管壁电势、电流对地层横向电阻导数的微分方程(称Jacobi矩阵微分方程)及边界条件;利用Jacobi矩阵微分方程边值问题导出了过套管电阻率测井反演地层参数的Jacobi矩阵系数的解析表示,利用Marquardt方法实现了过套管测井的地层电阻率反演;通过计算对Jacobi矩阵的特性进行了探讨,并获得了较快的计算速度(因为Jacobi矩阵是用解析解表示的),反演结果与地层模型取得了较好的逼近.本文实现了过套管电阻率测井地层参数的Jacobi系数矩阵的快速计算及地层电阻率反演,为进一步开展电阻率测井数据处理提供了理论依据和快速反演算法. 相似文献
15.
利用三维感应测井交叉分量之差响应曲线在边界附近的响应行为,建立了求解大斜度井中井斜角快速反演方法.响应曲线随真实垂直深度(TVD)变化的对数衰减率与随测量距离(MD)变化的对数衰减率与之比是井斜角的余弦函数,通过简单三角函数计算可快速反演井斜角.响应曲线形态规则,在层状各向异性地层中与垂直电导率无关,主要反映水平电导率反差造成的地层边界信息.因此本文的井斜角快速反演方法既适用于各向同性地层也适用于各向异性地层.我们考察了测井环境及仪器参数对响应曲线衰减率的影响以及衰减率估算偏差对井斜角反演精度的影响.数值模拟结果验证了该反演方法在大斜度井中的有效性和实用性. 相似文献
16.
可控源音频大地电磁数据的反演方法 总被引：18，自引：7，他引：18
从反演方程、构造目标函数和求解三方面对用于可控源音频大地电磁法(CSAMT)的实用反演方法中的四种进行了描述．水平层状地层CSAMT法资料的直接反演法首次尝试了一维空间的全资料CSAMT反演，效果较好，但该方法尚难应用于2D、3D复杂介质中；奥克姆反演方法既考虑了横向的光滑函数，又考虑了纵向的光滑函数，得到比较光滑的横向、纵向变化的背景电性结果，但有可能把一些小构造光滑掉．快速松驰反演算法和共轭梯度算法由于计算速度快，占内存少而被用于三维反演中，二者相比，快速松驰算法在求解雅可比矩阵时只做一次正演计算，在更新模型时解小型方程组，所以在速度上更胜一筹．在后三种算法中，由于复杂电性结构无解析解，正演计算都采用数值计算．数值计算的可靠性、速度影响着反演算法的有效性，这方面的研究也将是2D、3D复杂电性结构反演的研究方向之一. 相似文献
17.
声反射成像测井中常用的基于射线理论的Kirchhoff积分偏移算法和基于单程波理论的F-K偏移算法均可实现井旁缝洞反射体的快速偏移成像,但其仅适用于地层垂向变化较弱的速度场和高陡角裂缝的偏移成像,无法实现低角度反射体的准确偏移归位,产生偏移假象误导测井解释.逆时偏移基于全波动方程,可适应强垂向变化速度场,实现近似水平反射体的偏移成像.本文详细分析了将逆时偏移应用于声反射成像测井时存在的数据准备、时间采样间隔匹配和成像条件改进等若干问题,通过设置多组理论模型来说明算法对井旁不同反射体的识别能力.模拟资料和实际资料处理结果证实,较F-K偏移算法,逆时偏移算法成像精度更高、收敛性更好,可有效实现近似水平构造偏移归位.改进的归一化互相关成像条件可解决深部地层的远井壁成像衰减问题,降低测井解释的多解性.逆时偏移将成为声反射成像测井高精度偏移技术的发展方向. 相似文献
18.
水平层状介质中双侧向资料的全参数正则化迭代反演与应用 总被引：3，自引：1，他引：3
水平层状介质中的模型向量通常由地层原状电阻率、侵人带电阻率、侵入半径以及地层界面深度等参数组成,整个测井记录实质上是模型向量的非线性函数.本文将研究双侧向仪器Born近似表达式并建立用双侧向测井资料同时反演所有模型的全参数正则化迭代算法(RWPII).双侧向电极系上的电流和电位必须满足特定的聚焦条件,其正演模拟与Born近似表达式与非聚集型仪器完全不同.应用摄动原理,首先研究双侧向仪器响应的Born近似,给出电极系上的电流密度、电流和电位微小变化与电导率摄动之间的积分关系.然后,根据层状地层中电阻率分布特征以及半解析正演模拟技术,建立Fréchet导数矩阵的快速算法,该Fréchet导数矩阵包含了视电阻率对所有模型参数的偏导数.在此基础上,研究同时反演所有模型参数的正则化迭代算法,并给出正则化因子的具体计算方法.最后,应用理论模型和实际资料反演结果验证该反演算法的有效性. 相似文献
19.
套管外地层受异常地应力、油气开采的影响,在径向上表现出非均质性;采用声波测井可以对该非均质性进行探测,利用偶极子横波测井数据可以对横波速度的径向分布进行反演.本文建立了套管井外地层横波速度径向分层参考模型,采用修正的微扰法计算了该模型的偶极弯曲波频散曲线,建立了横波速度径向分布反演目标函数,采用高斯牛顿法和快速模拟退火法对目标函数进行了求解,得到了套管井外地层横波速度的径向分布.分析了偶极弯曲波频段、套管横波速度对反演结果的影响,对比了高斯牛顿法和快速模拟退火法对反演过程的影响.分析对比结果表明,采用偶极弯曲波激发强度较高的频段与采用全频段的反演结果相近;套管的横波速度准确度越高,反演结果越准确;高斯牛顿法和快速模拟退火法计算精度相同,都可以得到高精度的横波速度径向分布;快速模拟退火法的计算效率略低于高斯牛顿法,但其收敛性对初始值依赖更小,实际处理中应选择快速模拟退火法. 相似文献
20.
可控源音频大地电磁法在资源勘探等领域中发挥着重要的作用.我们把有限差分数值模拟方法用于可控源音频大地电磁三维正演,结合正则化反演方案和共轭梯度反演的思路,将反演中的雅可比矩阵计算问题转为求解两次"拟正演"问题,得到模型参数的更新步长,形成反演迭代,实现了可控源音频大地电磁三维共轭梯度反演算法.该反演算法可用于对有限长度电偶源激发下采集到的可控源音频大地电磁全区(近区、过渡区和远区)视电阻率和相位资料进行三维反演定量解释,获得地下三维模型的电阻率结构.理论模型合成数据的反演算例验证了所实现的可控源音频大地电磁三维共轭梯度反演算法的有效性和稳定性. 相似文献