Interpretation of transient electromagnetic soundings over three-dimensional structures for the central-loop configuration

Summary. An assessment is made of the bias of fitting constrained layered-earth models to transient electromagnetic data obtained over 3-D structures. In this assessment we use the central-loop configuration and show that accurate estimates of the depth of burial of 3-D structures can be obtained with layered-earth model fitting. However, layered-earth interpretations are not reliable for estimating depth extents and resistivities of 3-D structures. When layered earths are used for interpretation, it is advantageous in some cases to use data based on the magnetic field instead of the voltage. A magnetic-field definition of apparent resistivity, in contrast to a definition based on the voltage, eliminates apparent-resistivity overshoots and undershoots in the data. A resistivity undershoot in the data can produce an extraneous and misleading layer in an interpretation of a 3-D resistive structure. Due to 3-D effects, apparent-resistivity soundings (magnetic field and voltage) may rise so steeply at late times that it may not be possible to fit a sounding to a reasonable layered-earth model. Truncating such a sounding, over a buried conductor, allows for a reasonable layered-earth fit and an accurate estimate of the depth to the conductor. However, the resistivity of the conductor is overestimated.
Measurements of the horizontal field in the central-loop configuration can map 3-D structures, provided the sensor is located accurately at the centre of the transmitting loop. Horizontal-field calculations show that the transients peak on the flanks of a 3-D structure, but are depressed over the structure's centre. Weak transient responses flanked by two large transient responses, which are opposite in sign, locate the structure. The sign reversal is caused by a corresponding reversal in the currents that are channelled through or deflected away from conductive or resistive structures, respectively.     

Numerical estimation of REV and permeability tensor for fractured rock masses by composite element method

The Monte Carlo method is used to generate parent stochastic discrete fracture network, from which a series of fractured rock samples of different sizes and orientations are extracted. The fracture network combined with a regular grid forms composite element mesh of the fractured rock sample, in which each composite element is composed of sub‐elements incised by fracture segments. The composite element method (CEM) for the seepage is implemented to obtain the nodal hydraulic potential as well as the seepage flow rates through the fractured rock samples. The application of CEM enables a large quantity of stochastic tests for the fractured rock samples because the pre‐process is facilitated greatly. By changing the sizes and orientations of the samples, the analysis of the seepage characteristics is realized to evaluate the variation of the permeability components, the existence of the permeability tensor and the representative element volume. The feasibility and effectiveness are illustrated in a numerical example. Copyright © 2008 John Wiley & Sons, Ltd.     

Vertical variations of wave-induced radiation stress tensor

INTRODUcrIONThe concept of radiation stress was deve1oPed by tonguet--Higgins and Stewart (1964 ),who intreduced the definition of radiation stress as the excess mornentum due to the presence ofwaves, on the basis of time-averaged laws of Newtonian fluid mechanics and the assmption ofa unifOrm velocity distribution over depth. Subequently, the theory has been applied success-fully in the investigation of phenomena such as wave set-up and set--down (Bowen et al.,l968), longshore currents …     

滇西地区地壳上地幔电性结构与地壳构造活动的关系

本文根据滇西地区18个大地电磁测深点资料的数据处理和分析结果,对测区深部导电率在纵、横向上的变化特征进行了研究。结果表明:滇西地区深部电性为多层结构,大致可分四至五个电性结构层;深部电性结构横向变化大,明显受区域构造控制;该区上部地壳内普遍存在低阻层;上地幔高导层明显存在两个隆起区,一个以剑川—鹤庆为中心呈北北西向展布的隆起区,另一个以腾冲—潞西为轴呈南北向展布的隆起区。 本文还讨论了地壳上地幔电性结构与大地构造的关系,滇西北裂陷区盆地的形成,以及该区地震活动与深部构造的关系     

频率域电磁法中视电阻率全区定义

视电阻率定义在电法勘探中得到广泛的应用。目前在频率域电磁法中常常采用远区视电阻率定义。由于在频率域中收发距并非足够大,以致在低频段远区定义视电阻率往往不能反映地下电阻率值。我们定义全区电阻率及引入校正系数K,经二层三层模型及野外实测资料的计算,证明了方法的有效性。     

应用电性特征研究银川平原地下水水质空间分布规律

以蕴涵丰富地下水系统水循环信息的电性特征为线索,采用电测深法,在分析研究银川平原450 m深度内地下水系统的空间结构的基础上,对平原地下水水质空间分布规律进行了综合研究,为银川平原地下水资源可持续利用战略的制订提供了依据。     

华北地台地热勘探中利用大地电磁测深法(MT)圈定热田范围和探讨地热形成机制

利用大地电磁测深法（MT）反映地质体的电性差异,首先可以确定地质体的空间展布、地质结构和构连;而由于地质体电阻率受地热水影响而降低,进而据此圈定地热水的分布范围;并在此基础上可以进一步探讨地热形成机制。在华北地台地热勘探中MT法取得了较好的勘探效果。     

Simulations of 3D DC borehole resistivity measurements with a goal-oriented <Emphasis Type="Italic">hp</Emphasis> finite-element method. Part II: through-casing resistivity instruments

