Gas Layer Research Articles

SOFIA FEEDBACK Survey: The Pillars of Creation in [C ii] and Molecular Lines

We investigate the physical structure and conditions of photodissociation regions (PDRs) and molecular gas within the Pillars of Creation in the Eagle Nebula using SOFIA FEEDBACK observations of the [C ii] 158 μm line. These observations are velocity resolved to 0.5 km s−1 and are analyzed alongside a collection of complimentary data with similar spatial and spectral resolution: the [O i] 63 μm line, also observed with SOFIA, and rotational lines of CO, HCN, HCO+, CS, and N2H+. Using the superb spectral resolution of SOFIA, APEX, CARMA, and BIMA, we reveal the relationships between the warm PDR and cool molecular gas layers in context of the Pillars’ kinematic structure. We assemble a geometric picture of the Pillars and their surroundings informed by illumination patterns and kinematic relationships and derive physical conditions in the PDRs associated with the Pillars. We estimate an average molecular gas density cm−3 and an average atomic gas density n H ∼ 1.8 × 104 cm−3 and infer that the ionized, atomic, and molecular phases are in pressure equilibrium if the atomic gas is magnetically supported. We find pillar masses of 103, 78, 103, and 18 M ⊙ for P1a, P1b, P2, and P3, respectively, and evaporation times of ∼1–2 Myr. The dense clumps at the tops of the pillars are currently supported by the magnetic field. Our analysis suggests that ambipolar diffusion is rapid and these clumps are likely to collapse within their photoevaporation timescales.

Numerical investigation of production characteristics and interlayer interference during co-production of natural gas hydrate and shallow gas reservoir

The coexistent system of natural gas hydrate and shallow gas (NGH-SG) enhances stability compared to single resource reservoirs, rendering it a promising target for co-production and commercial exploitation. In this study, we established a reservoir model of the NGH-SG system based on actual geological conditions. Subsequently, a long-term production simulation was performed to contrast the production behaviors of three exploitation methods: co-production of NGH-SG (CP), single production of hydrate-bearing layer (SP-HBL), and single production of shallow gas layer (SP-SGL). The study revealed that the presence of free gas in SGL can significantly increase the overall gas production capacity, with SGL contributing over 85% of the total gas production during the entire co-production period. On the 5000th day, the cumulative gas production of CP exceeded that of SP-SGL by 1.07 times and surpassed SP-HBL by 9.60 times. Furthermore, the study extensively investigated the impact of interlayer permeability on production behaviors. Lower interlayer permeability was found to be more conducive to releasing the natural gas production capacity from the reservoir. Notably, the phenomenon of interlayer interference should be carefully considered during CP. The influence degree of interlayer interference during CP becomes increasingly pronounced with increasing interlayer permeability. Fundamental contributors to interlayer interference encompass heat transfer, pressure diffusion, and interlayer cross-flow. As interlayer permeability rises, mass transfer between the layers intensifies, thus promoting pressure diffusion. Interlayer pressure difference and permeability are critical factors that affect interlayer cross-flow. For NGH-SG reservoirs with higher interlayer permeability, the advantages of CP are more obvious during the early production period. The implementation of the findings of this study may improve production efficiency and contribute to the development of effective management strategies for offshore NGH-SG reservoirs.

A compact multi-planet system transiting HIP 29442 (TOI-469) discovered by TESS and ESPRESSO

