首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
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.  相似文献   

2.
In recent years, geophysical methods (shallow seismic, electromagnetic, resistivity, ground penetrating radar) have been increasingly applied to overburden investigations. Their effectiveness has been found to depend significantly on local geological conditions. Compared with advanced seismic techniques, EM methods are faster and hence more cost-effective, but they have not been considered sufficiently accurate. Analysis is carried out of data obtained with the multifrequency horizontal-loop method (HLEM) in northeastern Ontario, where the overburden consists of Quaternary glacial and glaciolacustrine sediments. Surveying along 1-6 km long profiles permitted recognition of bedrock inhomogeneities and selection of sites suitable for HLEM data interpretation using the layered model. Phasor diagrams and computer inversion based on the ridge regression technique were used to interpret HLEM soundings obtained at eight frequencies. Interpreted layer resistivities and thicknesses were correlated with the results of Rotasonic drilling at 70 sites. Relatively accurate estimates of overburden thickness (within 10%) could be obtained in about 80% of the cases. Nine examples of HLEM soundings are given and discussed: three each of one-, two- and three-layer situations. An appropriate interpretation model cannot be selected simply by minimizing the rms error or by analysing the parameter resolution matrix. Frequently, the most effective way of evaluating a solution is to consider whether resistivity values determined by inversion fit any of the ranges determined by statistical analyses of sediment resistivities. A previously published study of electrical properties of Quaternary sediments indicated that resistivities of clay, till and sand are stable within a fairly large area, such as the one under investigation. While the application of HLEM methods to mapping of Quaternary sediments can be considered a success, interpretation of EM data in regions covered by glacial sediments is more difficult than in weathered terrains, where near-surface layering is more predictable. The problem of equivalence causes non-uniqueness in interpretation. Thickness equivalence, which results in poor resistivity estimates, was found to affect areas convered by sand and till. Conductance equivalence caused poor resolution of thickness and resistivity for thin clay layers (less than 10 m).  相似文献   

3.
Airborne electromagnetic (AEM) methods are increasingly being used as tools of geological mapping, groundwater exploration and prospecting for coal and lignite. In such applications, quantitative interpretation is commonly based on the layered-earth model. A new approach, a damped least-squares inversion with singular value decomposition, is proposed for interpretation of time-domain, towed-bird AEM data. Studies using theoretical and field AEM data indicate that inversion techniques are dependable and provide fast converging solutions. An analysis has been made of the accuracy of model parameter determination, which depends on resistivity and thickness distribution. In the common case of conductive overburden, upper-layer resistivity and thickness are usually well determined, although situations exist where their separation becomes difficult. In the case of a resistive layer overlying a conductive basement, the layer thickness is the best-determined parameter. In both cases, estimates of basement resistivity are the least reliable. Field data obtained with the Chinese-made M-l AEM system in Dongling, Anhui Province, China, were processed using the described inversion algorithm. The survey area comprised fluvial Cenozoic clays and weathered Mesozoic sediments. Inversion of AEM data resulted in accurate depth-to-bedrock sections and realistic estimates of the resistivities of overburden and bedrock which agree with the results of drilling and resistivity sounding.  相似文献   

4.
Results of gravimetric surveys can be interpreted by comparing the Bouguer field values with master curves based on simplified geological models. It has been shown in a previous paper how this procedure can be transformed into routine processes which can be computerized. The application of this method has yielded useful results in detailed gravity surveys. The present paper discusses the application of the same interpretation principles to magnetic data. After some modifications, the method elaborated for the gravimetric data can be used for the interpretation of magnetometric survey results. Magnetometric-tectonic maps are obtained which show the structural picture by common geological symbols. In the case of faults, the dimensions of depth of burial and throw are indicated on the maps. The method is illustrated by an example where these procedures have been applied to gravimetric and magnetometric data of the same area. Two different maps are obtained: One shows the tectonics according to density contrasts and the other map depicts the tectonic situation on the basis of magnetic susceptibility contrasts.  相似文献   

