Numerical simulation of the dual laterolog for carbonate cave reservoirs and response characteristics

洞碳酸盐形成被异质描绘，它使电的木头预言困难。知道怎么使用双 laterolog 精确地识别当前是重要的并且份量上评估洞。使用数字模拟计算电的木头回答能为洞鉴定和评估提供一个理论基础。在这份报纸，基于双 laterolog 原则，我们首先学习球形的洞用有限元素方法建模的不同尺寸(女性)，决定在在包括地与不同充满的材料学习洞模型的木头回答以后充满的抵抗力和洞之间的一种关系，并且最后学习洞上的回答用页岩和页岩填满的双 laterolog 。数字结果为碳酸盐洞水库的鉴定和评估提供一个理论基础。

Applications of texture attribute analysis to seismic interpretation

The first generation coherence algorithm(namely C1 algorithm) is based on the statistical cross-correlation theory, which calculates the coherency of seismic data along both in-line and cross-line. The work, based on texture technique, makes full use of seismic information in different directions and the difference of multi-traces, and proposes a novel methodology named the texture coherence algorithm for seismic reservoir characterization, for short TEC algorithm. Besides, in-line and cross-line directions, it also calculates seismic coherency in 45° and 135° directions deviating from in-line. First, we clearly propose an optimization method and a criterion which structure graylevel co-occurrence matrix parameters in TEC algorithm. Furthermore, the matrix to measure the difference between multi-traces is constructed by texture technique, resulting in horizontal constraints of texture coherence attribute. Compared with the C1 algorithm, the TEC algorithm based on graylevel matrix is of the feature that is multi-direction information fusion and keeps the simplicity and high speed, even it is of multi-trace horizontal constraint, leading to significantly improved resolution. The practical application of the TEC algorithm shows that the TEC attribute is superior to both the C1 attribute and amplitude attribute in identifying faults and channels, and it is as successful as the third generation coherence.

Applications of fluid substitution effect analysis on seismic interpretation

The root mean square(RMS) difference of time-lapse seismic amplitudes is routinely used to identify the substituted fluid type in a reservoir during oil field production and recovery. By a time-lapse seismic method, we study the effects of fluid substitution in a physical model, which is an analogy of the three-dimensional inhomogeneous reservoir. For a weak inhomogeneous medium, gas/oil substitution results in positive anomalies in the reservoir layers, and negative anomalies below the bottom of the reservoir layers; while water/oil substitution causes only weak variations in the reservoir layers, but positive anomalies below the bottom of the reservoir layers. For the strong inhomogeneous medium, no matter what kind of fluid substitution(gas/oil or water/oil), there are significant anomalies in seismic amplitude difference attributes both in and below the reservoir layers. Therefore, for weak inhomogeneous media, such as tight sandstone or thin interbedded layers, the RMS amplitude difference attributes can be used to monitor fluid changes and predict the drilling direction; for inhomogeneous medium with karst carves or fractures, it is difficult to accurately determine the distribution of fluids with the RMS amplitude difference attributes.

A seismic texture coherence algorithm and its application

The first generation coherence algorithm （the C1 algorithm） that calculates the coherence of seismic data in-line and cross-line was developed using statistical cross-correlation theory, and it has the limitation that the technique can only be applied to horizons. Based on the texture technique, the texture coherence algorithm uses seismic information in different directions and differences among multiple traces. It can not only calculate seismic coherence in in-line and cross-line directions but also in all other directions. In this study, we suggested first an optimization method and a criterion for constructing the gray level co-occurrence matrix of the seismic texture coherence algorithm. Then the co-occurrence matrix was prepared to evaluate differences among multiple traces. Compared with the C1 algorithm, the seismic texture coherence algorithm suggested in this paper is better than the C1 in its information extraction and application. Furthermore, it implements the multi-direction information fusion and it, also has the advantage of simplicity and effectiveness, and improves the resolution of the seismic profile. Application of the method to field data shows that the texture coherence attribute is superior to that of C 1 and that it has merits in identification of faults and channels.

Deformation geometry and timing of the Wupoer thrust belt in the NE Pamir and its tectonic implications

