Prediction of Sedimentary Microfacies Distribution by Coupling Stochastic Modeling Method in Oil and Gas Energy Resource Exploitation

In view of the problem that a single modeling method cannot predict the distribution of microfacies, a new idea of coupling modeling method to comprehensively predict the distribution of sedimentary microfacies was proposed, breaking the tradition that different sedimentary microfacies used the same modeling method in the past. Because different sedimentary microfacies have different distribution characteristics and geometric shapes, it is more accurate to select different simulation methods for prediction. In this paper, the coupling modeling method was to establish the distribution of sedimentary microfacies with simple geometry through the point indicating process simulation, and then predict the microfacies with complex spatial distribution through the sequential indicator simulation method. Taking the DC block of Bohai basin as an example, a high-precision reservoir sedimentary microfacies model was established by the above coupling modeling method, and the model verification results showed that the sedimentary microfacies model had a high consistency with the underground. The coupling microfacies modeling method had higher accuracy and reliability than the traditional modeling method, which provided a new idea for the prediction of sedimentary microfacies.

Calcium sulfate whisker prepared by flue gas desulfurization gypsum:A physical–chemical coupling production process

In this paper,the solid waste desulfurization gypsum produced by coal-fired power plants was used as a raw material to prepare calcium sulfate whiskers with high application prospects.Calcium sulfate whiskers with uniform morphology and high aspect ratio can be prepared by hydrothermal method in sulfuric acid solution.A new process of desulfurization gypsum activated by high-energy grinding to reduce the reaction temperature and sulfuric acid concentration was developed.Through the comparison of product morphology,the best grinding time was determined to be 3.5 h.The mechanism of desulfurization gypsum through physical–chemical coupling to reduce energy consumption was clarified.The activation of desulfurization gypsum by grinding and the acidic environment provided by the sulfuric acid solution made the calcium sulfate solution reached rapid saturation and accelerated the nucleation rate.By calculating the conversion and crystallization rate of calcium sulfate whiskers,it was found that there were obvious"autocatalytic"kinetic characteristics during the crystallization process.

Mechanical criterion for coal and gas outburst:a perspective from multiphysics coupling

Although a series of hypotheses have been proposed,the mechanism underlying coal and gas outburst remains unclear.Given the low-index outbursts encountered in mining practice,we attempt to explore this mechanism using a multiphysics coupling model considering the effects of coal strength and gas mass transfer on failure.Based on force analysis of coal ahead of the heading face,a risk identification index C_(m)and a critical criterion(C_(m)≥1)of coal instability are proposed.According to this criterion,the driving force of an outburst consists of stress and gas pressure gradients along the heading direction of the roadway,whereas resistance depends on the shear and tensile strengths of the coal.The results show that outburst risk decreases slightly,followed by a rapid increase,with increasing vertical stress,whereas it decreases with increasing coal strength and increases with gas pressure monotonically.Using the response surface method,a coupled multi-factor model for the risk identification index is developed.The results indicate strong interactions among the controlling factors.Moreover,the critical values of the factors corresponding to outburst change depending on the environment of the coal seams,rather than being constants.As the buried depth of a coal seam increases,the critical values of gas pressure and coal strength decrease slightly,followed by a rapid increase.According to its controlling factors,outburst can be divided into stress-dominated,coal-strength-dominated,gas-pressure-dominated,and multi-factor compound types.Based on this classification,a classified control method is proposed to enable more targeted outburst prevention.

Heat-fluid-solid coupling model for gas-bearing coal seam and numerical modeling on gas drainage promotion by heat injection

