Imaging of groundwater resources in glacial deposits using high-resolution reflection seismics, Sweden

Two high-resolution reflection seismic profiles were acquired in the Heby area of eastern Sweden over glacial deposits for the purpose of mapping groundwater resources. The majority of shot points were located in clay resulting in good quality data along most of the profiles. On stacked and migrated sections, the uppermost clay is about 20 m thick and is characterized by its subhorizontal reflectivity. Sand/gravel deposits below it contain more dipping interfaces and have a chaotic reflectivity pattern. Depth to bedrock is interpreted to be 90 and 65 m on the respective profiles and occurs in about a 100-m-wide trough on both profiles. Reflections from the tops of sandy gravel zones generally have higher amplitudes. Clear reflections from a thin silt layer (20 cm thick) at about 10-m depth are observed on one of the profiles. Elastic finite difference modeling and the observation of this reflection in shot gathers show that the reflection is not an artifact of the acquisition nor the processing. The modeling also shows that there is no marked low-velocity waveguide in the near surface, but that an effective low Q zone may be present. Comparison with refraction profiling on the other profile shows that there is better agreement between the reflection seismic results and penetration tests than the refraction results with these tests. Both profiles allow the thickness of the overlying clay layers to be determined, as well as the thickness of the underlying sand/gravel deposits. This is important for estimating the amount of groundwater resources in an area.     

Mapping of quick clay by electrical resistivity tomography under structural constraint

Geotechnical projects usually rely on traditional sounding and drilling investigations. Drilling only provides point information and the geology needs to be interpolated between these points. Near surface geophysical methods can provide information to fill those gaps. Norwegian case studies are presented to illustrate how two-dimensional electrical resistivity tomography (ERT) can be used to accurately map the extent of quick clay deposits. Quick clay may be described as highly sensitive marine clay that changes from a relatively stiff condition to a liquid mass when disturbed. Quick clay slides present a geo-hazard and therefore layers of sensitive clay need to be mapped in detail. They are usually characterized by higher resistivity than non-sensitive clay and ERT is therefore a suitable approach to identify their occurrence. However, our experience shows that ERT cannot resolve this small resistivity contrast near large anomalies such as a bedrock interface. For this reason, a constrained inversion of ERT data was applied to delineate quick clay extent both vertically and laterally. As compared to the conventional unconstrained inversions, the constrained inversion models exhibit sharper resistivity contrasts and their resistivity values agree better with in situ measurements.     

The detection of King Xerxes' Canal by the use of shallow reflection and refraction seismics — Preliminary results

Shallow reflection and refraction seismic studies were carried out in Greece in the eastern leg of the Chalkidiki peninsula, in order to test the validity of reports in history books which describe a legendary canal built by the engineers of King Xerxes during the major Persian invasion of Europe through Greece in the 5th century B.C.   In the narrowest part of the Athos peninsula, where it is 2 km wide, an 85 m profile was topographically defined almost centrally between the two coastlines. The position of this profile was based on palaeogeographical, geomorphological and topographic studies and observations. A sledgehammer was used as the seismic source for the shallow target. Despite the presence of significant urban and coherent noise, a final stacked section was produced by a suitable choice of acquisition and processing parameters.   Both the reflection and refraction seismic studies illustrated clearly the existence of a channel-like structure of trapezoid cross-section, almost midway between the two opposite sides of the peninsula.     

Soil characterization using electrical resistivity tomography and geotechnical investigations

Electrical Resistivity Tomography (ERT) has been used in association with Standard Penetration Test (SPT) and Dynamic Cone Penetration Test (DCPT) for Geotechnical investigations at two sites, proposed for thermal power plants, in Uttar Pradesh (UP), India. SPT and DCPT tests were conducted at 28 points and two ERT profiles, each measuring 355 m long, were recorded using 72 electrodes deployed at 5 m spacing. Electrical characterization of subsurface soil was done using borehole data and grain size analysis of the soil samples collected from boreholes. The concept of electrical resistivity variation with soil strength related to the grain size distribution, cementation, porosity and saturation has been used to correlate the transverse resistance of soil with the number of blow counts (N-values) obtained from SPT and DCPT data. It was thus observed that the transverse resistance of soil column is linearly related with the number of blow counts (N-values) at these sites. The linear relationships are site-specific and the coefficients of linear relation are sensitive to the lithology of subsurface formation, which was verified by borehole data. The study demonstrates the usefulness of the ERT method in geotechnical investigations, which is economic, efficient and less time consuming in comparison to the other geotechnical methods, such as SPT and DCPT, used for the purpose.     

