甘肃井网的流量观测及震例分析

对甘肃省3口井流量的多年旬均值数据作了线性趋势分析处理,用1.5倍剩余标准差作为异常分辨控制线,研究了1985年6月24日宕昌5.0级地震、1987年1月8日迭部5.9级地震和1987年10月25日礼县5.1级地震前清水井、礼县井和成县井的异常形态和异常过程,给出了甘肃井网中流量井反映地震的能力明显高于水位井这个基本事实。     

Efficacy of an In‐Well Sonde to Determine Magnetic Susceptibility of Aquifer Sediment

Magnetite is a natural component of many aquifers. Abiotic degradation of chlorinated solvents by magnetite can be an important mechanism for natural attenuation of these contaminants. The quantity of magnetite in aquifer materials can be estimated by measuring the magnetic susceptibility of the materials. This is most commonly done by determining the magnetic susceptibility of core samples in an analytical laboratory using a magnetic susceptibility meter. Unfortunately, the cost of acquiring core samples often makes an evaluation of abiotic degradation by magnetite economically unrealistic. Downhole sondes (probes) are available for the determination of magnetic susceptibility. In this study, a downhole sonde was evaluated as an affordable alternative to acquiring and analyzing core samples. The sonde was introduced into 10 monitoring wells. The data from the sonde were then compared to data from core samples that were collected from the same elevation as the sonde data. The core samples analyzed in the laboratory were used as the standard against which the sonde data were compared. The downhole sonde reported values that were similar to values reported on core samples. At most wells, the means of the two measurements could not be distinguished at the 95% confidence interval. When the means could be distinguished, they still agreed within a factor of two.     

A brief review of kriging and its application to optimal interpolation and observation well selection

ABSTRACT

The current state of kriging in subsurface hydrology is critically reviewed. In an application to a region where boreholes already exist, methods of optimal location of additional observation wells for geophysical parameter investigation and optimal interpolation for the purpose of solving the inverse problem are investigated. The particular case of the location of wells for the measurements of transmissivity and hydraulic head in the Kennet Valley Chalk aquifer, UK, is examined. Results of interpolation of measured hydraulic conductivity values by kriging are compared with results from a standard graphical package for interpolation. Reference is also made to the distribution obtained by the inverse method (in which the conductivity distribution is obtained from the head distribution). On the basis of the application, the conditional simulation (in which the generated data are both consistent with field values and the field statistical structure) is deemed to be the best. It is also found that different methods of interpolation give widely different distributions in the case of hydraulic conductivity. It is suggested that the kriged map or conditional map of the transmissivity should serve as the basis for regional discretization to which corrections via the inverse model may be made.     

Ellipsoidal area mean gravity anomalies — precise computation of gravity anomaly reference fields for remove-compute-restore geoid determination

Gravity anomaly reference fields, required e.g. in remove-compute-restore (RCR) geoid computation, are obtained from global geopotential models (GGM) through harmonic synthesis. Usually, the gravity anomalies are computed as point values or area mean values in spherical approximation, or point values in ellipsoidal approximation. The present study proposes a method for computation of area mean gravity anomalies in ellipsoidal approximation (‘ellipsoidal area means’) by applying a simple ellipsoidal correction to area means in spherical approximation. Ellipsoidal area means offer better consistency with GGM quasigeoid heights. The method is numerically validated with ellipsoidal area mean gravity derived from very fine grids of gravity point values in ellipsoidal approximation. Signal strengths of (i) the ellipsoidal effect (i.e., difference ellipsoidal vs. spherical approximation), (ii) the area mean effect (i.e., difference area mean vs. point gravity) and (iii) the ellipsoidal area mean effect (i.e., differences between ellipsoidal area means and point gravity in spherical approximation) are investigated in test areas in New Zealand and the Himalaya mountains. The impact of both the area mean and the ellipsoidal effect on quasigeoid heights is in the order of several centimetres. The proposed new gravity data type not only allows more accurate RCR-based geoid computation, but may also be of some value for the GGM validation using terrestrial gravity anomalies that are available as area mean values.     

Remote measurement of river morphology via fusion of LiDAR topography and spectrally based bathymetry

