用梯度和感应测井联合确定地层电阻率各向异性

实际地层可能出现2种电阻率的宏观各向异性:微观的统计平均效应和裂缝走向等因素的构造效应.三分量感应测井方法可能是确定地层电阻率各向异性的最好方法,但是,该仪器目前在中国还没有应用.已有的数据一般是梯度和双感应曲线,梯度和双感应曲线单独使用都难以揭示各向异性,但是,将二者联合使用则有可能.在反演模型中考虑了电阻率各向异性的存在,给出了梯度和双感应曲线联合求地层电阻率各向异性的方法.通过直井水平层的二维人工模型表明,用感应方法可确定地层的水平电阻率和厚度.而梯度方法的视电阻率可以近似看作水平电阻率和纵向电阻率的几何平均值,梯度方法的视厚度是各向异性参数和真实厚度的乘积.因而不能用梯度方法单独确定各向异性参数和真实厚度2个参数中的任何1个.然而,联合梯度和双感应的方法则可能确定3个参数:各向异性参数,水平电阻率和地层厚度.人工数据的例子表明了这种联合的可行性.     

用梯度和感应测井联合确定地层电阻率各向异性(英文)

实际地层可能出现 2种电阻率的宏观各向异性 :微观的统计平均效应和裂缝走向等因素的构造效应。三分量感应测井方法可能是确定地层电阻率各向异性的最好方法 ,但是 ,该仪器目前在中国还没有应用。已有的数据一般是梯度和双感应曲线 ,梯度和双感应曲线单独使用都难以揭示各向异性 ,但是 ,将二者联合使用则有可能。在反演模型中考虑了电阻率各向异性的存在 ,给出了梯度和双感应曲线联合求地层电阻率各向异性的方法。通过直井水平层的二维人工模型表明 ,用感应方法可确定地层的水平电阻率和厚度。而梯度方法的视电阻率可以近似看作水平电阻率和纵向电阻率的几何平均值 ,梯度方法的视厚度是各向异性参数和真实厚度的乘积。因而不能用梯度方法单独确定各向异性参数和真实厚度 2个参数中的任何 1个。然而 ,联合梯度和双感应的方法则可能确定 3个参数 :各向异性参数 ,水平电阻率和地层厚度。人工数据的例子表明了这种联合的可行性     

数字滤波技术在电测井中的应用

本文采用数字滤波技术,研究了电阻率测井的理论曲线计算和R_a曲线的环境影响校正,并对影响数字滤波效果的因素进行了分析。结果表明,用数字滤波新技术进行环境校正,提高了校正质量,比传统的图版校正方法有明显优点。文中采用数字滤波方法绘制的R_a理论曲线,证明所推导的数学模型是正确的。     

Influence of Degree of Saturation in the Electric Resistivity–Hydraulic Conductivity Relationship

The relationship between aquifer hydraulic conductivity and aquifer resistivity, either measured on the ground surface by vertical electrical sounding (VES) or from resistivity logs, or measured in core samples have been published for different types of aquifers in different locations. Generally, these relationships are empirical and semi-empirical, and confined in few locations. This relation has a positive correlation in some studies and negative in others. So far, there is no potentially physical law controlling this relation, which is not completely understood. Electric current follows the path of least resistance, as does water. Within and around pores, the model of conduction of electricity is ionic and thus the resistivity of the medium is controlled more by porosity and water conductivity than by the resistivity of the rock matrix. Thus, at the pore level, the electrical path is similar to the hydraulic path and the resistivity should reflect hydraulic conductivity. We tried in this paper to study the effect of degree of groundwater saturation in the relation between hydraulic conductivity and bulk resistivity via a simple numerical analysis of Archie’s second law and a simplified Kozeny-Carmen equation. The study reached three characteristic non-linear relations between hydraulic conductivity and resistivity depending on the degree of saturation. These relations are: (1) An inverse power relation in fully saturated aquifers and when porosity equals water saturation, (2) An inverse polynomial relation in unsaturated aquifers, when water saturation is higher than 50%, higher than porosity, and (3) A direct polynomial relation in poorly saturated aquifers, when water saturation is lower than 50%, lower than porosity. These results are supported by some field scale relationships.     

