共查询到18条相似文献,搜索用时 156 毫秒
1.
2.
3.
提出了一种新的二维电阻率层析成像方案。与国内外已有成像方法相比,该方法无需计算雅可比矩阵,不必要作线性化近似,这不仅节省了计算时间,而且使成像迭代过程,稳定,可以对实际上广泛存在的大一对比度介质获得满意的成像效果,因而大大提高了电阻率层析成像技术的这用性。中给出的数值模拟结果显示了该方法的有效性。 相似文献
4.
对三维电阻率反演问题进行了深入研究,提供了一种利用地表观测数据实现三维反演的实用算法.该方法应用有限差分求正演解,并通过对粗糙度矩阵元素进行适当改进,使之适用于各种情况下粗糙度矩阵的求取,进而建立在模型的总粗糙度极小条件下的反演方程.对反演方程采用收敛速度快且稳定的最小二乘正交分解(LSQR)法进行迭代求解,在迭代求解过程中只需利用偏导数矩阵和其转置矩阵乘以一个向量的结果,回避了直接求偏导数矩阵的繁琐计算,节省了内存,加快了反演的计算速度.不同的计算实例表明上述方法是求解大规模三维电阻率反演问题的有效方法. 相似文献
5.
对于时间域航空电磁法二维和三维反演来说,最大的困难在于有效的算法和大的计算量需求.本文利用非线性共轭梯度法实现了时间域航空电磁法2.5维反演方法,着重解决了迭代反演过程中灵敏度矩阵计算、最佳迭代步长计算、初始模型选取等问题.在正演计算中,我们采用有限元法求解拉式傅氏域中的电磁场偏微分方程,再通过逆拉氏和逆傅氏变换高精度数值算法得到时间域电磁响应.在灵敏度矩阵计算中,采用了基于拉式傅氏双变换的伴随方程法,时间消耗只需计算两次正演,从而节约了大量计算时间.对于最佳步长计算,二次插值向后追踪法能够保证反演迭代的稳定性.设计两个理论模型,检验反演算法的有效性,并讨论了选择不同初始模型对反演结果的影响.模型算例表明:非线性共轭梯度方法应用于时间域航空电磁2.5维反演中稳定可靠,反演结果能够有效地反映地下真实电性结构.当选择的初始模型电阻率值与真实背景电阻率值接近时,能得到较好的反演结果,当初始模型电阻率远大于或远小于真实背景电阻率值时反演效果就会变差. 相似文献
6.
最近开发了一种针对二维大地电磁野外数据进行处理解释的新反演方法。该方法以加入阻抗相位信息的一维大地电磁连续介质曲线对比法为基础 ,把一维反演得到的电阻率和相位的数据集作为二维反演的初始模型 ,使用二维有限单元法做正演模拟。在程序的后继迭代中 ,深度方向上用一维反演修改模型的电阻率和深度值 ,沿测线方向由二维有限元作修改 ,反演结束可得到一个接近真实电性分布的电阻率数据集 ,并绘制成电阻率断面图。对模型的反演实验结果显示 ,该反演方法能够较真实地反映地下电性分布 ,而且避免了偏导数矩阵的计算 ,其原理简单 ,计算速度快 ,表明该反演方案是可行的 相似文献
7.
8.
提出了一种计算体积分方程的方法-改进型局域非线性迭代(MLNI),并用其对三维井间电磁场进行反演计算,该方法将井间大尺度散射体分为近场和远场区域两部分,它们的位置和尺寸均随场点位置的变化而改变,采用局域非线性近似计算近场区域的影响,将远场区域的影响作为外部激励源,采用迭代方法计算,该方法具有计算速度快,所需内存量少,收敛性好的优点,在反演中采用基于MLNI的非线性反演方法,利用这种方法可以将成像区域集中于一定范围内而不是整个三维空间,该反演方法由于考虑到了雅可比矩阵元素的非线性项,因而与传统的Born迭代反演方法相比更为精确和稳定,考虑到信息量和计算机内存的限制,第1次成像采用双重体元分割法进行较粗略的成像,然后缩小成像范围进行第2次较精确的成像,数值计算结果表明,MLNI是一种有效的计算井间大尺度异常体散射场的方法,将该方法用于反演过程能够得到较高分辨率的三维井间电导率图像。 相似文献
9.
过套管电阻率测井是通过测量套管壁电势实现测量地层的视电阻率,基于传输线方程理论,针对层状地层,给出了套管壁电势、电流对地层横向电阻导数的微分方程(称Jacobi矩阵微分方程)及边界条件;利用Jacobi矩阵微分方程边值问题导出了过套管电阻率测井反演地层参数的Jacobi矩阵系数的解析表示,利用Marquardt方法实现了过套管测井的地层电阻率反演;通过计算对Jacobi矩阵的特性进行了探讨,并获得了较快的计算速度(因为Jacobi矩阵是用解析解表示的),反演结果与地层模型取得了较好的逼近.本文实现了过套管电阻率测井地层参数的Jacobi系数矩阵的快速计算及地层电阻率反演,为进一步开展电阻率测井数据处理提供了理论依据和快速反演算法. 相似文献
10.
11.
A fast inversion technique for the interpretation of data from resistivity tomography surveys has been developed for operation on a microcomputer. This technique is based on the smoothness-constrained least-squares method and it produces a two-dimensional subsurface model from the apparent resistivity pseudosection. In the first iteration, a homogeneous earth model is used as the starting model for which the apparent resistivity partial derivative values can be calculated analytically. For subsequent iterations, a quasi-Newton method is used to estimate the partial derivatives which reduces the computer time and memory space required by about eight and twelve times, respectively, compared to the conventional least-squares method. Tests with a variety of computer models and data from field surveys show that this technique is insensitive to random noise and converges rapidly. This technique takes about one minute to invert a single data set on an 80486DX microcomputer. 相似文献
12.
13.
Resistivity anomaly imaging by probability tomography 总被引:10,自引:0,他引:10
Probability tomography is a new concept reflecting the inherently uncertain nature of any geophysical interpretation. The rationale of the new procedure is based on the fact that a measurable anomalous field, representing the response of a buried feature to a physical stimulation, can be approximated by a set of partial anomaly source contributions. These may be given a multiplicity of configurations to generate cumulative responses, which are all compatible with the observed data within the accuracy of measurement. The purpose of the new imaging procedure is the design of an occurrence probability space of elementary anomaly sources, located anywhere inside an explored underground volume. In geoelectrics, the decomposition is made within a regular resistivity lattice, using the Frechet derivatives of the electric potential weighted by resistivity difference coefficients. The typical tomography is a diffuse image of the resistivity difference probability pattern, that is quite different from the usual modelled geometry derived from standard inversion. 相似文献
14.
The boundary-element method is used to model the 2D terrain effect on the magnetotelluric (MT) field. Firstly, the boundary-value problem of a 2D magnetotelluric field is transformed into an integral equation problem by using Green's theorem. Then the boundary-element method is used to solve the integral equation and to obtain the MT field and its normal derivative on the terrain. From these values, the apparent resistivity can be calculated. Compared with the finite-element method, the boundary-element method is simpler in element division and the initial data preparation. The configuration of a terrain divided by the boundary-element method is more consistent with the practical terrain. The method proposed in this paper can be run on a microcomputer, so that it can be used in the field. 相似文献
15.
本文用边界单元法模拟二维地形对大地电磁场的影响.首先用格林公式将二维大地电磁的边值问题转变成积分方程,然后用边界单元法解积分方程,得到地形上的大地电磁场和它的法向导数,由此可计算电阻率.与有限单元法相比,本方法剖分后的地形与实际地形的拟合程度高,向计算机输入原始数据的工作简单,可以在微机上计算地形引起的视电阻率.本方法适合于在野外生产现场进行大地电磁法的地形改正.计算表明,地形对H_x型波的影响比对E_x型波严重得多;随着周期的增长,地形对H_x型波的影响变得稳定. 相似文献
16.
地震波场反演的BG-逆散射方法 总被引:5,自引:0,他引:5
本文讨论利用三维反射地震数据进行波场反演的一种方法,旨在取得高分辨率的地球模型.这种方法用Backus-Gilbert的理论构造波动方程非线性反问题的逐次线性化迭代格式,用逆散射原理导出泛函的Frechet导数,并用最佳折衷准则求解线性化后的方程组.根据迭代过程中不断提高分辨率的思想和减少计算成本的原则,设计了可供实用的反演算法流程. 相似文献
17.
本文讨论利用三维反射地震数据进行波场反演的一种方法,旨在取得高分辨率的地球模型.这种方法用Backus-Gilbert的理论构造波动方程非线性反问题的逐次线性化迭代格式,用逆散射原理导出泛函的Frechet导数,并用最佳折衷准则求解线性化后的方程组.根据迭代过程中不断提高分辨率的思想和减少计算成本的原则,设计了可供实用的反演算法流程. 相似文献