轻非水相液体污染过程的高密度电阻率成像法室内监测

为探讨高密度电阻率成像法监测多孔介质中轻非水相液体迁移过程的有效性,本文通过三维砂槽进行了非饱和带中轻非水相液体的污染试验,并利用高密度电阻率成像法进行了同步的动态监测.试验之后,将砂槽层层挖开,通过数码成像,获取了污染区域的实际范围与形状.结果表明,由高密度电阻率成像法圈定的污染区域在范围与形状上都与实际的结果比较接近,并可通过三维电阻率相对值的时间变化明显的看出轻非水相液体的污染过程.这说明利用高密度电阻率成像法对非饱和多孔介质中轻非水相液体的空间分布范围进行圈定并监测其迁移过程是完全可行的.     

Increasing the effectiveness of electrical resistivity tomography using γ11n configurations

A new array type, i.e., the γ_11n arrays, is introduced in this paper, in which the sequence of the current (C) and potential (P) electrodes is CPCP, and the distance between the last two electrodes is n times the distance between the first two ones and that of the second one and the third one. These arrays are called quasinull arrays because they are—according to their array and behaviour—between the traditional and null arrays. It is shown by numerical modelling that, in detecting small‐effect inhomogeneity, these configurations may be more effective than the traditional ones, including the optimized Stummer configuration. Certain γ_11n configurations—especially the γ₁₁₂, γ_113, and γ₁₁₄—produced better results both in horizontal and vertical resolution investigations. Based on the numerical studies, the γ_11n configurations seem to be very promising in problems where the anomalies are similar to the numerically investigated ones, i.e., they can detect and characterize, e.g., tunnels, caves, cables, tubes, abandoned riverbeds, or discontinuity, in a clay layer with greater efficacy than those of the traditional configurations. γ_11n measurements need less data than traditional configurations; therefore, the time demand of electrical resistivity tomography measurements can be shortened by their use.     

树干水分分布及运移的高密度电阻率成像法时空监测

在砍伐树干周围四个特定高度环形布置一定数量的金属电极,测定了对应横截面上各点视电阻率,再利用有限元反演计算,获得了树干各截面上电阻率的空间分布图像.由于树干内的水分含量和对应部分的导电特性密切相关,可以从得到的电阻率平面分布图上反映出树木内水分含量的空间分布和动态变化特征.结果表明,树干横截面上电阻率的空间分布呈现一系列的同心圆结构,且从核心部到外边缘的韧皮部电阻率有一个减小的趋势,说明截面内树木水分的含量沿径向向外有逐渐递增的趋势.时间变化过程中,各监测截面上的电阻率的分布和变化基本反映了截面水分的流逝过程,表明高密度电阻率成像法动态监测树干水分分布及变化是可行的.     

Characterization of seasonal local recharge using electrical resistivity tomography and magnetic resonance sounding

A groundwater recharge process of heterogeneous hard rock aquifer in the Moole Hole experimental watershed, south India, is being studied to understand the groundwater flow behaviour. Significant seasonal variations in groundwater level are observed in boreholes located at the outlet area indicating that the recharge process is probably taking place below intermittent streams. In order to localize groundwater recharge zones and to optimize implementation of boreholes, a geophysical survey was carried out during and after the 2004 monsoon across the outlet zone. Magnetic resonance soundings (MRS) have been performed to characterize the aquifer and measure groundwater level depletion. The results of MRS are consistent with the observation in boreholes, but it suffers from degraded lateral resolution. A better resolution of the regolith/bedrock interface is achieved using electrical resistivity tomography (ERT). ERT results are confirmed by resistivity logging in the boreholes. ERT surveys have been carried out twice—before and during the monsoon—across the stream area. The major feature of recharge is revealed below the stream with a decrease by 80% of the calculated resistivity. The time‐lapse ERT also shows unexpected variations at a depth of 20 m below the slopes that could have been interpreted as a consequence of a deep seasonal water flow. However, in this area time‐lapse ERT does not match with borehole data. Numerical modelling shows that in the presence of a shallow water infiltration, an inversion artefact may take place thus limiting the reliability of time‐lapse ERT. A combination of ERT with MRS provides valuable information on structure and aquifer properties respectively, giving a clue for a conceptual model of the recharge process: infiltration takes place in the conductive fractured‐fissured part of the bedrock underlying the stream and clayey material present on both sides slows down its lateral dissipation. Copyright © 2007 John Wiley & Sons, Ltd.     