Improving the absorbed gas to active desorption and seepage and delaying gas drainage attenuation are considered as key methods for increasing drainage efficiency and gas output.According to the solid mechanics theory,the nonlinear Darcy seepage theory and thermodynamics,the heat-fluid-solid coupling model for gassy coal has been improved.The numerical model was founded from the improved multi-field coupling model by COMSOL Multiphysics and gas drainage by borehole down the coal seam enhanced by heat injection was modelled.The results show that the heatfluid-solid model with adsorption effects for gassy coal was well simulated by the improved multi-field model.The mechanism of coal seam gas desorption seepage under the combined action of temperature,stress and adsorption can be well described.Gas desorption and seepage can be enhanced by heat injection into coal seams.The gas drainage rate was directly proportional to the temperature of injected heat in the scope of 30-150 ℃ and increasing in the whole modelleddrainage process (0-1000 d).The increased level was maximum in the initial drainage time and decreasing gradually along with drainage time.The increasing ratio of drainage rate was maximum when the temperature raised from 30 to 60 ℃.Although the drainage rate would increase along with increasing temperature,when exceeding 60 ℃,the increasing ratio of drainage rate with rising temperature would decrease.Gas drainage promotion was more effective in coal seams with lower permeability than with higher permeability.The coal seam temperature in a 5 m distance surrounding the heat injection borehole would rise to around 60 ℃ in 3 months.That was much less than the time of gas drainage in the coal mines in sites with low permeability coal seams.Therefore,it is valuable and feasible to inject heat into coal seams to promote gas drainage,and this has strong feasibility for coal seams with low permeability which are widespread in China.

Fluid-solid coupling model for studying wellbore instability in drilling of gas hydrate bearing sediments

As the oil or gas exploration and development activities in deep and ultra- deep waters become more and more, encountering gas hydrate bearing sediments （HBS） is almost inevitable. The variation in temperature and pressure can destabilize gas hydrate in nearby formation around the borehole, which may reduce the strength of the formation and result in wellbore instability. A non-isothermal, transient, two-phase, and fluid-solid coupling mathematical model is proposed to simulate the complex stability performance of a wellbore drilled in HBS. In the model, the phase transition of hydrate dissociation, the heat exchange between drilling fluid and formation, the change of mechanical and petrophysical properties, the gas-water two-phase seepage, and its interaction with rock deformation are considered. A finite element simulator is developed, and the impact of drilling mud on wellbore instability in HBS is simulated. Results indicate that the re- duction in pressure and the increase in temperature of the drilling fluid can accelerate hydrate decomposition and lead to mechanical properties getting worse tremendously. The cohesion decreases by 25% when the hydrate totally dissociates in HBS. This easily causes the wellbore instability accordingly. In the first two hours after the formation is drilled, the regions of hydrate dissociation and wellbore instability extend quickly. Then, with the soaking time of drilling fluid increasing, the regions enlarge little. Choosing the low temperature drilling fluid and increasing the drilling mud pressure appropriately can benefit the wellbore stability of HBS. The established model turns out to be an efficient tool in numerical studies of the hydrate dissociation behavior and wellbore stability of HBS.

Simulation of dilatancy-controlled gas migration processes in saturated bentonite using a coupled multiphase flow and elastoplastic H2 M model

Dilatancy-controlled gas flow in preferential pathways plays a key role in the safety analysis of radioactive waste repositories.This is particularly the case for bentonite,an often-preferred barrier material.Gas flow in preferential pathways is characterized by localization and spontaneous behavior,which is challenging to simulate in numerical models due to strong hydro-mechanical coupling.To analyze a laboratory experiment in the framework of the DECOVALEX-2023 project,this study introduced a new approach of combining continuous modelling methods with spatial material properties derived from material heterogeneities and experimental observations.The proposed model utilized hydro-mechanical spatial distributions,namely Young’s modulus and gas entry pressure,and elastoplasticity combined with a linear swelling model.A conceptual strain-dependent permeability approach simulated dilatancycontrolled gas flow based on hydro-mechanical coupling.To test the effectiveness of the presented approach,a gas injection test in a compacted,saturated bentonite sample was simulated using the opensource code OpenGeoSys 5.8 and compared with experimental observations.The presented methodology is capable of simulating localized gas flow in preferential pathways.The spatial distributions of Young’s modulus and gas entry pressure affect the swelling pressure,relative permeability and,in combination with the strain-dependent permeability model,also the intrinsic permeability.

