Sedimentary Environment and Mineralization of the Black Shale Polymetallic Layer in the Early Cambrian,SW China:Constraints from in situ LA-ICP-MS Analysis of Pyrite

A polymetallic layer is usually developed at the bottom of the early Cambrian black shale in Guizhou Province.The mineral that makes up the polymetallic layer is related to the sedimentary facies.To analyze the differentiation mechanism between polymetallic deposits(Ni-Mo and V),the Zhijin Gezhongwu profile located in the outer shelf and the Sansui Haishan V deposit located in the lower slope are selected to study the in situ sulfur isotopes and trace elements of pyrite.The results show that δ^(34)S values of pyrite vary widely from−7.8‰to 28‰in the Gezhongwu profile,while the δ^(34)S values are relatively uniform(from 27.8‰to 38.4‰)in the Haishan profile.The isotopic S composition is consistent with the transition that occurs in the sedimentary phase from the shelf to the deep sea on the transgressive Yangtze platform;this indicates that the δ^(34)SO_(4)^(2−)values in seawater must be differently distributed in depositional environments.The sulfur in the Ni-Mo layer is produced after the mixing of seawater and hydrothermal fluid,while the V layer mainly originates from seawater.Overall,the Ni-Mo and V deposits have been differentiated primarily on the basis of the combined effect of continental weathering and hydrothermal fluid.

Quantitative characterization of the brittleness of deep shales by integrating mineral content,elastic parameters,in situ stress conditions and logging analysis

Deep shale reservoirs(3500–4500 m)exhibit significantly different stress states than moderately deep shale reservoirs(2000–3500 m).As a result,the brittleness response mechanisms of deep shales are also different.It is urgent to investigate methods to evaluate the brittleness of deep shales to meet the increasingly urgent needs of deep shale gas development.In this paper,the quotient of Young’s modulus divided by Poisson’s ratio based on triaxial compression tests under in situ stress conditions is taken as SSBV(Static Standard Brittleness Value).A new and pragmatic technique is developed to determine the static brittleness index that considers elastic parameters,the mineral content,and the in situ stress conditions(BIEMS).The coefficient of determination between BIEMS and SSBV reaches 0.555 for experimental data and 0.805 for field data.This coefficient is higher than that of other brittleness indices when compared to SSBV.BIEMS can offer detailed insights into shale brittleness under various conditions,including different mineral compositions,depths,and stress states.This technique can provide a solid data-based foundation for the selection of‘sweet spots’for single-well engineering and the comparison of the brittleness of shale gas production layers in different areas.

Depositional Environment and Lithofacies Analyses of Eocene Lacustrine Shale in the Bohai Bay Basin:Insights from Mineralogy and Elemental Geochemistry

The effect of various depositional parameters including paleoclimate,paleosalinity and provenance,on the depositional mechanism of lacustrine shale is very important in reconstructing the depositional environment.The classification of shale lithofacies and the interpretation of shale depositional environment are key features used in shale oil and gas exploration and development activity.The lower 3rd member of the Eocene Shahejie Formation(Es_(3)^(x)shale)was selected for this study,as one of the main prospective intervals for shale oil exploration and development in the intracratonic Bohai Bay Basin.Mineralogically,it is composed of quartz(avg.9.6%),calcite(avg.58.5%),dolomite(avg.7%),pyrite(avg.3.3%)and clay minerals(avg.20%).An advanced methodology(thin-section petrography,total organic carbon and total organic sulfur contents analysis,X-ray diffraction(XRD),X-ray fluorescence(XRF),field-emission scanning electron microscopy(FE-SEM))was adopted to establish shale lithofacies and to interpret the depositional environment in the lacustrine basin.Six different types of lithofacies were recognized,based on mineral composition,total organic carbon(TOC)content and sedimentary structures.Various inorganic geochemical proxies(Rb/Sr,Ca/(Ca+Fe),Ti/Al,Al/Ca,Al/Ti,Zr/Rb)have been used to interpret and screen variations in depositional environmental parameters during the deposition of the Es_(3)^(x)shale.The experimental results indicate that the environment during the deposition of the Es_(3)^(x)shale was warm and humid with heightened salinities,moderate to limited detrital input,higher paleohydrodynamic settings and strong oxygen deficient(reducing)conditions.A comprehensive depositional model of the lacustrine shale was developed.The interpretations deduced from this research work are expected to not only expand the knowledge of shale lithofacies classification for lacustrine fine-grained rocks,but can also offer a theoretical foundation for lacustrine shale oil exploration and development.