5.
We investigate a novel way to introduce resistivity models deriving from airborne electromagnetic surveys into regional geological modelling. Standard geometrical geological modelling can be strengthened using geophysical data. Here, we propose to extract information contained in a resistivity model in the form of local slopes that constrain the modelling of geological interfaces. The proposed method is illustrated on an airborne electromagnetic survey conducted in the region of Courtenay in France. First, a resistivity contrast corresponding to the clay/chalk interface was interpreted confronting the electromagnetic soundings to boreholes. Slopes were then sampled on this geophysical model and jointly interpolated with the clay/chalk interface documented in boreholes using an implicit 3D potential‐field method. In order to evaluate this new joint geophysical–geological model, its accuracy was compared with that of both pure geological and pure geophysical models for various borehole configurations. The proposed joint modelling yields the most accurate clay/chalk interface whatever the number and location of boreholes taken into account for modelling and validation. Compared with standard geological modelling, the approach introduces in between boreholes geometrical information derived from geophysical results. Compared with conventional resistivity interpretation of the geophysical model, it reduces drift effects and honours the boreholes. The method therefore improves what is commonly obtained with geological or geophysical data separately, making it very attractive for robust 3D geological modelling of the subsurface.  相似文献   

6.
Since its development some thirty years ago, the airborne electromagnetic (AEM) method has been primarily used as a tool for mineral prospecting. However, advanced AEM systems are capable of other tasks, such as geological mapping and groundwater exploration. Excellent correlation between maps of apparent conductivity and geological maps was observed in several regions of Brazil where AEM surveys were performed. The degree of correlation seems to depend on the local climate. In humid and subhumid tropical regions, a weathered layer develops whose thickness and conductivity depend upon bedrock lithology. Therefore certain lithological types can be recognized from their conductivity signature; e.g., granites and Precambrian coarse clastic rocks are resistive, metavolcanic (particularly mafic) and volcanic rocks are conductive, Phanerozoic sediments are generally highly conductive. Two geophysical surveys are analyzed in the paper. The first was conducted with the time-domain, towed-bird AEM system in the Itapicuru greenstone belt in the state of Bahia. The apparent conductivity map correlated better with the local lithology than the magnetic map. Results of the AEM survey were successfully used to improve the regional geological map. A helicopter EM system was used in the second survey, which covered a portion of the Precambrian shield of Rio Grande do Sul. Also in this region ground checks confirmed the usefulness of conductivity surveys in geological mapping. The technique outlined in the paper holds great promise for countries of humid tropical climate, where few outcrops exist and access is often difficult. The tests performed in Brazil indicate that by executing AEM/aeromagnetic surveys during initial stages of mapping and exploration programs, time and expenditure required for geological field work can be considerably reduced. The resulting geological maps are more accurate and the inventory of mineral occurrences becomes more complete.  相似文献   

7.
The electrical properties of the weathered layer were investigated by means of Geonics VLF-EM 16/16R equipment in two areas of the Andhra Pradesh State. The resistivity and thickness of the weathered layer were found to be variable even over a small survey area. In the area underlain by Precambrian granite-gneiss, most of the recorded VLF-EM anomalies were caused by variations in the resistivity of the weathered layer. Changes in thickness were well reflected in the VLF-EMR curves. The second area was underlain by Cretaceous basalts and dolerites. Quantitative interpretation of the VLF-EMR data with a simple one-dimensional model yielded considerable detail about the weathered layer. For the granitic area, a prior estimate of at least one resistivity parameter of the ground is required. If this is not already available, a limited amount of direct-current resistivity surveying can provide the required information. A study of the EMR data from the basaltic area revealed the presence of a thin, highly conductive layer between the weathered layer and the bedrock. The parameters of this layer were found to be variable, making it necessary to use a set of diagrams for quantitative interpretation. Due to the presence of this highly conductive layer, the EMR data contain little information on the bedrock resistivity. Our field studies suggest that the VLF-EMR method can be used as a fast and inexpensive tool for mapping of the weathered layer in tropical regions with hard rock geology. Such mapping is of considerable importance because the weathered layer is an important source of groundwater.  相似文献   

