The production of natural gas hydrates will change the cementation strength, porosity, and effective stress in the stratum, which may lead to engineering and geological disasters. Sand production is a phenomenon where sand particles are carried out of the reservoir along with fluids during gas extraction, posing challenges to safe and sustainable production. This study explored the mechanism of fine particle migration in multiphase flow by a microscopic visualization test device. The device can inject a gas–liquid–solid phase at the same time and allow real-time observation. Experimental tests on fine particle migration of single- and two-phase fluid flow were carried out considering different conditions, i.e., fine particle concentration, fine particle size, fluid flow rate, and gas–liquid ratio. The results show that in single-phase fluid flow, the original gas will gradually dissolve in the liquid phase, and finally stay in the test device as bubbles, which can change the pore structures, resulting in the accumulation of fine particles at the gas–liquid interface. In two-phase fluid flow with mixed gas–water fluids, there are two flow modes of gas–liquid flow: mixed flow and separated flow. The interfacial tension at the gas–liquid interface can effectively migrate fine particles when the gas–liquid flows alternately and the sand production rate further increases as the gas–liquid ratio increases. In addition, changes in the concentration of fine particles, particle size, fluid flow rate, and the gas–liquid ratio will affect the migration of fine particles, leading to differences in the final sand production.

The degradation of soil bonding, which can be described by the evolution of bond degradation variables, is essential in the constitutive modeling of cemented soils. A degradation variable with a value of 0/1.0 indicates that the applied stress is completely sustained by bonded particles/unbounded grains. The discrete element method (DEM) was used for cemented soils to analyze the bond degradation evolution and to evaluate the degradation variables at the contact scale. Numerical cemented soil samples with different bonding strengths were first prepared using an advanced contact model (CM). Constant stress ratio compression, one-dimensional compression, conventional triaxial tests (CTTs), and true triaxial tests (TTTs) were then implemented for the numerical samples. After that, the numerical results were adopted to investigate the evolution of the bond degradation variables B_N and B₀. In the triaxial tests, B₀ evolves to be near to or larger than B_N due to shearing, which indicates that shearing increases the bearing rate of bond contacts. Finally, an approximate stress-path-independent bond degradation variable B_σ was developed. The evolution of B_σ with the equivalent plastic strain can be effectively described by an exponential function and a hyperbolic function.

The recognition of the geomechanical properties of methane hydrate-bearing soil (MHBS) is crucial to exploring energy resources. The paper presents the mechanical properties of a pore-filled MHBS at a critical state using the distinct element method (DEM). The pore-filled MHBS was simulated as cemented MH agglomerates to fill the soil pores at varying levels of methane hydration (MH) saturation. A group of triaxial compression (TC) tests were conducted, subjecting MHBS samples to varying effective confining pressures (ECPs). The mechanical behaviors of a pore-filled MHBS were analyzed, as it experienced significant strains leading to a critical state. The findings reveal that the proposed DEM successfully captures the qualitative geomechanical properties of MHBS. As MH saturation increases, the shear strength of MHBS generally rises. Moreover, higher ECPs result in increased shear strength and volumetric contraction. The peak shear strength of MHBS increases with rising MH saturation, while the residual deviator stress remains mainly unchanged at a critical state. There is a good correlation between fabric changes of the MHBS with variations in principal stresses and principal strains. With increasing axial strain, the coordination number (CN) and mechanical coordination number (MCN) increase to peak values as the values of MH saturation and ECPs increase, and reach a stable value at a larger axial strain.

The tunnels may traverse soft soil strata rich in hazardous gases such as methane, which poses several adverse effects on construction. These effects primarily include excessive deformation of the tunnel, leakage, or migration of hazardous gases causing soil deformation and settlement, thereby negatively impacting shield tunneling operations; and the ingress of hazardous gases into the tunnel, which could lead to fires or explosions upon encountering an ignition source, as well as endangering the health of construction personnel within the tunnel. Therefore, this paper focuses on the critical technical issues associated with the construction of the B Line tunnel in the presence of hazardous gases. Utilizing numerical simulations, the uneven settlement and reinforcement techniques are analyzed for shield tunneling through hazardous gas-bearing strata. Based on a comprehensive geological investigation report, a robust finite element model of the tunnel has been developed. The implementation of reinforcement measures and strata stabilization techniques has proven effective under identical leakage conditions. The maximum displacement of the tunnel strata was reduced from the initial 2.0–1.5 mm, resulting in a 25% decrease in overall strata displacement. The findings will provide scientific basis and technical support for similar construction projects.