Possible link between long-term and short-term water injections and earthquakes in salt mine and shale gas site in Changning,south Sichuan Basin, China

Late at night on 17 June 2019,a magnitude 6.0 earthquake struck Shuanghe Town and its surrounding area in Changning County,Sichuan,China,becoming the largest earthquake recorded within the southern Sichuan Basin.A series of earthquakes with magnitudes up to 5.6 occurred during a short period after the mainshock,and we thus refer to these earthquakes as the Changning M6 earthquake sequence(or swarm).The mainshock was located very close to a salt mine,into which for^3 decades fresh water had been extensively injected through several wells at a depth of 2.7–3 km.It was also near(within^15 km)the epicenter of the 18 December 2018 M5.7 Xingwen earthquake,which is thought to have been induced by shale gas hydraulic fracturing(HF),prompting questions about the possible involvement of industrial activities in the M6 sequence.Following previous studies,this paper focuses on the relationship between injection and seismicity in the Shuanghe salt field and its adjacent Shangluo shale gas block.Except for a period of serious water loss after the start of cross-well injection in 2005–2006,the frequency of earthquakes shows a slightly increasing tendency.Overall,there is a good correlation between the event rate in the Shuanghe area and the loss of injected water.More than 400 M≥3 earthquakes,including 40 M≥4 and 5 M≥5 events,had been observed by the end of August 2019.Meanwhile,in the Shangluo area,seismicity has increased during drilling and HF operations(mostly in vertical wells)since about 2009,and dramatically since the end of 2014,coincident with the start of systematic HF in the area.The event rate shows a progressively increasing background with some fluctuations,paralleling the increase in HF operations.More than 700 M≥3 earthquakes,including 10 M≥4 and 3 M≥5 in spatially and temporally clustered seismic events,are correlated closely with active fracturing platforms.Well-resolved centroid moment tensor results for M≥4 earthquakes were shown to occur at very shallow depths around shale formations with active HF,in agreement with some of the clusters,which occurred within the coverage area of temporary or new permanent monitoring stations and thus have been precisely located.After the Xingwen M5.7 earthquake,seismic activity in the salt well area increased significantly.The Xingwen earthquake may have created a unidirectional rupture to the NNW,with an end point close to the NW-trending fault of the Shuanghe earthquake.Thus,a fault in the Changning anticline might have terminated the fault rupture of the Xingwen earthquake,possibly giving the Xingwen earthquake a role in promoting the Changning M6 event.

Transport model for shale gas well leakage through the surrounding fractured zones of a longwall mine

The environmental risks associated with casing deformation in unconventional(shale)gas wells positioned in abutment pillars of longwall mines is a concern to many in the mining and gas well industry.With the recent interest in shale exploration and the proximity to longwall mining in Southwestern Pennsylvania,the risk to mine workers could be catastrophic as fractures in surrounding strata create pathways for transport of leaked gases.Hence,this research by the National Institute for Occupational Safety and Health(NIOSH)presents an analytical model of the gas transport through fractures in a low permeable stratum.The derived equations are used to conduct parametric studies of specific transport conditions to understand the influence of stratum geology,fracture lengths,and the leaked gas properties on subsurface transport.The results indicated that the prediction that the subsurface gas flux decreases with an increase in fracture length is specifically for a non-gassy stratum.The sub-transport trend could be significantly impacted by the stratum gas generation rate within specific fracture lengths,which emphasized the importance of the stratum geology.These findings provide new insights for improved understanding of subsurface gas transport to ensure mine safety.

Hydrophilic and strength-softening characteristics of calcareous shale in deep mines

To better understand the mechanism of the strength weakening process of soft rocks in deep mines after interacting with water, a self-developed experimental system, Intelligent Testing System for Water Absorption in Deep Soft Rocks （ITSWADSR）, is employed to analyze the hydrophilic behavior of deep calcareous shale sampled from Daqiang coal mine. Experimental results demonstrate that the relation between water absorption and time can be expressed by power functions, and the soakage rate decreases while the soakage increases with time. In order to quantitatively calculate the weight coefficients of the influential factors for water absorbing capacity of rocks, a series of testing methods are adopted, including scanning electron microscope （SEM）, X-ray diffraction and mercury injection test. It is demonstrated that the effective porosity has a positive correlation with the water absorbing capacity of rocks and the contents of illite and illite/smectite. The initial water content presents a negative correlation with the water absorption capacity of rocks. According to the absolute value of weight coefficients of various influential factors, the order of magnitude from high to low is captured： initial water content, illite, illite/smectite formation （S=5%）, and the effective porosity. After water absorption tests, uniaxial compressive strength （UCS） tests were performed on rock specimens allowing a linear relationship between the UCS and the water content of rock to be established, indicating that the strength of calcareous shale decreases linearly with the increasing water content.