Perspectives on monitoring resistivity changes with telluric signals at Parkfield, California: 1988–1999

An array designed to detect changes of resistivity in the upper 10 km of the crust has been in continuous operation since 1988. Long period (T=300–7200 s) telluric signals are recorded digitally on 5–16 km long dipoles. Telluric transfer functions are computed between field dipoles and two reference dipoles and then decomposed to examine daily fluctuations of those functions. Based on analysis of these daily fluctuations, stabilities of 0.1–0.3% are achieved. No intermediate term resistivity changes are seen associated with any of the 5 M_b>4.0 earthquakes that have occurred since 1989. Short term fluctuations are seen associated with most of these earthquakes, however. Minimum bounds on resistivity changes that could have caused the fluctuations in the telluric responses range from 1.7 to 13.0%. While these are within the ranges seen in other field studies, the statistical significance of these fluctuations is low because they also occur at many times that cannot be correlated to earthquakes. Comparison of the one potential resistivity change associated with a strain change results in a strain sensitivity that is unrealistically high. This possible resistivity change coincides spatially with a temporary variation in seismic travel time in May–June 1989, however. It is likely that many of the fluctuations in the telluric responses are random variations except for the one in 1989, and recording through the next M6 earthquake will be needed in order to detect a statistically significant resistivity change.     

Comparison of surface seismic sources at the CO2SINK site, Ketzin, Germany

In 2004 three seismic surface sources (VIBSIST, accelerated weight drop and MiniVib) were tested in a pilot study at the Ketzin test site, Germany, a study site for geological storage of CO₂ (EU project CO₂SINK). The main objectives of this pilot study were to 1) evaluate the response of the Ketzin site to reflection seismics, especially at the planned injection depth, 2) test different acquisition parameters and 3) use the results to guide the planning of the 3D survey. As part of these objectives, we emphasize the source performance comparison in this study. The sources were tested along two perpendicular lines of 2.4 km length each. Data were acquired by shooting at all stations (source and receiver spacing of 20 m) on both lines, allowing common‐midpoint stacked sections to be produced. The sources' signal characteristics based on signal‐to‐noise ratio, signal penetration and frequency content of raw shot records were analysed and stacked sections were compared. The results show that all three surface sources are suitable for reflection seismic studies down to a depth of about 1 km and provide enough bandwidth for resolving the geological targets at the site, i.e., the Weser and Stuttgart Formations. Near surface conditions, especially a thick weathering layer present in this particular area, strongly influence the data quality, as indicated by the difference in reflectivity and signal‐to‐noise ratio of the two common‐midpoint lines. The stacked sections of the MiniVib source show the highest frequency signals down to about 500 ms traveltime (approximately 500 m depth) but also the shallowest signal penetration depth. The VIBSIST source generates signals with the highest signal‐to‐noise ratio and greatest signal penetration depth of the tested sources. In particular, reflections below 900 ms (approximately 1 km depth) are best imaged by the VIBSIST source. The weight drop performance lies in between these two sources and might be recommended as an appropriate source for a 3D survey at this site because of the shorter production time compared to the VIBSIST and MiniVib sources.     

Pilot‐scale field validation of the long electrode electrical resistivity tomography method

