scholarly journals Simulation of DC dual-laterolog measurements in complex formations: A Fourier-series approach with nonorthogonal coordinates and self-adapting finite elements

Geophysics ◽  
2009 ◽  
Vol 74 (1) ◽  
pp. E31-E43 ◽  
Author(s):  
Myung Jin Nam ◽  
David Pardo ◽  
Carlos Torres-Verdín

Dual laterolog (DLL) makes use of a galvanic conduction principle to focus electrical currents into rock formations, thereby minimizing shoulder and borehole effects in the measurement of formation resistivity. The tool includes two separate focusing systems: deep-sensing (LLd) and shallow-sensing modes (LLs). Laterolog current-focusing systems were designed for operation primarily in vertical boreholes penetrating horizontal layers; only recently their design has been revised for operation in deviated wells in the presence of electrical anisotropy. We simulated three-dimensional (3D) DLL measurements in dipping, invaded, and electrically anisotropic formations and appraised the corresponding effects on apparent resistivity logs. Simulations were performed by combining the use of a Fourier-series expansion in a nonorthogonal system of coordinates with an existing 2D goal-oriented, higher-order, and self-adaptive finite-element method. This numerical algorithm yields accurate solutions in limited CPU time because only a few Fourier modes are needed to simulate practical applications. For the calculation of focused currents, we introduced an embedded postprocessing method that incorporates a synthetic focusing principle to compute current intensities at each iterative step of optimal mesh refinements. Our numerical method accurately simulates 3D DLL measurements in rock formations that exhibit extreme contrasts of electrical resistivity. Simulations indicate that LLs resistivity logs are more sensitive to both invaded and anisotropic layers than LLd resistivity logs. In deviated wells, shoulder-bed effects on apparent resistivity logs increase with an increase of dip angle, and are emphasized across thin conductive layers. Electrical anisotropy effects on apparent resistivity logs increase substantially with dip angle.

Geophysics ◽  
2010 ◽  
Vol 75 (3) ◽  
pp. F83-F95 ◽  
Author(s):  
Myung Jin Nam ◽  
David Pardo ◽  
Carlos Torres-Verdín

Borehole triaxial induction instruments were designed to diagnose and measure rock electrical conductivity parallel and perpendicular to the bedding plane. Experience has shown that the interpretation of triaxial induction measurements often requires numerical modeling for a proper diagnosis of rock electrical conductivity anisotropy in the presence of geometric effects such as dipping wells, layer boundaries, and invasion. We introduce a new algorithm to simulate triaxial induction measurements that combines a Fourier series expansion in a nonorthogonal system of coordinates with a 2D goal-oriented, self-adaptive, high-order [Formula: see text] finite-element method. This procedure enables the accurate and reliable simulation of triaxial induction measurements across reservoir rock formations with extreme contrasts of electrical conductivity while reducing the 3D computational complexity associated with deviated wells. Numerical results indicate that borehole dip effects on triaxial induction measurements are larger than on standard coaxial induction measurements. The sensitivity of triaxial induction measurements to transversely isotropic rock formations decreases with increasing dip angle.


2021 ◽  
Author(s):  
Alvin Wong

This research developed a numerical method that solves complicated fluid flow problems without requiring end-user expertise with the solver. This method is capable of obtaining a spatially accurate solution in the same time or better as a skilled user with a conventional solver. An explicit preconditioned multigrid solver was used in this research with a multistage relaxation method. The prosposed method utilizies a database with optimized relaxation method parameters for different local flow and mesh conditions. The parameters are optimized for the relaxation such that the error modes in a complex Fourier series expansion of the residual can be quickly reduced. The convergence time and iteration count of this method was compared against the same solver using default input values, as well as a pre-optimized solver, to simulate a skilled user for various geometries. Improvements in both comparisons were demonstrated.


2021 ◽  
Author(s):  
Alvin Wong

This research developed a numerical method that solves complicated fluid flow problems without requiring end-user expertise with the solver. This method is capable of obtaining a spatially accurate solution in the same time or better as a skilled user with a conventional solver. An explicit preconditioned multigrid solver was used in this research with a multistage relaxation method. The prosposed method utilizies a database with optimized relaxation method parameters for different local flow and mesh conditions. The parameters are optimized for the relaxation such that the error modes in a complex Fourier series expansion of the residual can be quickly reduced. The convergence time and iteration count of this method was compared against the same solver using default input values, as well as a pre-optimized solver, to simulate a skilled user for various geometries. Improvements in both comparisons were demonstrated.


2017 ◽  
Vol 84 (9) ◽  
Author(s):  
Siyuan Bao ◽  
Shuodao Wang ◽  
Bo Wang

A modified Fourier–Ritz approach is developed in this study to analyze the free in-plane vibration of orthotropic annular sector plates with general boundary conditions. In this approach, two auxiliary sine functions are added to the standard Fourier cosine series to obtain a robust function set. The introduction of a logarithmic radial variable simplifies the expressions of total energy and the Lagrangian function. The improved Fourier expansion based on the new variable eliminates all the potential discontinuities of the original displacement function and its derivatives in the entire domain and effectively improves the convergence of the results. The radial and circumferential displacements are formulated with the modified Fourier series expansion, and the arbitrary boundary conditions are simulated by the artificial boundary spring technique. The number of terms in the truncated Fourier series and the appropriate value of the boundary spring retraining stiffness are discussed. The developed Ritz procedure is used to obtain accurate solution with adequately smooth displacement field in the entire solution domain. Numerical examples involving plates with various boundary conditions demonstrate the robustness, precision, and versatility of this method. The method developed here is found to be computationally economic compared with the previous method that does not adopt the logarithmic radial variable.


1995 ◽  
Vol 02 (04) ◽  
pp. 489-494 ◽  
Author(s):  
E.E. MOLA ◽  
A.G. APPIGNANESSI ◽  
J.L. VICENTE ◽  
L. VAZQUEZ ◽  
R.C. SALVAREZZA ◽  
...  

The model for the angular orientational energy (AOE) has been extended to hexagonal submonolayer domains of Ag electrodeposited at a constant overpotential on a C(0001) surface. These domains which are characterized by an epitaxy angle θ=15±5° and an Ag−Ag distance d Ag−Ag =0.330± 0.016 nm, can be considered as precursors of 3D Ag crystal formation, according to a Volmer-Weber type mechanism. Calculations are based upon a simple Hamiltonian evaluated by introducing the concept of the commensurable unit cell. A Fourier series expansion for the substrate potential was used. Results from the model predict the existence of a commensurable cell in agreement with the experimental data derived from STM imaging.


Sign in / Sign up

Export Citation Format

Share Document