Context. One of the goals of the Echelle Spectrograph for Rocky Exoplanets and Stable Spectroscopic Observations (ESPRESSO) Guaranteed Time Observations (GTO) consortium is the precise characterisation of a selected sample of planetary systems discovered by TESS. One such target is the K0V star HIP 29442 (TOI-469), already known to host a validated sub-Neptune companion TOI-469.01, which we followed-up with ESPRESSO. Aims. We aim to verify the planetary nature of TOI-469.01 by obtaining precise mass, radius, and ephemeris, and constraining its bulk physical structure and composition. Methods. Following a Bayesian approach, we modelled radial velocity and photometric time series to measure the dynamical mass, radius, and ephemeris, and to characterise the internal structure and composition of TOI-469.01. Results. We confirmed the planetary nature of TOI-469.01 (now renamed HIP 29442 b), and thanks to the ESPRESSO radial velocities we discovered two additional close-in companions. Through an in-depth analysis of the TESS light curve, we could also detect their low signal-to-noise transit signals. We characterised the additional companions, and conclude that HIP 29442 is a compact multi-planet system. The three planets have orbital periods Porb,b = 13.63083 ± 0.00003, Porb,c = 3.53796 ± 0.00003, and Porb,d = 6.42975−0.00010+0.00009 days, and we measured their masses with high precision: mp,b = 9.6 ± 0.8 M⊕, mp,c = 4.5 ± 0.3 M⊕, and mp,d = 5.1 ± 0.4 M⊕. We measured radii and bulk densities of all the planets (the 3σ confidence intervals are shown in parentheses): Rp,b = 3.48−0.08(−0.28)+0.07(+0.19) R⊕ and ρp,b = 1.3 ± 0.2(0.3)g cm−3; Rp,c = 1.58−0.11(−0.34)+0.10(+0.30) R⊕ and ρp,c = 6.3−1.3(−2.7)+1.7(+6.0)g cm−3; Rp,d = 1.37 ± 0.11(−0.43)(+0.32) R⊕ and ρp,d = 11.0−2.4(−6.3)+3.4(+21.0)g cm−3. Due to noisy light curves, we used the more conservative 3σ confidence intervals for the radii as input to the interior structure modelling. We find that HIP 29442 b appears as a typical sub-Neptune, likely surrounded by a gas layer of pure H-He with amass of 0.27−0.17+0.24 M⊕ and a thickness of 1.4 ± 0.5 R⊕. For the innermost companions HIP 29442 c and HIP 29442 d, the model supports an Earth-like composition. Conclusions. The compact multi-planet system orbiting HIP 29442 offers the opportunity to study simultaneously planets straddling the gap in the observed radius distribution of close-in small-size exoplanets. High-precision photometric follow-up is required to obtain more accurate and precise radius measurements, especially for planets c and d. This, together with our determined high-precision masses, will provide the accurate and precise bulk structure of the planets, and enable an accurate investigation of the system’s evolution.

Numerical study of Richtmyer–Meshkov instability of light fluid layer with reshock

The Richtmyer–Meshkov instability of a light fluid layer driven by reflected shock is investigated numerically and theoretically to reveal the feature of the interfacial evolution of light fluid layer under reshock. Cases with single- and double-perturbation helium gas layers are investigated to study the evolution and merging of interface, and comparisons with diverse layer thicknesses and initial perturbation amplitude are conducted to explore the coupling effects of the interface. For the single-perturbed case, the amplitude variation of the left interface exhibits a distinct inflection point as impacted by the reflected shock, and the growth rate in the reflected stage is noticeably larger than that in the incident stage. During the merging process of the interfaces, the displacement difference between the two interfaces and the amplitude growth of the interface play a dominant role before and after the reshock, respectively. For the double-perturbed cases, the head-on collision of the two interfaces' spike occurs when the two interfaces are initially anti-phase, and the spike “catching up with” the bubble occurs when the two interfaces are initially in phase. While the initial fluid layer is very thin, the merging of interfaces accelerates and the interface-coupling effect increases. A modified model has been proposed to predict the amplitude growth of the interface after reshock, which agrees well with the numerical results. The distribution and development of vorticity are similar for the studied cases with different initial amplitudes and fluid layer thicknesses.

Contribution to modeling hydrogen permeation and thickness optimization in blow molded plastic liners for on-board compressed hydrogen tanks

Blow molded plastic liners (BMPLs) for on-board Type IV compressed gaseous hydrogen (CGH2) storage tanks consist of a high-molecular-weight polymer wall that serves as a permeation layer for hydrogen gas. This paper presents the latest numerical model integrated in the Blow View software for predicting the hydrogen permeation in a BMPL. The prediction model is based on Fick's solution-diffusion laws, with pressure and temperature dependent transport properties, for steady/unsteady state permeation regimes through a polymeric wall. An overview of key simulation results for an industrial BMPL for a CGH2 tank, as well as the adopted methodology addressing the challenge to optimize the wall thickness for light weighting and permeability performance, are also presented.

Laser-driven hierarchical "gas-needles" for programmable and high-precision proximity transfer printing of microchips.

