地质导向钻井关键技术及方法

On the basis of analyzing the principle and realization of geo-steering drilling system, the key technologies and methods in it are systematically studied in this paper. In order to recognize lithology, distinguish stratum and track reservoirs, the techniques of MWD and data process about natural gamma, resistivity, inductive density and porosity are researched. The methods for pre-processing and standardizing MWD data and for converting geological data in directional and horizontal drilling are discussed, consequently the methods of data conversion between MD and TVD and those of formation description and adjacent well contrast are proposed. Researching the method of identifying sub-layer yields the techniques of single well explanation, multi-well evaluation and oil reservoir description. Using the extremum and variance clustering analysis realizes logging phase analysis and stratum subdivision and explanation, which provides a theoretical method and lays a technical basis for tracing oil reservoirs and achieving geo-steering drilling. Researching the technique for exploring the reservoir top with a holdup section provides a planning method of wellpath control scheme to trace oil and gas reservoir dynamically, which solves the problem of how to control well trajectory on condition that the layer’s TVD is uncertain. The control scheme and planning method of well path for meeting the demands of target hitting, soft landing and continuous steering respectively provide the technological guarantee to land safely and drill successfully for horizontal, extended-reach and multi-target wells. The integrative design and control technologies are researched based on geology, reservoir and drilling considering reservoir disclosing ratio as a primary index, and the methods for planning and control optimum wellpath under multi-target restriction, thus which lets the target wellpath lie the favorite position in oil reservoir during the process of geo-steering drilling. The BHA (bottomhole assembly) mechanical model is discussed using the finite element method, and the BHA design methods are given on the basis of mechanical analyses according to the shape of well trajectory and the characteristics of BHA’s structure and deformation. The methods for predicting the deflection rate of bent housing motors and designing their assemblies are proposed based on the principle of minimum potential energy, which can clearly show the relation between the BHA’s structure parameters and deflection rate, especially the key factors’ effect to the deflection rate. Moreover, the interaction model between bit and formation is discussed through the process of equivalent formation and equivalent bit considering the formation anisotropy and bit anisotropy on the basis of analyzing the influence factors of well trajectory. Accordingly, the inherence relationship among well trajectory, formation, bit and drilling direction is revealed, which lays the theory basis and technique for predicting and controlling well trajectory.

地幔对流理论及其在解释地球物理观测的应用

