Inversão gravimétrica 3D e interpretação integrada dos métodos CSAMT, gravimétrico e magnetométrico no estudo do kimberlito Régis, MG, Brasil

Os diamantes são minerais raros e de alto valor econômico. Estes minerais se formam em condições mantélicas, numa profundidade aproximada de 150 a 200 km e ascendem à superfície englobada em magmas alcalinos e carbonatíticos, nos kimberlitos e lamproitos (depósitos primários). No presente trabalho objetivou-se o estudo da geometria do Kimberlito Régis, o qual pertence à província alcalina do Brasil Meridional, localizado no estado de Minas Gerais, na cidade de Carmo do Paranaíba. Para este estudo foi realizada uma inversão gravimétrica 3D e uma modelagem utilizando dados de gravimetria e os resultados dos estudos de magnetometria (Menezes e La Terra, 2011) e CSAMT (La Terra e Menezes, 2012). As anomalias Bouguer e magnética foram então modeladas a partir de um resultado CSAMT, apresentando como resultado um conduto vulcânico com cerca de 2800 metros de profundidade e 200 metros de largura. Já o resultado da inversão mapeou um corpo com forma de cone invertido, apresentando contrastes de densidade negativos. A metodologia de integração de métodos geofísicos mostrou-se eficiente para estudos exploratórios de kimberlitos, apresentando rapidez e baixo custo quando comparados com os métodos tradicionais.

有源电磁扩散场克希霍夫积分偏移

As active electromagnetic method, field data of CSAMT method follow the equation of diffusion. Propagting in solid earth media, diffusion EM signal has strong attenuation and dispersion, otherwise seismic wave shows weak attenuation and dispersion, therefore the resolution power of CSAMT method is not better than seismic reflection method. However, there is consistence and similarity between EM signal and seismic wave in wave equation, we can apply Kirchhoff integral migration technique, a proven one in seismic method in time domain, to carry out seduo-seismic processing for CSAMT signal in frequency domain so that the attenuation and dispersion could be made compensated in some extent, and the resolution power and interpretation precision of active EM wave could be improved. Satisfying passive homogeneous Helmholtz quation, we proceed with Green theorem and combine the active inhomogenous Helmholtz quation, the Kirchhoff integral formula could be derived. Given practical problems, if we only consider the surface integral value, and assume that the intergral value in other interface is zero, combined with Green theorem in uniform half space, the expression could be simplified, and we can obtain frequency-domain Kirchhoff integral formula in surface, which is also called downward continuation of EM field in frequency domain. With image conditions and energy compensation considered, in order to get image conditions in time domain Fourier inverse transformation in frequency domain can be performed, so we can formulate the active Kirchhoff integral migration expression. At first, we construct relative stratified model, with different frequency series taken into account, then we change the distances between transmitter and reciever, the EM response can be obtained. Analyzing the EM properties, we can clarify near and far zone that can instruct us to carry out transmitter layout in practical application. Combined with field data surveyed in far zone, We perform Kirchhoff integral migration and compare the results with model to interpret. Secondly, with far field EM data, we apply TM mode to get EM response of given 2D model, then apply Kirchhoff integral migration on modelling data and interpret the results.

大型反倾岩体变形破坏机制及稳定性分析———以金沙江宽谷坝址龙蟠边坡为例