Micro-transfer printing (μTP) techniques are essential for advanced electronics. However, current contact/noncontact μTP techniques fail to simultaneously achieve high selectivity and transfer accuracy. Here, a laser projection proximity transfer (LaserPPT) technique is presented, which assembles the microchips in an approach-and-release manner, combining high-precision parallelism with individual chip control. An embedded carbon layer with a thin gas layer is generated by an ultraviolet laser, followed by absorbing heat from the infrared laser, to enable the sequential expansion of hierarchical "gas-needles." The level 1 large gas-needle with a substantially growing height can reduce the gap between the microchip and the receiver. Then, the level 2 small gas-needles enable the gentle release of a chip. Therefore, the LaserPPT can obtain a strong adhesion modulation (~1000 times), excellent size scalability (<100 micrometers), and high transfer accuracy of ~4 micrometers. Last, the assembly of a micro-light-emitting diode display demonstrates the capabilities for deterministic assembly of microarrays.

Preparation and Evaluation of Composite Hydrogel for Reducing the Leakage Rate of Lost Circulation.

Fractured reservoirs are widely distributed and rich in hydrocarbon resources. When encountering fractured reservoirs during the drilling process, it is often accompanied by formation losses characterized by high leak-off rates, causing severe damage to the reservoir and hindering the detection of oil and gas layers, which is not conducive to the accurate and efficient development of the reservoirs. Conventional plugging materials have poor retention effects in the fractures, resulting in the low stability of the sealing layer. The treatment of malignant lost circulation in fractured formations is an urgent problem to be solved in drilling engineering. This article focuses on improving the success rate of formation plugging through the combined use of multiple plugging materials and develops a composite hydrogel that can effectively reduce leakage rates. This hydrogel is mainly composed of materials such as polyvinyl alcohol, borax, and sodium silicate. It has good temperature resistance, maintains good gel strength at 60 °C, and can maintain long-term performance stability under simulated geological water conditions with salinity of 12,500 mg/L. For immersion corrosion by water or gasoline, the amount of corrosion is small and its fundamental performance remains largely unchanged. Through indoor simulation of a leak formation scenario, the hydrogel demonstrates commendable sealing pressure-bearing capacity. In terms of delaying fluid leakage, mixing the hydrogel with cement slurry at a ratio of 1:1 can delay the leakage rate of the cement slurry by a factor of 5.29.

Model of cross-flow interference index and application for multi-layer commingled production in Sulige tight sandstone gas reservoir, Ordos Basin, China

Multi-layer commingled production is the main feature of gas well development in Sulige tight sandstone gas reservoir, Ordos Basin. Identifying potential interference between layers and devising methods for their characterization are crucial considerations for optimizing the development of gas reservoirs. To address these issues, we designed a physical simulation experiment process and interlayer commingled mining schemes, implementing various interlayer combination modes. The results show that a common occurrence in the process of multi-layer commingled production of tight gas and water layers, whether it involves gas layers production alone or simultaneous production of gas and water layers. This phenomenon involves the crossflow of gas and water between layers, resulting in interlayer interference and a subsequent reduction in gas reservoir recovery. Based on these observations, the concept of an interlayer interference index in multi-layer commingled production in tight sandstone gas reservoirs is proposed. The interference index model is obtained by fitting the multiple linear regression method, showcasing its correlation with the physical properties of the reservoir. High water saturation and a significant permeability ratio of the water layer to the gas layer (exceeding the critical value of 1) can result in the early occurrence of interlayer interference and yield a higher interference index. Furthermore, based on the interference index model, a novel method for productivity evaluation of gas wells in tight gas reservoirs is established. The calculations demonstrate that the interference index curve effectively characterizes the interlayer interference performance of gas wells. The productivity and production performance predictions derived from this model align closely with historical production data, affirming the model's effectiveness and accuracy. Therefore, the interference index model emerges as a valuable tool for predicting the productivity and production performance of gas wells in Sulige tight sandstone gas reservoirs. The research results have important theoretical guidance and practical significance for the efficient development of Sulige tight sandstone gas reservoirs.

THE CONTRIBUTION OF THE RELAXATION OF THE HEAT FLOW TO THE CHARACTERISTICS OF THE FUNDAMENTAL HARMONICS OF A NONLINEAR PHOTOACOUSTIC SIGNAL IN CONDENSED MEDIA

