共查询到20条相似文献，搜索用时 125 毫秒
1.
本文基于垂直磁场分量研究了SOTEM数据的一维OCCAM反演方法,并将其应用于理论三维数据及野外实测数据的反演.对于大部分一维模型,OCCAM反演可取得较好的反演效果,且反演结果不依赖于偏移距;噪声对SOTEM数据的OCCAM反演具有较大影响,但当信号含噪水平不超过5%时,反演结果仍具有较好的准确性;若浅层存在较厚的低阻层,OCCAM反演结果对下部地层的分辨能力下降,仅能获得具有平均效应的电阻率.将一维算法应用于SOTEM三维数据的反演,会产生较大的误差,尤其是在异常体边缘地带影响最为严重.该影响程度与异常体和背景电阻率之间的差异有关,对于大多数电性近似呈连续变化的真实大地而言,一维OCCAM反演算法仍可获得较好的效果.最后通过陕西某煤田深部富水性调查的实测SOTEM数据反演验证了本文的研究成果. 相似文献
2.
本文提出了磁总场异常垂直梯度三维相关成像方法,用于成像地下等效磁源分布.它首先将地下待成像空间划分为三维规则网格,然后直接计算每个网格节点磁偶板子在观测面理论磁总场异常垂直梯度与实测磁总场异常垂直梯度的互相关,其相关系数值表征等效磁偶极子分布(即磁偶极子发生的概率).理论长方体组合模型数据和实际某矿区磁测资料试验结果表明本文方法计算得到的相关系数值能基本反映地下的磁源分布,且分辨率明显高于磁总场异常三维相关成像的分辨率,也高于基于熵滤波分离异常的磁总场异常三维相关成像的分辨率. 相似文献
3.
David Beamish 《Journal of Applied Geophysics》1994,32(4)
Very low frequency electromagnetic (EM) methods using VLF transmitters have found many applications in subsurface geophysical investigations. Surface measurements involving both the vertical component of the magnetic field (VLF-EM or VLF-Z) and of the apparent resistivity (VLF-R) are increasingly common. Although extensive VLF data sets have been successfully used for mapping purposes, modelling and interpretation techniques which asess the third (i.e. depth) dimension appear limited.Given a profile of VLF-R measurements the main purpose of the present study is to demonstrate an automatic method for the construction of a resistivity cross-section. The technique used is one of a new generation of regularised inversion methods. These techniques attempt to overcome the problem of equivalence/non-uniqueness in EM sounding data by constructing the resistivity distribution with the minimum amount of structure that fits the data.VLF data represent a special case of plane-wave EM sounding in that they conform, in practice, to a single-frequency technique. This fact imposes a limitation in the amount of vertical resolution that we can expect using such data. In the case of two-dimensional modelling and inversion, resolution through the cross-section is a resultant attribute from both vertical and lateral resistivity gradients within the subsurface. In order to provide insight into the practical application of regularised inversion techniques to VLF data, both synthetic and field examples are considered. Both sets of examples are primarily concerned with VLF data applied to near-surface fault mapping where the main aim is to assess the location, dip and depth extent of conductive subsurface features. 相似文献
4.
本文提出了能提高异常体分辨能力,同时得到绝对电导率的地面磁电阻率数据三维反演方法.磁电阻率响应用准直流的低频磁场代替;数值模拟由频率域电场满足的Helmholtz方程出发,采用三维交错网格有限差分法;长直导线源作为发射源,其中源的计算包含在背景场中;结合地面磁电阻率数据各分量的特点,选择y分量进行反演研究;反演采用三维非线性共轭梯度反演技术,为了提高异常体的深度分辨能力,进行迭代重构反演;用印模法对初始模型进行重构,采用的是辅模型在浅部,元模型在深部的组合方式.从合成数据和实际数据的反演结果可以得到以下的认识:(1)由频率域麦克斯韦方程组出发,低频磁场数据反演可以直接得到电导率,而不是相对电导率之比;(2)采用印模法组合初始模型,进行迭代重构反演,可以提高地面磁电阻率数据反演对异常体的分辨能力,确定埋深位置,同时不会丧失对于浅部异常体的分辨能力;(3)在结合印模法的地面磁电阻率数据三维反演中,深部异常体的分辨能力受地表不均匀导电体影响较小;(4)确定印模深度可以采用上一次重构反演结束时的模型变化量,通过相邻两次重构反演结束时的模型变化量之差来确定迭代重构是否终止.因为静磁场与重力场在数学上的相似性,本文的反演方法可以被运用到重力场等位场的地面数据的反演中. 相似文献
5.
A code for 3-D resistivity modelling and inversion of vertical electrical soundings has been developed based on the finite-element technique and regularisation method. Synthetic data were used to test the effectiveness of the code and to examine the resolving capability of the Schlumberger array in investigating 3-D resistivity distributions. The code was applied to experimental data set constituted by 35 Schlumberger soundings collected near the Cairo city in order to study the subsurface resistivity distribution. The results have shown that valuable imaging of the subsurface resistivity distribution can be constructed even when the vertical electrical soundings are acquired in a sparse field data set. 相似文献
6.
J. van der Kruk J.A.C. Meekes P.M. van den Berg J.T. Fokkema 《Geophysical Prospecting》2000,48(6):1033-1052
Apparent resistivity is a useful concept for initial quickscan interpretation and quality checks in the field, because it represents the resistivity properties of the subsurface better than the raw data. For frequency‐domain soundings several apparent‐resistivity definitions exist. One definition uses an asymptote for the field of a magnetic dipole in a homogeneous half‐space and is useful only for low induction numbers. Another definition uses only the amplitude information of the total magnetic field, although this results in a non‐unique apparent resistivity. To overcome this non‐uniqueness, a complex derivation using two different source–receiver configurations and several magnetic field values for different frequencies or different offsets is derived in another definition. Using the latter theory, in practice, this means that a wide range of measurements have to be carried out, while commercial systems are not able to measure this wide range. In this paper, an apparent‐resistivity concept is applied beyond the low‐induction zone, for which the use of different source–receiver configurations is not needed. This apparent‐resistivity concept was formerly used to interpret the electromagnetic transients that are associated with the turn‐off of the transmitter current. The concept uses both amplitude and phase information and can be applied for a wide range of frequencies and offsets, resulting in a unique apparent resistivity for each individual (offset, frequency) combination. It is based on the projection of the electromagnetic field data on to the curve of the field of a magnetic dipole on a homogeneous half‐space and implemented using a non‐linear optimization scheme. This results in a fast and efficient estimation of apparent resistivity versus frequency or offset for electromagnetic sounding, and also gives a new perspective on electromagnetic profiling. Numerical results and two case studies are presented. In each case study the results are found to be comparable with those from other existing exploration systems, such as EM31 and EM34. They are obtained with a slight increase of effort in the field but contain more information, especially about the vertical resistivity distribution of the subsurface. 相似文献
7.
Geophysical signatures of uranium mineralization and its subsurface validation at Beldih,Purulia District,West Bengal,India: a case study 
下载免费PDF全文
下载免费PDF全文 Animesh Mandal W.K. Mohanty S.P. Sharma Arkoprovo Biswas Joydip Sen A. K. Bhatt 《Geophysical Prospecting》2015,63(3):713-726
The Beldih open cast mine of the South Purulia Shear Zone in Eastern India is well known for apatite deposits associated with Nb–rare‐earth‐element–uranium mineralization within steeply dipping, altered ferruginous kaolinite and quartz–magnetite–apatite rocks with E–W strikes at the contact of altered mafic–ultramafic and granite/quartzite rocks. A detailed geophysical study using gravity, magnetic, and gradient resistivity profiling surveys has been carried out over ～1 km2 area surrounding the Beldih mine to investigate further the dip, depth, lateral extension, and associated geophysical signatures of the uranium mineralization in the environs of South Purulia Shear Zone. The high‐to‐low transition zone on the northern part and high‐to‐low anomaly patches on the southeastern and southwestern parts of the Bouguer, reduced‐to‐pole magnetic, and trend‐surface‐separated residual gravity–magnetic anomaly maps indicate the possibility of highly altered zone(s) on the northern, southeastern, and southwestern parts of the Beldih mine. The gradient resistivity survey on either side of the mine has also revealed the correlation of low‐resistivity anomalies with low‐gravity and moderately high magnetic anomalies. In particular, the anomalies and modeled subsurface features along profile P6 perfectly match with subsurface geology and uranium mineralization at depth. Two‐dimensional and three‐dimensional residual gravity models along P6 depict the presence of highly altered vertical sheet of low‐density material up to a depth of ～200 m. The drilling results along the same profile confirm the continuation of uranium mineralization zone for the low‐density material. This not only validates the findings of the gravity model but also establishes the geophysical signatures for uranium mineralization as low‐gravity, moderate‐to‐high magnetic, and low‐resistivity values in this region. This study enhances the scope of further integrated geophysical investigations along the South Purulia Shear Zone to delineate suitable target areas for uranium exploration. 相似文献
8.
Interpretation of a single geophysical data set is not sufficient to get complete subsurface information. Cooperative or joint inversion of geophysical data sets is the preferred method for most case studies. In the present study, we present the results of the cooperative inversion approach of direct current resistivity and gravity data. The algorithm uses fuzzy c-means clustering to determine the petrophysical relationship between density and resistivity to obtain the similarity. Synthetic data set has demonstrated that the cooperative inversion approach can produce more reliable and better resistivity and density models of the subsurface as compared to those obtained through individual inversions. To utilize the presented cooperative inversion algorithm, the number of geologic units (number of clusters) in the study region must be known a priori. As a field study, the cooperative inversion approach was used to identify the extension of uranium-bearing target rock around the Beldih open cast mine. We noted the inconsistencies in both resistivity and density models obtained from the individual inversions. However, the presented cooperative inversion approach was able to produce similar resistivity and density models while maintaining the same error level of the respective individual inversions. We have considered four geologic units in the presented cooperative inversion as a field case study. We have also compared our cooperative results with drilled borehole and found to be a reliable tool to differentiate between the target rocks (kaolinite and quartz–magnetite–apatite rocks) and the ultramafic rock (host rock quartzite/alkaline granite). However, this study is subject to certain limitations such as the inability to differentiate between closely spaced kaolinite and quartz–magnetite–apatite rocks. 相似文献
9.
Integrated geophysical surveys using vertical electrical sounding (VES), very low frequency (VLF) EM, radiation counting,
total magnetic field and self-potential (SP) measurements are carried out to characterize the geothermal area around a hot
spring in the Nayagarh district, Orissa, India that lies in the East Indian geothermal province. The study was performed to
delineate the fracture pattern, contaminated groundwater movement and possible heating source. VES interpretations suggest
a three- to four-layer structure in the area. Resistivity survey near the hot spring suggests that weathered and fractured
formations constitute the main aquifer system and extend to 60 m depth. Current flow measured at various electrode separations
normalized by the applied voltage suggests that fractures extend to a greater depth. Detailed VLF study shows that fractures
extend beyond 70 m depth. VLF anomaly has also very good correlation with the total magnetic field measured along the same
profiles. Study results suggest that a gridded pattern of VLF survey could map the underground conductive fracture zones that
can identify the movement of contaminated groundwater flow. Therefore, precautionary measures can be taken to check further
contamination by delineating subsurface conducting structures. Self potential (SP) measured over the hot spring does not show
a large anomaly in favor of the presence of a sulphide mineral body. A small positive (5–15mV) SP anomaly is measured which
may be streaming potential due to subsurface fluid flow. A high radiation is measured about four kilometers from the hot spring,
suggesting possible radiogenic heating. However, the exact nature of the heating source and its depth is not known in the
area. Deep resistivity followed by a magneto-telluric survey could reveal the deeper structures. 相似文献
10.
We develop a two‐dimensional full waveform inversion approach for the simultaneous determination of S‐wave velocity and density models from SH ‐ and Love‐wave data. We illustrate the advantages of the SH/Love full waveform inversion with a simple synthetic example and demonstrate the method's applicability to a near‐surface dataset, recorded in the village ?achtice in Northwestern Slovakia. Goal of the survey was to map remains of historical building foundations in a highly heterogeneous subsurface. The seismic survey comprises two parallel SH‐profiles with maximum offsets of 24 m and covers a frequency range from 5 Hz to 80 Hz with high signal‐to‐noise ratio well suited for full waveform inversion. Using the Wiechert–Herglotz method, we determined a one‐dimensional gradient velocity model as a starting model for full waveform inversion. The two‐dimensional waveform inversion approach uses the global correlation norm as objective function in combination with a sequential inversion of low‐pass filtered field data. This mitigates the non‐linearity of the multi‐parameter inverse problem. Test computations show that the influence of visco‐elastic effects on the waveform inversion result is rather small. Further tests using a mono‐parameter shear modulus inversion reveal that the inversion of the density model has no significant impact on the final data fit. The final full waveform inversion S‐wave velocity and density models show a prominent low‐velocity weathering layer. Below this layer, the subsurface is highly heterogeneous. Minimum anomaly sizes correspond to approximately half of the dominant Love‐wavelength. The results demonstrate the ability of two‐dimensional SH waveform inversion to image shallow small‐scale soil structure. However, they do not show any evidence of foundation walls. 相似文献
11.
利用有限元2.5维可控源音频大地电磁法（简称CSAMT）数值模拟方法,对100Ωm均匀半空间介质中有限长度的电偶极源产生的电场、磁场及视电阻率、相位特征进行了数值模拟,研究了场的空间变化规律. 在一个象限中,场的特征存在双叶现象,当收发距大于4个趋肤深度时,电阻率较接近介质真实的电阻率,这些结果为观测系统和收发距的选择提供了依据. 波数域场的特征表明,低波数对源的贡献占较大的比例. 有限元法2.5维CSAMT数值模拟的优势在于能较准确地获得复杂介质结构的波场特征. 本文结合直立异常体、倾斜异常体及断陷模型对CSAMT电阻率、相位剖面特征及频率曲线特征的可靠性进行了研究. 数值模拟结果直观地给出了异常体的剖面异常形态. 通过对比研究异常体的剖面异常形态和半空间场的特征进一步说明本文方法和软件在模拟复杂介质结构场特征时是可靠的. 这为认识观测数据,指导反演解释提供了较好的依据. 相似文献
12.
A horizontal transmitter loop (vertical magnetic dipole) is used for frequency electromagnetic (FEM) soundings. The frequency ranges from approximately 6 Hz to about 4000 Hz. The vertical and radial magnetic field components are measured for 20 frequencies per decade several hundred meters from the transmitter loop. A very small bandwidth is selected for amplification using a reference signal. An Apple computer is used for data acquisition. A computer program for Marquardt inversion optimizes the parameters for the n-layer case: the resistivities and thicknesses of individual beds and a correction factor for the primary magnetic field. Interpretation of each component individually yields practically the same parameters. Examples from the field are given with interpretation; comparison with dc resistivity measurements is provided. The ratio of vertical/radial magnetic field components vs. frequency can be transformed simply into apparent resistivity vs. apparent depth. This can be done in the field to obtain an estimation of the depth of the layer boundaries. FEM results are compared with Schlumberger d.c. sounding obtained at the same site. 相似文献
13.
We describe a least-squares inversion approach to estimating the subsurface resistivity structure from cross-hole or borehole-to-surface electromagnetic data. It is assumed that the resistivity distribution is symmetric about the axis of a borehole and that vertical magnetic dipoles are located on the borehole axis. The receivers are placed either in another borehole or on the earth's surface. The inversion scheme uses the finite-element and smoothness-constrained least-squares methods. The computational effort required to obtain partial derivatives is reduced considerably by using the reciprocity principle. Numerical simulations show that the reconstructions are generally in good agreement with the true structures when the assumption of an axisymmetric earth structure holds. An example involving the breakdown of this assumption, which can be obtained by interchanging the source and receiver boreholes, suggests that the inversion result may also be useful for locating a general 3D anomaly although artifacts are present. 相似文献
14.
《Journal of Applied Geophysics》2006,58(2):158-175
The interpretation of airborne VLF data represents an important aspect of geophysical mapping of the upper few hundred meters of the Earth's crust, especially in areas with crystalline rocks. We have examined the ability of the single frequency VLF method to provide quantitative subsurface resistivity information using two generic models and standard airborne parameters with a flight altitude of 70 m and a frequency of 16 kHz. The models are long thin conductor (10 m thick, 10 Ω m resistivity and 1 km long) and a wider buried conductive dike (100 Ω m resistivity and 500 m wide). Using standard regularized inversion it turned out that for both models the conductivity of the conductors are underestimated and the vertical resolution is rather poor. The lateral positions of the minimum of the resistivity distributions coincide well with the true positions of the shallow conductors. For deeper conductors the position of the minimum resistivity moves from the edges of the conductor into the conductor. The depth to the minimum of the resistivity anomalies correlates well with the true depth to the top of the conductors although the latter is always smaller than the former.Interpretation of field airborne data collected at 70 m flight height resolved both small scale and large scale near surface conductors (conductance ∼1 S). Deeper conductors show up in the VLF data as very long wavelength anomalies that are particularly powerful in delineating the lateral boundaries of the conductors. Many of the VLF anomalies in the Stockholm area are dominated by these deep conductor responses with some near surface conductors superimposed. The deep conductors often follow topographic lows coinciding with metasediments. We interpret the frequent absence of near surface responses at 70 m flight height as a result of weak coupling between the primary VLF wave and the small scale (in all three dimensions) near-surface conductors.Radio magnetotelluric (RMT) ground measurements were carried out along a short profile coinciding with part of an airborne profile. Using data at 9 frequencies (14–250 kHz) small scale conductors in the upper few tens of meters, not identified from the airborne data, could be well resolved. Large scale deeper conductors could be identified by both methods at nearly the same positions. 相似文献
15.
Laust B. Pedersen Lena Persson Mehrdad Bastani Sören Byström 《Journal of Applied Geophysics》2009,67(3):250-258
Airborne VLF data are routinely collected by The Geological Survey of Sweden (SGU) as part of its bedrock mapping programme. In this paper we demonstrate that the novel Tensor VLF technique developed at Uppsala University and SGU can provide useful qualitative and quantitative information about the electrical conductivity distribution in the upper few hundred meters. Single transmitter scalar VLF maps emphasize those conductive structures that have dominant strikes in the direction of the transmitter. The tensor tipper (essentially the vertical magnetic field from currents along the strike direction) calculated from multiple transmitters is dependent only upon the underlying conductivity structure. Transformation of the tipper into the peaker (the horizontal divergence) has proven to enhance the lateral resolution while the transformation to the apparent resistivity can be used to discriminate different rock types. Two case histories from the application of VLF data are presented in this study. Two dimensional structures can be quantitatively modelled by modern inversion methods developed originally for deep electromagnetic MT soundings. Direct inversion of the real and imaginary parts of the tipper provides more quantitative information about the subsurface resistivity distribution. 相似文献
16.
复杂地形、地质条件的大地电磁数据解释容易出现假象,采用三维正演技术模拟地形和地表不均匀体的背景响应,对实测数据阻抗相位不变量进行校正,实现更准确的定性分析;对三维异常体模型的合成数据进行一维、二维多参数反演试算,以确定地形剧变区选择反演技术的最佳方案.合成数据的试反演结果显示一维反演水平切片假异常较多,二维反演能压制测向假异常,但不能压制走向的假异常,水平切片多出现测向条带.使用本文提出的阻抗相位不变量校正法扣除地形、地表背景响应,结合一维、二维反演,能使实际资料解释成果更加可靠. 相似文献
17.
地表起伏地形在野外矿产资源勘察中不可避免,其对直流电阻率法勘探影响巨大.近年来,电阻率三维正演取得诸多进展,特别是应用非结构网格我们能够进行任意复杂地形和几何模型的电阻率三维数值模拟,但面向实际应用的起伏地形下电阻率三维反演依然困难.本文基于非结构化四面体网格,并考虑到应用GPS/GNSS时,区域地球物理调查中可非规则布设测网的实际特点,实现了任意地形(平坦或起伏)条件下、任意布设的偶极-偶极视电阻率数据的不完全Gauss-Newton三维反演.合成数据的反演结果表明了方法的有效性,可应用于复杂野外环境下的三维电法勘探. 相似文献
18.
可控源音频大地电磁法(CSAMT)一直沿用大地电磁法(MT)的办法,通过计算电场分量与磁场分量的比值,求取卡尼亚视电阻率.而CSAMT场源已知,电场分量和磁场分量都与地下电阻率存在一定的关系,可以单独采用CSAMT电场分量或者磁场分量提取地下介质的视电阻率.本文通过分析电场分量与磁场分量的数据特性,提出利用CSAMT电场单分量数据进行视电阻率的计算,用改进的广义逆矩阵反演方法,使初始模型中的地电层数等于频道个数,克服了以往反演计算中层数较少的问题;实现全场区电场分量视电阻率曲线的拟合反演.同时对单分量视相位计算方法进行分析,结合山西大同地区积水采空区探测及数据解释结果,论证本文提出的单分量解释方法的有效性. 相似文献
19.
J. DUROUX 《Geophysical Prospecting》1967,15(4):564-583
The electromagnetic field radiated from a magnetic dipole lying on the ground is considered, in the extremely low frequency range (DC to 20 000 Hertz). Theoretical and experimental data are given on the characteristics of the surface wave (vanishing wave) generated at the air-ground interface, in the case of an homogeneous subsurface. The case of a subsurface with electrical resistivity varying with depth is considered. It is shown how the above-mentioned characteristics may be applied in the quantitative investigation of the electrical resistivity of the subsurface as a function of depth, in a method using the measurements on the ground of all the components of the radiated field: horizontal electrical components, vertical and horizontal magnetic components. 相似文献
20.
建立客观的磁力ΔT异常,获得可信的磁力ΔZ⊥异常、居里面成果, 在此基础上开展综合地质地球物理研究,对于解决南海资源等相关的地质、地球物理问题意义重大.本文汇集我国历年在南海取得的磁力实测资料,重新处理、校正、整合,特别是追溯南海周边的国际地磁台的日变资料进行日变改正,形成的南海磁力ΔT异常,其资料品质得到了极大提高.采用全磁纬变倾角化极技术,获得磁力ΔZ⊥异常,其与实测证实的磁性海山、海岛有非常好的匹配关系,建立起可以直接用于反演处理解释的南海磁力ΔZ⊥异常.选择磁力ΔZ⊥异常的小波多尺度分解四阶逼近,采用单一磁性界面反演技术进行深度反演,获得居里面又称磁性底界面.建立了认识、解释南海的磁场基础体系,证实了南海基性岩浆岩分布区与磁力ΔZ⊥异常的正磁异常高密切相关,区域性的磁力ΔZ⊥异常负磁异常区可能与居里面深度相关. 相似文献