瞬变电磁晚期场资料的一维反演

本文研讨了由垂直磁偶源和水平电偶源所激发的瞬变电磁晚期场的垂直磁场时间偏导数及其视电阻率的性态特征。探讨了其一维反演的某些问题，研制出了对该时间域电磁响应函数的一维实用化反演程序，表明具有较好的使用效果。     

Transient electromagnetic response due to a pair of horizontal conducting sheets in a loop field

The transient electromagnetic response due to a pair of horizontal conducting sheets induced by a loop field is obtained. The response characteristics in the three different regions (above the upper sheet, between the two sheets and below the lower sheet) are briefly discussed. In the first and second regions the composite response is characterised by an overburden response asymptote during early times and a target response asymptote during late times. The transition zones are characterised by shoulder-like bends. In the third region, the composite response increases in direct proportion to the induction number, eventually merging with the target response for large values of the induction number.     

Transient deformation following the 30 January 1997 dike intrusion at Kīlauea volcano,Hawai’i

On 30 January 1997 an intrusion on Kīlauea volcano opened a new fissure within the East Rift Zone (ERZ) at Nāpau Crater, 3 km uprift from the ongoing eruptions at Pu’u ’Ō’ō. The fissure eruption lasted 22 h and opened a 5.1 km long, nearly vertical dike 1.9 m, extending from the surface to a depth of 2.4 km (Owen et al. 2000b). During the eruption, the lava pond at Pu’u ’Ō’ō drained, and eruptions ceased there. Pu’u ’Ō’ō eventually refilled in late February and eruptions resumed there on 28 March 1997. Continuous GPS data show a large transient following the 30 January 1997 dike intrusion. After lengthening 40 cm during the initial eruption, the baseline between two stations spanning the ERZ lengthened an additional 10 cm over the following 6 months. A coastal station KAEP also exhibited transient deformation, as it continued to move southward (5 cm) over the same 6-month period. The baseline between two stations spanning Kīlauea’s summit caldera contracted sharply during the eruption, but gradually recovered to slightly longer than its previous length 2 months after the intrusion. We use the extended network inversion filter (McGuire and Segall 2003) to invert continuous GPS data for volume change of a spherical pressure source under Kīlauea’s summit, opening distribution on a nearly vertical dike in the ERZ and potential slip on a decollement 9 km beneath the south flank. Following the 30 January intrusion, rift extension continued below the initial dike intrusion for the duration of the transient. Decollement slip, regardless of its assumed depth, is not required to fit the data. The modeled transient summit reinflation and rift opening patterns under Nāpau crater coincide with changes in observed behavior of Pu’u ’Ō’ō’s lava pond. Rift opening accelerated while Pu’u ’Ō’ō eruptions paused and began to decelerate after the lava pond reappeared nearly a month after the Nāpau eruption. The transient deformation is interpreted as resulting from shallow accommodation of the new dike volume.     

On the meridional circulation and balance of momentum in the Southern Ocean of POP

 The circulation of the Southern Ocean is studied in the eddy-resolving model POP (Parallel Ocean Program) by an analysis of zonally integrated balances. The TEM formalism (Transformed Eulerian Mean) is extended to include topography and continental boundaries, thus deviations from a zonally integrated state involve transient and standing eddies. The meridional circulation is presented in terms of the Eulerian, eddy-induced, and residual streamfunctions. It is shown that the splitting of the meridional circulation into Ekman and geostrophic transports and the component induced by subgrid and Reynolds stresses is identical to a particular form of the zonally integrated balance of zonal momentum. In this balance, the eddy-induced streamfunctions represent the interfacial form stresses by transient and standing eddies and the residual streamfunction represents the acceleration of the zonal current by density fluxes in a zonally integrated frame. The latter acceleration term is directly related to the surface flux of density and interior fluxes due to the resolved and unresolved eddies. The eddy-induced circulation is extremely vigorous in POP. In the upper ocean a shallow circulation, reversed in comparison to the Deacon cell and mainly due to standing eddies, appears to the north of Drake Passage latitudes, and in the Drake Passage belt of latitudes a deep-reaching cell is induced by transient eddies. In the resulting residual circulation the Deacon cell is largely cancelled and the residual advection of the zonal mean potential density is balanced by diapycnal eddy and subgrid fluxes which are strong in the upper few hundred meters but small in the ocean interior. The balance of zonal momentum is consistent with other eddy-resolving models; a new aspect is the clear identification of density effects in the zonally integrated balance. We show that the wind stress and the stress induced by the residual circulation drive the eastward current, whereas both eddy species result in a braking. Finally, we extend the Johnson–Bryden model of zonal transport to incorporate all relevant terms from the zonal momentum balance. It is shown that wind stress and induction by the residual circulation carry an eastward transport while bottom form stress and the stress induced by standing eddies yield westward components of transport. Received: 26 June 2001 / Accepted: 2 November 2001     