This study developed and evaluated a hybrid approach to remote measurement of river morphology that combines LiDAR topography with spectrally based bathymetry. Comparison of filtered LiDAR point clouds with surveyed cross‐sections indicated that subtle features on low‐relief floodplains were accurately resolved by LiDAR but that submerged areas could not be detected due to strong absorption of near‐infrared laser pulses by water. The reduced number of returns made the active channel evident in a LiDAR point density map. A second dataset suggested that pulse intensity also could be used to discriminate land from water via a threshold‐based masking procedure. Fusion of LiDAR and optical data required accurate co‐registration of images to the LiDAR, and we developed an object‐oriented procedure for achieving this alignment. Information on flow depths was derived by correlating pixel values with field measurements of depth. Highly turbid conditions dictated a positive relation between green band radiance and flow depth and contributed to under‐prediction of pool depths. Water surface elevations extracted from the LiDAR along the water's edge were used to produce a continuous water surface that preserved along‐channel variations in slope. Subtracting local flow depths from this surface yielded estimates of the bed elevation that were then combined with LiDAR topography for exposed areas to create a composite representation of the riverine terrain. The accuracy of this terrain model was assessed via comparison with detailed field surveys. A map of elevation residuals showed that the greatest errors were associated with underestimation of pool depths and failure to capture cross‐stream differences in water surface elevation. Nevertheless, fusion of LiDAR and passive optical image data provided an efficient means of characterizing river morphology that would not have been possible if either dataset had been used in isolation. Copyright © 2011 John Wiley & Sons, Ltd.     

Test of statistical means for the extrapolation of soil depth point information using overlays of spatial environmental data and bootstrapping techniques

Hydrological modelling depends highly on the accuracy and uncertainty of model input parameters such as soil properties. Since most of these data are field surveyed, geostatistical techniques such as kriging, classification and regression trees or more sophisticated soil‐landscape models need to be applied to interpolate point information to the area. Most of the existing interpolation techniques require a random or regular distribution of points within the study area but are not adequate to satisfactorily interpolate soil catena or transect data. The soil landscape model presented in this study is predicting soil information from transect or catena point data using a statistical mean (arithmetic, geometric and harmonic mean) to calculate the soil information based on class means of merged spatial explanatory variables. A data set of 226 soil depth measurements covering a range of 0–6·5 m was used to test the model. The point data were sampled along four transects in the Stubbetorp catchment, SE‐Sweden. We overlaid a geomorphology map (8 classes) with digital elevation model‐derived topographic index maps (2–9 classes) to estimate the range of error the model produces with changing sample size and input maps. The accuracy of the soil depth predictions was estimated with the root mean square error (RMSE) based on a testing and training data set. RMSE ranged generally between 0·73 and 0·83 m ± 0·013 m depending on the amount of classes the merged layers had, but were smallest for a map combination with a low number of classes predicted with the harmonic mean (RMSE = 0·46 m). The results show that the prediction accuracy of this method depends on the number of point values in the sample, the value range of the measured attribute and the initial correlations between point values and explanatory variables, but suggests that the model approach is in general scale invariant. Copyright © 2009 John Wiley & Sons, Ltd.     

Spatial variability of in situ microbial activity: biotracer tests

Biotracer tests have been proposed as a means by which to characterize the in situ biodegradation potential for field-scale systems. In this study, field experiments were conducted at two sites to evaluate the utility of the biotracer method for characterizing the spatial variability of microbial activity. The first site is a mixed waste-contaminated surficial aquifer in Utah, and the second site is a chlorinated solvent-contaminated regional aquifer in Tucson, Arizona. Mass recovery of the biotracer decreased approximately linearly with increasing residence time for the Tucson site. Similar behavior was observed at the Utah site, except in the region adjacent to the injection zone, where percent recoveries were much lower than those predicted using a correlation determined using data collected downgradient of the injection zone. First-order biodegradation rate coefficients obtained from model calibration of the tracer data varied between 0.2 and 0.5/day for the Tucson site. For the Utah site, the values varied between 0.1 and 0.6/day downgradient of the injection wells, and between 0.7 and 2.6/day near the injection wells. Considering the large range over which biodegradation rate coefficients can vary, the rate coefficient exhibited relatively minimal spatial variability (factor of 2.5) for the Tucson site. Conversely, the spatial variability of the rate coefficient was an order of magnitude greater for the Utah site. These differences in variability are consistent with conditions associated with the respective sites. For example, the greater microbial activity observed in the vicinity of the injection wells for the Utah site is consistent with the biomass distribution determined from analysis of core samples, which shows larger bacterial cell densities for the region near the injection wells. These results illustrate the utility of biotracer tests for in situ characterization of microbial activity (e.g., biodegradation potential), including evaluation of potential spatial variability.     

