首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 9 毫秒
1.
受计算机硬件水平限制,大地电磁(MT)三维反演难以在实际中推广应用,MT数据解释仍以二维反演为主.地质构造具有一定的二维性特征,但不同的地质构造走向往往不同、MT测线无法与每个构造的走向垂直,这使得实际MT数据并不满足严格的二维反演条件.因此,有必要开展大地电磁测深二维反演对"准二维"地质构造的适应性研究.本文设计三个理论地电模型,通过三维正演计算获得各测点大地电磁响应,以此模拟实际观测数据.对MT数据进行相位张量分析,结果表明设计的地电模型主要表现为二维性特征.利用非线性共轭梯度(NLCG)反演算法对理论MT数据进行二维反演研究,重点讨论了测线方向、电性主轴旋转策略及反演模式选择,对反演结果的影响.对比分析反演结果,得到如下认识:1.测线方向对反演结果影响较小;2.电性主轴旋转角度对反演结果影响较大;3.TE+TM联合模式及单独TM模式的反演效果较好.研究结果表明:当一条测线下方在横向与纵向上存在走向不同的多个异常体时,对整条剖面分测点、分频段进行电性主轴旋转,反演所得结果最可靠.  相似文献   

2.
应用大地电磁测深法,对秦皇岛-唐山沿海地区进行地热资源勘查.在该测区布设6条测线,获得103个测点的野外资料.采用先进的数据处理和反演方法,得到大地电磁测深剖面二维反演电阻率模型,根据当地区域地质资料,对地层、断裂构造进行推断解释,并结合华北地区地热资源类型及其与地质构造特征的关系,圈定该测区地热资源远景区.  相似文献   

3.
大地电磁的多尺度反演   总被引:17,自引:6,他引:11       下载免费PDF全文
对于迭代方式的参数化反演方法,如何使反演结果稳定地收敛到整体极小仍是目前大地电磁(MT)反演中急需解决的问题.本文利用小波变换理论中的多尺度分析方法将大地电磁反问题分解为依赖于尺度变量的反问题序列,然后按尺度从大到小的次序依次求解,求解过程中前一个尺度反问题的解作为下一个尺度反问题的初始模型,直到来出对应于尺度为0的原反问题的解为止.该方法称为多尺度反演方法.数值试验和实际资料的反演结果表明,该方法可有效改善传统广义逆反演方法易陷入局部极小的弊端.  相似文献   

4.
大地电磁野外实测数据月前大多为二维剖面数据.如何反演这些二维剖面数据获得较为接近实际地电情况的结果,是多数大地电磁工作者关心的问题.我们通过对理论模型的三维响应进行分析和对合成数据及实测资料的反演结果进行对比研究,讨论了利用三维反演的方法来获得大地电磁二维剖面附近三维电阻率结构的可行性.结果表明:可用三维反演的方法来解...  相似文献   

5.
稀疏测线大地电磁资料三维反演研究:合成算例(英文)   总被引:2,自引:1,他引:1  
受勘探成本和工区环境等因素的影响,当前大多数大地电磁实际工作采取布置稀疏测线采集数据和使用二维反演方法解释这些稀疏测线数据的方式。然而,二维反演方法在解释三维地电构造数据时存在局限性,有时甚至做出错误的地质解释。本文尝试了使用三维反演方法对大地电磁稀疏测线数据进行反演解释。使用大地电磁全信息资料三维共轭梯度反演程序对理论模型合成稀疏测线数据进行了三维反演。结果表明:这种反演方案是可行与有效的。同时,我们发现在不同数据的三维反演结果中,四个张量阻抗元素和两个倾子数据同时反演的结果相对更为准确,更接近理论模型。  相似文献   

6.
大地电磁的人工鱼群最优化约束反演   总被引:3,自引:2,他引:1       下载免费PDF全文
大地电磁的反演问题是非线性,如果采用线性反演方法容易陷入局部极小,使得反演结果非唯一性严重.本文将人工鱼群算法引入到地球物理反演之中,提出了非线性的大地电磁人工鱼群最优化反演.该方法不需要进行偏导数的求取,可以对反演的范围进行约束,以减小反演结果的非唯一性.同时我们对搜索步长进行了改进,给出适用于大地电磁反演的人工鱼群参数.大量的理论数据试算表明,人工鱼群反演算法能够较好地寻找到全局最优解.实测数据的处理结果表明,该方法可以用来处理实际资料,并且能够取得很好的应用效果.  相似文献   