The main structure and tectonic features of the Chernorud-Mukhor site on the western shore of Lake Baikal from TEM and SP measurements

Transient electromagnetic (TEM), self-potential (SP) and geoelectrical mapping measurements were carried out at the Chernorud-Mukhor site in the Priolkhonje area on the western shore of the Lake Baikal. All measurements were made along several profiles across the main strike of the regional Primorsky fault. TEM measurements were carried out in a time range from a few tens of microseconds to several tens of milliseconds. The most important result of the 1D modelling of TEM soundings is the discovery of nearly horizontal boundaries that divide high resistive overlying and well conducting underlying rocks. The resistivity of the former is in the range from 100 Ωm to 1000 Ωm, while the resistivity of the latter varies from less than 1 Ωm to several tens of Ωm. This good conductive zone could also be verified by geoelectrical mapping using Schlumberger array (AB/2=100 m). Due to high conductivity of the underlying rocks only the upper boundary of the conductive layer could be determined by TEM soundings. A regional SP anomaly with amplitude of about −450 mV has also been observed above the low resistivity zone indicating the electron nature of its conductance. Geologically, the conductive zone is represented by a graphite-bearing layer within the region of archean rocks. Since that layer extends over a large area, it may be used as a key in studying structures and tectonics of the Priolkhonje area. A 1D TEM geoelectric section shows a wide, gently sloping syncline as a probable base structure of the Chernorud-Mukhor site. Neotectonic faults divide the syncline into vertically displaced blocks that form a wide complicated graben with a total amplitude of about 250 m.     

A study of transient variations in the Earth’s electromagnetic field at equatorial electrojet latitudes in western Africa (Mali and the Ivory Coast)

In the framework of the French-Ivorian participation to the IEEY, a network of 10 electromagnetic stations were installed at African longitudes. The aim of this experiment was twofold: firstly, to study the magnetic signature of the equatorial electrojet on the one hand, and secondly, to characterize the induced electric field variations on the other hand. The first results of the magnetic field investigations were presented by Doumouya and coworkers. Those of the electric field experiment will be discussed in this study. The electromagnetic experiment will be described. The analysis of the electromagnetic transient variations was conducted in accordance with the classical distinction between quiet and disturbed magnetic situations. A morphological analysis of the recordings is given, taking into consideration successively quiet and disturbed magnetic situations, with the results interpreted in terms of the characterization of external and internal sources. Particular attention was paid to the effects of the source characteristics on the induced field of internal origin, and to the bias they may consequently cause to the results of electromagnetic probing of the Earth; the source effect in electromagnetic induction studies. During quiet magnetic situations, our results demonstrated the existence of two different sources. One of these, the S_R^E source, was responsible for most of the magnetic diurnal variation and corresponded to the well-known magnetic signature of the equatorial electrojet. The other source (the S_R*^E source) was responsible for most of the electric diurnal variation, and was also likely to be an ionospheric source. Electric and magnetic diurnal variations are therefore related to different ionospheric sources, and interpreting the electric diurnal variation as induced by the magnetic field diurnal variation is not relevant. Furthermore, the magnetotelluric probing of the upper mantle at dip equator latitudes with the electromagnetic diurnal variation is consequently impossible to perform. In the case of irregular variations, the source effect related to the equatorial electrojet is also discussed. A Gaussian model of equatorial electrojet was considered, and apparent resistivities were computed for two models of stratified Earth corresponding to the average resistive structure of the two tectonic provinces crossed by the profile: a sedimentary basin and a cratonic shield. The apparent resistivity curves were found to depend significantly on both the model used and the distance to the center of the electrojet. These numerical results confirm the existence of a daytime source effect related to the equatorial electrojet. Furthermore, we show that the results account for the observed differences between daytime and night-time apparent resistivity curves. In particular, it was shown that electromagnetic probing of the Earth using the classical Cagniard-Tikhonov magnetotelluric method is impossible with daytime recordings made at dip latitude stations.     