A validation experiment, carried out in a scaled field setting, was attempted for the long electrode electrical resistivity tomography method in order to demonstrate the performance of the technique in imaging a simple buried target. The experiment was an approximately 1/17 scale mock‐up of a region encompassing a buried nuclear waste tank on the Hanford site. The target of focus was constructed by manually forming a simulated plume within the vadose zone using a tank waste simulant. The long electrode results were compared to results from conventional point electrodes on the surface and buried within the survey domain. Using a pole‐pole array, both point and long electrode imaging techniques identified the lateral extents of the pre‐formed plume with reasonable fidelity but the long electrode method was handicapped in reconstructing vertical boundaries. The pole‐dipole and dipole‐dipole arrays were also tested with the long electrode method and were shown to have the least favourable target properties, including the position of the reconstructed plume relative to the known plume and the intensity of false positive targets. The poor performance of the pole‐dipole and dipole‐dipole arrays was attributed to an inexhaustive and non‐optimal coverage of data at key electrodes, as well as an increased noise for electrode combinations with high geometric factors. However, when comparing the model resolution matrix among the different acquisition strategies, the pole‐dipole and dipole‐dipole arrays using long electrodes were shown to have significantly higher average and maximum values within the matrix than any pole‐pole array. The model resolution describes how well the inversion model resolves the subsurface. Given the model resolution performance of the pole‐dipole and dipole‐dipole arrays, it may be worth investing in tools to understand the optimum subset of randomly distributed electrode pairs to produce maximum performance from the inversion model.     

Feasibility of utilizing wavelet phase to map the CO2 plume at the Ketzin pilot site,Germany

Spectral decomposition is a powerful tool that can provide geological details dependent upon discrete frequencies. Complex spectral decomposition using inversion strategies differs from conventional spectral decomposition methods in that it produces not only frequency information but also wavelet phase information. This method was applied to a time‐lapse three‐dimensional seismic dataset in order to test the feasibility of using wavelet phase changes to detect and map injected carbon dioxide within the reservoir at the Ketzin carbon dioxide storage site, Germany. Simplified zero‐offset forward modelling was used to help verify the effectiveness of this technique and to better understand the wavelet phase response from the highly heterogeneous storage reservoir and carbon dioxide plume. Ambient noise and signal‐to‐noise ratios were calculated from the raw data to determine the extracted wavelet phase. Strong noise caused by rainfall and the assumed spatial distribution of sandstone channels in the reservoir could be correlated with phase anomalies. Qualitative and quantitative results indicate that the wavelet phase extracted by the complex spectral decomposition technique has great potential as a practical and feasible tool for carbon dioxide detection at the Ketzin pilot site.     

Joint interpretation of hydrological and geophysical data: electrical resistivity tomography results from a process hydrological research site in the Black Forest Mountains,Germany

The use of electrical resistivity tomography (ERT; non‐intrusive geophysical technique) was assessed to identify the hydrogeological conditions at a surface water/groundwater test site in the southern Black Forest, Germany. A total of 111 ERT transects were measured, which adopted electrode spacings from 0·5 to 5 m as well as using either Wenner or dipole‐dipole electrode arrays. The resulting two‐dimensional (2D) electrical resistivity distributions are related to the structure and water content of the subsurface. The images were interpreted with respect to previous classical hillslope hydrological investigations within the same research basin using both tracer methods and groundwater level observations. A raster‐grid survey provided a quasi 3D resistivity pattern of the floodplain. Strong structural heterogeneity of the subsurface could be demonstrated, and (non)connectivities between surface and subsurface bodies were mapped. Through the spatial identification of likely flow pathways and source areas of runoff, the deep groundwater within the steeper valley slope seems to be much more connected to runoff generation processes within the valley floodplain than commonly credited in such environmental circumstances. Further, there appears to be no direct link between subsurface water‐bodies adjacent to the stream channel. Deep groundwater sources are also able to contribute towards streamflow from exfiltration at the edge of the floodplain as well as through the saturated areas overlying the floodplain itself. Such exfiltrated water then moves towards the stream as channelized surface flow. These findings support previous tracer investigations which showed that groundwater largely dominates the storm hydrograph of the stream, but the source areas of this component were unclear without geophysical measurements. The work highlighted the importance of using information from previous, complementary hydrochemical and hydrometric research campaigns to better interpret the ERT measurements. On the other hand, the ERT can provide a better spatial understanding of existing hydrochemical and hydrometric data. Copyright © 2009 John Wiley & Sons, Ltd.     

