Formation Electrical Anisotropy Derived From Induction-Log Measurements in a Horizontal Well

2001 ◽  
Vol 4 (06) ◽  
pp. 483-488
Author(s):  
Jon Bang ◽  
Arne Solstad ◽  
Svein Mjaaland

Summary An existing theory describes how electrical anisotropy in the formationaffects the response of resistivity logging tools. We have related this theory to the processing of logging while drilling (LWD) induction logs and are thus able to calculate the anisotropic resistivities directly from the logs. The method has been demonstrated by application to logs from a horizontal well section. Anisotropy ratios of 2 to 5, and occasionally higher values, were obtained for this formation. We also addressed the accuracy of these numbers by using independent sets of input logs. The results indicate that the logs are influenced by factors like invasion, in addition to the anisotropy. Our approach provides a fast and efficient computer algorithm. The output is calculated at the depths of the input logs; hence, the resulting anisotropy becomes a depth-dependent formation property. Introduction Electrical anisotropy has gained considerable attention in recent years. If present in the formation, neglection of this property when interpreting resistivity logs may lead to erroneous saturation estimates and may thus have great consequences upon development and production strategies and the overall economic situation. Electrical anisotropy denotes that the resistivity shows directional dependence. In sedimentary formations, it is commonly assumed that the anisotropy is caused by the deposition process, which yields different small-scale (grain and pore-size scale) structural properties in the vertical and horizontal directions. Anisotropy may also occur on a lithology scale[i.e., as a result of thin layers (compared to the extension of the electricfield) having individual isotropic properties]. Because the effect is determined by the sedimentary structure, a formation can be expected to show anisotropy in several properties, such as electric, acoustic, and fluid-flow resistance (permeability) properties, simultaneously. A common way of describing anisotropy is to distinguish between the vertical direction and directions in the horizontal plane. In this paper, we shall denote the resistivities in these directions by RV andRH, respectively. However, the terms "vertical" and"horizontal" refer to the original deposition process and may no longer correspond to the actual orientation of the formation owing to small- or large-scale geological activity. For dipping beds, it is common practice to assume one resistivity (R H) in the bedding plane and one (RV) in the direction normal to the bed, unless evidence of intrabed disturbances suggests other orientations of the anisotropy. Numerous publications have addressed the influence of electrical anisotropy on resistivity logs. Among the effects that have been studied are anisotropy in dipping and thinly laminated formations1–3 and in crossbedded formations.4 Effort has been put on theoretical tool response modeling and simulation 5–7 and on anisotropy corrections to logs.8,9 From field cases, anisotropy ratios(RV/R H) up to the order of 5 to 10 have been reported.7,8,10 In this paper, we demonstrate a method for calculating the electrical anisotropy directly from well logs, based on the theory developed by Hagiwara.6 The method has been implemented and applied to log data from a horizontal North Sea well. Theory Hagiwara6 has analyzed the resistivity log's response in anisotropic formations. According to this reference, two different measurements are sufficient to determine the anisotropy unambiguously, as long as the anisotropy orientation is known. The measurements may differ with respect to one or more of the following:antenna spacing (which is a prerequisite for phase- and attenuation-derived resistivity),frequency, ordeviation angle between tool axis and anisotropy orientation. In our work, we consider the LWD induction response. For this instrument class, Hagiwara shows that the complex voltage V recorded by one transmitter-receiver pair of electrodes isEquation 1 where i=the imaginary unit (i=-11/2) and L=the antenna spacing. Further,Equation 2 where a2= RH/RV is the anisotropy ratio between horizontal and vertical resistivitiesRH and RV, and ?=the deviation of tool direction from the R V direction. Notice the interpretation of the terms "vertical" and "horizontal," as discussed in the introduction. The wave number k is defined byEquation 3 where ?=the measurement angular frequency, µ=the magnetic permeability, andeH=the horizontal dielectric constant. In this study, we used the free space magnetic permeability µ=µ0=4p×10–7 N/A, and approximated eH from the logged resistivity through an empirical relation. Both these approximations are considered to have negligible influence on the results.

2021 ◽  
Author(s):  
Yessica Fransisca ◽  
Karinka Adiandra ◽  
Vinda Manurung ◽  
Laila Warkhaida ◽  
M. Aidil Arham ◽  
...  