Heat, chloride, and specific conductance as ground water tracers near streams

Commonly measured water quality parameters were compared to heat as tracers of stream water exchange with ground water. Temperature, specific conductance, and chloride were sampled at various frequencies in the stream and adjacent wells over a 2-year period. Strong seasonal variations in stream water were observed for temperature and specific conductance. In observation wells where the temperature response correlated to stream water, chloride and specific conductance values were similar to stream water values as well, indicating significant stream water exchange with ground water. At sites where ground water temperature fluctuations were negligible, chloride and/or specific conductance values did not correlate to stream water values, indicating that ground water was not significantly influenced by exchange with stream water. Best-fit simulation modeling was performed at two sites to derive temperature-based estimates of hydraulic conductivities of the alluvial sediments between the stream and wells. These estimates were used in solute transport simulations for a comparison of measured and simulated values for chloride and specific conductance. Simulation results showed that hydraulic conductivities vary seasonally and annually. This variability was a result of seasonal changes in temperature-dependent hydraulic conductivity and scouring or clogging of the streambed. Specific conductance fits were good, while chloride data were difficult to fit due to the infrequent (quarterly) stream water chloride measurements during the study period. Combined analyses of temperature, chloride, and specific conductance led to improved quantification of the spatial and temporal variability of stream water exchange with shallow ground water in an alluvial system.     

SUBSURFACE GRAVITY MEASUREMENTS AND DENSITY MAPPING IN THE MAILARAM COPPER MINES,ANDHRA PRADESH,INDIA*

Gravity measurements were made in the Mailaram copper mines, Andhra Pradesh. The observations were distributed between the two shafts situated about 220 m apart and in the three levels up to a maximum depth of 100 m. Assuming a normal free-air gradient, average densities for the three layers were determined as 2.631, 2.604, and 2.823 g cm^-3, respectively. Upon incorporating the weighted mean density values from measurements on samples, the free-air gradients were found to be 0.315 mGal m^-1 for the second layer (i.e. between the first and second levels) and 0.2978 mGal m^-1 for the third layer (i.e. between the second and third levels). The density variation map obtained from the gravity data, the deduced anomalies, the weighted mean density values from measurements on rock samples, and the varying free-air coefficients all suggest correspondence with the concentration of ore lodes.     

COMPARISON OF CORE ANALYSIS AND DRAWDOWN-TEST RESULTS FROM A WATER-BEARING UPPER PENNSYLVANIAN SANDSTONE OF CENTRAL OKLAHOMA / Comparaison de l'analyse d'échantillons et des résultats de rabattements pour une couche perméable de grès du Pensylvanien Supérieur dans l'Oklahoma Central

Abstract

Drawdown-test data from a brine-saturated consolidated sandstone reservoir in central Oklahoma show that permeability can be determined accurately from laboratory measurement of small core samples. Such confirmation is significant because permeability values determined by core analyses are frequently considered not representative, at least quantitatively, of reservoir conditions. Yet such analyses are often the only way to obtain the permeability value necessary to predict the performance of a reservoir.

The field value for permeability was computed by applying the Theis typecurve solution of the non-equilibrium formula to drawdown data from 64 wells in an 8-acre area surrounding a single discharging well. The laboratory value was based on complete cores from 39 of these wells. For objective analysis, the gross reservoir thickness was defined as that entire interval in which flow was confiined, regardless of the permeability of the included rock strata. This thickness was measured directly from the cores. Each plug was saturated with reservoir water and permeability was measured with a constant-head discharging apparatus.

Both laboratory and drawdown-test methods show the reservoir to have a permeability of about 0.4 darcy and a mean transmissibility (permeability x thickness) of about 3.4 darcy-ft.     

Estimation of areal specific yield in sands using the central limit theorem

