共查询到18条相似文献,搜索用时 343 毫秒
1.
2.
3.
4.
采用有限元直接迭代算法实现了线源频率域测深电磁响应的二维正演计算. 首先给出了线源正演问题的有限元直接迭代格式,然后由迭代法进行求解. 在处理奇异源问题上,采用向内递推的组合网格技巧,在源点附近可进行局部加密,并实现粗细网格的对接,从而较好地解决了奇异源附近的计算问题. 还提出一种迭代求取全区视电阻率的方法,避免了远近区的划分. 通过对均匀半空间、层状介质和二维模型电磁响应的计算,获得了与大地电磁测深相似的视电阻率曲线,验证了算法的正确性;通过对计算结果的分析,在理论上说明了线源频率域近区测深的可行性. 相似文献
5.
利用有限单元法对水下激发极化测深进行数值模拟以研究其异常特征和适用条件.水下的电阻率测深和激电测深均可对水底地质体进行近距离探测,提高探测分辨率.而水底地形对视电阻率测深的影响很大.水底极化体有限元正演模拟的结果表明,水底地形对激发极化测深结果没有影响.采用水下激发极化法进行地质体探测是可行、有效的.通过计算不同水体和水底岩石电阻率对视极化率幅值的影响发现,水体电阻率是影响水下极化率测深的主要因素.由于海水的电阻率过低,观测视极化率异常微弱而不适合开展水下激电工作. 相似文献
6.
复电阻率法二维有限元数值模拟 总被引:9,自引:2,他引:9
伴随着复电阻率法的广泛应用,发展精确和快速的正演和反演算法成为复电阻率法研究的重点.本文采用基于三角单元剖分的有限单元法进行了复电阻率二维数值模拟研究.为了提高计算速度,对无穷远边界进行了近似处理.整个正演计算过程分为两步,首先采用有限单元法计算四个不同频率的视复电阻率数据,然后对前一步得到的视复电阻率数据采用递推算法计算视Cole-Cole参数.采用这种正演算法与一维正演的结果进行了对比,验证了本文方法的正确性.设计了两个二维极化模型,数值模拟结果表明视复电阻率和Cole-Cole视参数等值线断面图对于异常目标体都有比较明显的反映. 相似文献
7.
通过时间-频率转换关系,TEM数据可以转换成平面波场测深数据,从而可以对TEM资料进行拟平面波场处理解释.在对瞬变电磁视电阻率数据向平面波场测深视电阻率数据转换时,发现由于瞬变电磁使用晚期计算公式及装置问题,使测深曲线早期数据发生畸变.文中建立了视电阻率曲线进入晚期所满足的关系式,从理论上给出不同情况下瞬变电磁测深视电阻率曲线进入晚期的临界点.以瞬变电磁、大地电磁、CSAMT为例,对大量的模型进行正演计算,对计算结果进行对比分析,建立不同地表电性结构、不同时间延迟情况下,瞬变电磁早期数据误差的校正量板. 相似文献
8.
工程电法勘探中,常需要探测埋深不大、规模较小、分布复杂的目标体,对探测方法的精度要求较高.由于三维高密度电阻率法数据采集密度大,可对目标体进行多方位观测,是目前工程探测中常用的方法之一.本文对三维高密度电阻率E SCAN法的视电阻率异常特征进行研究.设计几种典型地质体模型,利用有限单元法进行正演计算,并分析其异常特征和分布规律.算例结果表明,E SCAN法对低阻体的分辨能力强于高阻体.无需对所有电极进行扫描供电观测,即可分辨地质体电性特征及水平位置,但观测结果难以分辨地质体纵向延伸. 相似文献
9.
工程电法勘探中,常需要探测埋深不大、规模较小、分布复杂的目标体,对探测方法的精度要求较高.由于三维高密度电阻率法数据采集密度大,可对目标体进行多方位观测,是目前工程探测中常用的方法之一.本文对三维高密度电阻率ESCAN法的视电阻率异常特征进行研究.设计几种典型地质体模型,利用有限单元法进行正演计算,并分析其异常特征和分布规律.算例结果表明,ESCAN法对低阻体的分辨能力强于高阻体.无需对所有电极进行扫描供电观测,即可分辨地质体电性特征及水平位置,但观测结果难以分辨地质体纵向延伸. 相似文献
10.
对常见的隐伏地质构造体进行大地电磁场正演模拟,并研究和总结了构造体对大地电磁场的响应特征.由麦克斯韦方程组导出大地电磁场的边界条件,应用有限单元法求解变分问题并得出刚度矩阵,根据线性方程组计算出相应的场值,再由改进的视电阻率定义式求出视电阻率.模拟结果显示:改进的视电阻率定义式相对于传统的Cagniard定义式能有效地提高视电阻率曲线的收敛速度和逼近程度;隐伏地质构造体对大地电磁场的响应具有一定的规律性,通过分析两种极化模式下的正演结果,可以有效地识别地质构造体的存在和分布情况.通过实例分析,得出隐伏地质构造正演模拟与实际观测资料的互相验证关系,为实践应用中的反演工作提供了重要的指导意义. 相似文献
11.
研讨了频率域电磁法中不同源装置的大地电磁测深、线源频率电磁测深和偶极源频率电磁测深阻抗视电阻率的源效应影响特征。在唯象分析的基础上,提出了几种电磁测深法阻抗视电阻率的相互换算法──源效应校正法(大地电磁测深二维TE极化视电阻率和其它两种电磁法的赤道装置二维阻抗视电阻车)。模型试验表明,利用这一源效应校正法可以由大地电磁二维视电阻率近似地计算出线源频率电磁二维阻抗视电阻率。这一方法被尝试应用于由线源频率电磁二维阻抗视电阻率估算偶极源频率电磁二维阻抗视电阻率。 相似文献
12.
Induced polarization and resistivity model studies over thin dykes for varying resistivity contrasts, depth of burial, and dip angles show striking parallelism of the surface apparent resistivity contours with the boundary of the body. This effect may be utilized for the estimation of strike length of the body. Results show that intermediate values of the electrode spacing is satisfactory for detection under widely varying conditions of resistivity contrast and depth. The percentage frequency effect (P.F.E.) and the metal factor (M.F.) responses are found to be more sensitive to the variation in the depth of burial than the resistivity responses. Pseudosections for P.F.E. and M.F. are concentrated much closer to the body than the resistivity pseudo-sections. 相似文献
13.
E. F. SCURTU 《Geophysical Prospecting》1972,20(3):605-625
The use of the electrical image method in calculating the potential of a source in the presence of a buried conductive sphere proves to be difficult because of the great number of terms making up the potential formula. The disadvantages contained by this method are avoided here by resorting to the theory of graphs. In view of this aim the “propagation” of the electrical images due to the sphere has been represented by an infinite graph to which an infinite matrix has been attached. Taking this as a basis several FORTRAN IV programs for the computation of resistivity pseudo-sections for various electrode configurations have been coded, and a great deal of such pseudo-sections have already been calculated with an IBM-360 computer. Using these results, recommendations regarding both the methodology of the geoelectrical prospecting of spherical shaped conductors and the most suitable way of building resistivity pseudo-sections are developed. 相似文献
14.
15.
The advantages of the Wenner tripotential method (Carpenter 1955) for apparent resistivity profiling are described and two new data processing techniques introduced as an aid to the interpretation of apparent resistivity sections (pseudo-sections). These techniques were developed from model data computed using a two-dimensional finite difference method. Oscillatory components present in anomalies on tripotential profiles and related to electrode spacing are shown to be effectively removed by linear filtering that also simplifies their form and aids recognition. Furthermore, it is shown that the ratio of the beta- and gamma-apparent resistivities is a good indicator of resistivity variation, and is particularly sensitive to lateral change. Model data indicates that, over a wide range of conditions, enough subsurface information can be obtained by inspection of tripotential resistivity and ratio profiles, and from space sections to make possible a useful—and sometimes semi-quantitative—interpretation. A rationale for the general interpretation of tripotential data is developed. Field data are described from an area of weathered granite basement in Nigeria. A model of the subsurface is developed using parameters derived from the processed observations. The observed and calculated apparent resistivity space sections are very similar. 相似文献
16.
E. MUNDRY 《Geophysical Prospecting》1984,32(2):336-346
For helicopter-borne electromagnetic systems, the distance between the transmitting and the receiving coils is small compared with the altitude above ground. For this case, a major simplification can be made for the calculation of model curves. Some two-layer curves for the interpretation of frequency measurements are presented. A very simple procedure is demonstrated for the conversion of the relative secondary field into apparent resistivity and apparent distance for the mapping of airborne electromagnetic data. Furthermore, an approximation is described for the determination of the thickness and the resistivity of a layer lying on a perfectly conducting half-space. 相似文献
17.
THE DISPERSION RELATION, OF EARTH''S ELECTROMAGNETIC RESPONSE FUNCTIONS AND THE JOINT INVERSION OF DATA 下载免费PDF全文
On the basis of the dispersion relations of MT field, the necessity and applied prospects of the joint inversions using a pair of MT response functions which are correlative with the dispersion relations, are infered. A filter coefficient algorithm is made, with which the corresponding impedance phase data can be estimated using a set of apparent resistivities. The tests for the observed MT data show that when comparing the impedance phase estimated using the dispersion relation with the ob served phase, it can be checked whether the dispersion relation between observed apparent resistivity and phase data is satisfied or not, and that the use of the phase data corrected using the dispersion relation in the joint inversion is advantageous to obtain more confident results. It is shown that joint inversions are more advantageous than single parameter inversions, and that in the most case the joint inversion using the apparent resistivities of impedance real and imaginary parts is more advantageous than the jointinversion using the normal apparent resistivity and impedance phase. The existence of the dipersion relations between the ratio apparent resistivity and corresponding impedance phase of the orthogonal electric and magnetic field horizontal Components in the frequency EM sounding with horizontal electric dipole(FEMS) are discussed, the better effect of the joint inversion using the pair of EM response functions is obtained. The problems on the one-dimensional joint inversion for the MT and FEMS apparent resistivities, for which the observed frequency bands partly overlape each other, are studied. It is shown that this joint inversion is applicable and effective:the joint inversions of the practical data for two kinds of EM methods at two sites give the results well corresponding to the drilling data. The simulated MT inversions for the data of two kinds of EM methods are made, and more confident results also are obtained. 相似文献
18.
文中利用"十.五"期间静海地震台布设的秒采样大地电磁网(N-MT)仪器系统和FHDZ-M15地磁组合系统的同步观测数据,实现了固定台站大地电磁测深(MT)观测。采用Robust技术对静海地震台2008年1—6月记录的天然电磁场观测数据进行了MT处理,得到了张量阻抗资料。经过对周期为6~600s范围的18个相同周期的视电阻率及阻抗相位随时间的变化进行分析,其中有多个周期的视电阻率及阻抗相位的时序曲线在2008年3月11日卢龙ML4.2地震前出现异常变化。研究使用的资料周期对应深度为5~50km,属华北地震震源深度范围,其结果反映了震源区地层地震前后震源深度的电性变化。MT方法比现行的地电观测方法具有更直接反映震源区地层的电性变化和更深的勘探深度的优势 相似文献