With the improving of mantle convection theory, the developing of computing method and increasing of the measurement data, we can numerically simulate more clearly about the effects on some geophysical observed phenomenons such as the global heat flow and global lithospheric stress field in the Earth's surface caused by mantle convection, which is the primary mechanism for the transport of heat from the Earth's deep interior to its surface and the underlying force mechanism of dynamics in the Earth.Chapter 1 reviews the historical background and present research state of mantle convection theory.In Chapter 2, the basic conception of thermal convection and the basic theory about mantle flow.The effects on generation and distribution of global lithospheric stres s field induced by mantle flow are the subject of Chapter 3. Mantle convection causes normal stress and tangential stresses at the bottom of the lithosphere, and then the sublithospheric stress field induces the lithospheric deformation as sixrface force and results in the stress field within the lithosphere. The simulation shows that the agreement between predictions and observations is good in most regions. Most of subduction zones and continental collisions are under compressive. While ocean ridges, such as the east Pacific ridge, the Atlantic ridge and the east African rift valley, are under tensile. And most of the hotspots preferentially occur in regions where calculated stress is tensile. The calculated directions of the most compressive principal horizontal stress are largely in accord with that of the observation except for some regions such as the NW-Pacifie subduction zone and Qinghai-Tibet Plateau, in which the directions of the most compressive principal horizontal stress are different. It shows that the mantel flow plays an important role in causing or affecting the large-scale stress field within the lithosphere.The global heat flow simulation based on a kinematic model of mantle convection is given in Chapter 4. Mantle convection velocities are calculated based on the internal loading theory at first, the velocity field is used as the input to solve the thermal problem. Results show that calculated depth derivatives of the near surface temperature are closely correlated to the observed surface heat flow pattern. Higher heat flow values around midocean ridge systems can be reproduced very well. The predicted average temperature as a function of function of depth reveals that there are two thermal boundary layers, one is close to the surface and another is close to the core-mantle boundary, the rest of the mantle is nearly isothermal. Although, in most of the mantle, advection dominates the heat transfer, the conductive heat transfer is still locally important in the boundary layers and plays an important role for the surface heat flow pattern. The existence of surface plates is responsible for the long wavelength surface heat flow pattern.In Chapter 5, the effects on present-day crustal movement in the China Mainland resulted from the mantle convection are introduced. Using a dynamic method, we present a quantitative model for the present-day crustal movement in China. We consider not only the effect of the India-Eurasia collision, the gravitational potential energy difference of the Tibet Plateau, but also the contribution of the shear traction on the bottom of the lithosphere induced by the global mantle convection. The comparison between our results and the velocity field obtained from the GPS observation shows that our model satisfactorily reproduces the general picture of crustal deformation in China. Numerical modeling results reveal that the stress field on the base of the lithosphere induced by the mantle flow is probably a considerable factor that causes the movement and deformation of the lithosphere in continental China with its eflfcet focuing on the Eastern China A numerical research on the small-scale convection with variable viscosity in the upper mantle is introduced in Chapter 6. Based on a two-dimensional model, small-scale convection in the mantle-lithosphere system with variable viscosity is researched by using of finite element method. Variation of viscosity in exponential form with temperature is considered in this paper The results show that if viscosity is strongly temperature-dependent, the upper part of the system does not take a share in the convection and a stagnant lid, which is identified as lithosphere, is formed on the top of system because of low temperature and high viscosity. The calculated surface heat flow, topography and gravity anomaly are associated well with the convection pattern, namely, the regions with high heat flow and uplift correspond to the upwelling flow, and vice versa.In Chapter 7, we give a brief of future research subject: The inversion of lateral density heterogeneity in the mantle by minimizing the viscous dissipation.

Energy spectrum and persistent currents in finite-width mesoscopic ring with radial potential barrier threading a magnetic flux through its hole

The energy spectrum and the persistent currents are calculated for a finite-width mesoscopic annulus with radial potential barrier, threading a magnetic flux through the hole of the ring. Owing to the presence of tunneling barrier, the coupling effect leads to the splitting of each radial energy subband of individual concentrical rings into two one. Thus, total currents and currents carried by single high-lying eigenstate as a function of magnetic flux exhibit complicated patterns. However, periodicity and antisymmetry of current curves in the flux still preserve.

Binding Energy and Spin-Orbit Splitting of a Hydrogenic Donor Impurity in AlGaN/GaN Triangle-Shaped Potential Quantum Well

In the framework of effective-mass envelope function theory, including the effect of Rashba spin-orbit coupling, the binding energy E-b and spin-orbit split energy Gamma of the ground state of a hydrogenic donor impurity in AlGaN/GaN triangle-shaped potential heterointerface are calculated. We find that with the electric field of the heterojunction increasing, (1) the effective width of quantum well (W) over bar decreases and (2) the binding energy increases monotonously, and in the mean time, (3) the spin-orbit split energy Gamma decreases drastically. (4) The maximum of Gamma is 1.22 meV when the electric field of heterointerface is 1 MV/cm.

DETERMINATION OF THE SUBBAND ENERGY IN THE V-SHAPED POTENTIAL WELL OF DELTA-DOPED GAAS BY DEEP-LEVEL TRANSIENT SPECTROSCOPY

Using deep level transient spectroscopy (DLTS) the conduction-subband energy levels in a V-shaped potential well induced by Si-delta doping in GaAs were determined. Self-consistent calculation gives four subbands in the well below the Fermi level. Experimentally, two DLTS peaks due to electron emission from these subbands were observed. Another two subbands with low electron concentration are believed to be merged into the adjacent DLTS peak. A good agreement between self-consistent calculation and experiment was obtained. (C) 1994 American Institute of Physics.

