The capillary pressure curve is a crucial basis for studying the pore structure and multiphase flow characteristics in oil and gas reservoirs. Due to the loose and unconsolidated nature of the clayey-silt sediment of natural gas hydrate reservoirs in the South China Sea, conventional methods such as mercury intrusion and centrifugation struggle to obtain capillary pressure curves for these sediments. In this study, X-ray diffraction analysis, scanning electron microscopy, nitrogen adsorption, and water-gas contact angle measurements are utilized to characterize the mineral composition, pore structure, pore size distribution, and wettability of the clayey-silt sediment. Subsequently, the filter paper method from soil mechanics is employed to determine the capillary pressure curve for the clayey-silt samples. The results indicate that the capillary pressure curve obtained through the filter paper method exhibits a saturation range of 18.39%-80.31% and a capillary pressure range of 19.04 to 46,481.42 kPa. It exhibits a distinct two-stage characteristic, where capillary pressure changes rapidly with water saturation below 61.05% and slowly above 61.05%. The pore radius calculated from the capillary pressure curve ranges from 2.41 nm to 5.91 μm. This alignment with the pore ranges obtains from nitrogen adsorption and Scanning Electron Microscopy confirms the accuracy of the obtained capillary pressure curve. Furthermore, in comparison with a literature capillary pressure curve obtained through centrifugation, the paper filtration method covers a broader range, providing better representation of capillary pressure in the multiscale pores of clayey-silt samples.

Natural gas hydrate reservoirs in the northern South China Sea primarily comprise clayey silt, making exploitation more challenging relative to sandy reservoirs in other countries and regions. This paper provides an overview of the latest research developments in the exploitation mechanism covering the past five years, focusing on hydrate phase transition, multiphase flow in the decomposition zone, the seepage regulation of reservoir stimulation zone, and production capacity simulation, all of which are relevant to the previously conducted two rounds of hydrate trial production in offshore areas of China. The results indicate that the phase transition of clayey-silt hydrate remains in a dynamic equilibrium, with the decomposition efficiency mainly controlled by the coupling of heat and flow and high heat consumption during decomposition. The decomposition zone exhibits strong hydrophilicity, easy adsorption, and sudden permeability changes. A temperature drop is present that is concentrated near the wellbore, and once a water lock has formed, the gas-phase flow capacity significantly decreases, leading to potential secondary hydrate formation. To enhance permeability and increase production, it is imperative to implement reservoir and temperature field reconstruction based on initial formation alterations, which will further optimize and improve the transport capacity of the reservoir.

Simulation of elastic wave propagation is an important method for oil and gas exploration. Accuracy and efficiency of elastic wave simulation in complex geological environment are always the focus issue. In order to improve the accuracy and efficiency in numerical modeling, a staggered grid fourth-order finite-difference scheme of modeling elastic wave in frequency-domain is developed, which can provide stable numerical solution with fewer number of grid points per wavelength. The method is implemented on first-order velocity-stress equation and a parsimonious spatial staggered-grid with fourth-order approximation of the first-order derivative operator. Numerical tests show that the accuracy of the fourth-order staggered-grid stencil is superior to that of the mixed-grid and other conventional finite difference stencils, especially in terms of shear-wave phase velocity. Measures of mass averaging acceleration and optimization of finite difference coefficients are taken to improve the accuracy of numerical results. Meanwhile, the numerical accuracy of the finite difference scheme can be further improved by enlarging the mass averaging area at the price of expanding the bandwidth of the impedance matrix that results in the reduction of the number of grid points to 3 per shear wavelength and computer storage requirement in simulation of practical models. In our scheme, the phase velocities of compressional and shear wave are insensitive to Poisson's ratio does not occur conventional finite difference scheme in most cases, and also the elastic wave modeling can degenerate to acoustic case automatically when the medium is pure fluid or gas. Furthermore, the staggered grid scheme developed in this study is suitable for wave propagation modeling in media with coupling fluid-solid interfaces that are not resolved for previous finite difference method.