Abstract This paper describes the combination of strategies deployed to optimize horizontal well placement in a 40 ft thick isotropic sand with very low resistivity contrast compared to an underlying anisotropic shale in Semoga field. These strategies were developed due to previously unsuccessful attempts to drill a horizontal well with multiple side-tracks that was finally drilled and completed as a high-inclined well. To maximize reservoir contact of the subject horizontal well, a new methodology on well placement was developed by applying lessons learned, taking into account the additional challenges within this well. The first approach was to conduct a thorough analysis on the previous inclined well to evaluate each formation layer’s anisotropy ratio to be used in an effective geosteering model that could better simulate the real time environment. Correct selections of geosteering tools based on comprehensive pre-well modelling was considered to ensure on-target landing section to facilitate an effective lateral section. A comprehensive geosteering pre-well model was constructed to guide real-time operations. In the subject horizontal well, landing strategy was analysed in four stages of anisotropy ratio. The lateral section strategy focused on how to cater for the expected fault and maintain the trajectory to maximize reservoir exposure. Execution of the geosteering operations resulted in 100% reservoir contact. By monitoring the behaviour of shale anisotropy ratio from resistivity measurements and gamma ray at-bit data while drilling, the subject well was precisely landed at 11.5 ft TVD below the top of target sand. In the lateral section, wellbore trajectory intersected two faults exhibiting greater associated throw compared to the seismic estimate. Resistivity geo-signal and azimuthal resistivity responses were used to maintain the wellbore attitude inside the target reservoir. In this case history well with a low resistivity contrast environment, this methodology successfully enabled efficient operations to land the well precisely at the target with minimum borehole tortuosity. This was achieved by reducing geological uncertainty due to anomalous resistivity data responding to shale electrical anisotropy. Recognition of these electromagnetic resistivity values also played an important role in identifying the overlain anisotropic shale layer, hence avoiding reservoir exit. This workflow also helped in benchmarking future horizontal well placement operations in Semoga Field. Technical Categories: Geosteering and Well Placement, Reservoir Engineering, Low resistivity Low Contrast Reservoir Evaluation, Real-Time Operations, Case Studies


2002 ◽  
Vol 7 (2) ◽  
pp. 45-52
Author(s):  
L. Jakučionis ◽  
V. Kleiza

Electrical properties of conductive thin films, that are produced by vacuum evaporation on the dielectric substrates, and which properties depend on their thickness, usually are anisotropic i.e. they have uniaxial anisotropy. If the condensate grow on dielectric substrates on which plane electrical field E is created the transverse voltage U⊥ appears on the boundary of the film in the direction perpendicular to E. Transverse voltage U⊥ depends on the angle γ between the applied magnetic field H and axis of light magnetisation. When electric field E is applied to continuous or grid layers, U⊥ and resistance R of layers are changed by changing γ. It means that value of U⊥ is the measure of anisotropy magnitude. Increasing voltage U0 , which is created by E, U⊥ increases to certain magnitude and later decreases. The anisotropy of continuous thin layers is excited by inequality of conductivity tensor components σ0 ≠ σ⊥. The reason of anisotropy is explained by the model which shows that properties of grain boundaries are defined by unequal probability of transient of charge carrier.


2002 ◽  
Vol 737 ◽  
Author(s):  
Theodore I. Kamins ◽  
Gilberto Medeiros-Ribeiro ◽  
Douglas A. A. Ohlberg ◽  
R. Stanley Williams

ABSTRACTWhen Ge is deposited epitaxially on Si, the strain energy from the lattice mismatch causes the Ge in layers thicker than about four monolayers to form distinctive, three-dimensional islands. The shape of the islands is determined by the energies of the surface facets, facet edges, and interfaces. When phosphorus is added during the deposition, the surface energies change, modifying the island shapes and sizes, as well as the deposition process. When phosphine is introduced to the germane/hydrogen ambient during Ge deposition, the deposition rate decreases because of competitive adsorption. The steady-state deposition rate is not reached for thin layers. The deposited, doped layers contain three different island shapes, as do undoped layers; however, the island size for each shape is smaller for the doped layers than for the corresponding undoped layers. The intermediate-size islands are the most significant; the intermediate-size doped islands are of the same family as the undoped, multifaceted “dome” structures, but are considerably smaller. The largest doped islands appear to be related to the defective “superdomes” discussed for undoped islands. The distribution between the different island shapes depends on the phosphine partial pressure. At higher partial pressures, the smaller structures are absent. Phosphorus appears to act as a mild surfactant, suppressing small islands.


