共查询到20条相似文献，搜索用时 813 毫秒
1.
对于迭代方式的参数化反演方法,如何使反演结果稳定地收敛到整体极小仍是目前大地电磁（MT）反演中急需解决的问题．本文利用小波变换理论中的多尺度分析方法将大地电磁反问题分解为依赖于尺度变量的反问题序列,然后按尺度从大到小的次序依次求解,求解过程中前一个尺度反问题的解作为下一个尺度反问题的初始模型,直到来出对应于尺度为0的原反问题的解为止．该方法称为多尺度反演方法．数值试验和实际资料的反演结果表明,该方法可有效改善传统广义逆反演方法易陷入局部极小的弊端． 相似文献
2.
研究地壳形变的理论和方法主要有两大类：地球物理模型方法和几何方法．由于基准误差、观测误差和模型误差等,两种方法得到的结果往往存在差异．为合理利用几何信息和物理信息,控制几何观测及物理模型误差对形变参数估计的影响,并平衡两类信息对形变参数的贡献,本文提出一种利用自适应滤波综合估计形变参数的方法．采用抗差等价权控制几何观测异常误差的影响,引入自适应因子平衡几何观测和地球物理模型信息对形变模型参数估计的贡献,利用高精度IGS站速度确定局部形变的基准．利用一实测GPS监测网进行计算,结果表明该混合估计策略可充分利用局部重复几何观测信息减弱地球物理模型信息带来的形变系统误差,提高了形变参数解算精度． 相似文献
3.
针对Mogi模型垂直位移与水平位移联合反演中的病态问题,改进火山形变总体最小二乘(Total Least Squares,TLS)联合反演的虚拟观测法,并使用方差分量估计(Variance Components Estimation,VCE)方法确定病态问题的正则化参数.将附有先验信息的参数作为观测方程,与垂直位移和水平位移的观测方程联合解算,推导了三类观测方程联合反演的求解公式及基于总体最小二乘方差分量估计确定正则化参数的表达式,给出了算法的迭代流程.通过算例实验,研究了总体最小二乘联合反演的虚拟观测法在火山Mogi模型形变反演中的应用;算例结果表明,三类数据的联合平差及方差分量估计方法可以确定权比因子并得到修正后的压力源参数,具有一定的实际参考价值. 相似文献
4.
5.
重复水准测量提供了点位高程变化的重要几何信息,但水准测量本身又受地球重力场的影响,在水准观测值中同时含有重力场的信息．重复绝对重力测量或相对重力测量提供了丰富的局部重力场的信息,由于点的重力值随高程变化而变化,因此,它们同时也含有点位高程变化的信息．经典平差方法将两类观测值分别解算,将互补的信息作为干扰加以改正,或作为噪声处理,这不仅丢失了有用的信息,而且有损于成果的精度．
本文在整体大地测量框架内构制了这两类观测值的联合解算模型,给出了解算方法,研究了形变场的拟合和预测问题,并用实例进行了验证,得到了一些有益的结果．研究表明:联合解算可同时给出地面点的高程及其变化率和局部重力场的参数及其变化率,有利于信息的提取;联合解算的精度优于单独解算的结果;秩亏抗差解和拟稳抗差解分别优于相应的非抗差解;拟稳抗差解算方法自动选择拟稳点,有效地抵制了粗差和形变异常对起算基准的干扰;对于多期复测的资料,使用动态拟稳抗差估计,不仅能动态地寻找可靠的拟稳点,而且能保证动态变化参数免受或少受各历元观测粗差的干扰,可明显地提高动态参数的验后精度;当点位丰富时,建立形变模型并据此内插未测点的变形和推测某个时刻地壳的形变状态也是可行的． 相似文献
6.
从地球物理反演的基本概念出发，认为地球物理反演是 对实测数据的地球物理属性的理解或解释，多方法的地球物理数据联合反演是一种多传感器 的数据融合. 本文分析了地球物理数据的模糊特性，采用基于语义的模糊化方法，使不同物 理意义和尺度的特征数据及测区的地质和地球物理背景成为一体，结合地球物理专家解释的 方法，利用基于模糊逻辑系统的神经网络实现了融合. 该方法充分利用了各种地球物理探测 数据的全部信息，避免了线性反演的复杂计算；其数据融合的观点,为解决地球物理联合反 演问题提供了新的思路. 通过模拟实验和应用实例验证，该方法是有效的. 相似文献
7.
8.
9.
用迭代模拟方法进行地球物理反演时必然涉及到用调整模型参数来拟合观测资料的问题。地震和位场模型响应均会受到调整岩石性质参数的影响。这种“联合反演”的目的是要得到一个模型,使它与有效的地表和钻孔地球物理资料相符。尽管地球物理资料的反演一 相似文献
10.
《世界地震译丛》2015,(3)
地震震源机制反演是少数确定地震发生深度主应力方向的可用手段之一。为进行这种反演已提出了多种方法。对于包括Jacques Angelier提出的一种方法在内的3种最常用的方法,本文讨论了相关的物理假设和误差确定,对其中的一种方法进行了改进并将这4种方法用于估计上莱茵地堑的应力场。地震数据采用锡伦茨1980年7月15日MW4.4地震后12小时布设的流动台网的记录。尽管根据不同方法确定的主应力方向的差别达28°(最小主应力的方向范围为从N051°E,俯角27°到N090°E,俯角为20°),但每种方法90%置信水平为11°~27°。不同的方法得到的表征主应力相对差的R值相当离散(自0.7~0.3,90%置信水平下的有0.2置信范围)。另外三种方法均给出了一些不好解释的震源机制。本文认为这是因为没有考虑到应力场的不均匀性。本文的结论是:如果不采用事件独立性假设,深部应力场的解算分辨率就较低。当采用该假设时,应力场的分辨率大约为15°。不同震源机制反演方法得到的最小主应力方向与距锡伦茨约20km的巴塞尔5km深处钻孔崩塌观测的方向(N054°E±14°)差4°~36°。与钻孔崩塌观测符合最好的解,是必要物理假设数目最少(3个)的解。 相似文献
11.
12.
为解决地球物理反演中多解性的问题,综合多种地球物理信息的联合反演受到了广泛的关注.本文依据不同地球物理响应可能由相同异常体引起,而不同地球物理分布参数之间存在相关性等特点,提出了一种基于局部Pearson相关系数约束的联合反演方法.该方法假设每个局部区域模型参数的分布具有线性相关特性,在拟合不同类型观测数据时,对局部模型参数施加相关性约束,进行联合反演以减少多解性.本文采用交替迭代联合反演流程,改善了同一目标函数下联合反演收敛性和速度问题.基于新的联合反演方法和流程,我们测试了三维大地电磁和重力仿真数据的联合反演.结果表明,本文提出的基于局部相关性约束的联合反演方法,能充分利用大地电磁和重力观测数据信息,有效改善单一地球物理反演收敛性和多解性的问题,反演效果得到明显提升. 相似文献
13.
Iteratively re‐weighted and refined least squares algorithm for robust inversion of geophysical data 
下载免费PDF全文
下载免费PDF全文 A robust metric of data misfit such as the ?1‐norm is required for geophysical parameter estimation when the data are contaminated by erratic noise. Recently, the iteratively re‐weighted and refined least‐squares algorithm was introduced for efficient solution of geophysical inverse problems in the presence of additive Gaussian noise in the data. We extend the algorithm in two practically important directions to make it applicable to data with non‐Gaussian noise and to make its regularisation parameter tuning more efficient and automatic. The regularisation parameter in iteratively reweighted and refined least‐squares algorithm varies with iteration, allowing the efficient solution of constrained problems. A technique is proposed based on the secant method for root finding to concentrate on finding a solution that satisfies the constraint, either fitting to a target misfit (if a bound on the noise is available) or having a target size (if a bound on the solution is available). This technique leads to an automatic update of the regularisation parameter at each and every iteration. We further propose a simple and efficient scheme that tunes the regularisation parameter without requiring target bounds. This is of great importance for the field data inversion where there is no information about the size of the noise and the solution. Numerical examples from non‐stationary seismic deconvolution and velocity‐stack inversion show that the proposed algorithm is efficient, stable, and robust and outperforms the conventional and state‐of‐the‐art methods. 相似文献
14.
微地震监测难以拾取准确初至,为了提高反演定位精度和减小多解性,研究了微地震贝叶斯差分进化反演方法.从分析讨论理论模型反演残差及其协方差分布特征出发,结合对比不加噪音和加入不同程度的噪音后残差协方差极小点位置移动、分布梯度变化特征,提出了先验信息解估计方法.针对后验估计中,由于难以获得先验信息解的方差估计致使无法计算加权系数问题,通过分析残差变化特征和解的变化关系,研究了利用残差求取加权系数的方法.为了加快寻优速度,讨论了差分进化反演方法,在变异操作方面使用差分策略,即利用种群中个体间的差分向量对个体进行扰动,实现个体变异,充分有效利用群体分布特性,提高算法的搜索能力,避免遗传算法中变异方式的不足.通过理论模型测试本方法的反演效果,并且和搜索方法反结果进行比较.测试结果证明本反演方法,对于不同程度初至干扰,反演结果向准确解逼近程度比搜索方法要好得多,实际资料的反演结果也好于搜索方法. 相似文献
15.
16.
地球物理反问题线性化处理之后, 各种反演算法归结为对病态线性方程组的求解. 为了快速准确地计算出地球物理参数, 本文提出了一种全新的基于LSQR算法的混合差分进化算法(Hybrid Differential Evolution Algorithm, HDE). 该算法利用LSQR算法给出DE算法的初始种群, 提高DE算法的计算速度和稳定性. 在不同噪声水平下, 对四种正则化方法Tikhonov、TSVD、LSQR和HDE的反演结果进行详细比较. 理论模型和实际数据反演的结果都表明: 改进的HDE算法应用于地球物理反问题的求解是成功的: 反演结果与原设定模型具有较高的相关性, 在稳定性和准确性上较常规的反演算法都具有一定的优势; 而且不需要给定正则化参数, 具有更强的实用性. 相似文献
17.
18.
2D多尺度混合优化地球物理反演方法及其应用(英文) 总被引：1，自引：0，他引：1
局部优化和全局优化方法广泛应用到地球物理反演,但是两者各有其优缺点。将两类方法结合起来可以取长补短。将退火遗传算法（SAGA）和单纯形算法相结合,得到了一种高效、健全的2D非线性混合地震走时反演方法。首先,利用SAGA进行大范围的全局搜索,然后由单纯形方法进行快速局部搜索。为了降低层析成像的多解性,我们采用了多尺度逐次逼近的技巧。把速度场划分为不同的空间尺度,定义网格节点上的速度作为待反演参数,采用双三次样条函数模型参数化,正问题采用有限差分走时计算方法,反问题采用多尺度混合反演方法。一个低速度异常体的数值模拟试验和抗走时扰动试验表明该方法是有效和健全的。我们将该方法应用到青藏高原东北缘阿尼玛卿rlet,Meyer,Marr,缝合带东段上部地壳速度结构研究中。数字模型试验和实际资料的应用表明了方法的有效性和健全性。 相似文献
19.
Raghu K. Chunduru Mrinal K. Sen Paul L. Stoffa R. Nagendra 《Geophysical Prospecting》1995,43(8):979-1003
The inversion of resistivity profiling data involves estimation of the spatial distribution of resistivities and thicknesses of rock layers from the apparent resistivity data values measured in the field as a function of electrode separation. The drawbacks of using traditional curve-matching techniques to solve this inverse problem have been overcome by iterative linear techniques but these require good starting models even if the shape of the causative body is asssumed known. In spite of the recent developments in inversion techniques, no robust method exists for the inversion of resistivity profiling data for the simple model of dikes and spheres which are the classical models of geophysical prospecting. We apply three different non-linear inversion schemes to invert synthetic resistivity profiling data for the classical models embedded in a uniform matrix of contrasting resistivity. The three non-linear algorithms used are called the Metropolis simulated annealing (SA), very fast simulated annealing (VFSA) and a genetic algorithm (GA). We compare the performance of the three algorithms using synthetic data for an outcropping vertical dike model. Although all three methods were successful in obtaining optimal solutions for arbitrary starting models, VFSA proved to be computationally the most efficient. 相似文献
20.
岩石物性约束是地球物理联合反演重要的实现方式.在以往联合反演研究中,这种约束通常以点对点的映射方式出现,要求经验性的岩石物性关联是相对精确的,无形中提高了岩石物性约束应用于联合反演的条件.基于此,本文结合全局寻优的模拟退火算法,提出了宽范围岩石物性约束技术.该技术将岩石物性约束与模拟退火模型扰动、模型接受准则相结合,实现了空间对空间的映射,将电阻率和速度耦合在一起.文章以MT和地震正则化同步联合反演为例,检验了不同的先验信息条件下的宽范围岩石物性约束技术的适用性.模型试验表明:新技术容易实现,容错性高,可降低"不精确"岩石物性关联信息引入带来的风险,提高岩石物性约束应用的灵活性. 相似文献