高分辨深反射地震探测采集处理关键技术综述

深反射地震勘探以其可靠性强和精度高的特点,一直作为深地探测的主要技术之一.其中分辨率又作为评价深反射地震勘探质量的关键指标之一,在理论研究及实际生产中都有着重要的意义.文章将深反射地震勘探中影响分辨率的主要因素划分为采集、地质、处理三方面,采集因素主要包括子波激发、空间采样率、接收及观测系统设计;地质因素则主要包括深部波阻抗变化、深部地质体尺度、浅部构造、吸收衰减等,针对以上因素人们提出了静校正、振幅恢复、反褶积、偏移等处理方法来提高深反射地震勘探的分辨率.本文阐述了影响深反射地震勘探分辨率的主要因素,并对提高深反射地震勘探分辨率的主要方法进行了总结与展望.     

Background velocity estimation using non-linear optimization for reflection tomography and migration misfit

We show that it is possible to estimate the background velocity for prestack depth migration in 2D laterally varying media using a non-linear optimization technique called very fast simulated annealing (VFSA). We use cubic splines in the velocity model parametrization and make use of either successive pairs of shot gathers or several constant-offset sections as input data for the inversion. A Kirchhoff summation scheme based on first-arrival traveltimes is used to migrate/model the input data during the velocity analysis. We evaluate and compare two different measures of error. The first is defined in the recorded data or (x,t) domain and is based on a reflection-tomography criterion. The second is defined in the migrated data or (x,z) domain and is based on a migration-misfit criterion. Depth relaxation is used to improve the convergence and quality of the velocity analysis while simultaneously reducing the computational cost. Further, we show that by coarse sampling in the offset domain the method is still robust. Our non-linear optimization approach to migration velocity analysis is evaluated for both synthetic and real seismic data. For the velocity-analysis method based on the reflection-tomography criterion, traveltimes do not have to be picked. Similarly, the migration-misfit criterion does not require that depth images be manually compared. Interpreter intervention is required only to restrict the search space used in the velocity-analysis problem. Extension of the proposed schemes to 3D models is straightforward but practical only for the fastest available computers.     

2D data modelling by electrical resistivity tomography for complex subsurface geology

A new tool for two‐dimensional apparent‐resistivity data modelling and inversion is presented. The study is developed according to the idea that the best way to deal with ill‐posedness of geoelectrical inverse problems lies in constructing algorithms which allow a flexible control of the physical and mathematical elements involved in the resolution. The forward problem is solved through a finite‐difference algorithm, whose main features are a versatile user‐defined discretization of the domain and a new approach to the solution of the inverse Fourier transform. The inversion procedure is based on an iterative smoothness‐constrained least‐squares algorithm. As mentioned, the code is constructed to ensure flexibility in resolution. This is first achieved by starting the inversion from an arbitrarily defined model. In our approach, a Jacobian matrix is calculated at each iteration, using a generalization of Cohn's network sensitivity theorem. Another versatile feature is the issue of introducing a priori information about the solution. Regions of the domain can be constrained to vary between two limits (the lower and upper bounds) by using inequality constraints. A second possibility is to include the starting model in the objective function used to determine an improved estimate of the unknown parameters and to constrain the solution to the above model. Furthermore, the possibility either of defining a discretization of the domain that exactly fits the underground structures or of refining the mesh of the grid certainly leads to more accurate solutions. Control on the mathematical elements in the inversion algorithm is also allowed. The smoothness matrix can be modified in order to penalize roughness in any one direction. An empirical way of assigning the regularization parameter (damping) is defined, but the user can also decide to assign it manually at each iteration. An appropriate tool was constructed with the purpose of handling the inversion results, for example to correct reconstructed models and to check the effects of such changes on the calculated apparent resistivity. Tests on synthetic and real data, in particular in handling indeterminate cases, show that the flexible approach is a good way to build a detailed picture of the prospected area.     

3-D inversion of borehole-to-surface electrical data using a back-propagation neural network

The “fluid-flow tomography”, an advanced technique for geoelectrical survey based on the conventional mise-à-la-masse measurement, has been developed by Exploration Geophysics Laboratory at the Kyushu University. This technique is proposed to monitor fluid-flow behavior during water injection and production in a geothermal field. However data processing of this technique is very costly. In this light, this paper will discuss the solution to cost reduction by applying a neural network in the data processing. A case study in the Takigami geothermal field in Japan will be used to illustrate this. The achieved neural network in this case study is three-layered and feed-forward. The most successful learning algorithm in this network is the Resilient Propagation (RPROP). Consequently, the study advances the pragmatism of the “fluid-flow tomography” technique which can be widely used for geothermal fields. Accuracy of the solution is then verified by using root mean square (RMS) misfit error as an indicator.     

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

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

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.     