Correlation of T-log with E-log in coal-bearing formations

The profile of temperature gradient versus depth (T-log) has been found to be very useful for correlation with electrical resistivity log (E-log) in coal-bearing formations. A positive correlation between electrical resistivity and thermal resistivity is observed in coal sections whereas a negative correlation is found in sandstone/shale beds, thus helping in coal prospect evaluation. T logs have been used to correct the location of coal bed which had apparently been misinterpreted by the E-log. Hole to hole correlation of T-log and E-log is found to be excellent and it is observed that thermal resistivity characteristics of given formations remain fairly uniform. A rough estimation of coal grade is possible from the detailed study of the T-logs. Abrupt changes of temperature gradient as also its reversals have been observed in burnt coal sections.     

Simple Electrical Logging Technique for Base Metal Exploration

—It is well known that electrical logs of boreholes can play a signi|fi|cant role in base metal exploration in identifying mineralised zones, especially when there is core loss or the borehole diameter is small or if drilling is by percussion. However, electrical logging is not widely utilised because of the additional burden on finances and time.¶A simple electrical logging technique, based on a pole-pole (hole-to-surface) configuration with one borehole electrode and nearly akin to the single point method, for S.P., resistivity and I.P. parameters, is presented. It is shown that it has the resolution of the single point method and the penetration of a very long normal sonde which is helpful for detection. Besides these features, the main advantage of this technique is that it can be easily carried out using ground I.P. (time domain) equipment.¶The electrical logs obtained in different base metal belts in Rajasthan, India, employing this technique and using ground I.P. equipment in connection with mise-à-la-masse surveys, illustrate the above features. A comparison with logs recorded by means of multi-electrode drill hole I.P. equipment of Scintrex, Canada, substantiates the same.     

On the automatic interpretation of direct current resistivity soundings

A method for the automatic inversion of resistivity soundings is presented. The procedure consists of two main stages. First, application of linear filters which transforms the apparent resistivity curve into the kernel function, and vice versa. In the second stage the first and second derivatives of the kernel function are calculated and used in a second-order modified Newton-Raphson iterative fitting procedure. The model obtained is optimal in the least squares sense. The method has been tried on some field examples and produced resistivity models which show a good agreement with the geological well logs.     

INVERSION OF SEISMOGRAMS AND PSEUDO VELOCITY LOGS*

Pseudo velocity logs can be obtained by seismogram inversion, using true amplitude processing and detailed investigation of move-out velocities. The precision of the results depends on the quality of the seismic data and on the possibility of deconvolving without increasing the noise. An investigation is made of the deformation of pseudo logs due to seismic signal variations and to imperfections of deconvolution. Both marine and land examples are shown, in some cases with adjustment on well logs. When the dips are large, time sections must be migrated and pseudo velocity logs must be computed from migrated sections. Comparison of sonic logs with pseudo velocity logs obtained in the same area is usually good enough to obtain information on lithological parameter variations by adjustment of pseudo velocity logs on sonic logs. Even when no well is available, pseudo velocity logs can give some indications on the nature of sediments between seismic horizons.     

Synthetic focusing and simulation of dual laterolog tool in axisymmetric subsurface models

The laterolog tools for electrical resistivity measurement are largely employed by the petroleum industry especially in boreholes filled with salt water-base mud. This work develops a synthetic focusing for many electrical tools, such as LLS, LLD and LL7. This is done by applying the principle of superposition to each normal arrangement. A finite element method (FEM) is used to simulate the tool response in axisymmetric subsurface models crossed vertically by the well with constant radius. Also, the logs are corrected for the borehole effects that are, in some cases, important. Despite the application of the synthetic focusing has been accomplished in two-dimensional models, it is possible to apply it in more complex simulations: three-dimensional models with deviated boreholes (or dipping beds) as well as anisotropic systems.     

Geomorphology of the Alluvial Sediments and Bedrock in an Intermontane Basin: Application of Variogram Modeling to Electrical Resistivity Soundings