A new approach is demonstrated that permits a reliable estimate of specific yield using published values of the van Genuchten water retention parameters and effective grain sizes and the measured effective grain sizes of soil samples. The specific yield distribution of the soil texture was computed using the published values of the van Genuchten parameters. The specific yield values and the published values of effective grain sizes were then used to construct a specific yield–effective grain size curve, which estimates the ‘point’ specific yield of the soil samples. Applying the central limit theorem, the point specific yields could be transformed into an ‘areal’ specific yield for a study area. Compared with other commonly used approaches, the present procedure requires relatively low computational efforts and readily obtainable data. It is cost effective and does not depend on soil texture classification. More importantly, it incorporates the depth to water table and the variations in grain sizes inherent in natural soil conditions in the estimation. The approach developed was applied for estimating the specific yield of an unconfined sandy aquifer created by land reclamation in the equatorial region. The values obtained were compared with field measurements and the typical ranges of specific yield from the literature. Instead of a single estimate of the specific yield, the method yields a confidence interval with a high confidence level of 95% and with a narrower range than the typical ranges from the literature. In addition, the estimated values are close to the field measurements; hence, the procedure provides a cost‐effective alternative to field measurement. The applicability of the present approach could be extended to sites with heterogeneity in the horizontal direction. Nevertheless, the applicability of the present approach for layered soil profiles requires further evaluations. Copyright © 2006 John Wiley & Sons, Ltd.     

Comparison of Methods for Sampling Dissolved Nitrogen in a Fractured Carbonate-Rock Aquifer

As part of an agricultural non-point-source study in the Conestoga River head waters area in Pennsylvania, different methods for collecting ground water samples from a fractured carbonate-rock aquifer were compared. Samples were collected from seven wells that had been cased to bedrock and drilled as open holes to the first significant water-bearing zone. All samples were analyzed for specific conductance, dissolved oxygen, and dissolved-nitrogen species. Water samples collected by a point sampler without pumping the well were compared to samples collected by a submersible pump and by a point sampler after pumping the well. Samples collected by using a point sampler, adjacent to major water-bearing zones in an open borehole without pumping the well, were not statistically different from samples collected from the pump discharge or from point samples collected adjacent to major water-bearing zones after pumping the well. Samples collected by using a point sampler without pumping the well at depths other than those adjacent to the water-bearing zones did not give the same results as the other methods, especially when the water samples were collected from within the well casings. It was concluded that, for the wells at this site, sampling adjacent to major water-bearing zones by using a point sampler without pumping the well provides samples that are as representative of aquifer conditions as samples collected from the pump discharge after reaching constant temperature and specific conductance, and by using a point sampler after pumping the well.     

Microbial and geochemical quality of shallow well water in high-density areas in Mzuzu City in Malawi

In Malawi, shallow wells constitute the most important water sources for domestic purposes. However, increasing human population coupled with poor sanitation and infrastructure is undermining the quality of shallow well water. An assessment of microbial and geochemical quality of shallow well water in high-density areas of Zolozolo, Ching’ambo and Chiputula in Mzuzu City, Northern Malawi, has been carried out. The study aimed at characterising domestic water sources, identifying possible sources of water contamination and determining levels of microbial and chemical contamination. Arc-view GIS was used to map the water sources. A questionnaire survey was carried out to elicit information on characteristics of drinking water sources. Water samples were collected from quasi-randomly selected shallow wells and analysed for microbial and chemical parameters using standard methods. HCA, performed using R-programme, was used to group sampled sites according to their bio-physicochemical characteristics. Compliance of the water with MBS/WHO water quality guidelines was determined. The WQI was computed to turn multifaceted data obtained from laboratory analyses into simple information that is comprehensible and useable by the public to assess overall quality of water at a specific water points. The GW-chart was used to show hydrogeochemical water types from each sampled site. Microbial analysis revealed that water from 96.3% of shallow wells recorded faecal coliforms ranging from 129 to 920 cfu per 100 ml which were significantly higher than the Malawi Standards and WHO thresholds. In general, shallow well water is of low mineralisation (EC range 80–500 μS cm⁻¹), with hydrogeochemical facies dominated by Ca–HCO₃, which evolves to Ca–Cl water type. The shallow well water registered a WQI range of 50.16–66.04%, with a medium WQ rating. This suggested that the water obtained from the shallow wells is unsuitable for direct human consumption. It was observed that 100% of the shallow wells were at risk of pollution from onsite sanitation because of their proximity to sanitary facilities. It was strongly recommended that onsite treatment interventions have to be mobilised and initiated to protect the households from further possible consequences of using the water.     