Dynamic Coupling Correlation of Gas Film in Dry Gas Seal with Spiral Groove

In working state, the dynamic performance of dry gas seal, generated by the rotating end face with spiral grooves, is determined by the open force of gas film and leakage flow rate. Generally, the open force and the leakage flow rate can be obtained by finite element method, computational fluid dynamics method and experimental measurement method. However, it will take much time to carry out the above measurements and calculations. In this paper, the approximate model of parallel grooves based on the narrow groove theory is used to establish the dynamic equations of the gas film for the purpose of obtaining the dynamic parameters of gas film. The nonlinear differential equations of gas film model are solved by Runge-Kutta method and shooting method. The numerical values of the pressure profiles, leakage flux and opening force on the seal surface are integrated, and then compared to experimental data for the reliability of the numerical simulation. The results show that the numerical simulation curves are in good agreement with experimental values. Furthermore, the opening force and the leakage flux are proved to be strongly correlated with the operating parameters. Then, the function-coupling method is introduced to analyze the numerical results to obtain the correlation formulae of the opening force and leakage flux respectively with the operating parameters, i.e., the inlet pressure and the rotating speed. This study intends to provide an effective way to predict the aerodynamic performance for designing and optimizing the groove styles in dry gas seal rapidly and accurately.

Optimization method of refracturing timing for old shale gas wells

Based on the elastic theory of porous media,embedded discrete fracture model and finite volume method,and considering the micro-seepage mechanism of shale gas,a fully coupled seepage-geomechanical model suitable for fractured shale gas reservoirs is established,the optimization method of refracturing timing is proposed,and the influencing factors of refracturing timing are analyzed based on the data from shale gas well in Fuling of Sichuan Basin.The results show that due to the depletion of formation pressure,the percentage of the maximum horizontal principal stress reversal area in the total area increases and then decreases with time.The closer the area is to the hydraulic fracture,the shorter the time for the peak of the stress reversal area percentage curve to appear,and the shorter the time for the final zero return(to the initial state).The optimum time of refracturing is affected by matrix permeability,initial stress difference and natural fracture approach angle.The larger the matrix permeability and initial stress difference is,the shorter the time for stress reversal area percentage curve to reach peak and return to the initial state,and the earlier the time to take refracturing measures.The larger the natural fracture approach angle is,the more difficult it is for stress reversal to occur near the fracture,and the earlier the optimum refracturing time is.The more likely the stress reversal occurs at the far end of the artificial fracture,the later the optimal time of refracturing is.Reservoirs with low matrix permeability have a rapid decrease in single well productivity.To ensure economic efficiency,measures such as shut-in or gas injection can be taken to restore the stress,and refracturing can be implemented in advance.

Optimizing the Diameter of Plugging Balls in Deep Shale Gas Wells

Deep shale gas reserves that have been fractured typically have many relatively close perforation holes. Due to theproximity of each fracture during the formation of the fracture network, there is significant stress interference,which results in uneven fracture propagation. It is common practice to use “balls” to temporarily plug fractureopenings in order to lessen liquid intake and achieve uniform propagation in each cluster. In this study, a diameteroptimization model is introduced for these plugging balls based on a multi-cluster fracture propagationmodel and a perforation dynamic abrasion model. This approach relies on proper consideration of the multiphasenature of the considered problem and the interaction force between the involved fluid and solid phases. Accordingly,it can take into account the behavior of the gradually changing hole diameter due to proppant continuousperforation erosion. Moreover, it can provide useful information about the fluid-dynamic behavior of the consideredsystem before and after plugging. It is shown that when the diameter of the temporary plugging ball is1.2 times that of the perforation hole, the perforation holes of each cluster can be effectively blocked.

Spin-orbit coupling adjusting topological superfluid of mass-imbalanced Fermi gas

