Geologic characteristics,exploration and production progress of shale oil and gas in the United States:An overview

We present a systematic summary of the geological characteristics,exploration and development history and current state of shale oil and gas in the United States.The hydrocarbon-rich shales in the major shale basins of the United States are mainly developed in six geological periods:Middle Ordovician,Middle-Late Devonian,Early Carboniferous(Middle-Late Mississippi),Early Permian,Late Jurassic,and Late Cretaceous(Cenomanian-Turonian).Depositional environments for these shales include intra-cratonic basins,foreland basins,and passive continental margins.Paleozoic hydrocarbon-rich shales are mainly developed in six basins,including the Appalachian Basin(Utica and Marcellus shales),Anadarko Basin(Woodford Shale),Williston Basin(Bakken Shale),Arkoma Basin(Fayetteville Shale),Fort Worth Basin(Barnett Shale),and the Wolfcamp and Leonardian Spraberry/Bone Springs shale plays of the Permian Basin.The Mesozoic hydrocarbon-rich shales are mainly developed on the margins of the Gulf of Mexico Basin(Haynesville and Eagle Ford)or in various Rocky Mountain basins(Niobrara Formation,mainly in the Denver and Powder River basins).The detailed analysis of shale plays reveals that the shales are different in facies and mineral components,and"shale reservoirs"are often not shale at all.The United States is abundant in shale oil and gas,with the in-place resources exceeding 0.246×10^(12)t and 290×10^(12)m^(3),respectively.Before the emergence of horizontal well hydraulic fracturing technology to kick off the"shale revolution",the United States had experienced two decades of exploration and production practices,as well as theory and technology development.In 2007-2023,shale oil and gas production in the United States increased from approximately 11.2×10^(4)tons of oil equivalent per day(toe/d)to over 300.0×10^(4)toe/d.In 2017,the shale oil and gas production exceeded the conventional oil and gas production in the country.In 2023,the contribution from shale plays to the total U.S.oil and gas production remained above 60%.The development of shale oil and gas has largely been driven by improvements in drilling and completion technologies,with much of the recent effort focused on“cube development”or“co-development”.Other efforts to improve productivity and efficiency include refracturing,enhanced oil recovery,and drilling of“U-shaped”wells.Given the significant resources base and continued technological improvements,shale oil and gas production will continue to contribute significant volumes to total U.S.hydrocarbon production.

Experimental study on the size effect on the equation of state of concretes under shock loading

Adopting the classical theory of hydrocodes,the constitutive relations of concretes are separated into an equation of state(EoS)which describes the volumetric behavior of concrete material and a strength model which depicts the shear properties of concrete.The experiments on the EoS of concrete is always challenging due to the technical difficulties and equipment limitations,especially for the specimen size effect on the EoS.Although some researchers investigate the shock properties of concretes by fly-plate impact tests,the specimens used in their tests are usually in one size.In this paper,the fly-plate impact tests on concrete specimens with different sizes are performed to investigate the size effect on the shock properties of concrete materials.The mechanical background of the size effect on the shock properties are revealed,which is related to the lateral rarefaction effect and the deviatoric stress produced in the specimen.According to the tests results,the modified EoS considering the size effect on the shock properties of concrete are proposed,which the bulk modulus of concrete is unpredicted by up to 20% if size effects are not accounted for.

The pressure compensation technology of deep-sea sampling based on the real gas state equation

Compressed gas is usually used for the pressure compensation of the deep-sea pressure-maintaining sampler.The pressure and volume of the recovered fluid sample are highly related to the precharged gas. To better understand the behavior of the gas under high pressure, we present a new real gas state equation based on the compression factor Z which was derived from experimental data. Then theoretical calculation method of the pressure and volume of the sample was introduced based on this empirical gas state equation. Finally, the proposed calculation method was well verified by the high-pressure vessel experiment of the sampler under 115 MPa.

A New Alpha-function for the Peng-Robinson Equation of State： Application to Natural Gas

Cubic equations of state(EOSs) are simple and easy at calculation. One way of improving the accuracy of a cubic EOS is through the modification of temperature-dependent energy parameter by using alpha-function.The industrial applications of natural gas are very wide and as a result, prediction of thermodynamic properties and phase behavior of natural gas is an important part of design for such processes. In this work we develop a newα-function for the Peng-Robinson(PR) EOS with the parameters optimized especially for natural gas components.The parameters are generalized as a linear function of acentric factor. The results are compared to the predictions from original PR EOS and other α-functions in literature. It is shown that the new α-function presents a good accuracy with the average deviation of 1.42% for natural gas components.

