Interpretation models and charts of production profiles in horizontal wells

Stratified flow is common for the gravity segregation and flow regimes are very complex because of borehole inclination, therefore, all the conventional production logging tools cannot be effectively applied in horizontal wells, thus significantly increasing the difficulties in log interpretation. In this paper, firstly, the overseas progress in updated integration tools for horizontal wells and production profile interpretation methods has been discussed in brief. Secondly, by means of theory study and experimental simulations, we have obtained the production profile interpretation model and experimental interpretation charts, which have been calibrated by the improved downhole technology and optimization methods. Finally, we have interpreted X-well with the production profile interpretation software designed by us, and it proves that the methods are useful for the production profile interpretation in horizontal wells.

     

Uniform head in horizontal and vertical wells

The steady-state head within a fully penetrating well may be estimated by evaluating the Thiem equation at the radius of the well. A method is presented here to extend results from the Thiem equation to horizontal wells and to partially penetrating wells. The particular model used in this investigation is based upon the analytic element method; it accurately reproduces a boundary condition of uniform head along the cylindrical surface at the perforated face of the well. This model is exercised over a representative range of parameters including the well's length, radius, and pumping rate, and the aquifer's hydraulic conductivity and thickness. Results are presented in a set of figures and tables that compare the well's drawdown to the drawdown that would have been obtained using the Thiem solution with the same pumping rate and radius. A methodology is presented to estimate the head within a horizontal or partially penetrating well by adding a correction term to results that can be readily obtained from computer models of vertical fully penetrating wells. This approach may also be used to contrast the differences in head between horizontal and vertical wells of various lengths, radii, and placement elevations.     

Electromagnetic holographic sensitivity field of two-phase flow in horizontal wells

Electromagnetic holographic data are characterized by two modes, suggesting that image reconstruction requires a dual-mode sensitivity field as well. We analyze an electromagnetic holographic field based on tomography theory and Radon inverse transform to derive the expression of the electromagnetic holographic sensitivity field (EMHSF). Then, we apply the EMHSF calculated by using finite-element methods to flow simulations and holographic imaging. The results suggest that the EMHSF based on the partial derivative of radius of the complex electric potential φ is closely linked to the Radon inverse transform and encompasses the sensitivities of the amplitude and phase data. The flow images obtained with inversion using EMHSF better agree with the actual flow patterns. The EMHSF overcomes the limitations of traditional single-mode sensitivity fields.     

Interpretation of vertical gravity and magnetic profiles

Summary Methods have been developed by quantitatively interpreting gravity and magnetic data in vertical bore-holes or shafts for geological bodies that can be approximated by spheres and cylinders. Master diagrams that can be directly used by the interpreter are provided.     

Stable pumping rates for horizontal wells in bank filtration systems

We consider common bank filtration systems and develop an explicit analytic solution representing steady, two-dimensional, groundwater flow to a horizontal well near a river in an unconfined aquifer. For the boundary-value problem investigated, we find that a unique solution exists for all negative well discharges. For positive discharges, a maximum value exists which corresponds to the formation of a cusp on the free surface. For positive discharges less than the maximum, the solution is not unique, consisting of two alternate configurations of the free surface. One solution includes a stable free surface on a single-valued physical plane, while the alternate solution includes a looped free surface lying on two sheets of a Riemann surface. Imposing a stability condition on the free surface results in a unique solution to the problem. We use the solution to investigate the behavior of the free surface under varying well discharges to identify stable pumping rates and predict well yield. In particular, we examine the well yield and the stability of the free surface when the head in the horizontal well is maintained at the top of the well screen. This condition is shown to produce a stable free surface for a wide range of well radii; the stability is independent of the hydraulic conductivity of the aquifer, the location of the well, or the presence of a skin resistance at the well.     

Permeability profiles in granular aquifers using flowmeters in direct-push wells

Numerical hydrogeological models should ideally be based on the spatial distribution of hydraulic conductivity (K), a property rarely defined on the basis of sufficient data due to the lack of efficient characterization methods. Electromagnetic borehole flowmeter measurements during pumping in uncased wells can effectively provide a continuous vertical distribution of K in consolidated rocks. However, relatively few studies have used the flowmeter in screened wells penetrating unconsolidated aquifers, and tests conducted in gravel-packed wells have shown that flowmeter data may yield misleading results. This paper describes the practical application of flowmeter profiles in direct-push wells to measure K and delineate hydrofacies in heterogeneous unconsolidated aquifers having low-to-moderate K (10(-6) to 10(-4) m/s). The effect of direct-push well installation on K measurements in unconsolidated deposits is first assessed based on the previous work indicating that such installations minimize disturbance to the aquifer fabric. The installation and development of long-screen wells are then used in a case study validating K profiles from flowmeter tests at high-resolution intervals (15 cm) with K profiles derived from multilevel slug tests between packers at identical intervals. For 119 intervals tested in five different wells, the difference in log K values obtained from the two methods is consistently below 10%. Finally, a graphical approach to the interpretation of flowmeter profiles is proposed to delineate intervals corresponding to distinct hydrofacies, thus providing a method whereby both the scale and magnitude of K contrasts in heterogeneous unconsolidated aquifers may be represented.     