Interpretation of shallow refraction seismic data by reflection/refraction tomography

A new algorithm for tomographic inversion of traveltimes of reflected and refracted seismic waves is developed. The inversion gives interface configurations and velocity distributions in layers. The important features of the algorithm are: (a) the inclusion of shot time delays in the list of unknown parameters; (b) the regularization is applied in such a way that the most probable model is characterized by the similarity of neighbouring interfaces. As the problem under consideration is non-linear, several iterations are necessary in order to obtain the final model. In the case of a very inexact initial model, a 'layer-by-layer' inversion strategy is recommended as a first inversion step. The inversion program is supplied with a user interface, thanks to which one can: (a) pick interactively and identify seismic traveltimes; (b) build and edit depth/velocity models; and (c) display calculated traveltime curves and compare them with picked traveltimes as well as with the original seismic sections. The efficiency of the inversion software developed is illustrated by a numerical example and a field example in which shallow seismic data are considered. Application to wide-aperture reflection/refraction profiling (WARRP) data is also possible.     

Enhancing hardrock seismic images: Reprocessing of high resolution seismic reflection data from Vihanti,Finland

Seismic reflection data were acquired across volcanic hosted massive sulfides (VHMS) of Vihanti in order to improve the understanding of the regional geological setting. Commercially processed seismic data from Vihanti are of good quality, but reprocessing can be used to extract additional information about geological structures. Especially, careful velocity analysis influences the quality of seismic images. Differentiating reflections caused by fractures from those caused by lithological contacts is very important for exploration and geological modeling. Reflections from fracture zones known from drilling stack with lower velocity (~ 5100 m/s) compared to typical stacking velocities of the Vihanti area (> 5500 m/s). The reprocessing also indicated that fracture zones are better imaged with low frequencies due to the better overall continuity of the fault zones at scales of hundreds of meters rather than at shorter seismic wavelengths.In full stacks, long offset data can mask structure close to the surface. More detailed seismic images of the shallow subsurface emerged by preferentially stacking short offset data wherever acquisition and processing lines lay close together and were nearly straight. Long offset data remains valuable for imaging deeper structures as well as dipping reflectors. Cross-dip-analysis revealed a bright diffractor located near the base of the Vihanti volcanic basin at 1.5 km depth. The seismic data allow a geological interpretation in which the Vihanti structure has developed through significant thrust faulting and displacement of the lithological contacts. Gentle folds that were formed prior to faulting are visible as undulating reflectivity in seismic sections. The reprocessed seismic section indicates a potential deep extension of the ore-hosting altered volcanic and calc-silicate rocks previously unexplored.     

垃圾填埋场渗滤液污染地下含水层及修复过程的三维动态监测实验

为实现垃圾填埋场渗滤液在地下含水层中动态扩散及污染层修复过程的实时监测,利用自行设计的三维电学观测系统开展了室内相关监测实验.实验表明渗滤液在含水层中的扩散过程会引起不同时期实测电剖面上低阻异常区的动态变化.对比分析这种变化特征可确定渗滤液扩散区的污染程度、扩散速度及扩散方向.在注水修复过程中,扩散区污染物含量的降低会引起对应区域实测视电阻率值的升高和异常区范围的变化.实验结果对于实现垃圾渗滤液污染地下含水层现状调查及动态监测具有重要意义.     

井间高分辨率纵横波层析成像研究井间油藏

采用视速度-偏振井间波场分离法分离得到纵波和横波波场,使用跨孔地震走时层析成像级联方法计算井间纵波和横波速度结构以及泊松比空间分布.综合纵横波速度、泊松比和测井参数,分析两井间岩性和储层特性.数值模型试验表明该方法垂直分辨率达到主频波长的1/4,并能准确分辨倾斜断层和垂直断层.处理相同纵波资料结果表明本文方法分辨率远远高于以往的井间走时层析成像方法,比井间声波层析成像分辨率更高.使用该方法处理垦71区井间资料,得到同台同源纵波和横波速度结构和泊松比分布.泊松比的低/高变化与测井参数指示的砂/泥岩层基本吻合.根据泊松比分布,参考岩石地震物理学测试数据,可以区分井间介质中的泥岩层、砂岩层,以及砂岩层中的饱水和饱油区,确定储层的连通性,圈定井间的剩余油藏.     