8.
In seismic hazard assessment studies, near-surface geophysical prospecting is increasingly used to help in identifying active faults and to constrain seismic deformation parameters (fault offset, recurrence…). In the current study, electrical tomography was used (1) to locate the Trévaresse slow reverse fault (Provence, SE France), the source of the most important seismic event recorded in the 20th century in France and (2) to test electrical tomography's ability to detect faults and tectonic deformation in a complex geological environment. Our study is supported by a detailed geologic and geomorphic field study. In the investigated area, where the geology is made of a succession of folded marl and limestone layers, electrical images are strongly influenced by electrical macro-anisotropy which plays a critical role in this survey, depending on the investigation scale. When resolution (i.e. the smallest target that can be resolved) is larger than the layer thickness (at depth and/or for large electrode spacing), inversion is unable to retrieve the real structure. To successfully interpret these electrical images, in terms of tectonically induced dip variations, we constrain field data interpretation by comparing measured and inverted electrical data with synthetic modelling results and geological data. The synthetic models were hypothetical coarse geological structures determined on the basis of geological/geomorphologic information (e.g. outcropping layers). For shallow depths, when the resolution is smaller or equal to the anisotropy scale, an electrical tomogram yields straightforward interpretation of the subsurface structure. Coupled with numerical modelling, electrical results successfully helped in positioning a trench for fault characterization.  相似文献   

9.
A sequence of glacial and alluvial deposits overlying the Cretaceous Chalk in Eastern England was characterised using two geophysical techniques: electrical resistivity imaging and electromagnetic (EM) induction. Extensive geological data were available from trenching and boreholes. Synthetic modelling of the electrical resistivity imaging technique was undertaken to identify its limitations and to optimise survey design. The EM induction method provided a quick and cost-effective reconnaissance technique for identifying large-scale lateral variation in lithology, and for siting resistivity profiles and further boreholes. The resistivity imaging technique provided detailed information on the vertical continuity of permeable units, and was able to identify permeable pathways through the sequence. Certain limitations in detecting thin sand or gravel layers underlying electrically conductive clay were seen in both the synthetic and field data. Nevertheless, the study shows that knowledge of these limitations allows interpretation for the purpose of groundwater vulnerability assessment, given that an appropriate amount of invasive investigation has been conducted.  相似文献   

10.
A sequence of glacial and alluvial deposits overlying the Cretaceous Chalk in Eastern England was characterised using two geophysical techniques: electrical resistivity imaging and electromagnetic (EM) induction. Extensive geological data were available from trenching and boreholes. Synthetic modelling of the electrical resistivity imaging technique was undertaken to identify its limitations and to optimise survey design. The EM induction method provided a quick and cost-effective reconnaissance technique for identifying large-scale lateral variation in lithology, and for siting resistivity profiles and further boreholes. The resistivity imaging technique provided detailed information on the vertical continuity of permeable units, and was able to identify permeable pathways through the sequence. Certain limitations in detecting thin sand or gravel layers underlying electrically conductive clay were seen in both the synthetic and field data. Nevertheless, the study shows that knowledge of these limitations allows interpretation for the purpose of groundwater vulnerability assessment, given that an appropriate amount of invasive investigation has been conducted.  相似文献   