PREDICTION OF THE PHASE BEHAVIOR OF GAS CONDENSATES USING PATEL-TEJA EQUATION OF STATE AND GAMMA DISTRIBUTION FUNCTION

The three-parameter Petal-Teja equation of state coupled with a characterization proceduref0r C7+-fraction based on gamma distribution function was employed to predict the phase behaviorof eight gas condensates.The lumping of the subdivided single carbon number（SCN）hydrocarbons inthe plus-fraction and the choice of empirical correlations for calculating the critical properties andacentric factor of SCN hydrocarbons were discussed.

Thermal equation of state for lattice Boltzmann gases

The Galilean invariance and the induced thermo-hydrodynamics of the lattice Boltzmann Bhatnagar-Gross-Krook model are proposed together with their rigorous theoretical background. From the viewpoint of group invariance, recovering the Galilean invariance for the isothermal lattice Boltzmann Bhatnagar-Gross-Krook equation （LBGKE） induces a new natural thermal-dynamical system, which is compatible with the elementary statistical thermodynamics.

A Modified Peng-Robinson Equation State for Prediction of Gas Adsorption Isotherm

This research proposes a modified two-dimensional Peng-Robinson equation model to predict adsorption isotherm in adsorbate-adsorbent systems. The parameters of the proposed model are calculated by using the optimization of experimental data for the different single gas adsorption systems at various temperatures. The experimental adsorption equilibrium data of adsorbate-adsorbent systems was compared with the calculated results in our proposed model and the two-dimensional Hill-deBoer equation model. The proposed model as indicated in the results shows a better prediction of the experimental results compared with two others.

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.

An insight into the role of the association equations of states in gas hydrate modeling:a review

Encouraged by the wide spectrum of novel applications of gas hydrates,e.g.,energy recovery,gas separation,gas storage,gas transportation,water desalination,and hydrogen hydrate as a green energy resource,as well as CO2 capturing,many scientists have focused their attention on investigating this important phenomenon.Of course,from an engineering viewpoint,the mathematical modeling of gas hydrates is of paramount importance,as anticipation of gas hydrate stability conditions is effective in the design and control of industrial processes.Overall,the thermodynamic modeling of gas hydrate can be tackled as an equilibration of three phases,i.e.,liquid,gas,and solid hydrate.The inseparable component in all hydrate systems,water,is highly polar and non-ideal,necessitating the use of more advanced equation of states(EoSs) that take into account more intermolecular forces for thermodynamic modeling of these systems.Motivated by the ever-increasing number of publications on this topic,this study aims to review the application of associating EoSs for the thermodynamic modeling of gas hydrates.Three most important hydrate-based models available in the literature including the van der Waals-Platteeuw(vdW-P) model,Chen-Guo model,and Klauda-Sandler model coupled with and SAFT EoSs were investigated and compared with cubic EoSs.It was concluded that the CPA and SAFT EoSs gave very accurate results for hydrate systems as they take into account the association interactions,which are very crucial in gas hydrate systems in which water,methanol,glycols,and other types of associating compounds are available.Moreover,it was concluded that the CPA EoS is easier to use than the SAFT-type EoSs and our suggestion for the gas hydrate systems is the CPA EoS.

THE EXISTENCE OF A BOUNDED INVARIANT REGION FOR COMPRESSIBLE EULER EQUATIONS IN DIFFERENT GAS STATES

In this article,by the mean-integral of the conserved quantity,we prove that the one-dimensional non-isentropic gas dynamic equations in an ideal gas state do not possess a bounded invariant region.Moreover,we obtain a necessary condition on the state equations for the existence of an invariant region for a non-isentropic process.Finally,we provide a mat hematical example showing that with a special state equation,a bounded invariant region for the non-isentropic process may exist.

CALCULATIONS OF THE PHASE EQUILIBRIA AND CRITICAL CURVES WITH THE EQUATION OF STATE FOR BINARY MIXTURES

A rational equation of state of the perturbation type with a repulsion and attraction term has been applied to reproduce critical curves of six different binary systems up to high temperatures and pressures. A square well potential for intermolecular interaction is used. With pairwise combination rules for these potentials three adjustable parameters are needed. The experimental critical point and phase equilibrium data are compared with the values predicted using the equation of state. Good agreement is obtained for the analysis of the critical pressure composition data and molar volumes.