Topological superfluid state is different from the normal superfluid one due to the excitation energy gap on the boundary.How to obtain the topological superfluid state by using spin-orbit coupling to control the s-waves paired mass-imbalanced Fermi gas is a recent novel topic.In this paper,we study the topological superfluid phase diagram of two-dimensional mass-imbalanced Fermi gas with Rashba spin-orbit coupling at zero temperature.We find that due to the competition among mass imbalance,pairing interaction and spin-orbit coupling,there is a double-well structure in the thermodynamic potential,which affects the properties of the ground state of the system.We comprehensively give the phase diagrams of the system on the plane of spin-orbit coupling and chemical potential,and the phase diagrams on the plane of the reduced mass ratio and two-body binding energy.This study not only points out the stable region of topolog-ical superfluid state of mass-imbalanced Fermi gas,but also provides a detailed theoretical basis for better observation of topological superfluid state in experiments.

Two-way coupling model for shock-induced laminar boundary-layer flows of a dusty gas

The present paper describes a numerical two-way coupling model for shock-induced laminar boundary-layer flows of a dust-laden gas and studies the transverse migration of fine particles under the action of Saffman lift force. The governing equations are formulated in the dilute twophase continuum framework with consideration of the finiteness of the particle Reynolds and Knudsen numbers. The full Lagrangian method is explored for calculating the dispersedphase flow fields （including the number density of particles） in the regions of intersecting particle trajectories. The computation results show a significant reaction of the particles on the two-phase boundary-layer structure when the mass loading ratio of particles takes finite values.

Topological superfluid in a two-dimensional polarized Fermi gas with spin-orbit coupling and adiabatic rotation

We study the properties of superfluid in a two-dimensional （2D） polarized Fermi gas with spin-orbit coupling and adiabatic rotation which are trapped in a harmonic potential. Due to the competition between polarization, spin-orbit coupling, and adiabatic rotation, the Fermi gas exhibits many intriguing phenomena. By using the Bardeen-Cooper-Schrieffer （BCS） mean-field method with local density approximation, we investigate the dependence of order parameter solution on the spin-orbit coupling strength and the rotation velocity. The energy spectra with different rotation velocities are studied in detail. Besides, the conditions for the zero-energy Majorana fermions in topological superfluid phase to be observed are obtained. By investigating distributions of number density, we find that the rotation has opposite effect on the distribution of number density with different spins, which leads to the enhancement of the polarization of Fermi gas. Here, we focus on the region of BCS pairing and ignore the Fulde-Ferrell-Larkin-Ovchinnikov state.

Spatial coupling relationships of gas hydrate formation in the Tibetan Plateau

At present, gas hydrates are known to occur in continental high latitude permafrost regions and deep sea sediments. For middle latitude permafrost regions of the Tibetan Plateau, further research is required to ascertain its potential development of gas hydrates. This paper reviewed pertinent literature on gas hydrates in the Tibetan Plateau. Both geological and ge- ographical data are synthesized to reveal the relationship between gas hydrate formation and petroleum geological evo- lution, Plateau uplift, formation of permafrost, and glacial processes. Previous studies indicate that numerous residual basins in the Plateau have been formed by original sedimentary basins accompanied by rapid uplift of the Plateau. Ex- tensive marine Mesozoic hydrocarbon source rocks in these basins could provide rich sources of materials forming gas hydrates in permafrost. Primary hydrocarbon-generating period in the Plateau is from late Jurassic to early Cretaceous, while secondary hydrocarbon generation, regionally or locally, occurs mainly in the Paleogene. Before rapid uplift of the Plateau, oil-gas reservoirs were continuously destroyed and assembled to form new reservoirs due to structural and thermal dynamics, forcing hydrocarbon migration. Since 3.4 Ma B.P., the Plateau has undergone strong uplift and extensive gla- ciation, periglacier processes prevailed, hydrocarbon gas again migrated, and free gas beneath ice sheets within sedi- mentary materials interacted with water, generating gas hydrates which were finally preserved under a cap formed by frozen layers through rapid cooling in the Plateau. Taken as a whole, it can be safely concluded that there is great temporal and spatial coupling relationships between evolution of the Tibetan Plateau and generation of gas hydrates.