A large number of catastrophic accidents were aroused by the instability and destruction of anti-dip rock masses in the worldwide engineering projects, such as hydropower station, mine, railways and so on. Problems in relation to deformation and failure about anti-dip rock slopes are significant for engineering geology research. This dissertation takes the Longpan slope in the Jinsha River as a case to study the deformation mechanism of large-scale anti-dip rock masses and the slope stability analysis method. The primary conclusions are as follows. The Dale Reach of Jinsha River, from Longpan to the debouchment of Chongjiang tributary, is located in the southeastern margin of the Qinghai-Tibet Plateau. Longpan slope is the right embankment of Dale dam, it is only 26 km to the Shigu and 18 km to Tiger Leaping Gorge. The areal geology tectonic structures here area are complicated and blurry. Base on the information of geophysical exploration (CSAMT and seismology) and engineering geological investigation, the perdue tectonic pattern of Dale Reach is put forward for the first time in this paper. Due to the reverse slip of Longpan fault and normal left-rotation of Baihanchang fault, the old faulted valley came into being. The thick riverbed sediments have layered characters of different components and corresponding causes, which attribute to the sedimentary environments according with the new tectonic movements such as periodic mountain uplifting in middle Pleistocene. Longpan slope consists of anti-dip alternate sandstone and slate stratums, and the deformable volume is 6.5×107m3 approximately. It was taken for an ancient landslide or toppling failure in the past so that Dale dam became a vexed question. Through the latest field surveying, displacement monitoring and rock masses deforming characters analyses, the geological mechanism is actually a deep-seated gravitational bending deformation. And then the discrete element method is used to simulate the deforming evolution process, the conclusion accords very well with the geo-mechanical patterns analyses. In addition strength reduction method based on DEM is introduced to evaluate the factor of safety of anti-dip rock slope, and in accordance with the expansion way of the shear yielding zones, the progressive shear failure mechanism of large-scale anti-dip rock masses is proposed for the first time. As an embankment or a close reservoir bank to the lower dam, the stability of Longpan slope especially whether or not resulting in sliding with high velocity and activating water waves is a key question for engineering design. In fact it is difficult to decide the unified slip surface of anti-dip rock slope for traditional methods. The author takes the shear yielding zones acquired form the discrete element strength reduction calculation as the potential sliding surface and then evaluates the change of excess pore pressure and factor of stability of the slope generated by rapid drawdown of ponded water. At the same time the dynamic response of the slope under seismic loading is simulated through DEM numerical modeling, the following results are obtained. Firstly the effective effect of seismic inertia force is resulting in accumulation of shear stresses. Secondly the discontinuous structures are crucial to wave transmission. Thirdly the ultimate dynamic response of slope system takes place at the initial period of seismic loading. Lastly but essentially the effect of earthquake load to bringing on deformation and failure of rock slope is the coupling effect of shear stresses and excess pore water pressure accumulation. In view of limitations in searching the critical slip surface of rock slope of the existing domestic and international software for limit equilibrium slope stability analyses, this article proposes a new method named GA-Sarma Algorithm for rock slope stability analyses. Just as its name implies, GA-Sarma Algorithm bases on Genetic Algorithm and Sarma method. GA-Sarma Algorithm assumes the morphology of slip surface to be a broken line with traceability to extend along the discontinuous surface structures, and the slice boundaries is consistent with rock mass discontinuities such as rock layers, faults, cracks, and so on. GA-Sarma Algorithm is revolutionary method that is suitable for global optimization of the critical slip surface for rock slopes. The topics and contents including in this dissertation are closely related to the difficulties in practice, the main conclusions have been authorized by the engineering design institute. The research work is very meaningful and useful for the engineering construction of Longpan hydropower station.

含激电信息的人工源频率域电磁法正演研究——大功率人工固定源极低频电磁法激电效应预研究

ABSTRACT Base on Dias-model, this paper study the forward modeling to IP response with MT souece and CSAMT source with the author’s program，that is to provide the theoretical basis for extracting the IP information from the far field, near field and transition field in extremely low frequency electromagnetic sounding with artificial strong fixed source signal. The outline of the paper is as follows. Firstly, review the history of the complex-resistivity-model, and analyses the bases for choosing Dias model; meanwhile, effects and responses of each parameter in Dias model are analyses. Afterward, study the forward modeling to 1D layered model with MT source, numerically simulate the IP effects of classic geoelectric objects; Find the clear anomaly and relationship between the peak value of the amplitude anomaly ratio and phase anomaly difference with and without IP parameters to be considered within a frequency range. On the basis of the modeling of MT, Study the CSAMT modeling with a dipole-source, obtain the anomaly responses and the relationship between anomaly and characters of object. Base on the infinite line source, study the 2D IP effect of geoelectric objects with 2D modeling. Calculate the response with different source distance, object depth and different wall rock resistivity systemically. Finally, conclude the viewpoints and give the discussion of the result, and point out the lack of this research also.

计算层状介质中电磁场的层矩阵法