11.
The many different existing models describing the spectral behaviour of the resistivity of geological materials at low frequency, combined with the lack of available field data, render the interpretation of complex resistivity (CR) data very difficult. With a recent interest in CR‐measurements for environmental applications and thanks to technological progress, the use of wide‐band frequency equipment seems promising, and it is expected to shed light on the different results among the published solutions to the electromagnetic (EM) coupling problem. We review the theory of EM‐coupling over a homogeneous half‐space with CR‐effects and study some aspects of the complex coupling function. We advocate the use of the CR‐based coupling function in the interpretation process, in order to obtain a better understanding of the physical processes involved in CR‐effects. Application of the model to real field data shows systematic good agreement in two simple cases, even over wide ranges of frequencies. Interpretation with a double Cole–Cole model is applied for comparison, and in spite of good fits to the data, large differences are observed in the interpreted low‐frequency dispersion. We conclude that the use of a second Cole–Cole model to describe EM‐coupling may corrupt the interpretation of the low‐frequency dispersion, even when only the normal range of frequencies (<100 Hz) is considered, and that the use of the actual EM‐coupling expression is essential when the goal is a better understanding of interaction between CR‐effects and EM‐coupling.  相似文献   

12.
Mapping of fault patterns is an important part of geophysical exploration. A computerized digital template analysis method is described which tests gravity maps for the effects of faults by comparing measured gravity data with calculated master curves. The interpreted gravity data are incorporated in a tectonic map using geological symbols and units, for ready use by the geologists. Tectonics can be investigated by gravity detailing if the smallest undulations on the Bouguer map are taken into consideration. Up to now, residual and derivative gravity maps have explained gravity effects by assuming spherical bodies which are almost unknown in geology. The method discussed here uses tectonic elements, such as fault-blocks and dikes as a basis of interpretation instead of the spherical bodies of the conventional interpretation methods. Gravity data can be easily and relatively cheaply obtained in the early phases of exploration by area wide spot coverage following lines of easy access such as roads, etc. Seismic studies, by contrast, have the disadvantage of being executed along predetermined profile lines. Only after sufficient detailing do seismic profiles permit-if at all -an areal mapping of faults. Thus a tectonic map is obtained only at the end of a geophysical survey instead of being available prior to the planning of costly seismic profiles. The use of gravimetric data and their interpretation by the suggested method provides tectonic detail maps in the early phases of geophysical studies. In addition, this article discusses a general geophysical interpretation method, using the investigation of faults by gravity as an example. Applications of this method for different gravity and magnetics problems as well as for combined interpretations are outlined. Detailed case histories will be published in later articles.  相似文献   

13.
The Millennium uranium deposit is located within the Athabasca Basin of northern Saskatchewan. The basement rocks, comprised primarily of paleo‐Proterozoic gneisses, are electrically resistive. However, the deposit is associated with highly conductive graphitic metasediments that are intercalated with the gneisses. An unconformity separates the basement rocks from the overlying, horizontally stratified, Proterozoic sandstones of the Athabasca Group (which are also highly resistive). The strike extents of the graphitic metasedimenary packages are extensive and therefore electromagnetic (EM) survey techniques are successful at identifying these zones but do not identify the specific locations where they are enriched in uranium. Through drilling it has been noted that hydrothermal processes associated with mineralization has altered the rocks in the vicinity of the deposits, which should in theory result in a resistivity low. A significant resistivity low has been mapped coincident with the Millennium deposit using ground resistivity survey techniques. However, a comparison of the airborne EM and ground resistivity results reveals that the two data sets have imaged different features. The resistive‐limit (on‐time) windows of the MEGATEM data show conductive features corresponding to lakes located to the west and south of the deposit. The late‐time windows show a feature to the east of the deposit, interpreted as being associated with the east‐dipping graphitic basement conductors (similar to that observed in historical ground EM data collected in this area). The early‐time TEMPEST windows (delay times less than 0.2 ms) show a broad resistivity low located at approximately the same location as where the alteration has been identified through drilling. Modelling the data is not easy but a response that decays prior to 0.3 ms is consistent with 500 Ωm material in the sandstone, a resistivity value close to the lower limit with respect to the hydrothermally altered Athabasca group sediments in this area. The MEGATEM system does not see a conductive zone over the alteration as clearly but the high signal‐to‐noise ratio in the late‐time MEGATEM data means that the conductive material at a greater depth is more coherently imaged.  相似文献   

