一种基于核磁共振测井计算低孔低渗气层孔隙度的新方法

川中地区须家河组发育大量的低孔低渗气藏,岩心分析结果表明,储层孔隙度为4%～8%,渗透率为0.1～1mD。对于这类低孔低渗气层,准确地计算储层孔隙度显得尤为重要。考虑到储层含气的影响,利用单一测井资料难以准确地计算储层孔隙度,本文从核磁共振测井基本原理和基于声波时差测井的岩石体积物理模型分析出发,提出了一种新的结合声波时差—核磁共振测井资料计算低孔低渗气层真实孔隙度的方法。实际资料应用表明,利用该方法计算的孔隙度与岩心分析结果吻合较好,计算结果可以真实地反映实际地层孔隙度。

Acoustic reflection well logging modeling using the frequency-domain finite-element method with a hybrid PML

In this paper, we propose a hybrid PML （H-PML） combining the normal absorption factor of convolutional PML （C-PML） with tangential absorption factor of Mutiaxial PML （M-PML）. The H-PML boundary conditions can better suppress the numerical instability in some extreme models, and the computational speed of finite-element method and the dynamic range are greatly increased using this HPML. We use the finite-element method with a hybrid PML to model the acoustic reflection of the interface when wireline and well logging while drilling （LWD）, in a formation with a reflector outside the borehole. The simulation results suggests that the PS- and SP- reflected waves arrive at the same time when the inclination between the well and the outer interface is zero, and the difference in arrival times increases with increasing dip angle. When there are fractures outside the well, the reflection signal is clearer in the subsequent reflection waves and may be used to identify the fractured zone. The difference between the dominant wavelength and the model scale shows that LWD reflection logging data are of higher resolution and quality than wireline acoustic reflection logging.

Ultrasonic attenuation estimation based on time-frequency analysis

The quality factor(or Q value)is an important parameter for characterizing the inelastic properties of rock.Achieving a Q value estimation with high accuracy and stability is still challenging.In this study,a new method for estimating ultrasonic attenuation using a spectral ratio based on an S transform(SR-ST)is presented to improve the stability and accuracy of Q estimation.The variable window of ST is used to solve the time window problem.We add two window factors to the Gaussian window function in the ST.The window factors can adjust the scale of the Gaussian window function to the ultrasonic signal,which reduces the calculation error attributed to the conventional Gaussian window function.Meanwhile,the frequency bandwidth selection rules for the linear regression of the amplitude ratio are given to further improve stability and accuracy.First,the feasibility and influencing factors of the SR-ST method are studied through numerical testing and standard sample experiments.Second,artificial samples with different Q values are used to study the adaptability and stability of the SR-ST method.Finally,a further comparison between the new method and the conventional spectral ratio method(SR)is conducted using rock field samples,again addressing stability and accuracy.The experimental results show that this method will yield an error of approximately 36%using the conventional Gaussian window function.This problem can be solved by adding the time window factors to the Gaussian window function.The frequency bandwidth selection rules and mean slope value of the amplitude ratio used in the SR-ST method can ensure that the maximum error of different Q values estimation(Q>15)is less than 10%.

Seismic physical modeling and quality factor

Accurate Q parameter is hard to be obtained, but there is great difference between Q measurements from different measurement methods in seismic physical modelling. The influence factors, stability and accuracy of different methods are analyzed through standard sample experiment and the seismic physical modelling. Based on this, we proposed an improved method for improving accuracy of pulse transmission method, in which the samples with similar acoustic properties to the test sample are selected as the reference samples. We assess the stability and accuracy of the pulse transmission, pulse transmission insertion, and reflection wave methods for obtaining the quality factor Q using standard and reference samples and seismic physical modeling. The results suggest that the Q-values obtained by the pulse transmission method are strongly affected by diffraction and the error is 50% or greater, whereas the relative error of the improved pulse transmission method is about 10%. By using a theoretical diffraction correction method and the improved measurement method, the differences among the Q-measuring methods can be limited to within 10%.