Journal of Natural Gas Science and Engineering最新文献

A spiral multilateral well network is a promising production method to enhance long-term gas recovery from prevalent fine-grained hydrate reservoirs. However, practical application is greatly restricted before the optimal well network parameters are determined and the mechanism behind a unique phenomenon in multilateral wells, namely inter-branch interference, is clear. In this study, we numerically optimized the well configuration and spacing when spiral multilateral wells were deployed in two typical fine-grained hydrate reservoirs, i.e., ultra-low permeability hydrate reservoirs (ULPHR, <1 mD) and low-permeability hydrate reservoirs (LPHR, >1 mD). The mechanism behind inter-branch interference was innovatively revealed. The results indicated that the number of spiral branches should be increased, and equidistant branches should be deployed uniformly in the lower ULPHR or throughout LPHR to enhance production efficiency. A wide spacing of spiral multilateral wells with long branches contributed to long-term productivity in fine-grained hydrate reservoirs with any permeability; however, narrow spacing was more favorable for short branches or short-term production. Our study found three inter-branch interference stages during gas production, namely, “no effect” stage, “positive” stage, and “negative” stage; all the three stages are controlled by reservoir permeability, production distance, and production time. Owing to the “positive” interference effect, longer equal-length branches resulted in superior long-term production enhancement in ULPHR, particularly for lengths greater than 30 m. Gas production from LPHR using only two optimal spiral multilateral wells exhibited high production performance similar to that of the sandy hydrate deposits in Japan, suggesting that the optimal spiral multilateral well network is promisingly suitable for commercial production in the future.

The dispersion and attenuation of seismic-wave propagation induced by ‘squirt flow’ effects in hydrocarbon-saturated reservoirs are significantly affected by their rock properties and fluid content. In this study, we analyse the frequency-dependent velocity, attenuation, and seismic responses when fractured porous rock is saturated with two immiscible fluids. First, when considering reservoir wettability, we calculate the effective fluid viscosity using a stable parameter, the capillary pressure, and a lattice Boltzmann model (LBM)-based relative permeability equation, which is a function of the saturation and viscosity ratio of the immiscible two-phase fluid. Then, we explore the frequency-dependent effects of fractured porous rocks saturated with two immiscible fluids under different cases of viscosity ratios and capillary pressure parameters by employing the Chapman model from the dynamic equivalent-medium theory. Then, we use a four-layer model to analyse the frequency-dependent seismic responses. The results show that the characteristics of frequency-dependent velocity and attenuation are both affected by the wettability, capillary pressure parameter, saturation, and viscosity ratio. The frequency-dependent features are greatly influenced by the capillary pressure parameter and viscosity ratio. For a larger viscosity ratio and lower capillary parameter, a dispersive effect can occur in the seismic frequency band. This indicates that the velocity dispersion anomalies are sensitive to wettability, capillary pressure parameter and viscosity ratio and should not be neglected. Synthetic seismic records demonstrate that the seismic reflection signatures, such as the waveform, amplitude, and reflective travel time, at the interfaces for saturated reservoirs are significantly affected by wettability and saturation. The numerical modeling helps to improve the wave propagation in rocks saturated by two immiscible fluids.

In this study, the non-intrusive EDFM (embedded discrete fracture model) method was presented to investigate the impact of different choke management strategies on well performance. Through the EDFM method, accurate simulation can be conducted to efficiently evaluate the fracture complexities. First, by implementing this powerful technology, a horizontal well with multi-stage hydraulic and natural fractures was set up, where the permeability can be distributed sequentially in each hydraulic fracture segment. Then various pressure drawdown scenarios from conservative to aggressive strategy were designed. The different levels of fracture closure can be properly modeled in each state. Additionally, pressure distribution for the matrix and fractures was depicted to provide straightforward insights for the choke management under two extreme strategies. Subsequently, a series of sensitivity studies were presented to evaluate the impacts of various factors on shale gas production, including fracture permeability modulus, fracture closure, and natural fractures network. The simulation results show that choke management can be simulated effectively by applying EDFM. After considering the fracture closure behavior and complex fracture networks, the conservative drawdown strategy can be addressed as the optimal strategy for the EUR, as it improves the cumulative gas production by maintaining the hydraulic fracture open through a steady pressure decline. The remained proppants enhance the fracture conductivity, thereby expanding its drainage influence towards larger zones of the reservoir. The influence of natural fractures, including the fracture length, fracture number, and fracture conductivity, are also studied. All these three variables play a significant impact on well performance. Consequently, the model becomes a valuable stencil to design fracture closure and complex fracture networks, which can be applied to optimize the choke management design for unconventional reservoirs.

