共查询到20条相似文献，搜索用时 46 毫秒
1.
2.
大地电阻率分布信息是影响磁共振地下水探测反演结果准确性的重要因素.在众多电磁法勘探技术中,瞬变电磁法具有高分辨率、高效率和大探测深度等优势,能准确探测地下几百米范围内的电阻率分布信息.因此磁共振与瞬变电磁联合解释方法具有重要意义.然而,利用单一测点拼接的磁共振与瞬变电磁联合解释方法进行模拟二维反演时存在解释结果不唯一,容易出现错误异常体等问题,尤其在复杂地质情况下,同一测线上相邻测点探测结果连续性差,解释结果偏离实际.基于此,本文提出磁共振与瞬变电磁横向约束联合反演方法（Laterally Constrained Inversion,简称LCI）,重点引入外推积分法（quadrature with extrapolation,简称QWE）,解决了传统正演过程中基于直接数值积分方法引起的求解效率低的问题,保证了联合反演方法的顺利实施,进而以相邻测点地下结构应具备连续性为依据,引入横向约束反演思想,通过在联合反演目标函数中加入相邻测点间各模型参数约束矩阵,提高磁共振解释结果准确性,加强探测剖面地质结构和含水模型连续性.经过理论模型证实,本文提出的LCI方法能有效提高传统一维反演结果的稳定性和唯一性.最后,对安徽黄山野外实际探测数据进行横向约束联合反演,验证了磁共振与瞬变电磁LCI联合反演方法的实用性.本文的研究成果将为磁共振与瞬变电磁空间约束联合反演奠定基础. 相似文献
3.
A set of geophysical data collected in an area in Iran are analyzed to check the validity of a geological map that was prepared in connection to a mineral prospecting project and also to image the spatial electrical resistivity distribution. The data set includes helicopter electromagnetic (HEM), airborne magnetic and ground electrical resistivity measurement. Occam approach was used to invert the HEM data to model the resistivity using a layered earth model with fixed thicknesses. The algorithm is based on a nonlinear inverse problem in a least-squares sense.The algorithm was tested on a part of an HEM dataset acquired with a DIGHEM helicopter EM system at Kalat-e-Reshm, Semnan in Iran. The area contains a resistive porphyry andesite that is covered by Eocene sedimentary units. The results are shown as resistivity sections and maps confirming the existence of an arc like resistive structure in the survey area. The resistive andesite seems to be thicker than it is indicated in the geological maps. The results are compared with the reduced to the pole (RTP) airborne magnetic anomaly field data as well as with two ground resistivity profiles. We found reasonable correlations between the HEM 1D resistivity models and 2D models from electrical resistivity tomography (ERT) inversions. A 3D visualization of the 1D models along all flight lines provided a useful tool for the study of spatial variations of the resistivity structure in the investigation area. 相似文献
4.
Bernhard Siemon 《Journal of Applied Geophysics》2009,67(3):206-218
For about three decades helicopter-borne electromagnetic (HEM) measurements have been used to reveal the resistivity distribution of the upper one hundred metres of the earth's subsurface. HEM systems record secondary fields, which are 3–6 orders of magnitude smaller than the transmitted primary fields. As both the primary fields and the secondary fields are present at the receivers, well-designed bucking coils are often used to reduce the primary fields at the receivers to a minimum. Remaining parts of the primary fields, the zero levels, are generally corrected by subtracting field values recorded at high altitudes (standard zero levelling) or estimated from resistivities of neighbouring lines or from resistivity maps (advanced zero levelling). These zero-levelling techniques enable the correction for long-term, quasi-linear instrumental drift. Short-term variations caused by temperature changes due to altitude variations, however, cannot be completely corrected by this procedure resulting in stripe patterns on thematic maps.Statistical methods and/or 2-D filter techniques called statistical levelling (tie-line levelling) and empirical levelling (microlevelling), respectively, used to correct stripe patterns in airborne geophysical data sets are, in general, not directly applicable to HEM data. Because HEM data levelling faces the problem that the parameter affected by zero-level errors, the secondary field, differs from the parameter generally levelled, the apparent resistivity. Furthermore, the dependency of the secondary field on both the resistivity of the subsurface and the sensor altitude is strongly nonlinear.A reasonable compromise is to microlevel both half-space parameters: apparent resistivity and apparent depth, followed by a recalculation of the secondary field components based on the half-space parameters levelled. Advantages and disadvantages of the diverse levelling techniques are discussed using a HEM data set obtained in a hilly region along the Saale River between the cities of Saalfeld and Jena in central Germany. It turns out from a comparison of apparent resistivity and apparent depth maps derived from levelled HEM data that manually advanced zero levelling of major level errors and automatic microlevelling of remaining minor level errors yield the best results. 相似文献
5.
We investigate helicopter electromagnetic (HEM) inversion schemes applied to synthetic and measured HEM sea ice profiling data. Direct HEM data-to-ice-thickness inversion is compared to three different formal, least squares layered earth inversion algorithms.By making several approximations, it is possible to directly invert a single channel measurement (i.e., the in-phase or quadrature component of a single frequency measurement) to obtain an estimate of sea ice thickness. Measurements from multiple input channels, however, can be used in a layered earth inversion to simultaneously recover several model parameters such as sea ice thickness, sea ice conductivity and sub-ice bathymetry. Synthetic data sets for a particular two-frequency HEM system showed that simple least squares inversion algorithms produce reliable estimates of sea ice thickness in cases where the ice is thicker than 3 m. These methods could also recover acceptable estimates of sea ice thickness when a thin, conductive, partially melted sea ice layer was present, and could determine shallow, sub-ice bathymetry in brackish water. As expected, 1D transformations and inversions of synthetic data for a three-dimensional pressure ridge keel structure contained artifacts, notably broadening of the apparent width of the keel.Prior to inverting a field data set acquired over rather thin (~ 0.5 m) Antarctic sea ice, we found it necessary to recalibrate the phase angle of the measurements using a phasor diagram-based method. Direct transformation of a single channel from the recalibrated data set produced more accurate estimates of sea-ice thickness than formal inversion of multi-channel data. We suggest that the least squares inversion methods are inferior in this situation because of the particular characteristics of the two-frequency HEM system used in this evaluation; the extreme differences in sensitivity of high and low frequency data components, the overall low sensitivity to sea ice conductivity (especially for thin ice), and the partially low signal-to noise ratios of the measurements. The data sets used in this study will be made available to the public to allow alternate inversion approaches to be applied and evaluated. It is suggested that inclusion of parameter bounds and other forms of regularization could help to improve the inversion results. 相似文献
6.
Sensitivity and inversion of marine electromagnetic data in a vertically anisotropic stratified earth 总被引：2，自引：0，他引：2
The controlled‐source electromagnetic (CSEM) and magnetotelluric method (MT) are two techniques that can be jointly used to explore the resistivity structure of the earth. Such methods have, in recent years, been applied in marine environments to the exploration and appraisal of hydrocarbons. In many situations the electric properties of the earth are anisotropic, with differences between resistivity in the vertical direction typically much higher than those in the horizontal direction. In cases such as this, the two modes of the time‐harmonic electromagnetic field are altered in different ways, implying that the sensitivity to the earth resistivity may vary significantly from one particular resistivity component (scalar, horizontal or vertical) to another, depending on the measurement configuration (range, azimuth, frequency or water depth). In this paper, we examine the sensitivity of the electromagnetic field to a vertically anisotropic earth for a typical set of configurations, compare inversion results of synthetic data characterizing a vertically anisotropic earth obtained using the isotropic and anisotropic assumptions and show that correctly accounting for anisotropy can prevent artefacts in inversion results. 相似文献
7.
Frequency-domain airborne electromagnetics is a proven geophysical exploration method. Presently, the interpretation is mainly based on resistivity—depth imaging and one-dimensional layered inversion; nevertheless, it is difficult to obtain satisfactory results for two- or three-dimensional complex earth structures using 1D methods. 3D forward modeling and inversion can be used but are hampered by computational limitations because of the large number of data. Thus, we developed a 2.5D frequency-domain airborne electromagnetic forward modeling and inversion algorithm. To eliminate the source singularities in the numerical simulations, we split the fields into primary and secondary fields. The primary fields are calculated using homogeneous or layered models with analytical solutions, and the secondary (scattered) fields are solved by the finite-element method. The linear system of equations is solved by using the large-scale sparse matrix parallel direct solver, which greatly improves the computational efficiency. The inversion algorithm was based on damping least-squares and singular value decomposition and combined the pseudo forward modeling and reciprocity principle to compute the Jacobian matrix. Synthetic and field data were used to test the effectiveness of the proposed method. 相似文献
8.
Piecewise 1D laterally constrained inversion of resistivity data 总被引：1，自引：0，他引：1
Esben Auken ers V. Christiansen Bo H. Jacobsen Nikolaj Foged Kurt I. Sørensen 《Geophysical Prospecting》2005,53(4):497-506
In a sedimentary environment, layered models are often capable of representing the actual geology more accurately than smooth minimum structure models. Furthermore, interval thicknesses and resistivities are often the parameters to which non‐geophysicist experts can relate and base decisions on when using them in waste site remediation, groundwater modelling and physical planning. We present a laterally constrained inversion scheme for continuous resistivity data based on a layered earth model (1D). All 1D data sets and models are inverted as one system, producing layered sections with lateral smooth transitions. The models are regularized through laterally equal constraints that tie interface depths and resistivities of adjacent layers. Prior information, e.g. originating from electric logs, migrates through the lateral constraints to the adjacent models, making resolution of equivalences possible to some extent. Information from areas with well‐resolved parameters will migrate through the constraints in a similar way to help resolve the poorly constrained parameters. The estimated model is complemented by a full sensitivity analysis of the model parameters, supporting quantitative evaluation of the inversion result. Examples from synthetic 2D models show that the model recognition of a sublayered 2D wedge model is improved using the laterally constrained inversion approach when compared with a section of combined 1D models and when compared with a 2D minimum structure inversion. Case histories with data from two different continuous DC systems support the conclusions drawn from the synthetic example. 相似文献
9.
The temporal variation in a soil moisture profile can be studied using resistivity sounding data acquired at different times. The layered earth model based estimation of soil moisture from apparent resistivity data is a two-step non-linear inversion. Firstly, the apparent resistivity data are inverted to derive the layer resistivity variations and thicknesses and, secondly, the moisture content is estimated from these layer resistivity variations using a calibration equation. The soil moisture–resistivity problem was studied using the one-dimensional formulation of resistivity problem. A generalized geoelectric earth model was considered to simulate the soil moisture distribution and its temporal variation in the unsaturated zone. An algorithm (RESMOS) for the interpretation of the apparent resistivity data in terms of soil moisture variations through this two-step inversion process is reported. 相似文献
10.
Airborne electromagnetic (AEM) surveys are currently being flown over populated areas and applied to detailed problems using high flight line densities. Interpretation information is supplied through a model of the subsurface resistivity distribution. Theoretical and survey data are used here to study the character and reliability of such models. Although the survey data were obtained using a fixed-wing system, the corresponding associations with helicopter, towed-bird systems are discussed. Both Fraser half-space and 1D inversion techniques are considered in relation to their ability to distinguish geological, cultural and environmental influences on the survey data. Fraser half-space modelling provides the dual interpretation parameters of apparent resistivity and apparent depth at each operational frequency. The apparent resistivity was found to be a remarkably stable parameter and appears robust to the presence of a variety of at-surface cultural features. Such features provide both incorrect altitude data and multidimensional influences. Their influences are observed most strongly in the joint estimate of apparent depth and this accounts for the stability of the apparent resistivity. Positive apparent depths, in the example data, result from underestimated altitude measurements. It is demonstrated that increasingly negative apparent depths are associated with increasing misfits between a 1D model and the data. Centroid depth calculations, which are a transform of the Fraser half-space parameters, provide an example of the detection of non-1D influences on data obtained above a populated area. 1D inversion of both theoretical and survey data is examined. The simplest use of the 1D inversion method is in providing an estimate of a half-space resistivity. This can be undertaken prior to multilayer inversion as an initial assessment. Underestimated altitude measurements also enter the problem and, in keeping with the Fraser pseudo-layer concept, an at-surface highly resistive layer of variable thickness can be usefully introduced as a constrained parameter. It is clearly difficult to ascribe levels of significance to a ‘measure’ of misfit contained in a negative apparent depth with the dimensions of metres. The reliability of 1D models is better assessed using a formal misfit parameter. With the misfit parameter in place, the example data suggest that the 1D inversion methods provide reliable apparent resistivity values with a higher resolution than the equivalent information from the Fraser half-space estimates. 相似文献
11.
Two-dimensional joint inversion of radiomagnetotelluric and direct current resistivity data 总被引：2，自引：0，他引：2
An algorithm for the two-dimensional (2D) joint inversion of radiomagnetotelluric and direct current resistivity data was developed. This algorithm can be used for the 2D inversion of apparent resistivity data sets collected by multi-electrode direct current resistivity systems for various classical electrode arrays (Wenner, Schlumberger, dipole-diplole, pole-dipole) and radiomagnetotelluric measurements jointly. We use a finite difference technique to solve the Helmoltz and Poisson equations for radiomagnetotelluric and direct current resistivity methods respectively. A regularized inversion with a smoothness constrained stabilizer was employed to invert both data sets. The radiomagnetotelluric method is not particularly sensitive when attempting to resolve near-surface resistivity blocks because it uses a limited range of frequencies. On the other hand, the direct current resistivity method can resolve these near-surface blocks with relatively greater accuracy. Initially, individual and joint inversions of synthetic radiomagnetotelluric and direct current resistivity data were compared and we demonstrated that the joint inversion result based on this synthetic data simulates the real model more accurately than the inversion results of each individual method. The developed 2D joint inversion algorithm was also applied on a field data set observed across an active fault located close to the city of Kerpen in Germany. The location and depth of this fault were successfully determined by the 2D joint inversion of the radiomagnetotelluric and direct current resistivity data. This inversion result from the field data further validated the synthetic data inversion results. 相似文献
12.
2.5维井间电磁反演在中国东部孤岛油田的应用 总被引：3，自引：1，他引：2
In this study, we present a practical technique of transforming cross-hole EM data into the inter-well resistivity distribution. The a priori information constraint is incorporated into an iterative regularized inversion procedure and a variable roughness is added into the inversion process. Finite element approximation based on a two and a half-dimensional （2.5D） model has been developed for the forward problem and the ＂pseudo-forward＂ problem needed for constructing the sensitivity matrix and synthetic data set. The regularized least-squares inversion scheme, constrained with the a priori information obtained from well logs, was adopted to reconstruct the inter-well resistivity profile from two synthetic electromagnetic data sets and field data acquired in the Gudao Oil Field, East China. The partial derivatives of the sensitivity matrix were computed by the adjoint equation based on the reciprocity principle. Inversion results of the synthetic and field data examples suggest that our method is robust and stable in the presence of random noise in the field data and can be used for cross-hole EM field data interpretation. 相似文献
13.
In order to interpret field data from small-loop electromagnetic (EM) instruments with fixed source–receiver separation, 1D inversion method is commonly used due to its efficiency with regard to computation costs. This application of 1D inversion is based on the assumption that small-offset broadband EM signals are insensitive to lateral resistivity variation. However, this assumption can be false when isolated conductive bodies such as man-made objects are embedded in the earth. Thus, we need to clarify the applicability of the 1D inversion method for small-loop EM data. In order to systematically analyze this conventional inversion approach, we developed a 2D EM inversion algorithm and verified this algorithm with a synthetic EM data set. 1D and 2D inversions were applied to synthetic and field EM data sets. The comparison of these inversion results shows that the resistivity distribution of the subsurface constructed by the 1D inversion approach can be distorted when the earth contains man-made objects, because they induce drastic variation of the resistivity distribution. By analyzing the integrated sensitivity of the small-loop EM method, we found that this pitfall of 1D inversion may be caused by the considerable sensitivity of the small-loop EM responses to lateral resistivity variation. However, the application of our 2D inversion algorithm to synthetic and field EM data sets demonstrate that the pitfall of 1D inversion due to man-made objects can be successfully alleviated. Thus, 2D EM inversion is strongly recommended for detecting conductive isolated bodies, such as man-made objects, whereas this approach may not always be essential for interpreting the EM field data. 相似文献
14.
基于飞艇的时间域地空电磁探测系统,具有勘探深度大、效率高、空间分辨率高、飞行控制容易等优势.但在低空飞行测量过程中,飞艇飞行高度、航迹、姿态等受风向、大气气流、地形、地面局部温度场变化等影响而发生变化,导致固定在艇囊前端的接收线圈发生运动,切割大地磁场,产生了电磁噪声、运动噪声、基线漂移等,从而影响电磁数据的电阻率成像质量.因此,研究地空电磁信号中多种噪声的去除方法,对数据的反演解释非常重要.由于地空电磁信号中有效信号频带与部分噪声频带相重叠,使用传统滤波或消噪方法具有一定局限性.因此,本文提出一种综合小波去噪法:根据地空电磁信号的特点,采用sym8小波基;基于小波多分辨率分析原理,利用小波高尺度近似分量估计基线漂移,以校正电磁数据中的基线;基于小波阈值收缩原理,采用5层小波分解、极小极大阈值配合硬收缩函数的消噪方法,来压制数据中的其余噪声.最后,通过异常环模型的理论响应和实测数据进行算法的验证,结果表明这种综合消噪法对多种噪声均有很好的抑制作用,是一种实用有效的时间域地空电磁数据消噪方法. 相似文献
15.
地电体对频率域电磁波激发源的响应为电磁感应和激电效应的综合响应.传统CSAMT法进行数据正反演时认为大地介质电阻率是与频率无关的实数,而实际上因为激电效应,地下可极化体的电阻率是一个与频率相关的复数.为推进二者总体响应研究,并扩展激电法的应用范围,同时提高电磁法勘探的精度,本文基于Dias模型,以复电阻率代替不考虑地电体极化效应的直流电阻率,对CSAMT场源一维层状模型进行了正演模拟,为提取CSAMT信号中所含激电信息提供理论基础.结果表明,考虑激电参数后,视电阻率及相位响应曲线出现明显异常（包括远场、过渡场、近场）;极化前后振幅比值异常峰值、相位差值异常峰值可直观体现激电异常;异常峰值与极化层层厚、埋深以及电阻率变化有连续的对应关系.认为从频率域电磁法信号中提取激电信息有乐观的前景. 相似文献
16.
电性源地空瞬变电磁法具有工作效率高、勘探深度大、采集信号信噪比高、适用于地形地质条件复杂地区等优点.但是,到目前为止,由于尚未建立起该方法的解释系统,大大制约了该方法的发展.本文旨在建立起完整的地空电磁探测系统,丰富整个探测系统的理论.本文围绕地空瞬变电磁法全域视电阻率定义、瞬变电磁虚拟波场的克希霍夫偏移成像、逆合成孔径成像方法三个科学问题进行了系统研究.提出了用磁场强度定义全域视电阻率的迭代算法,理论模型试验结果表明计算出的视电阻率曲线首支趋于第一层电阻率,尾支趋于最后一层电阻率,实现了全空域、全时域视电阻率的计算;在先前研究的基础上,实现了适合电性源地空装置的瞬变电磁虚拟波场的克希霍夫偏移成像;采用相关迭加技术,实现孔径内多测点数据合成,将传统的单点处理方式发展成为逐点推移多次覆盖的逆合成孔径处理方法.层状模型试验表明:(1)全域视电阻率能够光滑、完整、渐变地反映出模型的电性信息变化;(2)当改变三层模型中间层电阻率时,全域视电阻率曲线随着参数的改变分异明显,对电性层的识别容易且直观;(3)由于在电阻率计算中同时考虑了接收机高度、偏移距、时间等各参数的影响,全域视电阻率可实现全空域、全时域的视电阻率计算.含水采空区的复杂模型算例表明:(1)根据不同测线的全域视电阻率结果可以看出,在靠近采空区的位置,全域视电阻率断面可以清晰地反映出采空区的空间位置,随着测线离采空区越来越远,采空区异常越来越弱直至消失;(2)波场变换和偏移成像的结果显示存在两个电性差异较大的界面,上界面指示地表,由于空气和大地之间的电性差异较大,故该界面波场信号反映强烈,遍布整个区域,下界面异常信号则主要集中在中部,向外逐渐减弱,指示采空区;(3)逆合成孔径成像结果表明地表界面在合成前后没有变化,而采空区异常合成后范围明显变小,且异常边界清晰,指示的采空区位置与模型吻合很好.本文借助于逆合成孔径雷达成像的基本思想,建立了一套电性源瞬变电磁地空逆合成孔径成像方法.基于反函数思想结合迭代算法提出的电性源地空瞬变电磁法的全域视电阻率定义方法,实现了全空域、全时域的视电阻率计算;借鉴瞬变电磁拟地震偏移成像算法,实现了瞬变场的三维成像;借鉴逆合成孔径雷达的思想,提出电性源瞬变电磁地空逆合成孔径算法,进一步提高了成像的分辨率.采空区模型算例表明相关叠加合成确实具有增强有用信号、提高信噪比、提高分辨率的诸多优点,证实了瞬变电磁地空逆合成孔径成像方法的有效性. 相似文献
17.
采用Duhamel积分对斜阶跃电流激励的瞬变电磁场进行直接计算,应用叠加原理得到了多匝小回线全域瞬变电磁响应.研究了多匝小回线的早、晚延时临界条件以及等磁矩单匝大回线与多匝小回线的延时和关断时间特征,分析了电阻率、厚度等参数变化时斜阶跃电流对层状模型瞬变响应的影响特征,最后采用实测数据验证了关断效应的校正效果.结果显示,回线边长越小或表层电阻率越大,瞬变响应越早离开早期并进入晚期.多匝小回线的关断时间远大于等磁矩的单匝大回线,但由前者激发的瞬变响应进入晚期的时间却大幅提前,更有利于晚期条件下各种理论公式的使用.斜阶跃影响校正结果显示,增加多匝小回线的匝数可以扩展反演深度的范围,减少匝数则可以提高反演电阻率的精确性;实测数据处理结果显示校正方法产生了一定的效果.本研究补充了多匝小回线瞬变电磁理论内容,研究结果可以为复杂环境中小尺寸多匝线圈的工作参数选择提供参考. 相似文献
18.
Two different airborne electromagnetic methods were applied in the same area: the frequency-domain helicopter-borne electromagnetic (HEM) system operated by the Federal Institute for Geosciences and Natural Resources, Germany, and the time-domain SkyTEM system of the HydroGeophysics Group at the University of Aarhus, Denmark. For verification of and comparison with the airborne methods, ground-based transient electromagnetics and 2-D resistivity surveying were carried out. The target of investigation was the Cuxhaven valley in Northern Germany, which is a significant local groundwater reservoir. The course of this buried valley was revealed by drillings and the shape was determined by reflection seismics at several cross sections.We applied electrical and electromagnetic methods to investigate the structure of the valley filling consisting of gravel, sand, silt and clay. The HEM survey clearly outlines a shallow conductor at about 20m depth and a deeper conductor below 40m depth inside the valley. This is confirmed by 2-D resistivity surveying and a drilling. The thickness of the deeper conductor, however, is not revealed due to the limited investigation depth of the HEM system. The SkyTEM survey does not resolve the shallow conductor, but it outlines the thickness of the deeper clay layer inside the valley and reveals a conductive layer at about 180m depth outside the valley. The SkyTEM results are very consistent with ground-based transient electromagnetic soundings.Airborne electromagnetic surveying in general has the advantage of fast resistivity mapping with high lateral resolution. The HEM system is cost-efficient and fast, but the more expensive and slower SkyTEM system provides a higher depth of investigation. Ground-based geophysical surveys are often more accurate, but they are definitively slower than airborne surveys. It depends on targets of interest, time, budget, and manpower available by which a method or combination of methods will be chosen. A combination of different methods is useful to obtain a detailed understanding of the subsurface resistivity distribution. 相似文献
19.
Global optimization with very fast simulated annealing (VFSA) in association with joint inversion is performed for 1D earth structures. The inherent problems of equivalence and suppression in electromagnetic (EM) and direct current (DC) resistivity methods are studied. Synthetic phase data from multifrequency sounding using a horizontal coplanar coil system and synthetic apparent resistivity data from Schlumberger DC resistivity measurements are inverted individually and jointly over different types of layered earth structures. Noisy data are also inverted. The study reveals that global optimization of individual data sets cannot solve inherent equivalence or suppression problems. Joint inversion of EM and DC measurements can overcome the problem of equivalence very well. However, a suppression problem cannot be solved even after combination of data sets. This study reveals that the K-type earth structure is easiest to resolve while the A-type is the most difficult. We also conclude that the equivalence associated with a thin resistive layer can be resolved better than that for a thin conducting layer. 相似文献
20.
Parts of the consecutive apparent resistivity monitoring stations of China have recorded clear diurnal variations. The relative amplitudes of diurnal variations at these stations range from 1.3‰ to 5.8‰. The daily accuracies of apparent resistivity observation are better than 1‰, because the background electromagnetic noise is rather low at these stations. Therefore, the diurnal variations of apparent resistivity recorded at these stations are real phenomena. The diurnal variation shapes can be divided into two opposite types according to their characteristics. One type is that the apparent resistivity data decreases during the daytime but increases during the nighttime(Type 1). The other type is the apparent resistivity data increases during the daytime but decreases during the nighttime(Type 2). There is a correspondence between the diurnal and annual variation patterns of apparent resistivity. For the monitoring direction with diurnal variation of Type 1, the apparent resistivity decreases in summer and increases in winter. However, for the monitoring direction with diurnal variation of Type 2, the apparent resistivity increases in summer and decreases in winter.
We take an analysis on the mechanism of apparent resistivity diurnal variation, combining the influence factors of water-bearing medium's resistivity, the electric structure of stations, and the apparent resistivity sensitivity coefficient(SC)theory. Intuitively, diurnal variation of apparent resistivity is caused by diurnal variation of medium resistivity in the measured area. The diurnal variation of medium resistivity will inevitably be caused by the factors with diurnal variation. Among the possible factors, there is diural variation in earth tide and temperature.
Our analysis displays that apparent resistivity diurnal variation is not caused by the usually-believed earth tide, but by the ground temperature difference between daytime and nighttime. The earth tide strain is too small to cause remarkable effects on the apparent resistivity data. On the other hand, the daily tide strain has two peak-valley variations, and its phase and amplitude has a period of approximate 28 days. However, the apparent resistivity data do not show these corresponding features to earth tide. Furthermore, the detection range of current apparent resistivity stations is within a depth of several hundred meters. Within this depth range, the medium deformation caused by solid tide can be regarded as uniform change. Therefore, all monitoring directions and all stations will have the same pattern of diurnal variation.
In general, the temperature increases in the daytime but decreases in the nighttime. For most water-bearing rock and soil medium, its resistivity decreases as temperature increases and increases as temperature decreases. Diurnal temperature difference affects about 0.4m of soil depth. Therefore, resistivity of this surface thin soil layer decreases in the daytime while increases in the nighttime. Under layered medium model, SC of each layer represents its contribution to the apparent resistivity. For the stations with positive SC of surface layer, apparent resistivity decreases in the daytime but increases in the nighttime. While for the stations with negative SC of surface layer, apparent resistivity diurnal variations display the opposite shape. 相似文献