Based on the system of nonlinear heat conduction equations, which takes into account the relaxation of the heat flux for the gas layer, the sample, and the substrate, a system of equations is obtained for temperature fluctuations at the fundamental harmonic (FH) of a nonlinear photoacoustic (PA) signal. It is believed that all three layers in the PA-chamber have their own relaxation times, and the irradiated photon flux is modulated according to the harmonic law. The temperature dependence of thermophysical and optical parameters is represented by thermal coefficients. By solving the boundary problem, general expressions were obtained for the oscillatory component of the temperature perturbation in the exhaust gas, and then, according to a well-known scheme, for pressure fluctuations in the gas layer. The obtained expressions are analyzed. The frequency dependence of the amplitude and phase of this signal is numerically calculated for two possible substrates: the substrate is a thermal insulator and the substrate is an ideal conductor. It has been established that the frequency dependence of the PA-signal amplitude for these cases is antiphase. While, this dependence for the phase of the PA-signal is consistent in phases moreover, the values of the maximum and minimum coincide. The only significant difference is that in the region of zero modulation frequencies, the phase of the signal for a substrate made of a thermal insulator is positive, while for a substrate made of an ideal conductor is negative.

Scalable simulation of coupled adsorption and transport of methane in confined complex porous media with density preconditioning

The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. We analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.

Modeling of stratified two‐phase flows with non‐uniform evaporation based on the exact solution of convection equations

Within the framework of the Oberbeck–Boussinesq approximation, the exact solution of equations of thermoconcentration convection is studied, which has the group origin. The issue of applicability of the exact solution for describing steady‐state convective flows of a liquid and a co‐current gas flux under the conditions of inhomogeneous evaporation of the diffusive type in a flat horizontal channel is discussed. Algorithms for obtaining analytical representations of the required functions for various types of the boundary conditions for the temperature function on the outer channel wall are proposed. The influence of the external thermal load and boundary thermal conditions on the structure of the velocity and temperature fields, evaporation mass flow rate, and vapor content in the gas layer is investigated at an example of the HFE‐7100–nitrogen system. The solution correctly predicts hydrodynamical, temperature, and concentration parameters of convective regimes appearing in the two‐phase system. The main characteristics are compared with the characteristics of the system in the case of uniform evaporation with a constant intensity.

Sensible heat aided gas production from gas hydrate with an underlying water-rich shallow gas layer

Increasing the efficiency of gas production from marine gas hydrate still faces significant problems. Based on the geological survey of marine reservoirs, gas hydrates are generally associated with underlying free gas layers with high temperature and water saturation. Here the effect of this gas layer and its characteristics on the gas and water production from the hydrate layer are studied. It was found that a large temperature difference existed between the two layers during gas production due to the different heat demand and consumption in each layer. Slower gas production rate was observed in case of dense distribution of hydrates in the hydrate-bearing layer and the resulting sluggish gas permeation and enormous heat supply. Notably, this was not present when the free gas layer was more water-saturated. The resulting maximum gas production rate could be effectively enhanced by 41 % at most as well. This was attributed to the increased sensible heat available from the pore water and the interbedded heat exchange. It was thus indicated that a hydrate reservoir with an underlying water-rich layer could be more favorable for gas production in terms of heat supply.

Combined model for 15N, 13C, and spin-state chemistry in molecular clouds

We present a new gas-grain chemical model for the combined isotopic fractionation of carbon and nitrogen in molecular clouds. To this end, we have developed gas-phase and grain-surface chemical networks where the isotope chemistry of carbon and nitrogen is coupled with a time-dependent description of spin-state chemistry, which is important for nitrogen chemistry at low temperatures. We updated the rate coefficients of some isotopic exchange reactions considered previously in the literature, and here we present a set of new exchange reactions involving molecules substituted in 13C and 15N simultaneously. We applied the model to a series of zero-dimensional simulations representing a set of physical conditions across a prototypical prestellar core, exploring the deviations of the isotopic abundance ratios in the various molecules from the elemental isotopic ratios as a function of physical conditions and time. We find that the 12C/13C ratio can deviate from the elemental ratio to a high factor depending on the molecule, and that there are highly time-dependent variations in the ratios. The HCN/H13CN ratio, for example, can obtain values of less than ten depending on the simulation time. The 14N/15N ratios tend to remain close to the assumed elemental ratio within approximately 10%, with no clearly discernible trends for the various species as a function of the physical conditions. Abundance ratios between 13C-containing molecules and 13C+15N-containing molecules however show somewhat increased levels of fractionation as a result of the newly included exchange reactions, though they still remain within a few tens of percent of the elemental 14N/15N ratio. Our results imply the existence of gradients in isotopic abundance ratios across prestellar cores, suggesting that detailed simulations are required to interpret observations of isotopically substituted molecules correctly, especially given that the various isotopic forms of a given molecule do not necessarily trace the same gas layers.