As an important branch of electrical prospecting method, the artificial source frequency domain electromagnetism method has received more and more attention. But when conducts the fundamental research, people often isolated study some concrete method, so the research results of one method are very difficult to apply to another method directly. This article will possess the artificial source frequency domain EM method to an 1D model simply. It is stratified medium model, with an electric or magnetic source in or outside of it. Then take the horizontal electric dipole source as an example to introduce how to computing the EM field in stratified medium. Because layer matrix is the key of establishing equations, so we call it the layer-matrix method. The key of layer-matrix method is establishing equations by using layer matrixes in wavenumber(kx, ky, z) domain, then obtains the electromagnetic field value of wavenumber domain. After Fourier transform, we can get electromagnetic field of any position in spatial domain. The layer matrix technique theoretically can calculate electromagnetic field of any position for any source, is suitable for many kinds of electromagnetic method. After introduction of the layer matrix method, this article has done some CSAMT, MCSEM and Wireless Electro-Magnetic Method (WEM) modeling with layer matrix method separately. In CSAMT modeling, we get electromagnetic field dissemination characteristics considering wave number of the air, and obtain three-dimensional distribution characteristics of the electromagnetic field. In MCSEM modeling, we get electromagnetic field dissemination characteristics with and without considering the airwave, and obtain three-dimensional distribution characteristics of electromagnetic field. In WEM modeling, we get electromagnetic field’s difference between considering the ionosphere and not considering it, and recognize the ionosphere’s influence of electromagnetic field. With the layer matrix technique, we have got some new understandings of EM dissemination rules of different situations. All analysis results indicate that the layer-matrix technique is credible and effective, and are worthy of further thorough research and development.

胶东蓬家夼式金矿床矿床成因模型与找矿模型研究

The dynamic environments of mineralization in Mesozoic Jiaodong gold mine concentrated area can be devided into two types, compressive environment which related to intracontinental collision and extensional environment which related to intracontinental volcanic rift. The altered rock type (Jiaojia type) and quartz vein type (Linglong type) which related to the former one, were discovered for several years, and became the main types of gold deposits in recent years. A new type gold deposit, syn-detachment altered tectonic breccia type gold deposit, such as Pengjiakuang gold deposit and Songjiagou gold deposit has been discovered on the northeastern margin of Jiaolai Basin. In this paper, the new type of gold deposit has been studied in detail. The study area is located at the northeastern boundaries of Jiaolai Basin, and between the Taocun-Jimo Fault and Wji-Haiyang Fault, in the eastern part of the Jiaodong Block. Pengjiakuang gold deposit and Songjiagou gold deposit occur in a arc-shape detachment fault zone between conglomerate of Lower Cretaceous Laiyang Formation and metamorphic complex of Lower Proterozoic Jingshan Group. Regional geological studies show that Kunyuanshan and Queshan granite intrusions and Qingshanian volcanism were formed in different period of lithospheric thinning of East China in Mesozoic. Granite intrusions were formed in compressive environment, while Qingshanian volcanism were formed in extensional environment. They are all related to the detachment of Sulu Orogenic Belt and the sinistral motion of Tanlu Fault. The Pengjiakuang detachment systems which were formed in the the sinistral motion of Tanlu Fault are the important ore-controlling and ore-containing structure. The Pengjiakuang type gold deposit, controlled by detachment structure, was formed before Yanshanian volcanic period concerning with mixture of meteoric water and magmatic water found in fluid inclusions of gold ores. The minerogenetic epoch has been proposed in 90～120Ma. the host rocks have been extensively subjected to pyritization, silicification, sericitization and carbonatization. Individual ore-body has maximum length of 800m, oblique extension of 500～700m and gold grade of 1～43 * 10~(-6). Native gold is disseminated in silicified, phyllic or carbonatized tectonic breccia. Sulfur, carbon and lead isotope studies on gold ores and wall rocks show that the sulfur come from the metamorphic complex of Lower Proterozoic Jingshan Group, carbon comes from the marble in Jingshan Group, while a part of lead comes from the mantle. The mineralizing fluid is rich in Na~+ and Cl~-, but relatively impoverished in K~+ and F~-. According to the date from hydrogen and oxygen isotopic compositions (δ~(18)OH_2O = 0.59%～4.03%, δDH_2O = -89.5%～97.9%), the conclusion can be reached that the mineralizing fluid of Pengjiakuang gold deposit was a kind of mixed hydrothermal solution which was mainly composed of meteoric water and magmatic water. A genetical model has been formulated. Some apparent anomaly features which show low in the central part and high in the both sides corresponding to the gold-bearing structure, were sum up after analying a vast amount of date by prospecting the orebodies using gamma-ray spectrometer, electrogeochemical parameter technique, controlled source audio magnetic telluric (CSAMT) and shallow surface thermometry in Pengjiakuang gold deposit. The location forecasting problem of buried orebodies has been solved according to these features, and the successful rate is very high in well-drilling. The structural geological-geophysical-geochemical prospecting model has been formulated on the base of the study of geological, geophysical and geochemical characteristics of Pengjiakuang type gold deposit, and the optimum combinational process of geophysical and geochemical prospecting techniques has been summed up. A comparative study shows that the Pengjiakuang type gold deposit, the syn-detachment altered tectonic breccia type gold deposit, is different from Jiaojia type gold deposits and Linglong type gold deposits, in Jiaodong Block. In general, if formed under an extensional tectonic condition and located at detachment fault zone along the margin of Mesozoic Jiaolai basin, and the gold mineralization has also close genetic relationship with alkaline magamtism. Being a new type of gold deposit in Jiaodong gold mine concentrated area, it could be potential to explore in the same regions which processed the same ore-forming geological conditions and mineralization informations.