2000 ◽  
Vol 15 (8) ◽  
pp. 1702-1708
Author(s):  
Ruichao Zhang ◽  
Ren Xu

A novel two-step metalorganic chemical vapor deposition process was used in this study to prepare Sr1−xBaxNb2O6 (SBN) thin films. Two thin layers of single-phase SrNb2O6 and BaNb2O6 were deposited alternately on a silicon substrate, and the solid solution of SBN was obtained by high-temperature annealing. The stoichiometry control of the SrNb2O6 and the BaNb2O6 thin films was achieved through deposition process control, according to the evaporation characteristics of double metal alkoxide. The evaporation behavior of double metal alkoxide precursors SrNb2(1-OC4H9)12 and BaNb2(1-OC4H9)12 was studied, and the results were compared with the evaporation of single alkoxide Nb(1-OC4H9)5.


2016 ◽  
Author(s):  
Fabrice Chane Ming ◽  
Damien Vignelles ◽  
Fabrice Jegou ◽  
Gwenael Berthet ◽  
Jean-Batiste Renard ◽  
...  

Abstract. Coupled balloon-borne observations of Light Optical Aerosol Counter (LOAC), M10 meteorological global positioning system (GPS) sondes, ozonesondes and GPS radio occultation data, are examined to identify gravity-wave (GW) induced fluctuations on tracer gases and on the vertical distribution of stratospheric aerosol concentrations during the 2013 ChArMEx (Chemistry-Aerosol Mediterranean Experiment) campaign. Observations reveal signatures of GWs with short vertical wavelengths less than 4 km in dynamical parameters and tracer constituents which are also correlated with the presence of thin layers of strong local enhancements of aerosol concentrations in the upper troposphere and the lower stratosphere. In particular, this is evident from a case study above Ile du Levant (43.02 °N, 6.46 °E) on 26–29 July 2013. Observations show a strong activity of dominant mesoscale inertia GWs with horizontal and vertical wavelengths of 370–510 km and 2–3 km respectively, and periods of 10–13 h propagating southward at altitudes of 13–20 km and eastward above 20 km during 27–28 July which is also captured by the European Center for Medium range Weather Forecasting (ECMWF) analyses. Ray-tracing experiments indicate the jet-front system to be the source of observed GWs. Simulated vertical profiles of dynamical parameters with large stratospheric vertical wind maximum oscillations ± 40 mms−1 are produced for the dominant mesoscale GW using the simplified linear GW theory. Parcel advection method reveals signatures of GWs in the ozone mixing ratio and the specific humidity. Simulated vertical wind perturbations of the dominant GW and small-scale perturbations of aerosol concentration (aerosol size of 0.2–0.7 μm) are in phase in the lower stratosphere. Present results support the importance of vertical wind perturbations in the GW-aerosol relation. The observed mesoscale GW induces a strong modulation of the amplitude of tracer gases and the stratospheric aerosol background.


2012 ◽  
Vol 468 ◽  
pp. 57-69 ◽  
Author(s):  
KO Möller ◽  
M St John ◽  
A Temming ◽  
J Floeter ◽  
AF Sell ◽  
...  

1976 ◽  
Vol 71 ◽  
pp. 69-99 ◽  
Author(s):  
J. O. Stenflo

The observed properties of small-scale solar magnetic fields are reviewed. Most of the magnetic flux in the photosphere is in the form of strong fields of about 100–200 mT (1–2 kG), which have remarkably similar properties regardless of whether they occur in active or quiet regions. These fields are associated with strong atmospheric heating. Flux concentrations decay at a rate of about 107 Wb s-1, independent of the amount of flux in the decaying structure. The decay occurs by smaller flux fragments breaking loose from the larger ones, i.e. a transfer of magnetic flux from smaller to larger Fourier wave numbers, into the wave-number regime where ohmic diffusion becomes significant. This takes place in a time-scale much shorter than the length of the solar cycle.The field amplification occurs mainly below the solar surface, since very little magnetic flux appears in diffuse form in the photosphere, and the life-time of the smallest flux elements is very short. The observations further suggest that most of the magnetic flux in quiet regions is supplied directly from below the solar surface rather than being the result of turbulent diffusion of active-region magnetic fields.