Effect of combustion chamber diameter size on scramjet rotating detonation under high Mach number flow conditions

This study conducted 3D numerical simulations to investigate the impact of combustion chamber diameter on the combustion characteristics of liquid kerosene scramjet rotating detonation under high Mach number flow conditions (Ma6/28 km). After ignition, a local hotspot near the contact surface between the combustion product and fresh reactant facilitates the generation of new detonation wave heads, resulting in a co-directional multi-wave mode in the detonation combustion flow field. The higher total temperature of the incoming flow restricts the accumulation of a substantial fuel gas layer in the axial direction of the combustion chamber, resulting in a smaller fuel distribution area and a lower wave head. Increasing the inner diameter of the combustion chamber leads to an increase in the number of wave heads but a decrease in overall height. Specifically, when using diameters of 125 and 150 mm, we observed significant periodic low-frequency oscillations in the peak pressure of the detonation wave during stable propagation. The specific impulse of the fuel does not vary significantly across different combustion chamber diameters. However, when the inner diameter is 75 mm, periodic oscillations occur, which reduce thrust stability. These findings provide valuable insights into optimizing combustion chamber design and improving the efficiency and stability of liquid kerosene scramjet rotating detonation propulsion systems under high Mach number flow conditions.

Richtmyer–Meshkov instability of a single-mode heavy–light interface in cylindrical geometry

Richtmyer–Meshkov (RM) instability of a single-mode SF6–air interface subjected to a convergent shock is investigated experimentally. The convergent shock tube is specially designed with an opening tail to weaken the Rayleigh–Taylor effect and eliminate the reflected waves' effect. The gas layer scheme is used to create a heavy gas environment at the upstream side of the interface. Before phase inversion is finished, the amplitude reduction is accelerated, but the Bell–Plesset (BP) effect in this process is found to be negligible. After phase inversion is completed, the linear growth rate is generally predicted due to small amplitude and the weak BP effect. In nonlinear regime, an existing nonlinear model is revised based on the Padé approximation to give a better prediction of amplitude growth. The spike amplitude grows almost linearly, whereas the bubble amplitude gradually saturates and even reduces. For a heavy-light interface in convergent geometry, although both the spike and bubble amplitude growths are promoted by the BP effect, the spike growth is more promoted than the bubble. The BP effect enhances generation of the second-order harmonic, which results in saturation and reduction of the bubble amplitude. The discrepancy in the BP effect between light-heavy and heavy-light interfaces is qualitatively demonstrated for the first time.

REDUCED HEAT LOSS THROUGH WINDOW STRUCTURES

The results of theoretical and experimental studies of heat transfer through the window structures of the building are presented, the thermal characteristics of the windows are established and the level of heat loss realized through them is estimated. To determine the heat transfer resistance of two-chamber double-glazed units, numerical modeling of radiative-convective heat transfer is performed, taking into account the thermal conductivity of the gas layer between the glasses. Natural free convection of the gaseous medium in the double-glazed chambers occurs as a result of the temperature difference between the side glass surfaces of the chamber. Modeling is performed by numerically solving the system of fluid dynamics and energy equations for the air layer and glass. On the inner surfaces of the chambers, conditions of the fourth kind are set, which take into account the radiation and conductive components of the total heat flow coming from the glass surface. The results of simulation of heat transfer through double-glazed windows with ordinary glass indicate that about 60 % of heat is transferred by radiation. Therefore, an effective measure to reduce heat loss through windows is to reduce the radiation component of the total heat flow by applying a low-emissivity coating to the inner surfaces of the glass unit. This makes it possible to reduce the total heat flow (and, accordingly, heat loss to the environment) by 20–34 %, depending on the number of glass surfaces with this coating. Conducted independent experimental studies confirm the modeling results. In order to comply with the current requirements of state building regulations and to achieve the energy efficiency class of type C buildings (specific heat consumption for heating 1 sq. m of usable area is less than 75 kWh), for the climatic conditions of Kyiv, it is recommended to use double-chamber windows with two low-emission coatings with the formula of the glass unit: 4M1-10-i4M1-10-i4M1. Bibl. 22, Fig. 8.