Monitoring and volumetric estimation of injected CO2 using 4D seismic,petrophysical data,core measurements and well logging: a case study at Ketzin,Germany

More than 50 000 tons of CO₂ have been injected at Ketzin into the Stuttgart Formation, a saline aquifer, at approximately 620 m depth, as of summer 2011. We present here results from the 1^st repeat 3D seismic survey that was performed at the site in autumn 2009, after about 22 000 tons of CO₂ had been injected. We show here that rather complex time‐lapse signatures of this CO₂ can be clearly observed within a radius of about 300 m from the injection well. The highly irregular amplitude response within this radius is attributed to the heterogeneity of the injection reservoir. Time delays to a reflection below the injection level are also observed. Petrophysical measurements on core samples and geophysical logging of CO₂ saturation levels allow an estimate of the total amount of CO₂ visible in the seismic data to be made. These estimates are somewhat lower than the actual amount of CO₂ injected at the time of the survey and they are dependent upon the choice of a number of parameters. In spite of some uncertainty, the close agreement between the amount injected and the amount observed is encouraging for quantitative monitoring of a CO₂ storage site using seismic methods.     

Time-lapse processing of 2D seismic profiles with testing of static correction methods at the CO2 injection site Ketzin (Germany)

The Ketzin project provides an experimental pilot test site for the geological storage of CO₂. Seismic monitoring of the Ketzin site comprises 2D and 3D time-lapse experiments with baseline experiments in 2005. The first repeat 2D survey was acquired in 2009 after 22 kt of CO₂ had been injected into the Stuttgart Formation at approximately 630 m depth. Main objectives of the 2D seismic surveys were the imaging of geological structures, detection of injected CO₂, and comparison with the 3D surveys. Time-lapse processing highlighted the importance of detailed static corrections to account for travel time delays, which are attributed to different near-surface velocities during the survey periods. Compensation for these delays has been performed using both pre-stack static corrections and post-stack static corrections. The pre-stack method decomposes the travel time delays of baseline and repeat datasets in a surface consistent manner, while the latter cross-aligns baseline and repeat stacked sections along a reference horizon.Application of the static corrections improves the S/N ratio of the time-lapse sections significantly. Based on our results, it is recommended to apply a combination of both corrections when time-lapse processing faces considerable near-surface velocity changes. Processing of the datasets demonstrates that the decomposed solution of the pre-stack static corrections can be used for interpretation of changes in near-surface velocities. In particular, the long-wavelength part of the solution indicates an increase in soil moisture or a shallower groundwater table in the repeat survey.Comparison with the processing results of 2D and 3D surveys shows that both image the subsurface, but with local variations which are mainly associated to differences in the acquisition geometry and source types used. Interpretation of baseline and repeat stacks shows that no CO₂ related time-lapse signature is observable where the 2D lines allow monitoring of the reservoir. This finding is consistent with the time-lapse results of the 3D surveys, which show an increase in reflection amplitude centered around the injection well. To further investigate any potential CO₂ signature, an amplitude versus offset (AVO) analysis was performed. The time-lapse analysis of the AVO does not indicate the presence of CO₂, as expected, but shows signs of a pressure response in the repeat data.     

Mineshaft imaging using surface and crosshole 3D electrical resistivity tomography: A case history from the East Pennine Coalfield, UK

Hidden mineshafts located in urban areas are a significant problem across much of the industrialized world. Electrical resistivity tomography (ERT) is a technique that can detect and characterize hidden mine entries by exploiting resistivity contrasts between the shaft and surrounding materials, resulting from either compositional or structural differences. A case study is presented in which both surface and crosshole 3D ERT surveys are used to image a hidden backfilled mineshaft at a built environment site, situated on Carboniferous Lower Coal Measures strata in the UK.Backfilled shafts generally present the greatest challenge for detection using geophysical methods, as contrasts between the fill and bedrock are typically low compared to air or water-filled conditions. Nevertheless, the shaft in this case was identified by both the surface and crosshole 3D surveys. The shaft appeared as a strongly resistive anomaly relative to background materials, which we interpreted as resulting from the disturbed fabric of the fill materials rather than any significant compositional differences. The study highlighted the respective strengths and weaknesses of the surface and crosshole ERT methodologies for this type of problem. The surface survey, which covered a non-rectangular area to accommodate irregular boundaries and other physical obstructions, provided a relatively rapid means of investigating the study site. However, this method had a limited depth of investigation and was constrained in its coverage by the locations of buildings. By contrast, the 3D crosshole method was able to image the shaft to the level of the deepest borehole electrodes. Although crosshole ERT is too expensive to be used for large-scale mineshaft surveys, this study clearly demonstrates its suitability for targeted investigations where surface methods cannot be deployed, such as scanning beneath surface structures or in situations where it is essential for resolution to be maintained with depth.     