Sea gravity data in the Gulf of Naples. A contribution to delineating the structural pattern of the Phlegraean Volcanic District

Bathymetric and gravity surveys were carried out from 1988 to 1994, in the Gulf of Naples (Southern Italy) to offshore extend the already existing Bouguer anomaly map. In order to improve the knowledge of the structural setting beneath the active Neapolitan volcanoes (Vesuvio, Campi Flegrei and Ischia), 862 stations were surveyed within the isobath of 400 m; at the same time, and about 2000 on-land gravity values were also collected. A new Bouguer anomaly map spanning the whole volcanic region was drawn from the final data set. Gravity anomalies were referred to the new absolute gravity station in Naples and computed according to 1980 Geodetic Reference System. Finally, a density value of 2200 kg/m³ was used in the computation of the Bouguer and terrain effects. We carried out the inversion of the gravity anomalies adopting a 2.5-D modelling along selected profiles crossing the investigated area. The interpretative models were constrained to data obtained from deep wells and other geophysical investigations.     

Interpretation of Magnetic Anomalies and Estimation of Depth of Magnetic Crust in Slovakia

The magnetic map of Slovakia used in the paper was compiled as part of a project titled Atlas of Geophysical maps and profiles in 2001. The residual magnetic data were analyzed to produce Curie point estimates. To remove distortion of magnetic anomalies caused by the Earth’s magnetic field, reduction to pole transformation was applied to the magnetic anomalies using the magnetization angle of the induced magnetization. Anomalies reduced to the pole tend to be better correlated with tectonic structures. We applied a 3-km upward continuation to the residually compiled magnetic anomalies in order to remove effects of topography. The depth of magnetic dipoles was calculated by an azimuthally averaged power spectrum method for the entire area. Such estimates can be indicative of temperatures in the crust, since magnetic minerals lose their spontaneous magnetization according to Curie temperature of the dominant magnetic minerals in the rocks. The computed Curie point depths in the Slovakia region vary between 15.2 km and 20.9 km. Heat flow higher than 100 mWm⁻² occurs at the central volcanics and eastern part of Slovakia, where the Curie point depths values are shallow. The correlation between Curie point depths, heat flow and crust depth was investigated for two E-W cross sections. Heat flow and Curie point depth values are correlated with each other however, these values could not be correlated with crust depth. The Curie point isotherm, which separates magnetic and non-magnetic parts of the crust, is represented in two cross sections.     

Observing lingering hyporheic storage using electrical resistivity: variations around stream restoration structures,Crabby Creek,PA

Time‐lapse geophysical surveys can map lingering hyporheic storage by detecting changes in response to saline tracer. Tracer tests were conducted in Crabby Creek, an urban stream outside Philadelphia, to examine the influence of stream restoration structures and variable sediment thickness. We compared electrical resistivity surveys with extensive well sampling (57 wells) in two 13.5‐m‐long reaches, each with a step drop created by a J‐hook. The two step drops varied in tracer behaviour, based on both the well data and the geophysical data. The well data showed more variation in arrival time where the streambed sediment was thick and was more uniform where sediment was thin. The resistivity in the reach with thin sediment showed lingering tracer in the hyporheic zone both upstream and downstream from the J‐hook. In the second reach where the sediment was thicker, the lingering tracer in the hyporheic zone was more extensive downstream from the J‐hook. The contrasting results between the two reaches from both methods suggested that sediments influenced hyporheic exchange more than the step at this location. Resistivity inversion differed from well data in both reaches in that it showed evidence for tracer after well samples had returned to background, mapping lingering tracer either upstream or downstream of a step. We conclude that resistivity surveys may become an important tool for hyporheic zone characterization because they provide information on the extent of slow moving fluids in the hyporheic zone, which have the potential to enhance chemical reactions. Copyright © 2012 John Wiley & Sons, Ltd.     

基于表层岩土渗透对三峡井网井水位降雨影响的精准分析