Tracing of natural gas migration by light hydrocarbons:A case study of the Dongsheng gas field in the Ordos Basin,NW China

Based on the analysis of light hydrocarbon compositions of natural gas and regional comparison in combination with the chemical components and carbon isotopic compositions of methane,the indication of geochemical characteristics of light hydrocarbons on the migration features,dissolution and escape of natural gas from the Dongsheng gas field in the Ordos Basin is revealed,and the effect of migration on specific light hydrocarbon indexes is further discussed.The study indicates that,natural gas from the Lower Shihezi Formation(Pix)in the Dongsheng gas field displays higher iso-C5-7contents than n-C5-7contents,and the C6-7light hydrocarbons are composed of paraffins with extremely low aromatic contents(<0.4%),whereas the C7light hydrocarbons are dominated by methylcyclohexane,suggesting the characteristics of coal-derived gas with the influence by secondary alterations such as dissolution.The natural gas from the Dongsheng gas field has experienced free-phase migration from south to north and different degrees of dissolution after charging,and the gas in the Shiguhao area to the north of the Borjianghaizi fault has experienced apparent diffusion loss after accumulation.Long-distance migration in free phase results in the decrease of the relative contents of the methylcyclohexane in C7 light hydrocarbons and the toluene/n-heptane ratio,as well as the increase of the n-heptane/methylcyclohexane ratio and heptane values.The dissolution causes the increase of isoheptane values of the light hydrocarbons,whereas the diffusion loss of natural gas in the Shiguhao area results in the increase of n-C5-7contents compared to the iso-C5-7contents.

Normalized Ground State Solutions of Nonlinear Choquard Equations with Nonconstant Potential