Multi-Fidelity Simulation of Gas Turbine Overall Performance by Directly Coupling High-Fidelity Models of Multiple Rotating Components

Multi-fidelity simulations incorporate computational fluid dynamics(CFD) models into a thermodynamic model,enabling the simulation of the overall performance of an entire gas turbine with high-fidelity components.Traditional iterative coupled methods rely on characteristic maps,while fully coupled methods directly incorporate high-fidelity simulations.However,fully coupled methods face challenges in simulating rotating components,including weak convergence and complex implementation.To address these challenges,a fully coupled method with logarithmic transformations was developed to directly integrate high-fidelity CFD models of multiple rotating components.The developed fully coupled method was then applied to evaluate the overall performance of a KJ66 micro gas turbine across various off-design simulations.The developed fully coupled method was also compared with the traditional iterative coupled method.Furthermore,experimental data from ground tests were conducted to verify its effectiveness.The convergence history indicated that the proposed fully coupled method exhibited stable convergence,even under far-off-design simulations.The experimental verification demonstrated that the multi-fidelity simulation with the fully coupled method achieved high accuracy in off-design conditions.Further analysis revealed inherent differences in the coupling methods of CFD models between the developed fully coupled and traditional iterative coupled methods.These inherent differences provide valuable insights for reducing errors between the component-level model and CFD models in different coupling methods.The developed fully coupled method,introducing logarithmic transformations,offers more realistic support for the detailed and optimal design of high-fidelity rotating components within the overall performance platform of gas turbines.

Study on gas–liquid flow characteristics in stirred tank with dual-impeller based on CFD-PBM coupled model

Study on gas–liquid flow in stirred tank with two combinations of dual-impeller(six-bent-bladed turbine(6BT)+six-inclined-blade down-pumping turbine(6 ITD),the six-bent-bladed turbine(6BT)+six-inclinedblade up-pumping turbine(6ITU))was conducted using computational fluid dynamics(CFD)and population balance model(PBM)(CFD-PBM)coupled model.The local bubble size was captured by particle image velocimetry(PIV)measurement.The gas holdup,bubble size distribution and gas–liquid interfacial area were explored at different conditions through numerical simulation.The results showed that the 4 mm bubbles accounted for the largest proportion of 33%at the gas flow rates Q=0.76 m^(3)·h^(-1) and 22%at Q=1.52 m^(3)·h^(-1) for combined impeller of 6BT+6ITU,while the bubbles of 4.7 mm and 5.5 mm were the largest proportion for 6BT+6ITD combination,i.e.25%at Q=0.76 m^(3)·h^(-1) and 22%at Q=1.52 m^(3)·h^(-1),respectively,which indicated that 6BT+6ITU could reduce bubble size effectively and promote gas dispersion.In addition,the gas holdup around impellers was increased obviously with the speed compared with gas flow rate.So it was concluded that 6ITU impeller could be more conductive to the bubble dispersion with more uniform bubble size,which embodied the advantages of 6BT+6ITU combination in gas–liquid mixing.

Kinetics of oxidative coupling of methane:Bridging the gap between comprehension and description

The development of notions about the mechanism of the oxidative coupling of methane （OCM） over oxide catalysts and corresponding progress in its kinetic description are reviewed and discussed. The latter becomes essential at the stage of scaling up and optimization of the process in pilot and industrial reactors. It is demonstrated that the main achievements in the development of kinetic models can be reached by combining the approaches conventionally used in homogeneous gas-phase kinetics and in heterogeneous catalysis. In particular, some important features of the OCM process can be described if several elementary reactions of free radical species （formation and transformation） with surface active sites are included into the detailed scheme of methane oxidation in gas. However, some important features, such as a non-additive character of the reciprocal influence of methane and ethane in the case of their simultaneous presence in the reaction mixture, cannot yet be described and comprehended in the framework of schemes developed so far. Possible ways towards an advanced kinetic model, accounting the main principles of catalyst functioning （redox nature of active sites） and pathways of product formation （via free radicals） are traced.