大斜度井/水平井随钻方位电磁波测井资料实时反演方法

水平井钻进过程中,地层边界的准确预测对地质导向决策至关重要.但水平井测井环境复杂,随钻方位电磁波测井响应不直观,储层性质、地层界面等信息的精确、定量评价困难.本文提出了一种适用于大斜度井/水平井随钻方位电磁波测井资料的电阻率及地层界面信息的实时提取方法.该方法主要基于以下几个技术：（1）引入滑动开窗策略将高维随钻方位电磁波测井资料的反演问题转化为多个连续窗口的一维反演问题;（2）对反演参数进行多初值初始化,结合正则化Levenberg-Marquardt算法保证解的全局最优性;（3）充分利用随钻方位电磁波测井多种探测深度信息,采用分级反演策略提高反演速度.数值和实测资料反演结果表明：降维反演方法可实现对地层界面的实时、准确提取,且适用于多数复杂地层结构和任意井眼轨迹;利用浅探测地质信号可实现随钻方位电磁波测井资料的分级反演,最终反演速度可提升10倍以上.

     

Study on the simulation of acoustic logging measurements in horizontal and deviated wells

The conventional acoustic logging interpretation method, which is based on vertical wells that penetrate isotropic formations, is not suitable for horizontal and deviated wells penetrating anisotropic formations. This unsuitability is because during horizontal and deviated well drilling, cuttings will splash on the well wall or fall into the borehole bottom and form a thin bed of cuttings. In addition, the high velocity layers at different depths and intrinsic anisotropy may affect acoustic logging measurements. In this study, we examine how these factors affect the acoustic wave slowness measured in horizontal and deviated wells that are surrounded by an anisotropic medium using numerical simulation. We use the staggered-grid finite difference method in time domain (FDTD) combined with hybrid-PML. First, we acquire the acoustic slowness using a simulated array logging system, and then, we analyze how various factors affect acoustic slowness measurements and the differences between the effects of these factors. The factors considered are high-velocity layers, thin beds of cuttings, dipping angle, formation thickness, and anisotropy. The simulation results show that these factors affect acoustic wave slowness measurements differently. We observe that when the wavelength is much smaller than the distance between the borehole wall and high velocity layer, the true slowness of the formation could be acquired. When the wavelengths are of the same order (i.e., in the near-field scenarios), the geometrical acoustics theory is no longer applicable. Furthermore, when a thin bed of cuttings exists at the bottom of the borehole, Fermat's principle is still applicable, and true slowness can be acquired. In anisotropic formations, the measured slowness changes with increments in the dipping angle. Finally, for a measurement system with specific spacing, the slowness of a thin target layer can be acquired when the distance covered by the logging tool is sufficiently long. Based on systematical simulations with different dipping angles and anisotropy in homogenous TI media, slowness estimation charts are established to quantitatively determine the slowness at any dipping angle and for any value of the anisotropic ratio. Synthetic examples with different acoustic logging tools and different elastic parameters demonstrate that the acoustic slowness estimation method can be conveniently applied to horizontal and deviated wells in TI formations with high accuracy.     

利用水平井分段压裂微地震裂缝方位监测确定最大水平主应力方位

<正>水平井是石油钻井技术进步的产物,因其水平段在深度方向起伏很小,可以简化为地下几百米至几千米深度处的一个完全水平的空心管线,管线壁称为套管;在套管壁上射出若干簇孔,称为射孔;射孔使管线内外的流体可以进行交换,深度通常大于1m;射孔指向方向垂直于套管壁及水平段走向,人工压裂裂缝沿着最大水平主应力方向扩展,这是由人工裂缝方向判断最大水平主应力方向     

Interpretation of some gravity profiles from estern Sudan