RESONANT TUNNELING, EIGENVALUE AND ENERGY-BAND CALCULATION FOR POTENTIAL AND PERIODIC POTENTIAL STRUCTURES

Recursion formulae for the reflection and the transmission probability amplitudes and the eigenvalue equation for multistep potential structures are derived. Using the recursion relations, a dispersion equation for periodic potential structures is presented. Some numerical results for the transmission probability of a double barrier structure with scattering centers, the lifetime of the quasi-bound state in a single quantum well with an applied field, and the miniband of a periodic potential structure are presented.

Structures of energy spectrum and persistent currents in mesoscopic rings with radial potential barrier in magnetic field

The energy spectrum and the persistent currents are calculated for finite-width mesoscopic annular structures with radial potential barrier in the presence of a magnetic field. The introduction of the tunneling barrier leads to the creation of extra edge states around the barrier and the occurrence of oscillatory structures superimposed on the bulk Landau level plateaus in the energy spectrum. We found that the Fermi energy E-F increases with the number of electrons N emerging many kinks. The single eigenstate persistent current exhibits complicated structures with vortex-like texture, ''bifurcation'', and multiple ''furcation'' patterns as N is increased. The total currents versus N display wild fluctuations.

Effect of Surface Potential Barrier on the Electron Energy Distribution of NEA Photocathodes

By calculating the energy distribution of electrons reaching the photocathode surface and solving the Schrodinger equation that describes the behavior of an electron tunneling through the surface potential barrier,we obtain an equation to calculate the emitted electron energy distribution of transmission-mode NEA GaAs photocathodes. Accord- ing to the equation,we study the effect of cathode surface potential barrier on the electron energy distribution and find a significant effect of the barrier-Ⅰ thickness or end height,especially the thickness,on the quantum efficiency of the cath- ode. Barrier Ⅱ has an effect on the electron energy spread, and an increase in the vacuum level will lead to a narrower electron energy spread while sacrificing a certain amount of cathode quantum efficiency. The equation is also used to fit the measured electron energy distribution curve of the transmission-mode cathode and the parameters of the surface barri- er are obtained from the fitting. The theoretical curve is in good agreement with the experimental curve.

Isospin effects of the Skyrme potential and the momentum dependent interaction at intermediate energy heavy ion collisions

We improve the isospin dependent quantum molecular dynamical model by including isospin effects in the Skyrme potential and the momentum dependent interaction to obtain an isospin dependent Skyrme potential and an isospin dependent momentum interaction. We investigate the isospin effects of Skyrme potential and momentum dependent interaction on the isospin fractionation ratio and the dynamical mechanism in intermediate energy heavy ion collisions. It is found that the isospin dependent Skyrme potential and the isospin dependent momentum interaction produce some important isospin effects in the isospin fractionation ratio

Probing balance energy using momentum- and isospin-dependent potential

Using the momentum- and isospin-dependent Boltmann-Uehling-Uhlenbeck (BUU) model, we investigate the transverse flow and balance energy in two isotopic colliding systems Ca-48+Fe-58 and Cr-48+Ni-58 by adopting different symmetry potentials. By comparing the results between the two colliding systems, we find that the difference between the balance energies of two isotopic systems can be considered as a sensitive probe to the density dependence of symmetry energy.

Dynamical symmetry of screened Coulomb potential and isotropic harmonic oscillator

It is shown that for the screened Coulomb potential and isotropic harmonic oscillator, there exists an infinite number of closed orbits for suitable angular momentum values. At the aphelion (perihelion) points of classical orbits, an extended Runge-Lenz vector for the screened Coulomb potential and an extended quadrupole tensor for the screened isotropic harmonic oscillator are still conserved. For the screened two-dimensional (2D) Coulomb potential and isotropic harmonic oscillator, the dynamical symmetries SO3 and SU(2) are still preserved at the aphelion (perihelion) points of classical orbits, respectively. For the screened 3D Coulomb potential, the dynamical symmetry SO4 is also preserved at the aphelion (perihelion) points of classical orbits. But for the screened 3D isotropic harmonic oscillator, the dynamical symmetry SU(2) is only preserved at the aphelion (perihelion) points of classical orbits in the eigencoordinate system. For the screened Coulomb potential and isotropic harmonic oscillator, only the energy (but not angular momentum) raising and lowering operators can be constructed from a factorization of the radial Schrodinger equation.