The study describes a methodology used to integrate legacy resistivity data with limited geological data in order to build three-dimensional models of the near subsurface. Variogram analysis and inversion techniques more typically found in the petroleum industry are applied to a set of 1D resistivity data taken from electrical surveys conducted in the 1980s. Through careful integration with limited geological data collected from boreholes and outcrops, the resultant model can be visualized in three dimensions to depict alluvium layers as lithological and structural units within the bedrock. By tuning the variogram parameters to account for directionality, it is possible to visualize the individual lithofacies and geomorphological features in the subsurface. In this study, an electrical resistivity data set collected as part of a groundwater study in an area of the Peshawar basin in Pakistan has been re-examined. Additional lithological logs from boreholes throughout the area have been combined with local outcrop information to calibrate the data. Tectonic activity during the Himalayan orogeny has caused uplift in the area and generated significant faulting in the bedrock resulting in the formation of depressions which are identified by low resistivity values representing clays. Paleo-streams have reworked these clays which have been eroded and replaced by gravel–sand facies along paleo-channels. It is concluded that the sediments have been deposited as prograding fan-shaped bodies and lacustrine deposits with interlayered gravel–sand and clay–silt facies. The Naranji area aquifer system has thus been formed as a result of local tectonic activity with fluvial erosion and deposition and is characterized by coarse sediments with high electrical resistivities.     

A new look at multiphase invasion with applications to borehole resistivity interpretation

In well log interpretation, it is frequently necessary to correct logs for invasion. Invasion occurs in permeable formations when there is a radial differential pressure (RDP) between the borehole and formation. Other factors on which invasion depend include saturation, mobility, pressure (RDP) and capillary pressure, permeability and viscosity of fluids, and temperature transient effects associated with the mud filtrate injected into the formation. Thus, simulation of realistic invasion is not an easy task. This work reviews the famous Buckley–Leverett mathematical model in cylindrical coordinates appropriate for borehole geometries. The model predicts multiphase invasion in porous media when gravity, capillary pressure, and mud cake can be neglected. One application is to correct logging while drilling (LWD) and wireline resistivity logs for time-dependent invasion and formation temperature effects. This is important, for example, when there are possible large differences in formation and mud temperature. Modeling studies show these effects can be large enough to noticeably influence resistivity logs. However, after correction, difference in LWD and wireline logs arising from the time-dependent heat process are explained. Thus, the method, when coupled to a time-dependent heat flow model, and a response function formulation of resistivity, yields new insight into the influence of thermal and electrical transients in log interpretation.     

COMPUTATION OF THE ELECTRIC POTENTIAL IN THREE-DIMENSIONAL STRUCTURES*

A method has been developed for computation of the electrical DC potential of an arbitrary three-dimensional resistivity structure using a finite difference method. The threedimensionality is necessary for interpretation of geoelectrical soundings with controlled point sources over a laterally inhomogeneous medium. Lateral inhomogeneities should be considered in resistivity soundings with large layouts. The results obtained with the described method permit a more realistic representation of geological features. The resolution of the method is determined by the number of elements in the resistivity network. The problem of core memory space has been resolved by using random access disc files. The results computed using a Fortran program are in good agreement with analytically obtained solutions.     

The Method for Inferring a Buried Fault from Resistivity Tomograms and Its Typical Electrical Features

Electrical resistivity tomography (ERT) has been used to experimentally detect shallow buried faults in urban areas in the past a few years, with some progress and experience obtained. According to the results from Olympic Park, Beijing, Shandong Province, Gansu Province and Shanxi Province, we have generalized the method and procedure for inferring the discontinuity of electrical structures (DES) indicating a buried fault in urban areas from resistivity tomograms and its typical electrical features. In general, the layered feature of the electrical structure is first analyzed to preliminarily define whether or not a DES exists in the target area. Resistivity contours in resistivity tomograms are then analyzed from the deep to the shallow. If they extend upward from the deep to the shallow and shape into an integral dislocation, sharp flexure (convergence) or gradient zone, it is inferred that the DES exists, indicating a buried fault. Finally, horizontal tracing is be carried out to define the trend of the DES. The DES can be divided into three types-type AB, ABA and AC. In the present paper, the Zhangdian-Renhe fault system in Zibo city is used as an example to illustrate how to use the method to infer the location and spatial extension of a target fault. Geologic drilling holes are placed based on our research results, and the drilling logs testify that our results are correct. However, the method of this paper is not exclusive and inflexible. It is expected to provide reference and assistance for inferring the shallow buried faults in urban areas from resistivity tomograms in the future.     