Summary Gravimetrically determined stations in the province of Kordofan, western Sudan, are interpreted. Boreholes and surface geology are employed, and a Gier Algol print-plot programme has been used. The interpretation shows that the granite mass increases in average thickness from 2000 m in the west to about 20000 m in the east. A dyke-like pattern of gneisses and schists is observed.The paper is forwarded by ProfessorSvend Saxov, Laboratory of Applied geophysics, University of Aarhus, Aarhus, Denmark.     

大斜度井和水平井井地三维电阻率数值模拟和联合反演

井地电阻率成像法利用井套管作电流源向井下供入大功率直流电流,在地表测量由地下介质的电性变化形成的电位分布,通过反演可得到地下介质的电阻率分布.针对大斜度井和水平井开展井地三维电阻率数值模拟和反演研究,对油田注水及压裂效果监测具有重要意义.基于井地电阻率成像法原理,采用有限差分法和不完全切勒斯基共轭梯度法进行了三维正演模拟,研究了大斜度井和水平井的井地电位响应特征.提出了采用层状约束阻尼最小二乘法由浅到深地进行大斜度井和水平井的多层联合反演,并对实际水平井井地电位各个层段数据进行了三维反演.模拟结果表明,倾斜线源和水平线源会对地面电位响应产生明显影响,在反演中需要考虑线源形态.实际水平井井地电位反演成像表明,考虑倾斜线源或者水平线源的联合反演得到了准确的水平井三维注水层成像图,得出注水层的真电阻率分布,能够判断注水运移方向.

     

大斜度井和水平井井地三维电阻率数值模拟和联合反演

井地电阻率成像法利用井套管作电流源向井下供入大功率直流电流,在地表测量由地下介质的电性变化形成的电位分布,通过反演可得到地下介质的电阻率分布.针对大斜度井和水平井开展井地三维电阻率数值模拟和反演研究,对油田注水及压裂效果监测具有重要意义.基于井地电阻率成像法原理,采用有限差分法和不完全切勒斯基共轭梯度法进行了三维正演模拟,研究了大斜度井和水平井的井地电位响应特征.提出了采用层状约束阻尼最小二乘法由浅到深地进行大斜度井和水平井的多层联合反演,并对实际水平井井地电位各个层段数据进行了三维反演.模拟结果表明,倾斜线源和水平线源会对地面电位响应产生明显影响,在反演中需要考虑线源形态.实际水平井井地电位反演成像表明,考虑倾斜线源或者水平线源的联合反演得到了准确的水平井三维注水层成像图,得出注水层的真电阻率分布,能够判断注水运移方向.     

Interpretation of magnetic anomalies by horizontal and vertical derivatives of the analytic signal

Magnetic anomalies are often disturbed by the magnetization direction, so we can’t directly use the original magnetic anomaly to estimate the exact location and geometry of the source. The 2D analytic signal is insensitive to magnetization direction. In this paper, we present an automatic method based on the analytic signal horizontal and vertical derivatives to interpret the magnetic anomaly. We derive a linear equation using the analytic signal properties and we obtain the 2D magnetic body location parameters without giving a priori information. Then we compute the source structural index (expressing the geometry) by the estimated location parameters. The proposed method is demonstrated on synthetic magnetic anomalies with noise. For different models, the proposed technique can both successfully estimate the location parameters and the structural index of the sources and is insensitive to noise. Lastly, we apply it to real magnetic anomalies from China and obtain the distribution of unexploited iron ore. The inversion results are consistent with the parameters of known ore bodies.     

Comparison of longitudinal equilibrium profiles of estuaries in idealized and process-based models

The process-based morphodynamic model Delft3D-MOR and the idealized model of schuttelaars and De Swart (2000) are compared with each other. The differences between the two models in their mathematical-physical formulation as well as the boundary conditions are identified. Their effect on producing cross-sectionally averaged morphological equilibria of tidal inlets with arbitrary length and forced at the seaward boundary by a prescribed M₂ and M₄ sea- surface elevation is studied and an inventory is made of all relevant differences. The physical formulations in the source code of Delft3D-MOR are modified in various steps to resemble the formulations in the idealized model. The effect of each of the differences between the idealized and process-based model are studied by comparing the results of the idealized model to those of the adapted process-based model. The results of the idealized model can be qualitatively reproduced by the process-based model as long as the same morphological boundary condition is applied at the open sea end. This means that the simplifications concerning the mathematical formulation of the physical processes in the idealized model can be justified. Furthermore, it can be inferred that the morphological boundary condition at the open sea end is an essential element in controlling the behaviour of morphodynamic models for tidal inlets and estuaries.Responsible Editor: Jens Kappenberg     