The Pamir region, located to the northwest of the Tibetan Plateau, provides important information that can aid the understanding of the plateau＇s tectonic evolution. Here we present new findings on the deforma- tion geometry and timing of the Wupoer thrust belt at the northeastem margin of Pamir. Field investigations and interpretations of seismic profiles indicate that the eastern portion of the Wupoer thrust belt is dominated by an underlying foreland basin and an overlying piggy-back basin. A regional unconformity occurs between the Pliocene （N2） and the underlying Miocene （NI） or Paleogene （Pg） strata associated with two other local unconformities between Lower Pleistocene （Q1） and N2 and between Middle Pleistocene （Q2-4） and Q1 strata. Results of structural restorations suggest that compres- sional deformation was initiated during the latest Miocene to earliest Pliocene, contributing a total shortening magnitude of 48.6 km with a total shortening rate of 48.12%, most of which occurred in the period from the latest Miocene to earliest Pliocene. These results, com- bined with previous studies on the Kongur and Tarshkor- gan extensional system, suggest an interesting picture of strong piedmont compressional thrusting activity concur- rent with interorogen extensional rifting. Combining these results with previously published work on the lithospheric architecture of the Pamir, we propose that gravitational collapse drove the formation of simultaneous extensional and compressional structures with a weak, ductile middle crustal layer acting as a decollement along which both the extensional and compressional faults merged.

Flow pattern identification in oil wells by electromagnetic image logging

Petroleum production logging needs to determine the interpretation models first and flow pattern identification is the foundation, but traditional flow pattern identification methods have some limitations. In this paper, a new method of flow pattern identification in oil wells by electromagnetic image logging is proposed. First, the characteristics of gas-water and oil-water flow patterns in horizontal and vertical wellbores are picked up. Then, the continuous variation of the two phase flow pattern in the vertical and horizontal pipe space is discretized into continuous fluid distribution models in the pipeline section. Second, the electromagnetic flow image measurement responses of all the eight fluid distribution models are simulated and the characteristic vector of each response is analyzed in order to distinguish the fluid distribution models. Third, the time domain changes of the fluid distribution models in the pipeline section are used to identify the flow pattern. Finally, flow simulation experiments using electromagnetic flow image logging are operated and the experimental and simulated data are compared. The results show that the method can be used for flow pattern identification of actual electromagnetic image logging data.

基于矢量化反射率法的多波叠前联合AVA反演（英文）

目前,大多数多波叠前联合AVA(AVO)反演方法都是基于精确Zoeppritz方程及其近似公式。然而,这些方程仅仅反映了反射振幅与入射角和界面两侧弹性参数之间的关系,没有考虑几何扩散、吸收衰减、透射损失、多次波和纵横波不匹配等传播效应,导致不能精确地描述地震波真实传播特征。传统的AVA反演方法通常需要在反演之前消除这些传播效应,而实际中很难满足此要求。在一维地质模型的假设条件下,反射率法可以模拟地震波全波场响应。基于此,本文提出了一种基于矢量化反射率法的非线性多波叠前联合AVA反演方法,其中引入了快速非支配排序遗传算法优化处理多目标函数从而获得纵、横波速度和密度等多个参数。因为该方法能在未引入权重系数的情况下同时优化多个目标函数,所以具有较强的鲁棒性。模型测试验证了该叠前联合AVA反演方法的有效性。结果显示:当透射损失和层间多次波等传播效应没有被很好校正时,基于反射率法的反演结果优于基于精确Zoeppritz方程的反演结果。

Q full-waveform inversion based on the viscoacoustic equation

Presently, most full-waveform inversion methods are developed for elastic media and ignore the effect of attenuation. The calculation of the quality factor Q is based on velocity parameter inversion under the assumption of a given Q-model that is obtained by tomographic inversion. However, the resolution of the latter is low and cannot reflect the amplitude attenuation and phase distortion during wave propagation in viscoelastic media. Thus, a Q waveform inversion method is proposed. First, we use standard linear body theory to describe attenuation and then we derive the simplified viscoacoustic equation that characterizes amplitude attenuation and phase distortion. In comparison with conventional equations, the simplifi ed equation involves no memory variables and therefore requires less memory during computation. Moreover, the implementations of the attenuation compensation are easier. The adjoint equation and the corresponding gradient equation with respect to either L2-norm or the zero-lag cross-correlation objective function are then derived and the regularization strategy for overcoming the instability during numerical solution of the adjoint equation is proposed. The Q waveform inversion is developed using the limited-memory Broyden–Fletcher– Goldfarb–Shanno (L-BFGS) iteration method for known velocity. To alleviate the dependence of the waveform inversion on the initial model and overcome cycle skipping to some extent, we adopt multiscale analysis. Furthermore, anti-noise property and double-parameter inversion are assessed based on the results of numerical modeling.

Reservoir and development characteristics of the Da'anzhai tight oil in Sichuan Basin,SW China