APPLICATION OF BOREHOLE GEOPHYSICS AT AN EXPERIMENTAL WASTE STORAGE SITE*

A suite of electrical, radiation, and mechanical borehole probes were run in a 76-mm-diameter borehole drilled to a slant depth of 380 m in leptite and granite. The hole is located in Precambrian bedrock in central Sweden where a site is dedicated to in-situ experiments pertaining to the disposal of radioactive wastes. The challenge to borehole logging methods for such site investigations is to resolve geological features and fluid flow parameters in geological sites which are initially chosen for their homogeneity, low porosity, and minimal fracturing. The Stripa borehole is characterized by high electrical resistivity values in the 20–100 kΩm range, by acoustic velocities around 5800 m s^-1 (which is close to laboratory values on intact specimens), and by total porosity of around one volume percent. In this context, probe resolution was adequate to produce interpretable information on almost all of the logs. Two principal rock types were encountered in the hole: granite, of quartz monzonitic composition, and leptite. The granite and leptite intercepts are subdivided into units characterized by mafic mineral content, sulfide mineral content, and electrical and radiation properties. Iron-rich zones in the leptite are highly anomalous on the gamma-gamma and neutron logs; thin mafic zones in the granite can also be distinguished. Occurrences of a few percent pyrite are detected by the electrical, gamma-gamma, and neutron logs. Although overall porosity is quite low throughout the hole, analysis of the resistivity and neutron logs indicates the porosity increases by a few volume percent at fracture zones. The differential resistance and caliper probes detect borehole diameter roughness of less than 1 mm, helping to confirm acoustic waveform anomalies which are indicative of fracture zones. Compres-sional wave transit time and shear-wave interference patterns usually occur coincident with open fractures observed in core, the correlation being especially good at major fracture zones.     

Multidimensional Investigation of Bedrock Heterogeneity/Unconformities at a DNAPL‐Impacted Site

Organic solvent (i.e., dense nonaqueous phase liquid, DNAPL) migration in the subsurface is known to be extremely sensitive to geologic heterogeneity. There is often a focus on heterogeneity that results from changing depositional conditions over short spatial scales. Similar or even more extreme spatial heterogeneity can result postdeposition due to erosional processes. This study applies a synergistic approach based on a combination of high‐resolution lithologic logs of continuous cores, borehole geophysical logs, surface electrical resistivity, and seismic refraction tomography models to assess spatial heterogeneity in a shallow bedrock sequence subject to multiple unconformities and contaminated with a mixture of organic chemicals. The persistence of DNAPL in the source zone and an associated dissolved‐phase plume led to variable impacts on formation resistivity across the study site. Seismic refraction in combination with electrical resistivity tomography improved interpretation of highly irregular erosional boundaries by delineating sharp lateral transitions in lithologic composition near the source zone and across the dissolved‐phase plume. Electrical resistivity was effective at differentiating between clean and mud‐rich sandstones and their unconformable contact with an underlying dolostone. Geophysical measurements revealed eroded dolostone mounds encased by a network of younger mud‐rich sandstones channelized by clean semi‐lithified sand, all of which was buried beneath variable glacial drift. Our synergistic multidimensional approach resulted in the development of a detailed three‐dimensional shallow bedrock geospatial model, which has led to an improved understanding of DNAPL migration and contaminant plume heterogeneity.     

CHOOSING THE AVERAGING INTERVAL WHEN CALCULATING PRIMARY REFLECTION COEFFICIENTS FROM WELL LOGS1