(Ga, Gd, As) film growth on GaAs substrate by low-energy ion-beam deposit

(Ga, Gd, As) film was fabricated by the mass-analyzed dual ion-beam epitaxy system with the energy of 1000 eV at room temperature. There was no new peak found except GaAs substrate peaks (0 0 2) and (0 0 4) by X-ray diffraction. Rocking curves were measured for symmetric (0 0 4) reflections to further yield the lattice mismatch information by employing double-crystal X-ray diffraction. The element distributions vary so much due to the ion dose difference from AES depth profiles. The sample surface morphology indicates oxidizing layer roughness is also relative to the Gd ion dose, which leads to islandlike feature appearing on the high-dose sample. One sample shows ferromagnetic behavior at room temperature. (C) 2003 Elsevier B.V. All rights reserved.

Energy principle of ferroelectric ceramics and single domain mechanical model

Many physical experiments have shown that the domain switching in a ferroelectric material is a complicated evolution process of the domain wall with the variation of stress and electric field. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of ferroelectric ceramic and used to study the nonlinear constitutive behavior of ferroelectric body in this paper. The principle of stationary total energy is put forward in which the basic unknown quantities are the displacement u (i) , electric displacement D (i) and volume fraction rho (I) of the domain switching for the variant I. Mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion. On the basis of the domain switching criterion, a set of linear algebraic equations for the volume fraction rho (I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. Then a single domain mechanical model is proposed in this paper. The poled ferroelectric specimen is considered as a transversely isotropic single domain. By using the partial experimental results, the hardening relation between the driving force of domain switching and the volume fraction of domain switching can be calibrated. Then the electromechanical response can be calculated on the basis of the calibrated hardening relation. The results involve the electric butterfly shaped curves of axial strain versus axial electric field, the hysteresis loops of electric displacement versus electric filed and the evolution process of the domain switching in the ferroelectric specimens under uniaxial coupled stress and electric field loading. The present theoretic prediction agrees reasonably with the experimental results given by Lynch.

Quantum and variational monte-carlo interaction potentials for li2 (x1-sigma-g+)

Optimized trial functions are used in quantum Monte Carlo and variational Monte Carlo calculations of the Li₂(X ¹Σ⁺_g) potential curve. The trial functions used are a product of a Slater determinant of molecular orbitals multiplied by correlation functions of electron—nuclear and electron—electron separation. The parameters of the determinant and correlation functions are optimized simultaneously by reducing the deviations of the local energy E_L (E_L  Ψ⁻¹_THΨ_T, where Ψ_T denotes a trial function) over a fixed sample. At the equilibrium separation, the variational Monte Carlo and quantum Monte Carlo methods recover 68% and 98% of the correlation energy, respectively. At other points on the curves, these methods yield similar accuracies.

The mechanism of energy storage by shearing the solar magnetic field

In this paper, a complete set of MHD equations have been solved by numerical calculations in an attempt to study the dynamical evolutionary processes of the initial equilibrium configuration and to discuss the energy storage mechanism of the solar atmosphere by shearing the magnetic field. The initial equilibrium configuration with an arch bipolar potential field obtained from the numerical solution is similar to the configuration in the vicinity of typical solar flare before its eruption. From the magnetic induction equation in the set of MHD equations and dealing with the non-linear coupling effects between the flow field and magnetic field, the quantitative relationship has been derived for their dynamical evolution. Results show that plasma shear motion at the bottom of the solar atmosphere causes the magnetic field to shear; meanwhile the magnetic field energy is stored in local regions. With the increase of time the local magnetic energy increases and it may reach an order of 4×10^25 J during a day. Thus the local storage of magnetic energy is large enough to trigger a big solar flare and can be considered as the energy source of solar flares. The energy storage mechanism by shearing the magnetic field can well explain the slow changes in solar active regions.