The oil in the Jurassic Da'anzhai reservoirs in the Sichuan Basin is unconventional tight oil,which accumulated in or near source rocks,and did not experience extensive migration in a large-scale long distance.The first submember,second submember and third submember of Da'anzhai Member are dominated by shell limestone which is widely and continuously distributed,and are typical near-source lacustrine shell limestone tight reservoirs.Complex lithology,multiple types of reservoir space and complicated pore structure are developed in these reservoirs.The effective reservoir space mainly includes micro-pores and micro-fractures with strong fabric selectivity.The petrophysics experiment reveals that the average connected matrix porosity of tight oil reservoir in Jurassic Da'anzhai Member is about 2.13%,lower than that of other tight oil reservoirs but higher than the average effective porosity(0.97%)from previous single alcohol-saturated method.According to production performance data,the Da'anzhai shell limestone reservoir is not a simple porous or fractured reservoir,but has complex porethroat-fracture association or storage-seepage mode.Because the development of fossil shells controls the development of micro-fractures,fluids are difficult to enter into but easy to escape from the reservoirs.Although the pore-throat is fine,the sorting is poor and the displacement pressure is high,the movable fluid saturation and mercury ejection efficiency of the reservoir in the Da'anzhai Member is only slightly lower than that of some storage-seepage modes,and higher than that of Oil-bearing Group 7 of Yanchang Formation in the Ordos Basin.The reservoir in the Da'anzhai Member is one of the few tight oil reservoirs with high natural productivity.The tight oil in the shell limestone of the Da'anzhai Member has great development potential,but its extensive and effective development also has some challenges,such as high seepage resistance of matrix and ineffective single development mode.The development mode of the Da'anzhai tight oil should draw lessons from the Bakken Formation in North America and Oil-bearing Gourp 7 of Yanchang Formation in the Ordos Basin,and thus,effective development technologies based on volume fracturing and fine operation for shell limestone tight oil in the Da'anzhai Member in Sichuan Basin are developed to realize the development of profit and scale.

Seismic energy dispersion compensation by multi-scale morphology

Seismic energy decays while propagating subsurface, which may reduce the resolution of seismic data. This paper studies the method of seismic energy dispersion compensation which provides the basic principles for multi-scale morphology and the spectrum simulation method. These methods are applied in seismic energy compensation. First of all, the seismic data is decomposed into multiple scales and the effective frequency bandwidth is selectively broadened for some scales by using a spectrum simulation method. In this process, according to the amplitude spectrum of each scale, the best simulation range is selected to simulate the middle and low frequency components to ensure the authenticity of the simulation curve which is calculated by the median method, and the high frequency component is broadened. Finally, these scales are reconstructed with reasonable coefficients, and the compensated seismic data can be obtained. Examples are shown to illustrate the feasibility of the energy compensation method.

Numerical Simulation of a Multi-Frequency Resistivity Logging-While-Drilling Tool Using a Highly Accurate and Adaptive Higher-Order Finite Element Method

A novel,highly efficient and accurate adaptive higher-order finite element method(hp-FEM)is used to simulate a multi-frequency resistivity loggingwhile-drilling(LWD)tool response in a borehole environment.Presented in this study are the vector expression of Maxwell’s equations,three kinds of boundary conditions,stability weak formulation of Maxwell’s equations,and automatic hpadaptivity strategy.The new hp-FEM can select optimal refinement and calculation strategies based on the practical formation model and error estimation.Numerical experiments show that the new hp-FEM has an exponential convergence rate in terms of relative error in a user-prescribed quantity of interest against the degrees of freedom,which provides more accurate results than those obtained using the adaptive h-FEM.The numerical results illustrate the high efficiency and accuracy of the method at a given LWD tool structure and parameters in different physical models,which further confirm the accuracy of the results using the Hermes library(http://hpfem.org/hermes)with a multi-frequency resistivity LWD tool response in a borehole environment.

Quantitative description of steam channels after steam flooding

Steam channeling is one of the main barriers for EOR after steam flooding.In order to enhance the oil recovery in steam flooded reservoirs,steam channel volumes should be precisely known.In this paper,a set of methods has been established in order to study steam channeling quantitatively by using dynamic data.Firstly,steam channeling wells are identified through curves of watercut and temperature.Then,considering the hysteresis phenomenon,channeling relations are identified with the correlation coefficients between injection wells and production wells under different conditions.Lastly,an analytic model,in which steam condensation,pressure and temperature are considered,is established to calculate the steam channel volumes.A production well named L31615 in some block in Henan Oilfield in China is systematically analyzed by using the method established.The whole block is further analyzed,and the distribution map of steam channels and the steam channel volumes are obtained.The results show that steam channeling does not only occur inside a well group,but also sometimes occurs between wells in different well groups.The calculation of the steam channel volumes provides a theoretical basis for bringing a remedial action like plugging into operation.