Most seismic data is processed using a sample interval of 4 ms two-way time (twt). The study of the statistical properties of primary reflection coefficients showed that the power spectrum of primaries can change noticeably when the logs are averaged over blocks of 0.5, 1 and 2 ms twt (block-averaging). What is a suitable block-averaging interval for producing broadband synthetics, and in particular how should the power spectrum of primaries be constructed when it is to be used to correct 4 ms sampled deconvolved seismic data for the effects of coloured primary reflectivity? In this paper we show that for a typical sonic log, a block-averaging interval of 1 ms twt should satisfy some important requirements. Firstly, it is demonstrated that if the reflection coefficients in an interval are not too large the effect of all the reflection impulses can be represented by another much sparser set at intervals of Δt twt, The coefficient amplitudes are given by the differences in the logarithmic acoustic impedances, thus justifying block-averaging. However, a condition for this to hold up to the aliasing (Nyquist) frequency is that Δt takes a maximum value of about 1 ms twt. Secondly, an event on a log should be represented in the seismic data. For this the acoustic impedance contrast must have sufficient lateral extent or continuity. By making some tentative suggestions on the relation between continuity and bed-thickness, a bed-thickness requirement of 0.15 m or more is obtained. Combining this requirement with the maximum number of beds allowable in an interval in order that multiple reflections do not contribute significantly to the reflections in the interval, again suggests a value of about 1 ms for the block-averaging interval. With this in mind an experiment was performed on three sonic logs. The logs were block-averaged at 1 ms, and primary reflection coefficients calculated. These primaries were then anti-alias filtered and resampled to get a series of primaries at 4 ms, followed by ARMA spectrum fitting. The same logs were also block-averaged at 4 ms directly and primaries computed, followed by ARMA spectrum fitting. In all three cases the first approach gave the ARM A model spectrum with greatest dynamic range, which strongly suggests that direct 4 ms block-averaging introduces significant aliased energy into low frequencies of the primaries spectrum. The conclusion is that routine computation of broadband synthetics (primaries only or primaries plus multiples) should be carried out using a block-averaging interval of 1 ms twt, followed by anti-alias filtering and thinning to the desired final sample interval. In theory it would be advantageous to go to even finer intervals-say 0.5 ms-but in practice at this level the averaging of slowness imposed by the somic logging tool appears to attenuate high-wave number fluctuations, i.e. it interferes with the‘real’data. The 1ms choice is thus a reasonable compromise which will help minimize non-trivial aliasing effects and should give better matches to the seismic data.     

Constructing the Apparent Geological Model by Fusing Surface Resistivity Survey and Borehole Records

We constructed an apparent geological model with resistivity data from surface resistivity surveys. We developed a data fusion approach by integrating dense electrical resistivity measurements collected with Schlumberger arrays and wellbore logs. This approach includes an optimization algorithm and a geostatistic interpolation method. We first generated an apparent formation factor model from the surface resistivity measurements and groundwater resistivity records with an inverse distance method. We then converted the model into a geology model with the optimized judgment criteria from the algorithms relating the apparent formation factors to the borehole geology. We also employed a non-parametric bootstrap method to analyze the uncertainty of the predicted sediment types, and the predictive uncertainties of clay, gravel, and sand were less than 5%. Overall, our model is capable of capturing the spatial features of the sediment types. More importantly, this approach can be arranged in a self-updated sequence to enable adjustments to the model to accommodate newly collected core records or geophysical data. This approach yields a more detailed apparent geological model for use in future groundwater simulations, which is of benefit to multi-discipline studies.     

裸眼井声波全波测井中纵波和斯通利波的数值研究

本文利用射线展开理论、割线积分技术及留数定理,研究了裸眼井声波全波测井中纵波和斯通利波的性质。下面简要说明数值计算方法。 计算方法 考虑在无限大均匀地层中,有一充满流体的深井。在井轴上放一点状声源,则在井轴上离声源距离为z处的声压为     

双侧向测井曲线形状与地层侵入关系研究

通过对具有不同侵入特征地层的双侧向理论曲线的计算,研究双侧向测井曲线的形状与泥浆侵入地层状态的关系.发现双侧向视电阻率曲线在高侵剖面的地层界面处出现‘反冲’,曲线‘峰’和‘谷’均较尖锐,而在低侵剖面,深侧向和浅侧向视电阻率曲线在地层界面处出现‘过渡’趋势,且曲线形状较圆滑.研究表明：产生上述现象的原因是由于原状地层对电流的吸引和排斥,使主电极和屏蔽电极供出的电流增大或减小,因而屏流比变化,进而导致视电阻率的变化.     

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.