Hongliang Wu

《地球科学》(ISSN 1000-2383;CN42-1874/P) (原刊名《地球科学——中国地质大学学报》ISSN 1000-2383;CN42-1233/P)，创刊于1957年，是由教育部主管、中国地质大学主办的综合性地学学术理论期刊，月刊，2015年改名为《地球科学》,她以反映中国地质大学教学与科研最新的、高水平的基础地质、应用地质、资源与环境地质及地学工程技术研究成果为主要任务。

Unconventional reservoirs, such as shale and tight formations, have become increasingly vital contributors to oil and gas production. In these reservoirs, fractures serve as crucial spaces for fluid migration and storage, making their precise assessment essential. Array acoustic logging stands out as a pivotal method for evaluating fractures. To investigate the impact of fracture width, fracture-filling conditions, and acoustic frequency on compressional and shear waves, a three-dimensional variable mesh finite difference program was employed for acoustic logging numerical simulation. Firstly, numerical models representing fractured formations with varying fracture widths and distinct fluid-filling conditions were established, and array acoustic logging numerical simulations were conducted at different frequencies. Subsequently, the waveform data were processed to extract acoustic characteristic parameters, such as velocities and amplitude attenuations of compressional and shear waves. Finally, a quantitative analysis was conducted to examine the variation patterns of characteristic parameters of refracted compressional and shear waves in relation to fracture properties. The research results indicate that amplitude attenuation information derived from borehole wave modes is particularly sensitive to the changes in fracture properties. As fracture width increased, we observed a significant amplitude attenuation in both compressional and shear waves, proportional to the logarithm of the attenuation coefficients. Furthermore, when the fracture width was constant, gas-filled fractures exhibited more prominent amplitude attenuation than water-filled fractures, with shear wave attenuation being more sensitive to the filling material. Moreover, from a quantitative perspective, the analysis revealed that the attenuation coefficients of refracted compressional and shear waves exhibited an exponential variation with gas saturation. Notably, once fracture width and filling conditions were established, the amplitudes of compressional and shear waves at the dominant frequency of 40 kHz were significantly reduced compared to those at 8 kHz, accompanied by increased attenuation. Subsequent quantitative analysis revealed that, when the product of fracture width and dominant frequency remains constant, the corresponding attenuation coefficient ratios approach 1. This indicates that the attenuation process of acoustic propagation in fractured media follows the principle of acoustic similarity. The findings of this study provide reference for further research on fracture property evaluation methods based on array acoustic logging data.

Oil saturation index (OSI) serves as an important indicator for potential movable hydrocarbons evaluation of shale oil reservoirs, which is mainly obtained through rock pyrolysis experiments. A new method is proposed to evaluate the OSI of shale quantitatively by NMR logging. The OSI value can be accurately obtained through the experimental measurement of organic carbon content (TOC) and rock pyrolysis of shale samples, which can identify the development of mobile hydrocarbons. Subsequently, the mobile fluid porosity can be obtained based on NMR logging. In order to establish the relationship between OSI value and mobile fluid porosity, it is important to determine the T₂ cutoff value corresponding to the mobile fluid porosity. Take shale samples from the first member of the Qingshankou Formation (“Qing 1 Member”) as an example, based on 2D NMR experimental analysis in three different states (original, dried state at 105°C, saturated kerosene), the NMR T₂ cutoff value of movable fluid porosity in the shale of Qing 1 Member is clarified as 8 ms. Integrating rock pyrolysis and 2D NMR experiments, it suggests that the NMR bin porosity with T₂>8 ms has a good linear relationship with the OSI value obtained by pyrolysis analysis. The NMR bin porosity with T₂> 8 ms reflects the OSI value of shale effectively. The larger the NMR bin porosity with T₂>8 ms, the higher the mobile oil content of shale reservoir, which is consistent with the understanding of oil-bearing large pores in the Gulong Shale. The NMR bin porosity can continuously evaluate the vertical variation of the mobile hydrocarbon content. Compared with the traditional experimental measurement of finite depth points, this method has significant advantages, and can avoid the possibility of missing potentially movable oil layers.

Horizontal wells play an important role in expanding the drilling volume of reservoirs and oil production area, and are widely used in unconventional reservoirs. Fractures have a positive effect on reservoir permeability, but fractures can also cause accidents such as casing deformation and inter-well frac-hit. It is of great significance to identify and evaluate fractures intersected with horizontal wells in tight reservoirs. In this paper, a three-dimensional numerical model of horizontal wells and fractures in tight reservoirs is designed. The responses of monopole P-wave and dipole S-wave to fractures with different width, dip angle and filling medium are systematically studied, by using three-dimensional finite difference algorithm. The results show that when the fracture is filled with calcite, the amplitude attenuation of monopole P-wave and dipole S-wave has a monotonic exponential increase with the increase of fracture width and the decrease of fracture dip angle. In the real data processing, the amplitude attenuation of P- and S-waves can be used to jointly evaluate the fracture filled with calcite. When the fracture is filled with water, both P- and S-waves have prominent amplitude attenuation. P wave amplitude attenuation does not have a monotonic variation with the increase of fracture width but it has a monotonic increase with the decrease of fracture dip angle. S wave amplitude attenuation has a monotonic increase with the increase of fracture width and the decrease of fracture dip angle. The amplitude attenuation of P- and S- waves rises significantly when the fracture is filled with natural gas. This study is crucial for better understanding the response of P- and S-waves to fractures intersected with borehole in tight reservoirs, and it provides useful information for the inversion of fracture parameters by using P- and S-waves.

Unconventional oil and gas reservoirs have broad exploration and development prospects. Fracture parameters and effectiveness evaluation are two of the key tasks for the evaluation of these types of reservoirs. Array acoustic logging can be used for fracture evaluation to compensate for the deficiencies of the image logging fracture evaluation method. Therefore, to develop acoustic logging evaluation methods as well as nondestructive testing methods for fractures, experiments were conducted to study the shear wave transmission in fractured media. Experiment data demonstrate a good correlation between the shear wave attenuation coefficient and fracture width, and the shear wave attenuation coefficients rise logarithmically with the increase in the fracture width for all models with different porosities and distinct dip angles of fractures. The shear wave attenuation coefficient changes relatively faster with the fracture width when the fracture width is within 250 μm. In addition, the shear wave attenuation is affected by the core porosity and fracture dip angle. When the fracture width is constant, the shear wave attenuation caused by the 0° fracture is relatively larger and is obviously greater than that of the fractures at other angles, which is consistent with the existing experimental results. The results of this study can be used to guide further research on amplitude compensation methods for sonic signal transmission in fractured media and fracture evaluation methods.

To assess the productivity of oil- and gas-bearing carbonate reservoirs, acoustic waves are used to extract information on the reservoir properties. We performed shock-induced broadband linear acoustic wave experiments on fractured carbonate samples with different fracture apertures. Samples were prepared from a full-diameter carbonate rock core in which the aperture of a single fracture could be varied by means of spacers. Microseismograms were constructed from multiple wave experiments. We measured the wave speeds of different wave modes in the annulus along the cylindrical core sample and the transmission coefficient of the Stoneley wave over the fracture. We found that the transmission coefficient of the annulus Stoneley wave can directly be related to the fracture aperture and the fracture length. Comparison with a straightforward quasi-1D theory yielded best-fit values for fracture aperture and length that were in agreement with the actual values, for the fracture aperture determined by the size of the spacers. The values of the fracture length were not accurately predicted probably due to the underestimation of the fracture compliance in the experiment. The inversion method can be optimized for carbonate reservoir conditions to extract fracture data from Stoneley wave transmission and reflection coefficients.