In this paper,we mainly focus on the following Choquard equation-{△u-V(x)(I_(a*)|u|^(p))|u|^(p-2)u=λu,x∈R^(N),u∈H^(1)(R^(N))where N≥1,λ∈R will arise as a Lagrange multiplier,0<a<N and N+a/N<p<N+a+2/N Under appropriate hypotheses on V(x),we prove that the above Choquard equation has a normalized ground state solution by utilizing variational methods.

Desired Dynamic Equation for Primary Frequency Modulation Control of Gas Turbines

Gas turbines play core roles in clean energy supply and the construction of comprehensive energy systems.The control performance of primary frequency modulation of gas turbines has a great impact on the frequency control of the power grid.However,there are some control difficulties in the primary frequency modulation control of gas turbines,such as the coupling effect of the fuel control loop and speed control loop,slow tracking speed,and so on.To relieve the abovementioned difficulties,a control strategy based on the desired dynamic equation proportional integral(DDE-PI)is proposed in this paper.Based on the parameter stability region,a parameter tuning procedure is summarized.Simulation is carried out to address the ease of use and simplicity of the proposed tuning method.Finally,DDE-PI is applied to the primary frequency modulation system of an MS6001B heavy-duty gas turbine.The simulation results indicate that the gas turbine with the proposed strategy can obtain the best control performance with a strong ability to deal with system uncertainties.The proposed method shows good engineering application potential.

EQUATION OF STATE BASED HYDRATE MODEL FOR NATURAL GAS SYSTEMS CONTAINING BRINE AND POLAR INHIBITOR

A new thermodynamic model for gas hydrates was established by combining the modified Patel-Teja equation of state proposed for aqueous electrolyte systems and the simplified Holder -John multi -shell hydrate model. The new hydrate model is capable of predicting the hydrate formation/dissociation conditions of natural gas systems containing pure water/formation water (brine) and polar inhibitor without using activity coefficient model. Extensive test results indicate very encouraging results.

Control Factors and Diversities of Phase State of Oil and Gas Pools in the Kuqa Petroleum System

Based on the analysis of the hydrocarbon geochemical characteristics in the Kuqa petroleum system of the Tarim Basin, this study discusses the causes and controlling factors of the phase diversities and their differences in geochemical features. According to the characteristics and differences in oil and gas phase, the petroleum system can be divided into five categories： oil reservoir, wet gas reservoir, condensate gas-rich reservoir, condensate gas-poor reservoir and dry gas reservoir. The causes for the diversities in oil and gas phases include diversities of the sources of parent material, maturity of natural gas and the process of hydrocarbon accumulation of different hydrocarbon phases. On the whole, the Jurassic and Triassic terrestrial source rocks are the main sources for the hydrocarbon in the Kuqa Depression. The small differences in parent material may cause diversities in oil and gas amount, but the impact is small. The differences in oil and gas phase are mainly affected by maturity and the accumulation process, which closely relates with each other. Oil and gas at different thermal evolution stage can be captured in different accumulation process.

Strangeness to increase the density of finite nuclear systems in constraining the high-density nuclear equation of state

As the high-density nuclear equation of state(EOS) is not very well constrained, we suggest that the structural properties from the finite systems can be used to extract a more accurate constraint. By including the strangeness degrees of freedom, the hyperon or anti-kaon, the finite systems can then have a rather high-density core which is relevant to the nuclear EOS at high densities directly. It is found that the density dependence of the symmetry energy is sensitive to the properties of multi-K hypernuclei, while the high-density EOS of symmetric matter correlates sensitively to the properties of kaonic nuclei.

Basic quantities of the equation of state in isospin asymmetric nuclear matter

Based on the Hugenholtz-Van Hove theorem six basic quantities of the EoS in isospin asymmetric nuclear matter are expressed in terms of the nucleon kinetic energy t(k),the isospin symmetric and asymmetric parts of the single-nucleon potentials U_(0)(ρ,k)and U_(sym,i)(ρ,k).The six basic quantities include the quadratic symmetry energy E_(sym,2)(ρ),the quartic symmetry energy E_(sym,4)(ρ),their corresponding density slopes L_(2)(ρ)and L_(4)(ρ),and the incompressibility coefficients K_(2)(ρ)and K_(4)(ρ).By using four types of well-known effective nucleon-nucleon interaction models,namely the BGBD,MDI,Skyrme,and Gogny forces,the density-and isospin-dependent properties of these basic quantities are systematically calculated and their values at the saturation density q_(0)are explicitly given.The contributions to these quantities from t(k)U_(0)(ρ,k),and U_(sym,i)(ρ,k)are also analyzed at the norma nuclear density q_(0).It is clearly shown that the first-order asymmetric term U_(sym,1)(ρ,k)(also known as the symmetry potential in the Lane potential)plays a vital role in determining the density dependence of the quadratic symmetry energy E_(sym,2)(ρ).It is also shown that the contributions from the high-order asymmetric parts of the single-nucleon potentials(U_(sym,i)(ρ,k)with i>1)cannot be neglected in the calculations of the other five basic quantities Moreover,by analyzing the properties of asymmetric nuclear matter at the exact saturation densityρ_(sat)(δ),the corresponding quadratic incompressibility coefficient is found to have a simple empirical relation K_(sat,2)=K_(2)(ρ_(0))-4.14L_(2)(ρ_(0))

Representation of Phase Behavior of Ionic Liquids Using the Equation of State for Square-well Chain Fluids with Variable Range

An equation of state （EOS） for square-well chain fluids with variable range （SWCF-VR） developed based on statistical mechanics for chemical association was employed for the calculations of pressure-volume-temperature （pVT） and phase equilibrium of pure ionic liquids （ILs） and their mixtures. The new molecular parameters for 23 ILs were obtained by fitting their experimental density data over a wide temperature and pressure ranges. The mo- lecular parameters of ILs composed of homologous organic cation and an identical anion such as [Cxmim][NTf2] are good linear with respect to their molecular weight, indicating that the molecular parameters of homologous substances, subsequently p VT and vapor-liquid equilibria vapor-liquid equilibria （VLE） can be predicted using the generalized parameter when no experimental data were available. The new set of parameters were satisfactorily used for calculations of the property of solvent and ILs mixture and the solubility of gas in various ILs at low pressure only using one temperature-independent binary interaction parameter.

New Approach for Solving Master Equations in Quantum Optics and Quantum Statistics by Virtue of Thermo-Entangled State Representation

By introducing a fictitious mode to be a counterpart mode of the system mode under review we introduce the entangled state representation （η｜, which can arrange master equations of density operators p（t） in quantum statistics as state-vector evolution equations due to the elegant properties of （η｜. In this way many master equations （respectively describing damping oscillator, laser, phase sensitive, and phase diffusion processes with different initial density operators） can be concisely solved. Specially, for a damping process characteristic of the decay constant k we find that the matrix element of p（t） at time t in 〈η｜ representation is proportional to that of the initial po in the decayed entangled state （ηe^-kt｜ representation, accompanying with a Gaussian damping factor. Thus we have a new insight about the nature of the dissipative process. We also set up the so-called thermo-entangled state representation of density operators, ρ = f（d^2η/π）（η｜ρ〉D（η）, which is different from all the previous known representations.