Mineralogy and fracture development characteristics of marine shale-gas reservoirs: A case study of Lower Silurian strata in southeastern margin of Sichuan Basin, China

Mineral contents and fractures of shale from well Yuye-1 and outcrops were examined mainly based on systematic description of the cores and outcrops, and data from experimental analyses. The data enabled us to thoroughly explore the mineralogy and developmental features of shale of the Lower Silurian Longmaxi Formation in the study area. The results show that,the Lower Silurian Longmaxi Shale(SLS) in the southeastern margin of Sichuan Basin, China, is primarily characterized by a high content of brittle minerals and a relatively low content of clay minerals. The total content of brittle minerals is approximately 57%,including 27% quartz, 12.2% feldspar, 11.2% carbonate and 2.4% pyrite. The total content of clay minerals reaches 41.6%,composed of illite(23.8%), mixed-layer of illite and smectite(I/S)(10.8%) and chlorite(7.0%). The SLS accommodates the widespread development of various types of fractures, including tectonic fractures, diagenetic fractures, inter-layer fractures and slip fractures. The developmental level of the fracture in the SLS is mainly influenced by faults, lithology, mineral contents and total organic carbon content(TOC) in study area.

Physical Mechanism of Organic Matter-Mineral Interaction in Longmaxi Shale,Sichuan Basin,China

Objective Shale gas is as an important kind of unconventional natural gas,with a great resource potential,and its exploration and development has attracted much attention around the world.Organic matter（OM）pores are a common constituent in shales and form the dominant pore network of many shale gas systems.

Geochemistry, Petrology and Mineralogy of Coal Measure Shales in the Middle Jurassic Yanan Formation from Northeastern Ordos Basin, China: Implications for Shale Gas Accumulation

The Jurassic Yanan Formation is one of the most important coal-producing formations and hydrocarbon source rocks in the Ordos Basin, North China. To evaluate the shale gas potential of the Yanan shale, a total of 48 samples from north Ordos Basin were sampled, and their geochemical, petrological, mineralogical and pore characteristics were investigated. It was found that the shale samples are a suite of early mature source rock. The total organic carbon(TOC) content ranges from 0.33% to 24.12% and the hydrogen index(HI) ranges from 43.31 mg/g to 330.58 mg/g. The relationship between Tmax and HI indicates the organic matter is type Ⅱ-Ⅲ. This conclusion is also supported by the organic petrological examination results, which shows that the kerogen is mainly liptinite and vitrinite. Minerals in the samples are composed mainly of quartz, clay and feldspar, and the clay minerals are composed of prevailing kaolinite, illite/smectite, chlorite and a small amount of illite. Under scanning electron microscope, OM pores in the Yanan shale are scarce except pores come from the kerogen intrinsic texture or clay aggregates within the organic particles. As the weak compaction caused by shallow burial depth, interparticle pores and intraparticle pores are common, the hydrocarbon storage capacity of the Yanan shale was improved. According to evaluation, the Yanan shale is considered as a good shale gas reservoir, but its hydrocarbon potential is more dependent on biogenic and coal-derived gas as the thermogenic gas is limited by the lower thermal maturity.

Variation in Pore Space and Structure of Organic-rich Oilprone Shales from a Non-marine Basin:Constraints from Organic Matter and Minerals

