共查询到19条相似文献，搜索用时 171 毫秒
1.
2.
基于边界积分方程,提出了依据二极采样数据、但利用电位差的某种线性组合的"准三极"电阻率层析成像方法.研究表明,该方法不仅较好地解决了对高阻异常体的成像问题,而且使低阻异常体的成像质量获得了明显改善和提高,因此增强了电阻率层析成像技术的实用性. 相似文献
3.
4.
基于边界积分方程 ,提出了依据二极采样数据、但利用电位差的某种线性组合的“准三极”电阻率层析成像方法 .研究表明 ,该方法不仅较好地解决了对高阻异常体的成像问题 ,而且使低阻异常体的成像质量获得了明显改善和提高 ,因此增强了电阻率层析成像技术的实用性 . 相似文献
5.
用变分玻恩迭代方法重建二维非均匀介质结构 总被引：8，自引：1，他引：7
提出了用于二维轴对称非均匀介质结构的反演和成像的一种新的反演迭代方法──变分玻恩迭代方法（ＶＢＩＭ）．首先利用玻恩近似将非线性积分方程线性化,然后应用变分方法导出用于反演的电场积分方程．正演数据则利用高效的数值模式匹配方法获得．数值结果表明,ＶＢＩＭ与ＢＩＭ相比,其收敛速度、成像质量等均得到较大的改善。 相似文献
6.
基于非线性积分方程,利用井轴上阵列感应成像测井仪（AIT）的单频测量信号,将变形玻恩（Bom）迭代方法用于地球物理测井反演,用正则化方法求出地层电导率分布。文中首先给出对于未知电导率分布的非线性积分方程,然后用玻恩近似使非线性积分方程线性化,求出电导率分布。在迭代过程中,数值模式匹配法用于求解阵列感应的正演问题。数值计算表明,用井轴上 AIT 的单频测量信号得到的成像结果与真电导率剖面分布吻合较好。 相似文献
7.
系统地论述了用有限单元法研究复杂地形条件下三维直流电阻率的正演计算技术.首先给出了三维构造中点源电场的边值问题以及相应的变分问题;然后利用有限单元法求解变分问题,采用四面体单元对研究区域进行剖分,在单元中进行三线性函数插值,将变分方程化为线性代数方程组;最后,考虑到节约计算时间,利用对称超松弛顸条件共轭梯度迭代算法求解大型线性方程组,得到了各节点的电位值,进而计算出地表的视电阻率.通过理论模型的计算检验了算法的可行性之后,给出了几种常见纯地形异常的数值模拟结果和一个组合模型的计算结果,其研究工作为研究三维直流电阻率反演奠定了基础. 相似文献
8.
本文用边界单元法模拟二维地形对大地电磁场的影响.首先用格林公式将二维大地电磁的边值问题转变成积分方程,然后用边界单元法解积分方程,得到地形上的大地电磁场和它的法向导数,由此可计算电阻率.与有限单元法相比,本方法剖分后的地形与实际地形的拟合程度高,向计算机输入原始数据的工作简单,可以在微机上计算地形引起的视电阻率.本方法适合于在野外生产现场进行大地电磁法的地形改正.计算表明,地形对H_x型波的影响比对E_x型波严重得多;随着周期的增长,地形对H_x型波的影响变得稳定. 相似文献
9.
通过对波场的时间二阶积分运算以增强地震数据中的低频成分,提出了一种可有效减小对初始速度模型依赖性的地震数据全波形反演方法—时间二阶积分波场的全波形反演方法.根据散射理论中的散射波场传播方程,推导出时间二阶积分散射波场的传播方程,再利用一阶Born近似对时间二阶积分散射波场传播方程进行线性化.在时间二阶积分散射波场传播方程的基础上,利用散射波场反演地下散射源分布,再利用波场模拟的方法构建地下入射波场,然后根据时间二阶积分散射波场线性传播方程中散射波场与入射波场、速度扰动间的线性关系,应用类似偏移成像的公式得到速度扰动的估计,以此建立时间二阶积分波场的全波形迭代反演方法.最后把时间二阶积分波场的全波形反演结果作为常规全波形反演的初始模型可有效地减小地震波场全波形反演对初始模型的依赖性.应用于Marmousi模型的全频带合成数据和缺失4Hz以下频谱成分的缺低频合成数据验证所提出的全波形反演方法的正确性和有效性,数值试验显示缺失4Hz以下频谱成分数据的反演结果与全频带数据的反演结果没有明显差异. 相似文献
10.
提出了一种新的二维电阻率层析成像方案。与国内外已有成像方法相比,该方法无需计算雅可比矩阵,不必要作线性化近似,这不仅节省了计算时间,而且使成像迭代过程,稳定,可以对实际上广泛存在的大一对比度介质获得满意的成像效果,因而大大提高了电阻率层析成像技术的这用性。中给出的数值模拟结果显示了该方法的有效性。 相似文献
11.
以二维声波方程为模型,在时间域深入研究了全波形速度反演.全波形反演要解一个非线性的最小二乘问题,是一个极小化模拟数据与已知数据之间残量的过程.针对全波形反演易陷入局部极值的困难,本文提出了基于不同尺度的频率数据的"逐级反演"策略,即先基于低频尺度的波场信息进行反演,得出一个合理的初始模型,然后再利用其他不同尺度频率的波场进行反演,并且用前一尺度的迭代反演结果作为下一尺度反演的初始模型,这样逐级进行反演.文中详细阐述和推导了理论方法及公式,包括有限差分正演模拟、速度模型修正、梯度计算和算法描述,并以Marmousi复杂构造模型为例,进行了MPI并行全波形反演数值计算,得到了较好的反演结果,验证了方法的有效性和稳健性. 相似文献
12.
For wave propagation simulation in piecewise heterogeneous media, Gaussian-elimination-based full-waveform solutions to the generalized Lippmann–Schwinger integral equation (GLSIE) are highly accurate, but involved with extremely time-consuming computations because of the very large size of the resulting boundary–volume integral equation matrix to be inverted. Several flexible approximations to the GLSIE are scaled in an iterative way to adapt numerical solutions to the smoothness of heterogeneous media in terms of incident wavelengths, with a great saving of computing time and memory. Among various typical iterative schemes to the GLSIE matrix, the generalized minimal residual method (GMRES) is an efficient approach to reduce the computational intensity to some degree. The most efficient approximation can be obtained using a Born series, as an alternative iterative solution, to both the boundary-scattering and volume-scattering waves, leading to the Born-series approximation (BSA) scheme and the improved Born-series approximation (IBSA) scheme. These iteration schemes are validated by dimensionless frequency responses to a heterogeneous semicircular alluvial valley, and then applied to a heterogeneous multilayered model by calculating synthetic seismograms to evaluate approximation accuracies. Numerical experiments, compared with the full-waveform numerical solution, indicate that the convergence rates of these methods decrease gradually with increasing velocity perturbations. The comparison also shows that the BSA scheme has a faster convergence than the GMRES method for velocity perturbations less than 10 percent, but converges slowly and even hardly achieves convergence for velocity perturbations greater than 15 percent. The IBSA scheme gives a superior performance over the other methods, with the least iterations to achieve the necessary convergence. 相似文献
13.
GPS无线电掩星技术是崭新的、高效的地球大气层和电离层探测技术,但仍在发展和完善之中.本文详细推导了Abel积分和绝对TEC电离层反演方法,研究了如何解决Abel积分产生的上下限异常问题;用COSMIC发布的GPS原始数据进行了反演计算,将结果与地面电离层测高仪数据进行了比较,最后讨论了周跳对反演结果的影响问题.结果表明：（1）在较高轨道高度（约800 km）,Abel积分与绝对TEC方法的反演结果基本一致,都与电离层测高仪反演结果符合良好;在较低轨道高度（约500 km）,绝对TEC反演精度优于Abel积分反演精度;（2）绝对TEC反演的最大电子密度Nm较Abel积分法获得的结果更接近于电离层测高获得的峰值电子密度NmF2,绝对TEC反演法更加严密和有效;（3）周跳对绝对TEC反演结果的影响较Abel积分反演结果的影响更为敏感,但无论哪种方法,周跳对反演精度都造成严重损失.综合而言,绝对TEC反演法是更优的方法. 相似文献
14.
Zhdanov提出的拟线性(QL)近似方法有效解决了积分方程计算速度慢、占用内存多的缺陷,但因为三维频谱激电(SIP)换源频繁,每次换源都需要重新计算剖分单元的并矢格林函数和一次场,所以用该方法模拟三维SIP的速度仍然缓慢.本文根据一次场及格林函数的空间对称性,提出了一种适用于多源电磁法的快速QL近似正演方法.在此基础上,实现了Cole-Cole参数范围约束的三维SIP共轭梯度反演方法.理论试算结果表明:快速QL近似方法适用于三维SIP正演模拟,计算速度较QL近似方法有了很大提高;Cole-Cole参数范围约束的三维SIP共轭梯度反演方法对零频电阻率、极化率等参数的反演效果良好,将异常区域剖分成2000块时迭代一次只需约0.2s的时间,占用内存大约为140 MB. 相似文献
15.
16.
提出了一种计算体积分方程的方法－改进型局域非线性迭代（MLNI），并用其对三维井间电磁场进行反演计算，该方法将井间大尺度散射体分为近场和远场区域两部分，它们的位置和尺寸均随场点位置的变化而改变，采用局域非线性近似计算近场区域的影响，将远场区域的影响作为外部激励源，采用迭代方法计算，该方法具有计算速度快，所需内存量少，收敛性好的优点，在反演中采用基于MLNI的非线性反演方法，利用这种方法可以将成像区域集中于一定范围内而不是整个三维空间，该反演方法由于考虑到了雅可比矩阵元素的非线性项，因而与传统的Born迭代反演方法相比更为精确和稳定，考虑到信息量和计算机内存的限制，第1次成像采用双重体元分割法进行较粗略的成像，然后缩小成像范围进行第2次较精确的成像，数值计算结果表明，MLNI是一种有效的计算井间大尺度异常体散射场的方法，将该方法用于反演过程能够得到较高分辨率的三维井间电导率图像。 相似文献
17.
本文用三维有限差分法,在计算电导率异常的电磁响应时作了一次网格收敛试验.该收敛解和其它几种数值解法得到的结果作了比较.这些方法是Hohmann的积分方程法、Lee等的混合法和Gupta等的混合法.我们的结果基本上和Hohmann的以及Gupta等的结果一致.试验表明,有限差分法所用的计算机时间比SANGAM混合法的少.为了缩减迭代过程的计算时间,我们使用了一种新的技巧--变区域迭代法,这种方法的收敛速度和一般的（在全区域迭代的）迭代法的收敛速度相同. 相似文献
18.
Time-lapse monitoring is a powerful tool for observing dynamic changes in the subsurface. In particular it offers the potential for achieving inversion results with increased fidelity through the inclusion of complementary information from multiple time-steps. This inclusion of complementary information can reduce the need for spatial smoothing, without adding inversion artifacts to the resulting images. Commonly used time-lapse inversion methods include the ratio method, cascaded time-lapse inversion, difference inversion and differencing independent inversions. We introduce two additional methods in which both time-lapse data sets are inverted simultaneously. In the first, called temporally constrained time-lapse inversion, inversion of both datasets is done under a single optimization procedure and constraints are added to the regularization to ensure that the changes from one time to another are smooth. In the second method, called simultaneous time-lapse inversion, the inversions at time 1 and time 2 are performed simultaneously and constraints of smoothness and closeness to a reference model are applied to the difference image produced at each iteration, and subsequently, the constraints are updated at each iteration. Through both a numerical and a field example we compare the results of common time-lapse inversion methods as well as the introduced approaches. We found that of the commonly used time-lapse inversion methods the difference inversion method produced the best resolution of time-lapse changes and was the most robust in the presence of noise. However, we found that the alternative approach of simultaneous time-lapse inversion produced the best reconstruction of modeled EC changes in the numerical example and easily interpretable high resolution difference images in the field example. Moreover, there was less tailoring of regularization parameters with our simultaneous time-lapse approach, suggesting that it will lend itself well to an automated inversion code. 相似文献