7.
地下介质电性在很多地区呈现非一维性,大地电磁测深(MT)在实际勘探时常会遇到测量主轴(电极方向)和电性主轴(测线方向)不重合的情况,这时就有必要将测量主轴上的地电信息转换到电性主轴上来.本文以吉洪诺夫-卡尼尔模型为理论基础,结合湘、鄂西地区的实际大地电磁测深资料,通过理论推导,计算出经过电性主轴偏转后的视电阻率值.再经过静态校正、地形校正,最后反演得出二维电阻率剖面,解释结果符合研究区域的实际地质构造特征,表明采用电性主轴偏转计算处理大地电磁测深资料效果较好.  相似文献   

8.
长周期大地电磁(LMT)是基于常规MT理论发展起来的电磁测深技术,青藏高原东部岩石圈较厚、视电阻率较低,应用LMT方法能够弥补常规MT仪器对低频信号响应的不足,获得有用的深部信息。文中概述了LMT的数据采集、处理及反演解释,并介绍了将LMT应用于东喜马拉雅构造结及其周围地区完成的长周期观测剖面——下察隅—玉树LMT剖面的岩石圈结构探测实例。实际观测表明,在重叠频段内LMT结果与常规MT具有良好的一致性;LMT数据处理实现的阻抗和倾子估计,在深部信息上具有更高的分辨率。LMT剖面结果揭示了青藏高原东部普遍存在地下低阻体,为局部地区部分熔融和地下流体存在的可能性提供了依据。LMT测深较好地弥补了常规MT方法在巨厚岩石圈和遇到低阻体时探测深度的局限性。  相似文献   

9.
航空电磁探测数据量大,二维、三维反演算法复杂、计算速度慢,通常采用一维反演,利用层状模型拼接描述地下复杂结构,但航空电磁数据信噪比低,容易引起一维反演结果横向连续性差等问题.本文针对上述问题,基于一维反演算法,通过整合测线观测数据,建立了测线数据整体的目标函数,并根据Tikhonov正则化反演理论,引入包含空间粗糙度和先验信息的模型参数约束项,确定了拟二维整体反演的目标函数,推导了反演迭代方程组,利用超松弛共轭梯度算法,求得由于整条测线整体反演所致的大型稀疏矩阵的极小化解,实现了对整条测线数据同时反演的固定翼航空电磁数据的拟二维整体反演算法.在反演迭代过程中,正则化因子采用线性搜索自适应迭代的方法自动选取,提高了反演结果的稳定性.对比分析了仿真数据的一维反演与拟二维整体反演结果,得出拟二维整体反演算法横向连续性较好,对高导覆盖层下的导体分辨率优于一维反演,同时受高斯噪声的影响较小.最后,将直升机飞行实测噪声加入仿真数据中,拟二维整体反演结果平均相对误差较一维反演结果降低了31.6%,进一步验证了拟二维整体反演算法的有效性.  相似文献   

10.
大地电磁地震前兆初步分析系统   总被引:2,自引:0,他引:2  
介绍了大地电磁(MT)地震前兆初步分析系统的建立及主要功能,运用该系统对2003年度祁连山地震带8个大地电磁复测点进行了MT多参量综合对比、变化幅度的计算及异常前兆判定,取得较好效果。  相似文献   

11.
马尔科夫链蒙特卡洛方法(MCMC)是一种启发式的全局寻优算法,可以用来解决概率反演的问题.基于MCMC方法的反演不依赖于准确的初始模型,可以引入任意复杂的先验信息,通过对先验概率密度函数的采样来获得大量的后验概率分布样本,在寻找最优解的过程中可以跳出局部最优得到全局最优解.MCMC方法由于计算量巨大,应用难度较高,在地...  相似文献   

12.
地球物理抗差估计和广义逆方法   总被引:15,自引:3,他引:12       下载免费PDF全文
为进行观测数据误差较大的地球物理资料的反演,引人抗差估计,称作地球物理抗差估计;为使病态方程组的反演解更可靠,又结合广义逆方法进行了算法改进.首先介绍抗差最小二乘(RLS)的基本原理,然后推寻出适于广义逆反演方法的改进算式,最后举例加以讨论.分析表明,抗差估计可以有效地抑制地球物理观测异常的影响,得出正常模式下的较好估计值;而用改进后的算式和广义逆反演可以使反演解更加改善,不仅如此,改进后的方法还能直接给出解估计的可靠性评价.  相似文献   