Adiabatic wall temperature in the supersonic flow of moist air with spontaneous condensation

The expansion of moist air in supersonic nozzles is accompanied by spontaneous condensation of water vapor due to a sharp decrease in flow temperature and pressure. However, in the near-wall boundary layer due to the viscous dissipation, moist air temperature decreases slightly in comparison with initial one. This leads to the fact that the vapor remains overheated in the near-wall region and condensation does not spread to its entire depth. Droplets formed in the core flow as a result of condensation can penetrate into the near-wall region due to turbulent diffusion and evaporate due to higher temperature in the near-wall region. This process causes cooling of the near-wall gas layers and, consequently, a decrease in the adiabatic wall temperature. A decrease in the adiabatic wall temperature relative to the stagnation temperature leads to a decrease in aerodynamic heating and, consequently, in the heat flux from the gas to the body. In this paper, measurements of the adiabatic temperature of the supersonic nozzle wall for moist air flow with condensation are presented. Temperature measurements were made using an infrared thermal imager. Initial moisture content varied in the range of 1.5–20.2 g/(kg. dry air) and initial relative humidity in the range 20–90%. It is shown that in the region behind the zone of spontaneous condensation (“condensation shock”), adiabatic wall temperature depends on the initial moisture content and can be either greater or less than the value of a single-phase (dry air) flow with identical initial flow parameters. Estimates of the flow parameters behind the condensation region are obtained using one-dimensional diabatic flow equations and experimental static pressure distributions along the nozzle length. The correlation between the initial moisture content, the average droplet size and the temperature recovery factor is shown.

Review on hydrodynamic instabilities of a shocked gas layer

Review on hydrodynamic instabilities of a shocked gas layer

Near-Wall Cavitation Effect: A Molecular Dynamics Study.

Cavitation has been the subject of abundant studies, but the internal mechanism of cavitation is less well known. In this article, a microlevel near-wall model was established by using LAMMPS to present the process of cavitation effect. The results of molecular dynamics simulation revealed the fluctuation process of the liquid near the wall with the change in pressure. Molecular dynamics was also used to evaluate the void volume fraction and density distribution of the system. The results exhibited that the cavitation process can be divided into two stages: the initial cavitation stage and the rapid growth stage. Based on these results, the effects of wettability and initial system temperature on the near-wall cavitation effect were demonstrated. The results indicated that the hydrophobic near-wall forms a gas layer to weaken the density fluctuation, while the hydrophilic wall is opposite. Increasing the temperature could positively affect molecular motion and cavitation. This work provides a theoretical basis for further exploration of the cavitation effect.

Numerical study on natural gas hydrate production by hot water injection combined with depressurization

The low efficiency and unstable gas production of hydrate exploitation have become the difficulties restricting the commercial development of combustible ice. In this work, based on the existing geological conditions, this paper establishes a low-permeability hydrate heterogeneous reservoir model to analyze the effects of hot water injection on well and reservoir. The results show that the greater the wasted heat from the produced water and some injected water migrating to the gaseous hydrocarbon layer and the underlying at the middle and late stage of hydrate production as the hot water injection rate. Moreover, the gas and water production will suddenly appear surge phenomenon during the hydrate production due to the disappearance of the water lock effect affected by thermal-hydraulic-chemical multi-field coupling, which depends on the hot fluid properties. Therefore, the new technology of slow hot water injection after depressurization or different injection rates at different stages in the injection well is proposed in this paper. And it is also necessary to optimize the vertical position of the injection well or add a thin layer with low permeability at the bottom of the gaseous hydrocarbon layer. Finally, the hot water injection temperature should be controlled above 50 °C.