为研究三峡井网表层岩土渗透对井水位降雨的影响,采取井区表层岩土垂向渗透性测试方法试验,测得表层岩土垂向渗透性,并建立数学模型,用于降雨渗入补给分析。在此模型基础上,通过三峡井网8口井水位、气象三要素的对比观测资料对井水位日动态、月动态、年动态的影响进行精准分析与验证。结果表明:这种影响的特征是相当复杂的,同一个降雨过程在不同井上产生的影响特征不同,这一方面可能与各井的水文地质条件不同有关,另一方面可能还与各井点的降雨过程的差异也有关。     

Geoelectrical exploration for groundwater in Al Maha Forest,Ain Jouhra,Morocco

An area of about 30 km² located in Ain Jouhra, south of Rabat, Morocco, was the subject of a geoelectric resistivity investigation. The main goal of the investigation was the assessment of the groundwater potential of the uppermost aquifer. The aquifer conditions such as depth, thickness and boundaries were also investigated. The obtained apparent resistivity curves were first analysed qualitatively and classified using simple curve shapes. Thereafter, the data were converted to resistivity and thickness pairs semi‐quantitatively by means of master curves and then quantitatively by computer modelling using ATO and Winsev software (Zohdy, 1989; Zohdy and Bisdrof, 1989). Lithological control from the available single well with a stratigraphic log aided in the correlation of the resistivity values to different rock units. Three different AB‐spacing iso‐resistivity maps, an isopach map of the main groundwater‐bearing horizon, the depth to the aquifer substratum map and five geoelectric cross‐sections were constructed. The interpretation of these soundings indicates the presence of an unconfined to semi‐confined sandy aquifer with relatively important extent and varying thickness. The maximal thickness of the aquifer is recorded in the central part of the investigated area and is thinning southwards to pinch out farther to the south. Geophysical as well as field data indicate a hydraulic connection between the upper and deeper aquifers. Indeed, the two aquifers are separated from each other by a marly substratum that is indicated throughout the area by the lowest values of the interpreted true resistivity. The value of this resistivity varies laterally, most likely due to the lateral variation in the shale‐to‐sand ratio. The altitude of the substratum decreases towards the north, and increases southwards. Regarding the availability of the groundwater in the study area, zones with high potential are theoretically expected to occur in the central part where the transversal resistance is greatest. However, sufficient water supply and high flow rates from wells intended to produce restrictively from the most upper aquifer are not likely to exist. This conclusion, which seems to be very pessimistic, is evidenced from two real field and experimental observations. The first is the rapid fall of the level of Gharnoug lake, despite the ongoing feeding by three wells. Hence, the amount of water level drop cannot be accounted for by the evaporation alone. That means that the deeper aquifer is continuously draining the upper aquifer at a high flow rate. Very low rates are recorded in all the wells that penetrated only the upper aquifer, the exception being the well that reached deeper into the lower aquifer. The flow rate in this lower aquifer measured 18 litre s^?1. Copyright © 2007 John Wiley & Sons, Ltd.     

MAPPING MEAN AREAL PRECIPITATION

Abstract

A method is presented for using point mean precipitation data to estimate areal values in regions of high relief. Variation of precipitation with altitude is determined. Local anomalies from this relationship are mapped, and lines of equal anomaly are drawn. By use of the mean relation corrected for the local anomaly, the mean precipitation at any point can be determined and an isohyetal map drawn. A similar approach can be used to determine mean temperature for studies of snowmelt or of potential evapotranspiration.     

Ground‐shaking mapping for a scenario earthquake considering effects of geological conditions: a case study for the 1995 Hyogo‐ken Nanbu,Japan earthquake

In order to examine the applicability of ground‐shaking mapping techniques to a near‐field earthquake, a peak ground velocity map of the 1995 Hyogo‐ken Nanbu, Japan earthquake computed from seismic zoning methods that consider the effects of geological conditions is compared with the actual observed intensity map. When computing the ground‐shaking map, the site amplification at each site is calculated in terms of the average shear‐wave velocity of the ground estimated from the corresponding geomorphological conditions. This map shows a relatively good agreement with the observed intensity map. However, the computations provide smaller values for certain disastrous areas of the earthquake, where the effects on ground motion of a deep, irregular underground structure have been reported. The effect of such structures on site response is examined implementing 2D FEM analyses, thereby being also incorporated into the method. Results considering the effect of the irregular underground structure show better agreement with the observed intensity map. Copyright © 2002 John Wiley & Sons, Ltd.