Organic matter(OM)and minerals are major particle components of lacustrine shales.The influence of OM and minerals on pore space and structure in organic-rich oil-prone shales containing a large range of total organic carbon(TOC)contents is poorly understood.In this study,we investigated the variation in pore space and structure of low mature lacustrine shales in the Songliao Basin(NE China),based on a study of the mineralogy,petrography,geochemistry,and geophysical properties of shales.Different pore types make markedly different contributions to the mineral surface area(MSA)and pore volume(PV)of the shales.There exists a negative correlation between MSA/PV and TOC in mesopores(r^(2)=0.75/0.65)and macropores(r^(2)=0.74/0.68),and a positive correlation in micropores(r^(2)=0.59/0.64),which are associated with the variation of mineral and TOC contents.A positive relationship between the throat/pore ratio and TOC(r^(2)=0.82)shows an increase in throat radius and decrease in pore radius with increasing TOC content.This relationship is supported by the reduction in mean pore diameter(MPD)for large pores and increase in MPD for small pores.These variations are related to the decreased pores by quartz plus feldspar(Q+F)content,increased throats by clay minerals,an d enhanced pore-fill by OM.We propose that the variation of OM and minerals is a key control on the pore space and structure of low mature organic-rich oil-prone shales.

Uraniferous Black Shale and Related Uranium Mineralization Features in South China

Black shales are marine sediments with argillaceous, silty and siliceous compositions and high contents of organic materials, disseminated pyrite and uranium. Uraniferous black shale has uranium content of more than 20 ppm.

Wettability of different clay mineral surfaces in shale:Implications from molecular dynamics simulations

Shale contains a lot of clay minerals. Clay minerals mainly exist in nano- and micro-meter sized particles, and the pore structure is complex, which leads to its extremely complex wettability. The surface wettability of clay minerals significantly affects the oil and gas-bearing capacity of shale reservoirs. Therefore, studying the wettability of common clay minerals in shale at the nanoscale is of great significance for shale hydrocarbon exploration and development. In this study, the wetting behavior of water in n-hexane and toluene on different clay mineral surfaces at the nanoscale was systematically studied using Molecular dynamics (MD) simulation. And the influencing factors of wettability were analyzed. Through the analysis of the morphological changes of water, relative concentration of water, RDF and interaction energy, it is concluded that the following order of water wettability on the surfaces of clay minerals: montmorillonite > chlorite > kaolinite > illite. Through the analysis of interaction energy, it is concluded that the hydrophilicity of four clay minerals is stronger than that of lipophilicity. And the main interactions between water and oil and the mineral surfaces were van der Waals force and electrostatic force. In addition, the temperature, liquid hydrocarbon type, and mineralization of water affected the wettability of clay minerals. The concentration of water on the surfaces of montmorillonite, kaolinite, and illite decreased with increasing temperature, and the water wettability decreased. At 298 K, the hydrophilicity of the surfaces of the clay minerals in toluene follows the order montmorillonite > chlorite > kaolinite > illite. The higher the NaHCO3 concentration in water, the weaker the wettability of the clay mineral surfaces to water. By comparing the previous experimental results with the MD simulation results, similar wetting characteristics were obtained, and the reliability of the simulation results was verified. MD simulation was used to explore the water wetting of the surfaces of four clay minerals in a shale reservoir from the micro level. This makes up for the lack of experimental means for clarifying the flow and production mechanisms of shale oil and gas and effectively improves the evaluation technology of shale.

Influence of Bedding and Mineral Composition on Mechanical Properties and Its Implication for Hydraulic Fracturing of Shale Oil Reservoirs

The premise of hydraulic fracturing is to have an accurate and detailed understanding of the rock mechanical properties and fracture propagation law of shale reservoirs. In this paper,a comprehensive evaluation of the mechanical properties of the shale oil reservoir in the south of Songliao Basin is carried out. Based on the experiments and the in-situ stress analysis, the fracture propagation law of three types of shale reservoirs is obtained,and the suggestions for fracturing are put forward. The results have shown that the fracture propagation of pure shale and low mature reservoir is easy to open along the bedding plane under compression loading,which is greatly influenced by the bedding. Sand-bearing shale is slightly better,the fractures of which are not easy to open along the bedding plane. The mechanical experimental results show that all the samples have the characteristics of low compressive strength,low Young’s modulus and strong anisotropy,indicating that the shale oil reservoir is certain plastic,which is related to its high clay mineral content and controlled by the bedding development. Compared with pure shale and low mature shale,the sandbearing shale has less clay content and less developed bedding,which maybe the main reason for its slightly better brittleness. Overall,the expansion of hydraulic fracture is controlled by in-situ stress and bedding. Because of the development of bedding,it is easy to form horizontal fractures. Thus it is not suitable for horizontal well fracturing.Because of the high content of clay minerals,the applicability of conventional slick hydraulic fracturing fluid is poor. It is suggested to use vertical well or directional well to carry out volume fracturing. In this way,the effect of bedding can be effectively used to open and connect the bedding and form a larger fracture network.