1994 ◽  
Vol 20 ◽  
pp. 61-66 ◽  
Author(s):  
A.L. Veazey ◽  
M.O. Jeffries ◽  
K. Morris

The small-scale variability of physical properties and structural characteristics of multiple pairs of fast-ice cores obtained during the austral summer of 1991-92 at two Antarctic sites, McMurdo Sound (MCM) and Pine Island Bay (PIB), are examined and discussed with respect to the growth and decay of the sea ice. The ice at the MCM site was thicker than that at the PIB site and was covered by a somewhat thinner snowpack. While mean salinity and temperature of the ice at the two sites were similar, small-scale variations in both salinity and temperature were greater at PIB than at MCM. The ice sheet at MCM was a two-layer medium consisting of congelation ice overlying platelet ice. The ice from the PIB site, however, was composed of mainly frazil ice and layers of congelation ice with occasional thin layers of snow-ice at the surface of the cores. Crystal sub-structure measurements, c-axis orientation and brine-layer spacing from the MCM cores revealed that the congelation ice had moderately aligned, horizontally oriented c axes, suggesting that east-west currents off the southwest tip of Hut Point Peninsula control crystal-growth orientation.d others: Variability of physical and structural characteristics of Antarctic fast ice


2009 ◽  
Vol 27 (1) ◽  
pp. 417-425 ◽  
Author(s):  
N. V. Erkaev ◽  
V. S. Semenov ◽  
I. V. Kubyshkin ◽  
M. V. Kubyshkina ◽  
H. K. Biernat

Abstract. One-fluid ideal MHD model is applied for description of current sheet flapping disturbances appearing due to a gradient of the normal magnetic field component. The wave modes are studied which are associated to the flapping waves observed in the Earth's magnetotail current sheet. In a linear approximation, solutions are obtained for model profiles of the electric current and plasma densities across the current sheet, which are described by hyperbolic functions. The flapping eigenfrequency is found as a function of wave number. For the Earth's magnetotail conditions, the estimated wave group speed is of the order of a few tens kilometers per second. The current sheet can be stable or unstable in dependence on the direction of the gradient of the normal magnetic field component. The obtained dispersion function is used for calculation of the flapping wave disturbances, which are produced by the given initial Gaussian perturbation at the center of the current sheet and propagating towards the flanks. The propagating flapping pulse has a smooth leading front, and a small scale oscillating trailing front, because the short wave oscillations propagate much slower than the long wave ones.


2019 ◽  
Vol 84 (6) ◽  
pp. 1203-1219
Author(s):  
Anthony Radjen ◽  
Gabriele Gradoni ◽  
Richard Tew

Abstract Reflection and transmission phenomena associated with high-frequency linear wave incidence on irregular boundaries between adjacent acoustic or electromagnetic media, or upon the irregular free surface of a semi-infinite elastic solid, are studied in two dimensions. Here, an ‘irregular’ boundary is one for which small-scale undulations of an arbitrary profile are superimposed upon an underlying, smooth curve (which also has an arbitrary profile), with the length scale of the perturbation being prescribed in terms of a certain inverse power of the large wave-number of the incoming wave field. Whether or not the incident field has planar or cylindrical wave-fronts, the associated phase in both cases is linear in the wave-number, but the presence of the boundary irregularity implies the necessity of extra terms, involving fractional powers of the wave-number in the phase of the reflected and transmitted fields. It turns out that there is a unique perturbation scaling for which precisely one extra term in the phase is needed and hence for which a description in terms of a Friedlander–Keller ray expansion in the form as originally presented is appropriate, and these define a ‘distinguished’ class of perturbed boundaries and are the subject of the current paper.


Sign in / Sign up

Export Citation Format

Share Document