Analyzing crude oils from the Junggar Basin(NW China) using comprehensive two-dimensional gas chromatography coupled with time-of-flight mass spectrometry(GC×GC-TOFMS)

As a new technology of analyzing crude oils,comprehensive two-dimensional gas chromatography coupled with time-of-flight mass spectrometry(GC×GCTOFMS) has received much research attention.Here we present a case study in the Junggar Basin of NW China.Results show that the hydrocarbons,including saturates and aromatics,were all well-separated without large coelution,which cannot be realized by conventional one-dimensional GC-MS.The GC×GC technique is especially effective for analyzing aromatics and low-to-middlemolecular-weight hydrocarbons,such as diamondoids.The geochemical characteristics of crude oils in the study area were investigated through geochemical parameters extracted by GC×GC-TOFMS,improving upon the understanding obtained by GC-MS.Thus,the work here represents a new successful application of GC×GCTOFMS,showing its broad usefulness in petroleum geochemistry.

Coupled thermo-hydro-mechanical simulation of CO2 enhanced gas recovery with an extended equation of state module for TOUGH2MP-FLAC3D

As one of the most important ways to reduce the greenhouse gas emission,carbon dioxide（CO2）enhanced gas recovery（CO2-EGR） is attractive since the gas recovery can be enhanced simultaneously with CO2sequestration.Based on the existing equation of state（EOS） module of TOUGH2 MP,extEOS7C is developed to calculate the phase partition of H2O-CO2-CH4-NaCl mixtures accurately with consideration of dissolved NaCI and brine properties at high pressure and temperature conditions.Verifications show that it can be applied up to the pressure of 100 MPa and temperature of 150℃.The module was implemented in the linked simulator TOUGH2MP-FLAC3 D for the coupled hydro-mechanical simulations.A simplified three-dimensional（3D）1/4 model（2.2 km×1 km×1 km） which consists of the whole reservoir,caprock and baserock was generated based on the geological conditions of a gas field in the North German Basin.The simulation results show that,under an injection rate of 200,000 t/yr and production rate of 200,000 sm3/d,CO2breakthrough occurred in the case with the initial reservoir pressure of 5 MPa but did not occur in the case of 42 MPa.Under low pressure conditions,the pressure driven horizontal transport is the dominant process;while under high pressure conditions,the density driven vertical flow is dominant.Under the considered conditions,the CO2-EGR caused only small pressure changes.The largest pore pressure increase（2 MPa） and uplift（7 mm） occurred at the caprock bottom induced by only CO2injection.The caprock had still the primary stress state and its integrity was not affected.The formation water salinity and temperature variations of ±20℃ had small influences on the CO2-EGR process.In order to slow down the breakthrough,it is suggested that CO2-EGR should be carried out before the reservoir pressure drops below the critical pressure of CO2.

Experimental and theoretical study on the dynamic effective stress of loaded gassy coal during gas release

In the process of mining coalbed methane(CBM),an unsteady state often arises due to the rapid extraction,release and pressure relief of CBM.In this case,the effective stress of coal changes dynamically,affecting the stability of the gassy coal seam.In this paper,gas release tests of gassy coal under conventional triaxial compression were performed,and the dynamic effective stress(DES)during gas release was obtained indirectly based on a constitutive equation and deformation of coal.The results show that the maximum increases in DES caused by the release of free gas and adsorbed gas under the stress of 1.1 MPa were 0.811 and 5.418 MPa,respectively,which seriously affected the stress state of the coal.During the gas release,the free gas pressure and the adsorbed gas volume were the parameters that directly affected the DES and showed a positive linear relationship with the DES with an intercept of zero.The DES of the coal sample increased exponentially with time,which was determined by the contents of free and adsorbed gas.Based on the experimental results and theoretical analysis,an effective stress model was obtained for loaded gassy coal during gas release.The results of verification indicated accuracy greater than 99%.