Fungal diversity in major oil-shale mines in China

As an insufficiently utilized energy resource,oil shale is conducive to the formation of characteristic microbial communities due to its special geological origins.However,little is known about fungal diversity in oil shale.Polymerase chain reaction cloning was used to construct the fungal ribosomal deoxyribonucleic acid internal transcribed spacer（r DNA ITS）clone libraries of Huadian Mine in Jilin Province,Maoming Mine in Guangdong Province,and Fushun Mine in Liaoning Province.Pure culture and molecular identification were applied for the isolation of cultivable fungi in fresh oil shale of each mine.Results of clone libraries indicated that each mine had over 50% Ascomycota（58.4%–98.9%）and 1.1%–13.5%unidentified fungi.Fushun Mine and Huadian Mine had 5.9% and 28.1% Basidiomycota,respectively.Huadian Mine showed the highest fungal diversity,followed by Fushun Mine and Maoming Mine.Jaccard indexes showed that the similarities between any two of three fungal communities at the genus level were very low,indicating that fungi in each mine developed independently during the long geological adaptation and formed a community composition fitting the environment.In the fresh oil-shale samples of the three mines,cultivable fungal phyla were consistent with the results of clone libraries.Fifteen genera and several unidentified fungi were identified as Ascomycota and Basidiomycota using pure culture.Penicillium was the only genus found in all three mines.These findings contributed to gaining a clear understanding of current fungal resources in major oil-shale mines in China and provided useful information for relevant studies on isolation of indigenous fungi carrying functional genes from oil shale.

The new multistage water adsorption model of Longmaxi Formation shale considering the spatial configuration relationship between organic matter and clay minerals

The water adsorption by shale significantly affects shale gas content and its seepage capacity.However,the mechanism of water adsorption by shale is still unclear due to its strong heterogeneity and complicated pore structure.The relationship between the adsorbed water content at different relative humidities(RHs)and shale compositions,as well as shale pore structure and the spatial configuration relationship between organic matter(OM)and clay minerals,was investigated to clarify the controlling factors and mechanisms of water adsorption by Longmaxi Formation shale from the Southern Sichuan Basin in China.Consequently,the water adsorption process could be generally divided into three different stages from 0%RH to 99%RH.Furthermore,the Johnston’s clay mine ral interlayer pore structure model(JCM),the Freundlich model(FM)and the Dubinin-Astakhov model(DAM)were tested to fit the three water adsorption stages from low RH to high RH,respectively.The fitting results of the JCM and FM at lower RHs were far from good,while the fitting results of DAM at higher RHs were acceptable.Accordingly,two revised models(LRHM and MRHM)considering the spatial configuration relationship between OM and clay minerals were proposed for the two stages with lower RHs,and performed better fitting results indicating the pronounced effect of the spatial configuration relationship between OM and clay minerals on the water adsorption process of Longmaxi Formation shale.The outcomes of this study will contribute to clarifying the water distribution characteristics in the pore network of shale samples with variable water contents.

Mineralogy and Thermal Analysis of Natural Pozzolana Opal Shale with Nano-pores

Thermodynamic stability, microvoid distribution and phases transformation of natural pozzolana opal shale（POS） were studied systematically in this work. XRD analysis showed that opal-CT, including microcrystal cristobalite and tridymite, is a major component of POS. DTA and FT-IR indicated that there were many hydroxyl groups and acid sites on the surface of amorphous SiO2 materials. FE-SEM analysis exhibited amorphous SiO2 particles（opal-A） covering over stacking sequences microcrystal cristobalite and tridymite. Meanwhile, MIP analysis demonstrated that porosity and pore size distribution of POS remained uniform below 600 ℃. Because stable porous microstructure is a key factor in improving photocatalyst activity, POS is suited to preparing highly active supported.

A semi-analytical model for coupled flow in stress-sensitive multi-scale shale reservoirs with fractal characteristics