Mathematical models of groundwater flow with a horizontal drain

An analytical review is given to drainage theory, the foundation of which was laid down by researchers in our country in the first half of the 20th century. Later the author of this paper developed and thoroughly studied a complex of mathematical models describing two-dimensional groundwater flow in drained soils based on the boundary problems of the theory of analytical functions. These models were used to provide a hydrodynamic substantiation for the regularities in subsoil water dynamics that were established by ameliorators in relation to the leaching of salts from soils; to identify the specific features of groundwater flow with drainage in the presence in the flow of one or two unknown free boundaries; and to develop approaches to the analysis of such flows in a direct formulation.     

1683-1949年中国地磁测量、地磁图和地磁模型的总考评与研究

本文全面系统回顾考察了新中国成立前中国的地磁测量、地磁图和地磁模型研究概况:主要介绍中国人独立进行的地磁测量;评述了1932年山东半岛地磁图,1915.0和1936.0年中国(部分地区)地磁图;1945.0年东亚地磁图;1946.1年四川北碚地区地磁图以及1946.1年四川北碚地区地磁场泰勒多项式模型.并根据1936年地磁数据计算研究了1936.0年中国地磁场的球冠谐模型和曲面Spline模型,以及1900.0-1940.0年5个年代的中国地磁数据集、地磁模型和地磁图.对今后我国地磁测量研究工作给出了几点建议.     

低流量水平井油水两相分层流速度场数值模拟与实验研究

低流量水平井中,由于存在重力分异,油水两相的流型主要以分层流为主,流体速度场的分布沿重力方向发生变化,致使油田现场生产测井的流量和持水率测量及多相流分相流量解释非常困难.为了充分了解水平井流体流动特征提高生产测井解释精度,本文以低流量水平井油水两相分层流为研究对象,采用流体动态测量实验和数值模拟计算相结合的方法,研究其...     

Interpretation of self potential data for vertical and nearly vertical sheets of infinite horizontal extension

Summary Quantitative interpretation of SP anomalies due to vertical and nearly vertical, thin sheets of infinite horizontal extension has been developed in this paper. The methods proposed utilise selective SP values at some selected points on the profile thus enabling us to get the different parameters of the causative body directly. Also, for the case of vertical sheet, an interpretation technique has been formulated when the base station is within the zone of influence of the potential field of the causative body.     

Interpretation of ground VLF-EM data in terms of inclined sheet-like conductor models

The response of two-dimensional, inclined, sheet-like conductors with, low conductance values to plane wave electromagnetic fields in the very low frequency (VLF) range has been evaluated by using a numerical technique. The conductance values of the conductors considered are appropriate for those produced by water and/or clay-filled fracture and shear zones in the Precambrian crystalline rocks of the Canadian Shield. The surrounding host rock was assumed to be, resistive with resistivities in the 1–10 k.m range to reflect the high resistivities over the shield areas. No overburden was assumed in this analysis.The results of the computations are presented in the form of characteristic interpretation diagrams to interpret ground VLF data in the field, where facilities for direct numerical modelling may not be available. A method for interpreting ground VLF data using such characteristic diagrams has been proposed in this paper which requires a prior knowledge of the host rock resistivity and the inclination of the conductor. These two parameters may be derived from a VLF resistivity survey and from appropriate filtering of the VLF tilt angle response. The interpretation method was applied to a ground VLF anomaly obtained at a research site near Atikokan in NW Ontario, which yielded an interpretation compatible with information from geological mapping.Geological Survey of Canada Contribution No. 51888.     

Temperature profiles in the earth of importance to deep electrical conductivity models

Deep in the Earth, the electrical conductivity of geological material is extremely dependent on temperature. The knowledge of temperature is thus essential for any interpretation of magnetotelluric data in projecting lithospheric structural models. The measured values of the terrestrial heat flow, radiogenic heat production and thermal conductivity of rocks allow the extrapolation of surface observations to a greater depth and the calculation of the temperature field within the lithosphere. Various methods of deep temperature calculations are presented and discussed. Characteristic geotherms are proposed for major tectonic provinces of Europe and it is shown that the existing temperatures on the crust-upper mantle boundary may vary in a broad interval of 350–1,000°C. The present work is completed with a survey of the temperature dependence of electrical conductivity for selected crustal and upper mantle rocks within the interval 200–1,000°C. It is shown how the knowledge of the temperature field can be used in the evaluation of the deep electrical conductivity pattern by converting the conductivity-versustemperature data into the conductivity-versus-depth data.