13.
Fourier reconstruction with sparse inversion   总被引:2,自引:0,他引:2  
The problem of seismic data reconstruction is posed as an inverse problem where the objective is to obtain the Fourier coefficients that synthesize the signal. Once the coefficients have been found, they are used to reconstruct the data on a uniformly spaced grid. A non‐quadratic model weight function is included to stabilize the inversion and to provide the additional information required to interpolate through gaps. In the reconstruction of a non‐uniformly sampled trace, an image and a marine 3D VSP shot‐record, the method shows improved reconstruction in large gaps and is less sensitive to the spatial bandwidth used in the inversion compared to Fourier reconstruction without the non‐quadratic model weight function.  相似文献   

14.
非线性二次规划贝叶斯叠前反演   总被引:23,自引:11,他引:12       下载免费PDF全文
叠前反演的目的是基于弹性波理论从地震数据中获得地层参数的可靠估计,进而用于描述地层的流体和岩性特征.然而叠前反演问题都是高维的和非适定的,并且容易受各种噪声和采集过程中不确定因素的影响,因此,为了获得稳定可靠的解必需对反演过程加以合理的约束.本文提出了一种基于非线性二次规划的叠前三参数反演方法.首先基于贝叶斯参数估计理论,假设似然函数服从高斯分布,并使待反演的参数服从于改进的Cauchy分布,从而提高了反演结果的分辨率;其次用协方差矩阵来描述参数间的相关程度,进一步提高了反演结果的稳定性;最后将问题转化为一个非线性二次规划的求解问题,并在多种约束下得到问题的解.仿真实验和实际应用皆已表明,本文提出的反演方法运算速度快捷,既使在信噪比很低的情况下也可获得较好的反演结果,为储层的进一步识别提供更多的物性参数.  相似文献   

15.
Sound velocity inversion problem based on scattering theory is formulated in terms of a nonlinear integral equation associated with scattered field. Because of its nonlinearity, in practice, linearization algorisms (Born/single scattering approximation) are widely used to obtain an approximate inversion solution. However, the linearized strategy is not congruent with seismic wave propagation mechanics in strong perturbation (heterogeneous) medium. In order to partially dispense with the weak perturbation assumption of the Born approximation, we present a new approach from the following two steps: firstly, to handle the forward scattering by taking into account the second-order Born approximation, which is related to generalized Radon transform (GRT) about quadratic scattering potential; then to derive a nonlinear quadratic inversion formula by resorting to inverse GRT. In our formulation, there is a significant quadratic term regarding scattering potential, and it can provide an amplitude correction for inversion results beyond standard linear inversion. The numerical experiments demonstrate that the linear single scattering inversion is only good in amplitude for relative velocity perturbation ( \( \delta_{c}/c_{0} \) ) of background media up to 10 %, and its inversion errors are unacceptable for the perturbation beyond 10 %. In contrast, the quadratic inversion can give more accurate amplitude-preserved recovery for the perturbation up to 40 %. Our inversion scheme is able to manage double scattering effects by estimating a transmission factor from an integral over a small area, and therefore, only a small portion of computational time is added to the original linear migration/inversion process.  相似文献   

16.
水平井钻进过程中,地层边界的准确预测对地质导向决策至关重要.但水平井测井环境复杂,随钻方位电磁波测井响应不直观,储层性质、地层界面等信息的精确、定量评价困难.本文提出了一种适用于大斜度井/水平井随钻方位电磁波测井资料的电阻率及地层界面信息的实时提取方法.该方法主要基于以下几个技术:(1)引入滑动开窗策略将高维随钻方位电磁波测井资料的反演问题转化为多个连续窗口的一维反演问题;(2)对反演参数进行多初值初始化,结合正则化Levenberg-Marquardt算法保证解的全局最优性;(3)充分利用随钻方位电磁波测井多种探测深度信息,采用分级反演策略提高反演速度.数值和实测资料反演结果表明:降维反演方法可实现对地层界面的实时、准确提取,且适用于多数复杂地层结构和任意井眼轨迹;利用浅探测地质信号可实现随钻方位电磁波测井资料的分级反演,最终反演速度可提升10倍以上.  相似文献   