14.
The Yurihara oil and gas field is located on the southern edge of Akita Prefecture, northeastern Japan. In this area, drilling, surface geological surveys and many seismic surveys have been used to investigate the geological structure. Wells drilled into the Nishikurosawa Basalt Group (NBG) of Miocene age found oil and gas reservoirs at depths of 1.5–2 km. Oil and gas are now being produced commercially and further exploration is required in the surrounding areas. However, since the neighbouring areas are covered with young volcanic products from the Chokai volcano, and have a rough topography, the subsurface distribution of the NBG must be investigated using other methods in addition to seismic reflection. According to the well data, the resistivity of the NBG is comparatively higher than that of the overlying sedimentary formations, and therefore the magnetotelluric (MT) method is expected to be useful for the estimation of the distribution of the NBG. An MT survey was conducted along three survey lines in this area. Each line trended east–west, perpendicular to the regional geological strike, and was composed of about 25 measurement sites. Induction vectors evaluated from the magnetic field show that this area has a two-dimensional structure. The evaluated resistivity sections are in agreement with the log data. In conclusion, we were able to detect resistive layers (the NBG) below conductive layers. The results indicate that the NBG becomes gradually less resistive from north to south. In the centre of the northern line, an uplifted resistive area is interpreted as corresponding to the reservoir. By comparison with a seismic section, we prove the effectiveness of the integration of seismic and MT surveys for the investigation of the morphology and internal structure of the NBG. On other survey lines, the resistive uplifted zones are interpreted as possible prospective areas.  相似文献   

15.
航空电磁拟三维模型空间约束反演   总被引:1,自引:0,他引:1       下载免费PDF全文
为了克服时间域航空电磁数据单点反演结果中常见的电阻率或层厚度横向突变造成数据难以解释的问题,通过引入双向约束实现航空电磁拟三维空间约束反演.除考虑沿测线方向相邻测点之间的横向约束外,同时还考虑了垂直测线方向测点在空间上的相互约束.为此,首先设计拟三维模型中固定层厚和可变层厚两种空间约束反演方案,然后通过在目标函数中引入沿测线和垂直测线方向上的模型参数约束矩阵,并使用L-BFGS算法使目标函数最小化,获得最优拟三维模型空间反演解.基于理论模型和实测数据反演,对单点反演与两种空间约束反演方案的有效性进行比较,证明本文空间约束反演算法对于噪声的压制效果好,反演的界面连续光滑,同时内存需求和反演时间少,是一种快速有效的反演策略.  相似文献   

16.
The island of Socorro is located in the eastern Pacific Ocean, 650 km off the coast of Mexico. It is a rare example of an oceanic volcanic island whose above sea level volume is made up mostly of peralkaline trachytes and rhyolites, with subordinate mafic rocks. Subaerial volcanism started several hundred thousand years ago and continues until recent times. We present an investigation of surface and subsurface geology of the island, based on the first detailed extensive geophysical survey on the island. Acquired airborne magnetic and gamma-ray data were compared to existing geological information and supplemented with field investigations and satellite imagery. Magnetic data show a wide minimum in the central part of the island, possibly connected to a high-temperature zone in the deeper central portion of the volcano, likely to be due to a still hot magma body. The data also depict two parallel edges possibly suggesting the existence of a nested caldera. Analysis on upward continued magnetic data by recent imaging techniques highlighted two deep sources located around 5 km b.s.l., interpreted as feeding structures that are now filled with crystalline rocks. Gamma-ray data have been interpreted through integration with the geological survey results. Several previously known volcanic deposits have been identified based on radioelement distribution, and others have been redefined based on field evidence. A new succession of volcanic members is proposed, to be verified through more detailed geological mapping, geochemical analyses of rock samples and radiometric dating.  相似文献   