A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes due to stress sensitivity, which plays a crucial role in controlling pressure propagation and oil flow. This paper proposes a multi-scale coupled flow mathematical model of matrix nanopores, induced fractures, and hydraulic fractures. In this model, the micro-scale effects of shale oil flow in fractal nanopores, fractal induced fracture network, and stress sensitivity of multi-scale media are considered. We solved the model iteratively using Pedrosa transform, semi-analytic Segmented Bessel function, Laplace transform. The results of this model exhibit good agreement with the numerical solution and field production data, confirming the high accuracy of the model. As well, the influence of stress sensitivity on permeability, pressure and production is analyzed. It is shown that the permeability and production decrease significantly when induced fractures are weakly supported. Closed induced fractures can inhibit interporosity flow in the stimulated reservoir volume (SRV). It has been shown in sensitivity analysis that hydraulic fractures are beneficial to early production, and induced fractures in SRV are beneficial to middle production. The model can characterize multi-scale flow characteristics of shale oil, providing theoretical guidance for rapid productivity evaluation.

Effect of fracture fluid flowback on shale microfractures using CT scanning

The field data of shale fracturing demonstrate that the flowback performance of fracturing fluid is different from that of conventional reservoirs,where the flowback rate of shale fracturing fluid is lower than that of conventional reservoirs.At the early stage of flowback,there is no single-phase flow of the liquid phase in shale,but rather a gas-water two-phase flow,such that the single-phase flow model for tight oil and gas reservoirs is not applicable.In this study,pores and microfractures are extracted based on the experimental results of computed tomography(CT)scanning,and a spatial model of microfractures is established.Then,the influence of rough microfracture surfaces on the flow is corrected using the modified cubic law,which was modified by introducing the average deviation of the microfracture height as a roughness factor to consider the influence of microfracture surface roughness.The flow in the fracture network is simulated using the modified cubic law and the lattice Boltzmann method(LBM).The results obtained demonstrate that most of the fracturing fluid is retained in the shale microfractures,which explains the low fracturing fluid flowback rate in shale hydraulic fracturing.

Identification and evaluation of shale oil micromigration and its petroleum geological significance

Taking the Lower Permian Fengcheng Formation shale in Mahu Sag of Junggar Basin,NW China,as an example,core observation,test analysis,geological analysis and numerical simulation were applied to identify the shale oil micro-migration phenomenon.The hydrocarbon micro-migration in shale oil was quantitatively evaluated and verified by a self-created hydrocarbon expulsion potential method,and the petroleum geological significance of shale oil micro-migration evaluation was determined.Results show that significant micro-migration can be recognized between the organic-rich lamina and organic-poor lamina.The organic-rich lamina has strong hydrocarbon generation ability.The heavy components of hydrocarbon preferentially retained by kerogen swelling or adsorption,while the light components of hydrocarbon were migrated and accumulated to the interbedded felsic or carbonate organic-poor laminae as free oil.About 69% of the Fengcheng Formation shale samples in Well MY1 exhibit hydrocarbon charging phenomenon,while 31% of those exhibit hydrocarbon expulsion phenomenon.The reliability of the micro-migration evaluation results was verified by combining the group components based on the geochromatography effect,two-dimension nuclear magnetic resonance analysis,and the geochemical behavior of inorganic manganese elements in the process of hydrocarbon migration.Micro-migration is a bridge connecting the hydrocarbon accumulation elements in shale formations,which reflects the whole process of shale oil generation,expulsion and accumulation,and controls the content and composition of shale oil.The identification and evaluation of shale oil micro-migration will provide new perspectives for dynamically differential enrichment mechanism of shale oil and establishing a“multi-peak model in oil generation”of shale.

Analysis of CH_(4) and H_(2) Adsorption on Heterogeneous Shale Surfaces Using aMolecular Dynamics Approach

Determining the adsorption of shale gas on complex surfaces remains a challenge in molecular simulation studies.Difficulties essentially stem from the need to create a realistic shale structure model in terms of mineral heterogeneityand multiplicity.Moreover,precise characterization of the competitive adsorption of hydrogen andmethane in shale generally requires the experimental determination of the related adsorptive capacity.In thisstudy,the adsorption of adsorbates,methane(CH_(4)),and hydrogen(H_(2))on heterogeneous shale surface modelsof Kaolinite,Orthoclase,Muscovite,Mica,C_(60),and Butane has been simulated in the frame of a moleculardynamic’s numerical technique.The results show that these behaviors are influenced by pressure and potentialenergy.On increasing the pressure from 500 to 2000 psi,the sorption effect for CH_(4)significantly increasesbut shows a decline at a certain stage(if compared to H_(2)).The research findings also indicate that raw shalehas a higher capacity to adsorb CH_(4)compared to hydrogen.However,in shale,this difference is negligible.