Computation of synthetic seismograms in 2-D and 3-D media using the Gaussian beam method

Summary The Gaussian beam method is applied to vertical seismic profiling in 2-D and 3-D models. A simple approach to the computation of Gaussian beam seismograms in the vicinity of structural interfaces is proposed. The effects of (a) the radiation pattern of a point source, (b) noncausal attenuation, (c) transverse inhomogeneities in synthetic seismograms are studied on numerical examples.

---

m¶rt; a n unm ¶rt; mua u nuuau ¶rt; u m ¶rt; ¶rt;. nmm n¶rt;¶rt; am a mmu mm a¶rt;a. a nua ummuu a uam m (a) uu m umua, () nu nu u () nn ¶rt;¶rt;m.

     

High resolution 2-D maps of OI 630.0 nm thermospheric dayglow from equatorial latitudes

The first-ever high resolution 2-D maps of OI 630.0 nm dayglow obtained from equatorial latitudes clearly reveal the movement as a large-scale feature of the equatorial ionization anomaly (EIA). These also show the presence of wave-like features classified as gravity waves presumably originating at the crest of the EIA, similar to the equatorial electrojet acting as a source of these waves. These results are presented and discussed.     

Lysimetric evaluation of SEBAL using high resolution airborne imagery from BEAREX08

In this study, Surface Energy Balance Algorithm for Land (SEBAL) was evaluated for its ability to derive aerodynamic components and surface energy fluxes from very high resolution airborne remote sensing data acquired during the Bushland Evapotranspiration and Agricultural Remote Sensing Experiment 2008 (BEAREX08) in Texas, USA. Issues related to hot and cold pixel selection and the underlying assumptions of difference between air and surface temperature (dT) being linearly related to the surface temperature were also addressed. Estimated instantaneous evapotranspiration (ET) and other components of the surface energy balance were compared with measured data from four large precision weighing lysimeter fields, two each managed under irrigation and dryland conditions. Instantaneous ET was estimated with overall mean bias error and root mean square error (RMSE) of 0.13 and 0.15 mm h⁻¹ (23.8 and 28.2%) respectively, where relatively large RMSE was contributed by dryland field. Sensitivity analysis of the hot and cold pixel selection indicated that up to 20% of the variability in ET estimates could be attributed to differences in the surface energy balance and roughness properties of the anchor pixels. Adoption of an excess resistance to heat transfer parameter model into SEBAL significantly improved the instantaneous ET estimates.     

Effects of meteorological and hydrogeological factors on gross recharge percentage at unconfined sandy aquifers with an equatorial climate

A thorough understanding of rainfall recharge processes and their controlling factors is essential for management of groundwater systems. This study investigates the effects of various meteorological and hydrogeological factors on the gross recharge percentages, the rainfall–recharge relationships and the recharge threshold values for unconfined sandy aquifers under an equatorial climate. Among the meteorological factors investigated, rainfall intensity was found to have the most significant impact on the gross recharge rate. The effects of potential evaporation rate, relative humidity and air temperature on the gross recharge percentage were significant when the vadose zone thickness is larger than 2·5 m. The recharge threshold values were found to depend strongly on the vadose zone thickness. The rainfall–recharge relationships could generally be well defined by a normal–log relationship. The rainfall–recharge relationships derived here are applicable to yield estimates of gross recharge percentages for unconfined sandy aquifers under an equatorial climate, using rainfall intensity and vadose zone thickness as input variables. In this study, a theory was developed and validated to provide physical explanations for the observations, based on the residence time of the percolated rainwater within the vadose zone. Among the soil hydraulic parameters tested, porosity and saturated hydraulic conductivity were found to have the most pronounced effects on the gross recharge percentage. Utilizing the sensitivity results and the theory derived, an approach was developed for extending the application of the derived rainfall–recharge relationships to other sand textures. The approach was found to be capable of producing rough and fast estimations of gross recharge percentage for other sand textures. Copyright © 2007 John Wiley & Sons, Ltd.     

高分辨瞬时振幅和瞬时相位在薄互煤层识别中的应用

东北地区煤田中的煤层一般是一个典型的低速薄层,可采煤厚约1～8 m,这使得常规地震方法勘探煤层受到很大的限制.但由于煤层与围岩组成的是一个有规律的薄互层组合,所以通过研究薄互层的地震波波形、振幅、能量分配关系等特征来进行识别成为地震勘探煤层的一种可能.而高分辨三瞬技术能够提供精细化的波形、振幅等信息,因此本文将高分辨三...