17.
Archaeological indications near Mugardos (Ferrol, NW Spain) suggest the existence of a Roman settlement. In fact, in the area were found pavements, walls with north-south and east-west orientations and some structures that endured heating. These remains are covered by soil, more than 1 m thick, and lie over schists. In order to determine the archaeological potential of the area and to delimit future excavations a geophysical survey, consisting of a joint resistivity and magnetic survey, was planned and carried out. The square array of electrodes was used and the data are discussed as apparent square array resistivity maps and azimuthal inhomogeneity ratio (AIR) maps. The magnetic survey included total field measurements using sensor heights of 0.30 and 2.30 m above the ground, so that a magnetic gradient could be computed.A combined interpretation of both resistivity and magnetic data is discussed. Later excavations have confirmed the geophysical interpretation.  相似文献   

18.
Helicopter-borne frequency-domain electromagnetic (HEM) surveys are used for fast high-resolution, three-dimensional resistivity mapping. Standard interpretation tools are often based on layered earth inversion procedures which, in general, explain the HEM data sufficiently. As a HEM system is moved while measuring, noise on the data is a common problem. Generally, noisy data will be smoothed prior to inversion using appropriate low-pass filters and consequently information may be lost.For the first time the laterally constrained inversion (LCI) technique has been applied to HEM data combined with the automatic generation of dynamic starting models. The latter is important because it takes the penetration depth of the electromagnetic fields, which can heavily vary in survey areas with different geological settings, into account. The LCI technique, which has been applied to diverse airborne and ground geophysical data sets, has proven to be able to improve the HEM inversion results of layered earth structures. Although single-site 1-D inversion is generally faster and — in case of strong lateral resistivity variations — more flexible, LCI produces resistivity — depth sections which are nearly identical to those derived from noise-free data.The LCI results are compared with standard single-site Marquardt–Levenberg inversion procedures on the basis of synthetic data as well as field data. The model chosen for the generation of synthetic data represents a layered earth structure having an inhomogeneous top layer in order to study the influence of shallow resistivity variations on the resolution of deep horizontal conductors in one-dimensional inversion results. The field data example comprises a wide resistivity range in a sedimentary as well as hard-rock environment.If a sufficient resistivity contrast between air and subsurface exists, the LCI technique is also very useful in correcting for incorrect system altitude measurements by using the altitude as a constrained inversion parameter.  相似文献   

19.
Helicopter electromagnetic (HEM) systems are commonly used for conductivity mapping and the data are often interpreted using an isotropic horizontally layered earth model. However, in regions with distinct dipping stratification, it is useful to extend the model to a layered earth with general anisotropy by assigning each layer a symmetrical 3 × 3 resistivity tensor. The electromagnetic (EM) field is represented by two scalar potentials, which describe the poloidal and toroidal parts of the magnetic field. Via a 2D Fourier transform, we obtain two coupled ordinary differential equations in the vertical coordinate. To stabilize the numerical calculation, the wavenumber domain is divided into two parts associated with small and large wavenumbers. The EM field for small wavenumbers is continued from layer to layer with the continuity conditions. For large wavenumbers, the EM field behaves like a DC field and therefore cannot be sensed by airborne EM systems. Thus, the contribution from the large wavenumbers is simply ignored. The magnetic fields are calculated for the vertical coaxial (VCX), horizontal coplanar (HCP) and vertical coplanar (VCP) coil configurations for a helicopter EM system. The apparent resistivities defined from the VCX, VCP and HCP coil responses, when plotted in polar coordinates, clearly identify the principal anisotropic axes of an anisotropic earth. The field example from the Edwards Aquifer recharge area in Texas confirms that the polar plots of the apparent resistivities identify the principal anisotropic axes that coincide well with the direction of the underground structures.  相似文献   

20.
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.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号