Toxicity of Interactions of Zinc —Nickel,Copper —Nickel and Zinc —Nickel —Copper to a Freshwater Teleost,Lebistes reticulatus (PETERS)1

The lethal toxicity of mixtures of Zn²⁺ —Ni²⁺, Cu²⁺ —Ni²⁺ and Zn²⁺ —Cu²⁺ —Ni²⁺ to common guppy at 21£C in hard water (total hardness = 260 mg/l as CaCO₃) was studied under static bioassays test conditions with renewal of the test solutions every 24 h. The heavy metals were tested separately and in mixtures. The 48 h median lethal concentrations (LC₅₀) for individual salts were 75 mg/l Zn²⁺, 37 mg/l for Ni²⁺ and 2.5 mg/l for Cu²⁺. Concentrations were expressed in “toxic units” by taking them as proportions of LC₅₀ values. Experiments showed that in the Zn²⁺-Ni²⁺ mixture, when Ni²⁺ was more in proportion, the toxicity was more than additive. The 48 h LC₅₀ value and 95% confidence limits in the Ni²⁺-Cu²⁺ mixture were 0.684 (0.484 … 0.807) toxic units and the mixture produced more than the additive toxicity (synergism.). The LC₅₀ value and its 95% confidence limits in a Zn²⁺?Cu²⁺?Ni²⁺ mixture also suggested that the mixture was again strictly additive. The results indicate that heavy metallic mixtures would pose a greater toxicological danger to fish than the respective individual metals.     

Tracing solute infiltration using a combined method of dye tracer test and electrical resistivity tomography in an undisturbed forest soil profile

An accurate prediction of solute infiltration in a soil profile is important in the area of environmental science, groundwater and civil engineering. We examined the infiltration pattern and monitored the infiltration process using a combined method of dye tracer test and electrical resistivity tomography (ERT) in an undisturbed field soil (1 m × 1 m). A homogeneous matrix flow was observed in the surface soil (A horizon), but a preferential flow along macropores and residual rock structure was the dominant infiltration pattern in the subsurface soil. Saturated interflow along the slopping boundaries of A and C1 horizons and of an upper sandy layer and a lower thin clay layer in the C horizon was also observed. The result of ERT showed that matrix flow started first in A horizon and then the infiltration was followed by the preferential flows along the sloping interfaces and macropores. The ERT did not show as much detail as the dye‐stained image for the preferential flow. However, the area with the higher staining density where preferential flow was dominant showed a relatively lower electrical resistivity. The result of this study indicates that ERT can be applied for the monitoring of solute transportation in the vadose zone. Copyright © 2010 John Wiley & Sons, Ltd.     

Monitoring CO2 saturation using time‐lapse amplitude versus offset analysis of 3D seismic data from the Ketzin CO2 storage pilot site,Germany