17.
基于逆算子估计的AVO反演方法研究   总被引:1,自引:1,他引:0       下载免费PDF全文
传统反演算法以优化算法为主,而基于逆算子估计的AVO反演算法则利用了直接求逆的思路.算法的关键在于寻找存在逆函数的子域,进而可以在子域内直接求逆,这种解决反问题的思路不同于一般的优化类算法所采用的直接搜索解的方式,具有更高的效率.AVO反演利用了振幅随着偏移距的变化特征,反演的精度受到地震资料质量的影响,通过加入L1范数约束以及合理的初始模型有助于提高反演的稳定性以及准确度.模型测算和实际应用表明,基于逆算子估计的AVO反演方法具有较高的精确程度和可靠性.  相似文献   

18.
核磁共振与瞬变电磁三维联合解释方法   总被引:3,自引:1,他引:2       下载免费PDF全文
传统核磁共振地下含水量解释多采用基于均匀半空间或层状导电模型的一维反演,分层给出地下含水信息.然而,这些方法忽略了地下复杂电阻率分布信息对结果的影响,也不能很好地反映局部三维含水构造.本文从三维电介质中核磁共振响应的正演理论出发,提出首先利用瞬变电磁数据进行基于等效导电平面法的快速电阻率成像,然后将成像结果作为核磁共振三维反演的电性模型,进行联合解释.激发磁场的分布采用有限元法直接求解,通过引入伪δ源实现电流源的加载,并强加散度条件排除了三维磁场模拟中"弱解"的影响.针对核磁共振灵敏度矩阵的病态性和数据中存在的干扰信号,提出考虑罚项的非线性拟合目标函数,利用线性化方法进行核磁共振反演.模型数据表明该方法能较准确反映地下三维含水构造,实测算例进一步证明了方法的有效性.本研究将促使核磁共振方法在岩溶、裂隙水、孤立水体等复杂水文地质条件及隧道、矿井灾害水源探测等方面得到有效应用.  相似文献   

19.
基于遗传算法的CSAMT最小构造反演   总被引:15,自引:4,他引:11       下载免费PDF全文
利用遗传算法进行不考虑近场校正的全场资料CSAMT反演研究.遗传算法属于全局最优化方法,具有对初始模型依赖小,不易陷入局部极值的优点,然而,当未知数较多时,多解性仍是该方法的瓶颈.为了减小多层反演的多解性,在反演中引入最小构造约束,针对CSAMT的遗传算法反演问题定义了最小构造目标函数,经过模型试验找到了其具体表达式,并找到了适合CSAMT资料反演的拉格朗日乘子的最佳取值μ=0.5,实现了基于遗传算法的CSAMT最小构造反演.利用H、A、K、Q和HKH、KHA模型对方法进行了数值试验,在无噪和加入10%噪声情况下,反演结果与模型一致;加入20%噪声后,反演仍取得良好结果,与理论模型基本吻合.将该方法用于水平层状地层和横向变化地层的实测资料反演,结果与地质资料吻合.不同的计算实例表明了该方法的有效性.  相似文献   

20.
Two seismic source studies usingteleseismic data are performed by Simulated Annealing(SA), a non-linear inversion method. The Very FastSimulated Annealing (VFSA) algorithm is used and onlyteleseismic data are inverted. We have designed a fastand efficient way of multiple direct problemevaluation, which is based on pre-calculating theelementary Green's function. During the process we setthe values of the inversion control parameters(temperature, number of iterations) and modified thecooling schedule. In the current version, theinversion seeks for the point source mechanism, thedepth of the source, the scalar moment and the sourcetime function (STF). The method is applied to twoearthquakes: 18 Nov. 1992 in Greece, M 5.9 and 14 Sep.1995 in Mexico, M7.3. The calculation is performed ona simple 1D model of the structure. For the firstearthquake the inversion recovered the solution fairlywell; for the second the solution was less acceptable.However, we do not consider this fact to be a failureof the method, but a consequence of an inadequatemodel of the source and of the medium structure. Acasual attempt of reliability determination was alsoperformed; the obtained values of errors arereasonable, except for a few cases when the methodfailed.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号