We simulate direct current (DC) borehole resistivity measurements acquired in steel-cased deviated wells for the assessment of rock formation properties. The assumed data acquisition configuration considers one current (emitter) and three voltage (collector) electrodes that are utilized to measure the second difference of the electric potential along the well trajectory. We assume a homogeneous, 1.27-cm-thick steel casing with resistivity equal to 10^− 5 Ω· m. Simulations are performed with two different numerical methodologies. The first one is based on transferring two-dimensional (2D) axisymmetric optimal grids to a three-dimensional (3D) simulation software. The second one automatically produces optimal 3D grids yielded by a 3D self-adaptive goal-oriented algorithm. Both methodologies utilize high-order finite elements (FE) that are specially well-suited for problems with high-contrast coefficients and rapid spatial variations of the electric field, as it occurs in simulations that involve steel-cased wells. The method based on transferring 2D-optimal grids is efficient in terms of CPU time (few seconds per logging position). Unfortunately, it may produce inaccurate 3D simulations in deviated wells, even though the error remains below 1% for the axisymmetric (vertical) well. The method based on optimal 3D grids, although less efficient in terms of CPU time (few hours per logging position), produces more accurate results that are validated by a built-in a posteriori error estimator. This paper provides the first existing simulations of through-casing resistivity measurements in deviated wells. Simulated resistivity measurements indicate that, for a 30° deviated well, measurements in conductive layers 0.01 Ω· m) are similar to those obtained in vertical wells. However, in resistive layers (10,000 Ω· m), we observe 100% larger readings in the 30° deviated well. This difference becomes 3,000% for the case of a 60° deviated well. For this highly-deviated well, readings corresponding to the conductive formation layer are about 30% smaller in magnitude than those in a vertical well. Shoulder effects significantly vary in deviated wells.     

Rupture process of large earthquakes in the northern Mexico subduction zone

The Cocos plate subducts beneath North America at the Mexico trench. The northernmost segment of this trench, between the Orozco and Rivera fracture zones, has ruptured in a sequence of five large earthquakes from 1973 to 1985; the Jan. 30, 1973 Colima event (M _s 7.5) at the northern end of the segment near Rivera fracture zone; the Mar. 14, 1979 Petatlan event (M _s 7.6) at the southern end of the segment on the Orozco fracture zone; the Oct. 25, 1981 Playa Azul event (M _s 7.3) in the middle of the Michoacan gap; the Sept. 19, 1985 Michoacan mainshock (M _s 8.1); and the Sept. 21, 1985 Michoacan aftershock (M _s 7.6) that reruptured part of the Petatlan zone. Body wave inversion for the rupture process of these earthquakes finds the best: earthquake depth; focal mechanism; overall source time function; and seismic moment, for each earthquake. In addition, we have determined spatial concentrations of seismic moment release for the Colima earthquake, and the Michoacan mainshock and aftershock. These spatial concentrations of slip are interpreted as asperities; and the resultant asperity distribution for Mexico is compared to other subduction zones. The body wave inversion technique also determines theMoment Tensor Rate Functions; but there is no evidence for statistically significant changes in the moment tensor during rupture for any of the five earthquakes. An appendix describes theMoment Tensor Rate Functions methodology in detail.The systematic bias between global and regional determinations of epicentral locations in Mexico must be resolved to enable plotting of asperities with aftershocks and geographic features. We have spatially shifted all of our results to regional determinations of epicenters. The best point source depths for the five earthquakes are all above 30 km, consistent with the idea that the down-dip edge of the seismogenic plate interface in Mexico is shallow compared to other subduction zones. Consideration of uncertainties in the focal mechanisms allows us to state that all five earthquakes occurred on fault planes with the same strike (N65°W to N70°W) and dip (15±3°), except for the smaller Playa Azul event at the down-dip edge which has a steeper dip angle of 20 to 25°. However, the Petatlan earthquake does prefer a fault plane that is rotated to a more east-west orientation—one explanation may be that this earthquake is located near the crest of the subducting Orozco fracture zone. The slip vectors of all five earthquakes are similar and generally consistent with the NUVEL-predicted Cocos-North America convergence direction of N33°E for this segment. The most important deviation is the more northerly slip direction for the Petatlan earthquake. Also, the slip vectors from the Harvard CMT solutions for large and small events in this segment prefer an overall convergence direction of about N20°E to N25°E.All five earthquakes share a common feature in the rupture process: each earthquake has a small initial precursory arrival followed by a large pulse of moment release with a distinct onset. The delay time varies from 4 s for the Playa Azul event to 8 s for the Colima event. While there is some evidence of spatial concentration of moment release for each event, our overall asperity distribution for the northern Mexico segment consists of one clear asperity, in the epicentral region of the 1973 Colima earthquake, and then a scattering of diffuse and overlapping regions of high moment release for the remainder of the segment. This character is directly displayed in the overlapping of rupture zones between the 1979 Petatlan event and the 1985 Michoacan aftershock. This character of the asperity distribution is in contrast to the widely spaced distinct asperities in the northern Japan-Kuriles Islands subduction zone, but is somewhat similar to the asperity distributions found in the central Peru and Santa Cruz Islands subduction zones. Subduction of the Orozco fracture zone may strongly affect the seismogenic character as the overlapping rupture zones are located on the crest of the subducted fracture zone. There is also a distinct change in the physiography of the upper plate that coincides with the subducting fracture zone, and the Guerrero seismic gap to the south of the Petatlan earthquake is in the wake of the Orozco fracture zone. At the northern end, the Rivera fracture zone in the subducting plate and the Colima graben in the upper plate coincide with the northernmost extent of the Colima rupture zone.