瞬变电磁方法在探测煤矿浅层高阻采空区中的应用

过去在煤矿浅层高阻采空区探测中,普遍认为高密度电法为较好的方法,而瞬变电磁法通常被认为其对低阻反映较灵敏,对高阻反映不灵敏,而很少采用,但我国内蒙、山西等地,地表沟壑纵横,高差较大,不具备高密度电法的施工条件。因此,笔者给出一种在较小线框下进行的瞬变电磁方法,通过配合高发射频率和大发射电流,能够比较清晰的反映出地下高阻目标体,以此说明,针对地形起伏较大,不适合高密度电法施工的地方,小框瞬变电磁方法可以在一定程度上解决煤矿中不含水采空区的物探手段问题。文中参考煤矿不含水采空区埋深较浅的情况,设计了理论模型,正演计算得到的效果良好,对高阻的反映较为灵敏,并在鄂尔多斯地区某煤矿实际应用了该方法,得到的结果与开挖情况相符合,对煤矿生产起到了很好的指导作用,这表明该方法是切实有效的。     

瞬变电磁法在福建某铅锌矿井下勘查中的应用

老矿山深部、外围找矿是当前"以矿找矿"为重要方式。文章以福建某铅锌矿为例,通过坑道内开展小线框、大电流的瞬变电磁法获取电阻率异常,利用坑内钻进行验证,取得了良好的找矿效果,为老矿山勘探提供了新的有效的物探方法。     

用线性规划法求解瞬变电磁场反射系数

讨论了在水平层状介质中的反射系数的物理意义 ,推出时间域含有以反射系数序列为未知项的欠定方程组 ,提出用线性规划法求解瞬变电磁场反射系数序列 ,并给出求解思路和程序框图。对两个理论模型进行了处理 ,计算出了其反射系数。算例分析结果表明 ,算法可靠 ,程序可行。这一方法可用于瞬变电磁成像数值计算中。     

A mathematical model for predicting coupled heat and water movement in unsaturated soil

An investigation was conducted to develop a comprehensive moisture model for predicting non-isothermal moisture conditions in soils. An extensive literature review indicated that a model based on the Philip and de Vries equations for non-isothermal moisture movement and heat conduction would give the best results. By using numerical methods, the implicit finite difference approximations to the moisture movement and heat-transfer equations were programmed for computer solution of water content and temperature in the soil with time. Validation studies indicate that the moisture model can be used to predict accurately moisture conditions in the soil. The model was validated by using hydraulic data from laboratory studies conducted on soil columns compacted with AASHO A-3 and AASHO A-4 soil. The application of the moisture model to the study of non-isothermal moisture movement in the field is demonstrated. The influence of parameters such as water table depth, precipitation, and soil hydraulic properties on soil moisture content are shown by use of the moisture model. The model is shown to be applicable to a wide range of boundary conditions and that it predicts the moisture-temperature regime with time in soils utilizing climatic input data.