人工源频率域大地电磁法正反演研究

The CSAMT method is playing an important role in the exploration of geothermal and the pre-exploration in tunnel construction project recently. In order to instruct the interpretation technique for the field data, the forward method from ID to 3D and inversion method in ID and 2D are developed in this paper for the artificial source magnetotelluric in frequency domain. In general, the artificial source data are inverted only after the near field is corrected on the basis of the assumption of half-homogeneous space; however, this method is not suitable for the complex structure because the assumption is not valid any more. Recently the new idea about inversion scheme without near field correction is published in order to avoid the near field correction error. We try to discuss different inversion scheme in ID and 2D using the data without near field correction.The numerical integration method is used to do the forward modeling in ID CSAMT method o The infinite line source is used in the 2D finite-element forward modeling, where the near-field effect is occurred as in the CSAMT method because of using artificial source. The pseudo-delta function is used to modeling the source distribution, which reduces the singularity when solving the finite-element equations. The effect on the exploration area is discussed when anomalous body exists under the source or between the source and exploration area; A series of digital test show the 2D finite element method are correct, the results of modeling has important significant for CSAMT data interpretation. For 3D finite-element forward modeling, the finite-element equation is derived by Galerkin method and the divergence condition is add forcedly to the forward equation, the forward modeling result of the half homogeneous space model is correct.The new inversion idea without near field correction is followed to develop new inversion methods in ID and 2D in the paper. All of the inversion schemes use the data without near field correction, which avoid introducing errors caused by near field correction. The modified grid parameter method and the layer-by-layer inversion method are joined in the ID inversion scheme. The RRI method with artificial source are developed and finite-element inversion method are used in 2D inversion scheme. The inversion results using digital data and the field data are accordant to the model and the known geology data separately, which means the inversion without near field correction is accessible. The feasibility to invert the data only in exploration area is discussed when the anomalous body exists between the source and the exploration area.

Delineation of a quick clay zone at Smørgrav, Norway, with electromagnetic methods under geotechnical constraints

In many coastal areas of North America and Scandinavia, post-glacial clay sediments have emerged above sea level due to iso-static uplift. These clays are often destabilised by fresh water leaching and transformed to so-called quick clays as at the investigated area at Smørgrav, Norway. Slight mechanical disturbances of these materials may trigger landslides. Since the leaching increases the electrical resistivity of quick clay as compared to normal marine clay, the application of electromagnetic (EM) methods is of particular interest in the study of quick clay structures.

For the first time, single and joint inversions of direct-current resistivity (DCR), radiomagnetotelluric (RMT) and controlled-source audiomagnetotelluric (CSAMT) data were applied to delineate a zone of quick clay. The resulting 2-D models of electrical resistivity correlate excellently with previously published data from a ground conductivity metre and resistivity logs from two resistivity cone penetration tests (RCPT) into marine clay and quick clay. The RCPT log into the central part of the quick clay identifies the electrical resistivity of the quick clay structure to lie between 10 and 80 O m. In combination with the 2-D inversion models, it becomes possible to delineate the vertical and horizontal extent of the quick clay zone. As compared to the inversions of single data sets, the joint inversion model exhibits sharper resistivity contrasts and its resistivity values are more characteristic of the expected geology. In our preferred joint inversion model, there is a clear demarcation between dry soil, marine clay, quick clay and bedrock, which consists of alum shale and limestone.