Improper pairing of filler and polymer together with inappropriate filler loadings into polymer matrix may lead to structural defects such as large aggregations and interface void formations. Subsequently, the structural defects may sacrifice the selectivity of CO₂ over CH_4, which was unfavorable. In the current work, NH₂-MIL-125 (Ti) (MIL = Material Institute Lavoisier), which possesses NH₂-groups and theoretically capable of forming strong hydrogen bonding with F-groups of polyvinylidene fluoride (PVDF), was selected to spin hollow fiber mixed matrix membranes (HFMMMs). Besides, NH₂-MIL-125 (Ti) can interact better with CO₂ over CH₄ via quadrupole moment, and NH₂-groups also aid in CO₂ selectivity due to its high CO₂ adsorption capability. The HFMMMs were spun using a dry-wet spinning technique of filler loadings percentage ranging from 1 to 3 wt percent (wt%). The effect of filler and loadings percentage over HFMMMs properties, including contact angle, mechanical strength, thermal stability and cross-sectional morphology was investigated. The compatibility at interface of filler and polymer was observed to be good, and dispersion was observed to be acceptable up to 2 wt% filler loadings. However, apparent aggregation was observed beyond this point. The wt% of Ti, O, and N elements were found to increase from 0.72 to 2.05, 3.27 to 4.53, and 0.52 to 1.55, respectively, with increasing filler loading into HFMMMs. Subsequently, PVDF-2 membrane displayed the highest CO₂/CH₄ ideal selectivity with contact angle of 83.44 ± 1.45, ultimate tensile strength (UTS) of 1.33, 29.12 Young's Modulus, and 72.2% elongation at break. Therefore, optimizing loading percentage and selecting appropriate filler are considered practical methods to ensure good morphology and better hazardous CO₂ removal.

Understanding the gas-water two-phase flow behavior in a shale gas formation is important for reservoir simulation and production optimization. A molecular simulation study of gas-water flow in quartz and kerogen nano-slits (2–6 nm) at shale reservoir conditions (temperature: 313.15–393.15 K, pressure: 20–60 MPa) is reported in this work. The simulation results show that the existence of water in the hydrophilic quartz slits will form water films on the slit walls; while the presence of water in the hydrophobic kerogen slits will form water clusters in the central of the gas phase at low water saturation and a water layer at high water saturation. In both wetting conditions, water will take flow space and reduce gas flow path. However, water affects gas flow velocities in the two wetting types nano-slits in different ways due to the two opposite occupancies in the slits. The momentum transfer between water and methane molecules in the gas-water interface region plays an important role in the gas-water two-phase flow. The gas flow is more readily affected by water content in the quartz slit with an aperture greater than 2 nm. When the slit aperture is reduced to 2 nm, it is difficult to form a continuous gas or water phase, and the existence of water in both types of slits will reduce the velocity of methane. Increasing the temperature will accelerate the flow of methane and water because hydrogen bonds between water molecules as well as hydrogen bonds between water molecules and the walls are reduced. High pressure promotes the mixing of the methane and water molecules, resulting in the gas velocity decreasing in both quartz and kerogen slits. The flow mechanism of methane and water in nano-slits provide insights into theoretical models for shale gas production.

Drilling operations in oil and gas industry are associated with issues such as pipe sticking, poor wellbore cleaning and high fluid loss. Mitigation of such problems in water-based drilling mud (WBM) necessitates the application of nanotechnology in improving their filtration and rheological characteristics. In the present work, an attempt has been made to analyze the effect of nanofluid prepared using copper oxide (CuO) nanoparticles (NPs) dispersed in chia seed solution on WBM characteristics. Therefore, three samples of chia seed based nanofluids are synthesized using two-step method by varying the concentration of CuO nanoparticle from 0.2 wt% to 0.6 wt%. The resulting nanofluids are then mixed with WBM to prepare Nanofluid enhanced Water based Drilling Mud (NFWBM). The synthesized nanofluids are then characterized for their stability and thermal decomposition respectively using Scanning Electron Microscope (SEM) and Thermo-Gravimetric Analyzer (TGA). The NFWBMs are then analyzed for rheological and filtrate-loss properties at different temperatures of 30 °C, 50 °C, 70 °C and 90 °C. The hot roll aging process is carried out at 90 °C for 16 h maintaining the pressure at 0.1 MPa. The analysis projected a significant enhancement in the thermal stability of the WBM, with a reduction in viscosity of about 61.7% at 90 °C, which is critically observed to recover back to a significant extent of about 14% for chia based 0.4 wt% CuO nanofluid enhanced WBM and 19% for chia based 0.6 wt% CuO nanofluid enhanced WBM. Such improvement is observed in the rheological properties post hot rolling too. Further, the API fluid loss is observed to reduce from 7.2 ml to 6.8 ml, 6 ml, and 4.8 ml, respectively, before hot rolling, while the same reduced from 12.4 ml 11.4 ml, 10.2 ml, and 9.4 ml, respectively, for chia based 0.2 wt%, 0.4 wt%, and 0.6 wt% of CuO nanofluid enhanced water based drilling muds (NFWBMs). The present study aids in the development of novel and green additives for water-based muds to enhance their properties.