The injection of CO₂ at the Ketzin pilot site commenced in June 2008 and was terminated in August 2013 after 67 kT had been injected into a saline formation at a depth of 630–650 m. As part of the site monitoring program, four 3D surface seismic surveys have been acquired to date, one baseline and three repeats, of which two were conducted during the injection period, and one during the post‐injection phase. The surveys have provided the most comprehensive images of the spreading CO₂ plume within the reservoir layer. Both petrophysical experiments on core samples from the Ketzin reservoir and spectral decomposition of the 3D time‐lapse seismic data show that the reservoir pore pressure change due to CO₂ injection has a rather minor impact on the seismic amplitudes. Therefore, the observed amplitude anomaly is interpreted to be mainly due to CO₂ saturation. In this study, amplitude versus offset analysis has been applied to investigate the amplitude versus offset response from the top of the sandstone reservoir during the injection and post‐injection phases, and utilize it to obtain a more quantitative assessment of the CO₂ gaseous saturation changes. Based on the amplitude versus offset modelling, a prominent decrease in the intercept values imaged at the top of the reservoir around the injection well is indeed associated solely with the CO₂ saturation increase. Any change in the gradient values, which would, in case it was positive, be the only signature induced by the reservoir pressure variations, has not been observed. The amplitude versus offset intercept change is, therefore, entirely ascribed to CO₂ saturation and used for its quantitative assessment. The estimated CO₂ saturation values around the injection area in the range of 40%–60% are similar to those obtained earlier from pulsed neutron‐gamma logging. The highest values of 80% are found in the second seismic repeat in close vicinity to the injection and observation wells.     

Vertical seismic profiling using a daisy‐chained deployment of fibre‐optic cables in four wells simultaneously – Case study at the Ketzin carbon dioxide storage site

The geological storage of carbon dioxide is considered as one of the measures to reduce greenhouse gas emissions and to mitigate global warming. Operators of storage sites are required to demonstrate safe containment and stable behaviour of the storage complex that is achieved by geophysical and geochemical monitoring, combined with reservoir simulations. For site characterization, as well as for imaging the carbon dioxide plume in the reservoir complex and detecting potential leakage, surface and surface‐borehole time‐lapse seismic monitoring surveys are the most widespread and established tools. At the Ketzin pilot site for carbon dioxide storage, permanently installed fibre‐optic cables, initially deployed for distributed temperature sensing, were used as seismic receiver arrays, demonstrating their ability to provide high‐resolution images of the storage formation. A vertical seismic profiling experiment was acquired using 23 source point locations and the daisy‐chained deployment of a fibre‐optic cable in four wells as a receiver array. The data were used to generate a 3D vertical seismic profiling cube, complementing the large‐scale 3D surface seismic measurements by a high resolution image of the reservoir close to the injection well. Stacking long vibro‐sweeps at each source location resulted in vertical seismic profiling shot gathers characterized by a signal‐to‐noise ratio similar to gathers acquired using geophones. A detailed data analysis shows strong dependency of data quality on borehole conditions with significantly better signal‐to‐noise ratio in regions with good coupling conditions.     

Challenges of using electrical resistivity method to locate karst conduits—A field case in the Inner Bluegrass Region, Kentucky

Conduits serve as major pathways for groundwater flow in karst aquifers. Locating them from the surface, however, is one of the most challenging tasks in karst research. Geophysical methods are often deployed to help locate voids by mapping variations of physical properties of the subsurface. Conduits can cause significant contrasts of some physical properties that can be detected; other subsurface features such as water-bearing fractures often yield similar contrasts, which are difficult to distinguish from the effects of the conduits. This study used electrical resistivity method to search for an unmapped karst conduit that recharges Royal Spring in the Inner Bluegrass karst region, Kentucky, USA. Three types of resistivity techniques (surface 2D survey, quasi-3D survey, and time-lapse survey) were used to map and characterize resistivity anomalies. Some of the major anomalies were selected as drilling targets to verify the existence of the conduits. Drilling near an anomaly identified by an electrical resistivity profile resulted in successful penetration of a major water-filled conduit. The drilling results also suggest that, in this study area, low resistivity anomalies in general are associated with water-bearing features. However, differences in the anomaly signals between the water-filled conduit and other water-bearing features such as water-filled fracture zones were undistinguishable. The electrical resistivity method is useful in conduit detection by providing potential drilling targets. Knowledge of geology and hydrogeology about the site and professional judgment also played important roles in locating the major conduit.     

Evaluation of trenching,ground penetrating radar (GPR) and electrical resistivity tomography (ERT) for sinkhole characterization