The liquid-phase pipeline is the optimal choice for large-amount and long-distance ethane transportation. Selecting the optimal diameter is necessary for the economical design of the pipeline. However, the special critical temperature 32.2 °C and critical pressure 4.87 MPa of ethane makes it easy to become liquid-vapor phase change, which is not considered in the traditional natural gas or crude oil pipeline design. In this paper, a new mathematical model is built to calculate the optimal diameter of the ethane pipe by using the ‘pump station + pipeline’ unit as the research object. The model selects the lowest total pipeline construction and operation costs as the objective function, and the constraints include the ethane liquid-vapor phase change, the pipe maximum allowable stress, and pipe specifications. In particular, the liquid-vapor phase change constraint is added to the traditional model to avoid the ethane liquid-vapor phase change, which is obtained by quantitatively analyzing the variation of physical parameters of ethane close to the pressure-temperature phase boundary. The optimization model is solved by use of the genetic algorithm. Finally, optimal pipe diameters are calculated for the conditions of transmission capacity from 1000 t/d to 10,000 t/d. Comparisons of calculated pipe diameter with eight actual cases show that the results are feasible with the average and maximum deviations being less than 5% and 8%, respectively. The effects of pipe materials and electricity prices on the pipe diameter are analyzed. It is demonstrated that the pipe material has a negligible effect on the optimal diameter, whereas increasing the electricity price will lead to the increase of the optimal diameter in the case of large transmission volumes.

Large-scale slickwater fracturing is an important technical method for the effective development of shale gas, which generates complex fractures with fracture width of millimeters in the reservoir. It is known that the transport law of proppant in complex fractures is the premise for realizing effective propping. Taking the behavior of 70/140 mesh proppant particles commonly used in shale gas fracturing as the object，a numerical model based on the computational fluid dynamics-discrete element method (CFD-DEM) and the geometric model of the intersection of the main and secondary fractures are established. It is used to study the two-phase flow law under the conditions of different fracture widths and angles, pump displacements, and fluid viscosities. The results show that the proppant enters the secondary fracture in two ways: carried by the fluid in a suspended manner and rolling into the fracture from the sand bank surface. Particles suspended in the fracture can be transported to the distal end of the secondary fracture. Owing to the influence of the inertia force of particles, the particle flow rate entering the secondary fracture is much smaller than the fluid flow rate in the secondary fracture. As the included angle between the secondary and main fractures decreases, the fluid and particle flow rate increase, and particles can easily enter the secondary fracture. As the displacement, secondary fracture width and fracturing fluid viscosity increase, proppant particles are easier to enter secondary fractures. The absolute values of the main and secondary fracture widths become smaller, and the relative value remains the same, making it more difficult for proppant particles to enter the fractures.

An understanding of the mechanical properties of reservoir shale is of great significance for the efficient development of shale gas. Both marine and continental shale gas reservoirs in China have considerable development potential, but their different depositional environments may lead to substantial differences in their mechanical properties, which can result in production efficiency differences. In this study, nanoindentation, X-ray diffraction, backscattered electron imaging, and energy-dispersive X-ray spectroscopy were used to determine and analyze the mechanical properties and microtexture of marine and continental shale samples. The geogenesis of the microtexture of marine and continental shales and its influence on the mechanical properties were discussed. The results show that the elastic modulus of marine and continental shale samples are similar, but the hardness of latter is greater. The similar elastic modulus of the two shale samples may result from similar mineralogy. Due to differences in deposition and diagenesis, the marine shale sample forms a clay support matrix and the continental shale sample forms a rigid clastic support matrix, which results in lower hardness in the former and higher hardness in the latter. The low hardness of the shale with a clay support matrix indicates that it may be subject to more severe proppant embedment issue. The experimental results provide a useful reference for the development of these two types of shale gas reservoirs.

In the Czech Republic, coal-bearing siliciclastic sediments have been deposited during the Serpukhovian and Bashkirian (Carboniferous). Until now, no attention has been paid to inorganic geochemical assessment of the coals and associated non-coal rocks from the mixed shallow-marine to continental sediments (Ostrava Formation), and continental non-marine settings (Karviná Formation). Samples were collected from a 750 m deep coal exploration borehole at the ČSM Mine. The bulk parameters, total organic carbon TOC, total inorganic carbon TIC, total sulphur TS, major elements, trace elements, and REEs were measured on these samples, and their mineral associations have been investigated using microscopy combined with the principal component analysis (PCA). Common redox proxies V/Cr, U/Th, Ni/Co, Mo/U, and the ratio S/TOC have been tested on the samples to investigate their usefulness for studying anoxia. Research concludes that redox proxies such as U/Th, Ni/Co and V/Cr have been strongly influenced by the clastic input and carbonates, which it hinders for them to be reliable indicators of anoxia. On the basis of Eu anomaly and REEs distribution, the primary source of detrital elements comes from the parent rock, being governed more by physical than redox processes.