This paper explores the suitability and advantages of combining the trenching technique with geophysical surveys [ground penetrating radar (GPR) and electrical resistivity tomography (ERT)] for sinkhole characterization in a mantled karst area. The approach is applied to two active sinkholes concealed by anthropogenic deposits and formed by contrasting subsidence mechanisms; collapse and sagging. The ERT section acquired across the collapse sinkhole images the clayey fill of the depression as an obvious low resistivity area, showing the approximate location of the sinkhole edges. Spatially dense GPR surveys provide information on the position of the boundaries of the concealed subsidence structures and their three‐dimensional (3D) internal geometry, revealing the dominant subsidence mechanism. We illustrate the impact of several factors on the quality of the GPR data such as sinkhole size, nominal frequency of the antennas, antenna shielding, and the presence of backfilled excavations and above‐surface objects. Trenches provided detailed information on the subsurface structure of the sinkhole, subsidence magnitude, partitioning of the strain, and the position of the sinkhole edges, especially when they are deep enough and excavated across the central sector and perpendicular to the boundaries. The stratigraphic and structural relationships observed in the trench were then used to infer the spatial evolution of the sinkholes (e.g. enlargement), their kinematic behavior (episodic versus progressive), and to differentiate discrete subsidence events and their associated magnitude. Numerical dates were used to estimate average subsidence rates and the recurrence of subsidence events. Such integrated data sets may be used as an objective basis to forecast the future behavior of potentially damaging sinkholes and to assess the associated hazard and risk. Copyright © 2013 John Wiley & Sons, Ltd.     

Microearthquakes at the puhagan geothermal field, Philippines — A case of induced seismicity

The Puhagan area in Southern Negros is the only known Philippine geothermal field where there is a clear correlation between increased levels of local seismicity, and the development and early production phases of a geothermal power project. During commissioning of the Palinpinon I power plant in May 1983, a large increase in the microseismic event rate, occasionally exceeding 100 events per day, was noted. This seismic activity is characterized by swarms of events lasting from several hours to a month, separated by long periods of reduced activity. The largest events have local magnitudes of 2.4. Because the swarms appear to be triggered by both reinjection and production of fluids, it is difficult to relate them to a single triggering mechanism.An epicenter study was conducted during July to October 1983, using a simplified joint determination algorithm modified for a uniform velocity structure. The vast majority of the hypocenters occur in a narrow zone with a WNW lineation in the production sector of the field (correlating with a known fault trace) with very little activity in the reinjection sector (1 km to the north). First motions suggest activity has been induced on several non-parallel faults in the area, however, a majority of the events are consistent with normal faulting or oblique slip on steeply dipping NW-SE-trending planes.     

The electrical conductivity of Posidonia black shales — from magnetotelluric exploration to rock samples

We carried out a magnetotelluric field campaign in the South–East Lower Saxony Basin, Germany, with the main goal of testing this method for imaging regional Posidonia black shale sediments. Two‐dimensional inversion results of the magnetotelluric data show a series of conductive structures correlating with brine‐saturated sediments but also with deeper, anthracitic Westphalian/Namurian coals. None of these structures can be directly related with the Posidonia black shale, which appears to be generally resistive and therefore difficult to resolve with the magnetotelluric method. This assumption is supported by measurements of electrical resistivity on a set of Posidonia shale samples from the Hils syncline in the Lower Saxony basin. These rock samples were collected in shallow boreholes and show immature (0.53% Ro), oil (0.88% Ro), and gas (1.45% Ro) window thermal maturities. None of the black shale samples showed low electrical resistivity, particularly those with oil window maturity show resistivity exceeding 10⁴ Ωm. Moreover, we could not observe a direct correlation between maturity and electrical resistivity; the Harderode samples showed the highest resistivity, whereas the Haddessen samples showed the lowest. A similar trend has been seen for coals in different states of thermal maturation. Saturation of the samples with distilled and saline water solutions led to decreasing electrical resistivity. Moreover, a positive correlation of electrical resistivity with porosity is observed for the Wickensen and Harderode samples, which suggests that the electrical resistivity of the Posidonia black shale is mainly controlled by porosity.