Numerical Simulation on Production Trials by Using Depressurization for Typical Marine Hydrate Reservoirs:Well Type and Formation Dip

Natural gas hydrate has huge reserves and is widely distributed in marine environment.Its commercial development is of great significance for alleviating the contradiction between energy supply and demand.As an efficient research method,numerical simulation can provide valuable insights for the design and optimization of hydrate development.However,most of the current production models simplify the reservoir as a two-dimensional(2D)horizontal layered model,often ignoring the impact of formation dip angle.To improve the accuracy of production prediction and provide theoretical support for the optimization of production well design,two three-dimensional(3D)geological models with different dip angles based on the geological data from two typical sites are constructed.The vertical well,horizontal well and multilateral wells are deployed in these reservoirs with different permeabilities to perform production trial,and the sensitivity analysis of dip angles is also carried out.The short-term production behaviors in high and low permeability reservoirs with different dip angles are exhibited.The simulation results show that 1)the gas and water production behaviors for different well types in the two typical reservoirs show obviously different variation laws when the short-term depressurization is conducted in the inclined formation;2)the inclined formation will reduce the gas production and increase the water extraction,and the phenomena becomes pronounced as the dip angle increases,particularly in the low-permeability reservoirs;3)and the impact of formation dip on hydrate recovery does not change significantly with the variation of well type.

Saturation Estimation with Complex Electrical Conductivity for Hydrate-Bearing Clayey Sediments:An Experimental Study

Clays have considerable influence on the electrical properties of hydrate-bearing sediments.It is desirable to understand the electrical properties of hydrate-bearing clayey sediments and to build hydrate saturation(S_(h))models for reservoir evaluation and monitoring.The electrical properties of tetrahydrofuran-hydrate-bearing sediments with montmorillonite are characterized by complex conductivity at frequencies from 0.01 Hz to 1 kHz.The effects of clay and Sh on the complex conductivity were analyzed.A decrease and increase in electrical conductance result from the clay-swelling-induced blockage and ion migration in the electrical double layer(EDL),respectively.The quadrature conductivity increases with the clay content up to 10%because of the increased surface site density of counterions in EDL.Both the in-phase conductivity and quadrature conductivity decrease consistently with increasing Sh from 0.50 to 0.90.Three sets of models for Sh evaluation were developed.The model based on the Simandoux equation outperforms Archie’s formula,with a root-mean-square error(E_(RMS))of 1.8%and 3.9%,respectively,highlighting the clay effects on the in-phase conductivity.The fre-quency effect correlations based on in-phase and quadrature conductivities exhibit inferior performance(E_(RMS)=11.6%and 13.2%,re-spectively)due to the challenge of choosing an appropriate pair of frequencies and intrinsic uncertainties from two measurements.The second-order Cole-Cole formula can be used to fit the complex-conductivity spectra.One pair of inverted Cole-Cole parameters,i.e.,characteristic time and chargeability,is employed to predict S_(h) with an E_(RMS) of 5.05%and 9.05%,respectively.

Deep learning CNN-APSO-LSSVM hybrid fusion model for feature optimization and gas-bearing prediction

Conventional machine learning(CML)methods have been successfully applied for gas reservoir prediction.Their prediction accuracy largely depends on the quality of the sample data;therefore,feature optimization of the input samples is particularly important.Commonly used feature optimization methods increase the interpretability of gas reservoirs;however,their steps are cumbersome,and the selected features cannot sufficiently guide CML models to mine the intrinsic features of sample data efficiently.In contrast to CML methods,deep learning(DL)methods can directly extract the important features of targets from raw data.Therefore,this study proposes a feature optimization and gas-bearing prediction method based on a hybrid fusion model that combines a convolutional neural network(CNN)and an adaptive particle swarm optimization-least squares support vector machine(APSO-LSSVM).This model adopts an end-to-end algorithm structure to directly extract features from sensitive multicomponent seismic attributes,considerably simplifying the feature optimization.A CNN was used for feature optimization to highlight sensitive gas reservoir information.APSO-LSSVM was used to fully learn the relationship between the features extracted by the CNN to obtain the prediction results.The constructed hybrid fusion model improves gas-bearing prediction accuracy through two processes of feature optimization and intelligent prediction,giving full play to the advantages of DL and CML methods.The prediction results obtained are better than those of a single CNN model or APSO-LSSVM model.In the feature optimization process of multicomponent seismic attribute data,CNN has demonstrated better gas reservoir feature extraction capabilities than commonly used attribute optimization methods.In the prediction process,the APSO-LSSVM model can learn the gas reservoir characteristics better than the LSSVM model and has a higher prediction accuracy.The constructed CNN-APSO-LSSVM model had lower errors and a better fit on the test dataset than the other individual models.This method proves the effectiveness of DL technology for the feature extraction of gas reservoirs and provides a feasible way to combine DL and CML technologies to predict gas reservoirs.

Auto-Detection Method Using Convolution Neural Network for Bottom-Simulating Reflectors

In studies on gas hydrate,bottom-simulating reflectors(BSR)are used to determine the potential hydrate-bearing sedimentary layers.Usually,BSR detection is performed manually by experienced interpreters.Therefore,a method for implementing an auto-matic BSR detection process should be established.In this study,we develop a novel architecture for BSR characterization using the convolutional neural network(CNN)technique.We propose the use of Stokes’transform(ST)to obtain a time-frequency spectrum for the input of CNN.ST fully uses the frequency content of the seismic data,and a part of the 3D seismic data collected from the Blake Ridge is utilized to train the CNN.Synthetic seismic records with variable signal-to-noise ratios(SNR),as well as Blake Ridge seismic data,were used to validate the detection effect of the CNN.Results show that the CNN trained by this method exhibits excellent performance in noise-resistant testing and achieves an accuracy of more than 89% in field seismic data detection.

Deformation Characteristics of Hydrate-Bearing Sediments

The safe and efficient development of natural gas hydrate requires a deep understanding of the deformation behaviors of reservoirs.In this study,a series of triaxial shearing tests are carried out to investigate the deformation properties of hydrate-bearing sediments.Variations of volumetric and lateral strains versus hydrate saturation are analyzed comprehensively.Results indicate that the sediments with high hydrate saturation show dilative behaviors,which lead to strain-softening characteristics during shearing.The volumetric strain curves have a tendency to transform gradually from dilatation to compression with the increase in effective confining pressure.An easy prediction model is proposed to describe the relationship between volumetric and axial strains.The model coefficientβis the key dominating factor for the shape of volumetric strain curves and can be determined by the hydrate saturation and stress state.Moreover,a modified model is established for the calculation of lateral strain.The corresponding determination method is provided for the easy estimation of model coefficients for medium sand sediments containing hydrate.This study provides a theoretical and experimental reference for deformation estimation in natural gas hydrate development.

The suitability assessment on site selection for bottom seeding scallop culture based on analytic hierarchy process

Scallop culture is an important way of bottom-seeding marine ranching,which is of great significance to improve the current situation of fishery resources.However,there are some problems in site-selection evaluation of marine ranching,such as imperfect criteria system,complex structure,untargeted criteria quantification,etc.In addition,no site-selection evaluation method of bottom-seeding culture areas for scallops is available.Therefore,we established a hierarchy structure model according to the analytic hierarchy process(AHP)theory,in which social,physical,chemical,and biological environments are used as main criteria,and marine functional zonation,water depth,current,water temperature,salinity,substrate type,water quality,sediment quality,red tide,phytoplankton,and zooplankton are used as sub-criteria,on which a multi-parameter evaluation system is set up.Meanwhile,the dualism method,assignment method,and membership function method were used to quantify sub-criteria,and a quantitative evaluation for the entire criteria was added,including the evaluation and analysis of two types of unsuitable environmental situations.By overall consideration in scallop yield,quality,and marine ranching construction objectives,the weight of the main criteria could be determined.Five grades in the suitability corresponding to the evaluation result were divided,and the Python language was used to create an evaluation system for efficient calculation and intuitive presentation of the evaluation outcome.Eight marine cases were simulated based on existing survey data,and the results prove that the method is feasible for evaluating and analyzing the site selection of bottom-seeding culture areas for scallops under various environmental situations.The proposed evaluation method can be promoted for the site selection of bottom-seeding marine ranching.This study provided theoretical and methodological references for the site selection evaluation of other types of marine ranching.

A method for extracting the preseismic gravity anomalies over the Tibetan Plateau based on the maximum shear strain using GRACE data

The occurrence of earthquakes is closely related to the crustal geotectonic movement and the migration of mass,which consequently cause changes in gravity.The Gravity Recovery And Climate Experiment(GRACE)satellite data can be used to detect gravity changes associated with large earthquakes.However,previous GRACE satellite-based seismic gravity-change studies have focused more on coseismic gravity changes than on preseismic gravity changes.Moreover,the noise of the north–south stripe in GRACE data is difficult to eliminate,thereby resulting in the loss of some gravity information related to tectonic activities.To explore the preseismic gravity anomalies in a more refined way,we first propose a method of characterizing gravity variation based on the maximum shear strain of gravity,inspired by the concept of crustal strain.The offset index method is then adopted to describe the gravity anomalies,and the spatial and temporal characteristics of gravity anomalies before earthquakes are analyzed at the scales of the fault zone and plate,respectively.In this work,experiments are carried out on the Tibetan Plateau and its surrounding areas,and the following findings are obtained:First,from the observation scale of the fault zone,we detect the occurrence of large-area gravity anomalies near the epicenter,oftentimes about half a year before an earthquake,and these anomalies were distributed along the fault zone.Second,from the observation scale of the plate,we find that when an earthquake occurred on the Tibetan Plateau,a large number of gravity anomalies also occurred at the boundary of the Tibetan Plateau and the Indian Plate.Moreover,the aforementioned experiments confirm that the proposed method can successfully capture the preseismic gravity anomalies of large earthquakes with a magnitude of less than 8,which suggests a new idea for the application of gravity satellite data to earthquake research.

Dynamic Location Method for Shallow Ocean Bottom Nodes Using the Levenberg-Marquart Algorithm

Ocean bottom node(OBN)data acquisition is the main development direction of marine seismic exploration;it is widely promoted,especially in shallow sea environments.However,the OBN receivers may move several times because they are easily affected by tides,currents,and other factors in the shallow sea environment during long-term acquisition.If uncorrected,then the imaging quality of subsequent processing will be affected.The conventional secondary positioning does not consider the case of multiple movements of the receivers,and the accuracy of secondary positioning is insufficient.The first arrival wave of OBN seismic data in shallow ocean mainly comprises refracted waves.In this study,a nonlinear model is established in accordance with the propagation mechanism of a refracted wave and its relationship with the time interval curve to realize the accurate location of multiple receiver movements.In addition,the Levenberg-Marquart algorithm is used to reduce the influence of the first arrival pickup error and to automatically detect the receiver movements,identifying the accurate dynamic relocation of the receivers.The simulation and field data show that the proposed method can realize the dynamic location of multiple receiver movements,thereby improving the accuracy of seismic imaging and achieving high practical value.

Magmatic-hydrothermal Evolution and Mineralization Mechanisms of the Wangjiazhuang Cu(-Mo)Deposit in the Zouping Volcanic Basin,Shandong Province,China:Constraints from Fluid Inclusions

The Wangjiazhuang Cu(-Mo)deposit,located within the Zouping volcanic basin in western Shandong Province,China,is unique in this area for having an economic value.In order to expound the metallogenetic characteristics of this porphyry-like hydrothermal deposit,a detailed fluid inclusion study has been conducted,employing the techniques of representative sampling,fluid inclusion petrography,microthermometry,Raman spectroscopy,LA-ICP-MS analysis of single fluid inclusions,as well as cathode fluorescence spectrometer analysis of host mineral quartz.The deposit contains mainly two types of orebodies,i.e.veinlet-dissemination-stockwork orebodies in the K-Si alteration zone and pegmatiticquartz sulfide veins above them.In addition,minor breccia ore occurs locally.Four types of fluid inclusions in the deposit and altered quartz monzonite are identified:L-type one-or two-phase aqueous inclusions,V-type vapor-rich inclusions with V/L ratios greater than 50%-90%,D-type multiphase fluid inclusions containing daughter minerals or solids and S-type silicate-bearing fluid inclusions containing mainly muscovite and biotite.Ore petrography and fluid inclusion study has revealed a three-stage mineralization process,driven by magmatic-hydrothermal fluid activity,as follows.Initially,a hydrothermal fluid,separated from the parent magma,infiltrated into the quartz monzonite,resulting in its extensive K-Si alteration,as indicated by silicate-bearing fluid inclusions trapped in altered quartz monzonite.This is followed by the early mineralization,the formation of quartz veinlets and dissemination-stockwork ores.During the main mineralization stage,due to the participation and mixing of meteoric groundwater with magmatic-sourced hydrothermal fluid,the cooling and phase separation caused deposition of metals from the hydrothermal fluids.As a result,the pegmatitic-quartz sulfide-vein ores formed.In the late mineralization stage,decreasing fluid salinity led to the formation of L-type aqueous inclusions and chalcopyrite-sulfosalt ore.Coexistence of V-type and D-type inclusions and their similar homogenization temperatures with different homogenization modes suggest that phase separation or boiling of the ore-forming fluids took place during the early and the main mineralization stages.The formation P-T conditions of S-type inclusions and the early and the main mineralization stages were estimated as ca.156-182 MPa and 450-650℃,350-450℃,18-35 MPa and 280-380℃,8-15 MPa,respectively,based on the microthermometric data of the fluid inclusions formed at the individual stages.

First-Arrival Picking Method for Active Source Data with Ocean Bottom Seismometers Based on Spatial Waveform Variation Characteristics

The precision and reliability of first-arrival picking are crucial for determining the accuracy of geological structure inversion using active source ocean bottom seismometer(OBS)refraction data.Traditional methods for first-arrival picking based on sample points are characterized by theoretical errors,especially in low-sampling-frequency OBS data because the travel time of seismic waves is not an integer multiple of the sampling interval.In this paper,a first-arrival picking method that utilizes the spatial waveform variation characteristics of active source OBS data is presented.First,the distribution law of theoretical error is examined;adjacent traces exhibit variation characteristics in their waveforms.Second,a label cross-correlation superposition method for extracting highfrequency signals is presented to enhance the first-arrival picking precision.Results from synthetic and field data verify that the proposed approach is robust,successfully overcomes the limitations of low sampling frequency,and achieves precise outcomes that are comparable with those of high-sampling-frequency data.

A Marine or Continental Nature of the Deltas in the Early Cretaceous Lingshandao Formation--Evidences from Trace Elements

Objective A succession with well-developed soft-sediment deformation structures and well-exposed gravity-flow deposits in the Early Cretaceous Lingshandao Formation has caused considerable controversy on whether the nondisputed deltaic environment was marine or continental.This dispute cannot be solved by mere sedimentological analysis and,moreover,too few fossils are present to be decisive.

Numerical Simulation for Accuracy of Velocity Analysis in Small-Scale High-Resolution Marine Multichannel Seismic Technology

When used with large energy sparkers, marine multichannel small-scale high-resolution seismic detection technology has a high resolution, high-detection precision, a wide applicable range, and is very flexible. Positive results have been achieved in submarine geological research, particularly in the investigation of marine gas hydrates. However, the amount of traveltime difference information is reduced for the velocity analysis under conditions of a shorter spread length, thus leading to poorer focusing of the velocity spectrum energy group and a lower accuracy of the velocity analysis. It is thus currently debatable whether the velocity analysis accuracy of short-arrangement multichannel seismic detection technology is able to meet the requirements of practical application in natural gas hydrate exploration. Therefore, in this study the bottom boundary of gas hydrates(Bottom Simulating Reflector, BSR) is used to conduct numerical simulation to discuss the accuracy of the velocity analysis related to such technology. Results show that a higher dominant frequency and smaller sampling interval are not only able to improve the seismic resolution, but they also compensate for the defects of the short-arrangement, thereby improving the accuracy of the velocity analysis. In conclusion, the accuracy of the velocity analysis in this small-scale, high-resolution, multi-channel seismic detection technology meets the requirements of natural gas hydrate exploration.

Formation mechanism,experimental method,and property characterization of graindisplacing methane hydrates in marine sediment:A review

Grain-displacing hydrate deposits exist at many marine sites,which constitute an important part of methane hydrate resources worldwide.Attributed to the difficulties in acquiring field data and synthesizing experimental samples,the formation and property characterization of grain-displacing hydrate remains less understood and characterized than the pore-filling hydrate in current literature.This study reviews the formation mechanisms of grain-displacing hydrate from the perspective of geological accumulation and microscale sedimentary property.The experimental methods of synthesizing grain-displacing hydrate in the laboratory and the current knowledge on the property of grain-displacing hydrate sediment are also introduced.Shortcomings in current theories and suggestions for future study are proposed.The work is hoped to provide valuable insights for the research into the hydrate accumulation,geophysics,and hydrate exploitation targeted at the grain-displacing hydrate in the marine sediments.

A three-dimensional environmental monitoring system for the production of marine gas hydrates

1.Objectives As about 99% of the natural gas hydrate (NGH)resources are stored on the seafloor,it is inevitable that in its process of exploration and production,there will be geohazards such as methane leakage,seabed subsidence,ecosystem damage and even landslide.Therefore,a set of integrated systems for environmental monitoring in the process of NGH production has been proposed in the consideration of environmental disasters that may arise in the process of hydrate production.

Distinct BSRs and their implications for natural gas hydrate formation and distribution at the submarine Makran accretionary zone

To investigate the nature of gas hydrates in the Makran area,new high-resolution geophysical data were acquired between 2018-2019.The data collected comprise multibeam and two-dimensional multi-channel seismic reflection data.The multibeam bathymetry data show East-North-East(ENE)ridges,piggy-back basins,canyon and channel systems,and the morphology of the abyssal plain.Continuous and discontinuous bottom simulating reflectors(BSRs)occur in the piggy-back basins on most of the seismic profiles available.The BSRs cut the dipping layers with strong amplitude and reversed polarity.Discontinuous BSRs indicate a transition along a dipping high-permeable sand layers from gas-rich segment to the gas hydrate-bearing segment and sugge st alternating sediments of fine and relatively coarse grain size.Double BSRs are highly dynamic and attributed to slumps occurring in the study area.The BSRs induced by slumps are located both at deep and shallow depths,responding to the temperature or pressure variation.For the first time,BSRs are observed in the abyssal plain of the Makran area,being associated with anticline structures,which do not show large spatial continuity and are strongly conditioned by structural conditions such as anticlines and fluid migration pathways,including deep fault,gas chimney,and high-permeable sedimentary layer.Our results may help to assess the gas hydrate potential within the piggy-back basins and to determine the most promising target areas.Moreover,results about the abyssal plain BSR may help to locate hydrocarbon reservoirs in the deep ocean.

Differential tectonic evolution and formation mechanism of three subsags in Wenchang Sag of Pearl River Mouth Basin,South China Sea

Notable differences in the structural characteristics and evolution of three adjacent sub-sags,i.e.,the Wenchang sub-sags A,B,and C,on the downthrown side of the Zhu IlI South Fault in the Wenchang Sag,are significant as they affect the formation and distribution of the oil and gas in these three sub-sags.However,the differences in their tectonic evolutions and formation mechanisms have not yet been adequately explained.In this paper,stress analysis,equilibrium profiles,and paleogeomorphic restora-tion,are used to investigate the dynamic settings,formation mechanisms,and influencing factors of the structural deformation related to the formation of the Wenchang Sag based on interpretation of seismic data.The results of the stress analysis suggest clockwise deflection of the regional tensile stress direction from a WNW-ESE trend during the Early Paleocene to NW-SE and NNW-SSE trends during the Eocene,to a nearly N-S trend during the Oligocene,and finally to a NNE-SSW trend during the Miocene.This clockwise rotation of the regional tensile stress direction led to the formation of a dextral strike-slip stress component parallel to the NE-trending Zhu I South Fault.This strike-slip stress component formed a releasing bend in sub-sag A,and may be associated with the continuous subsidence of a thick sedimentary layer in sub-sag A.It also created a restraining bend in sub-sag B,which underwent multiple structural inversions during its extension and subsidence and has a relatively s mall sedimentary thick-ness.The double restraining bend in sub-sag C is considered to have been strongly uplifted and eroded in response to this strike-slip stress component.Four obvious structural inversions in sub-sag B are iden-tified in this paper.These structural inversions correspond to the last four regional tectonic movements.This interpretation suggests that the formation of the structural inversions was likely related to the strong tensile stress and the small intersection angle between the direction of the regional tensile stress and the pre-existing boundary fault.The rotation of the tensile stress direction was responsible for the strike-slip movement on the pre-existing boundary fault and the formation of the releasing bend and restraining bend,which controlled the structural evolutions of the sub-sags.This reasonably explains the differential tectonic evolution of these three sub-sags in the Wenchang Sag,and provides a crucial idea forstructuralanalysisof similarbasins.

Laboratory insights into the effects of methane hydrate on the anisotropic joint elastic-electrical properties in fractured sandstones

Fractured hydrate-bearing reservoirs show significantly anisotropic geophysical properties. The joint application of seismic and electromagnetic explorations is expected to accurately assess hydrate resources in the fractured reservoirs. However, the anisotropic joint elastic-electrical properties in such reservoirs that are the key to the successful application of the joint explorations, remain poorly understood. To obtain such knowledge, we designed and implemented dedicated laboratory experiments to study the anisotropic joint elastic-electrical properties in fractured artificial silica sandstones (with fracture density of about 6.2%, porosity of approximately 25.7%, and mean grainsize of 0.089 mm) with evolving methane hydrate. The experimental results showed that the anisotropic compressional wave velocities respectively increased and decreased with the forming and dissociating hydrate, and the variation in the increasing trend and the decreasing extent of the velocity perpendicular to the fractures were more significant than that parallel to the fractures, respectively. The experimental results also showed that the overall decreasing trend of the electrical conductivity parallel to the fractures was steeper than that perpendicular to the fractures during hydrate formation, and the general variations of the two conductivities with complex trend were similar during hydrate dissociation. The variations in the elastic and electrical anisotropic parameters with forming and dissociating hydrate were also found to be distinct. Interpretation of the experimental results suggested that the hydrate binding to the grains evolved to bridge the surfaces of fractures when saturation exceeded 10% during hydrate formation, and the bridging hydrate gradually evolved to floating in fractures during dissociation. The experimental results further showed that the anisotropic velocities and electrical conductivities were correlated with approximately consistent trends of different slopes during hydrate formation, and the joint elastic-electrical anisotropic parameters exhibited a sharp peak at the hydrate saturation of about 10%. The results suggested that the anisotropic joint properties can be employed not only to accurately estimate hydrate saturation but also possibly to identify hydrate distribution in the fractures.

Experimental study on solid particle migration and production behaviors during marine natural gas hydrate dissociation by depressurization

Sand production is one of the main obstacles restricting gas extraction efficiency and safety from marine natural gas hydrate(NGH)reservoirs.Particle migration within the NGH reservoir dominates sand production behaviors,while their relationships were rarely reported,severely constrains quantitative evaluation of sand production risks.This paper reports the optical observations of solid particle migration and production from micrometer to mesoscopic scales conditioned to gravel packing during depressurization-induced NGH dissociation for the first time.Theoretical evolutionary modes of sand migration are established based on experimental observations,and its implications on field NGH are comprehensively discussed.Five particle migration regimes of local borehole failure,continuous collapse,wormhole expansion,extensive slow deformation,and pore-wall fluidization are proved to occur during depressurization.The types of particle migration regimes and their transmission modes during depressurization are predominantly determined by initial hydrate saturation.In contrast,the depressurization mainly dominates the transmission rate of the particle migration regimes.Furthermore,both the cumulative mass and the medium grain size of the produced sand decrease linearly with increasing initial methane hydrate(MH)saturation.Discontinuous gas bubble emission,expansion,and explosion during MH dissociation delay sand migration into the wellbore.At the same time,continuous water flow is a requirement for sand production during hydrate dissociation by depressurization.The experiments enlighten us that a constitutive model that can illustrate visible particle migration regimes and their transmission modes is urgently needed to bridge numerical simulation and field applications.Optimizing wellbore layout positions or special reservoir treatment shall be important for mitigating sand production tendency during NGH exploitation.

Geochemical,mineralogical,and petrological analyses for the interpretation of the sedimentary environment of the Middle-Late Ordovician Majiagou Formation(northern China)as a tool for more effective gas exploration

Core samples from the deeply buried Ordovician Majiagou Formation below the Huainan Coalfield(E China) have been investigated for their carbonate types,major and trace elements(including rare earth elements) and C and O isotopes,The objective was to get a better insight into the possible occurrences of gas(and possibly oil) derived from Carboniferous coals.It was found that the carbonates are dolomites with strongly varying amounts of CaO and MgO.The low concentrations of SiO_(2),Al_(2)O_(3),and Fe_(2)O_(3)indicate deposition in a normal marine environment with little terrigenous input,The Na_(2)O/K_(2)O,Fe/Mn and Sr/Ba ratios,as well as the Ga values indicate mainly a marine salinity and a hot and humid climate.The slight depletion of Ce and Eu,the depletion of heavy rare earth elements(HREE) and the enrichment of light rare earth elements(LREE) indicate deposition in a reducing environment.It thus appears that the Majiagou Formation below the Huainan Coalfield closely resembles that in the eastern part of the Ordos Basin,where several gas reservoirs are present,so that the Majiagou Formation under the Huainan Coalfield represents a promising target for hydrocarbon exploration.

Evolution of nC_(16)H_(34)-water-mineral systems in thermal capsules and geological implications for deeply-buried hydrocarbon reservoirs

Organic-inorganic interactions between hydrocarbons and most minerals in deeply buried reservoirs remain unclear.In this study,gold capsules and fused silica capillary capsules(FSCCs)with different com-binations of nC_(16)H_(34),water(distilled water,CaCl_(2) water)and minerals(quartz,feldspar,calcite,kaolinite,smectite,and illite)were heated at 340℃ for 3-10 d,to investigate the evolution and reaction pathways of the organic-inorganic interactions in different hot systems.After heating,minerals exhibited little alteration in the anhydrous systems.Mineral alterations,how-ever,occurred obviously in the hydrous systems.Different inorganic components affected nC_(16)H_(34) degra-dation differently.Overall,water promoted the free-radical thermal-cracking reaction and step oxidation reaction but suppressed the free-radical cross-linking reaction.The impact of CaCl_(2) water on the nC_(16)H_(34) degradation was weaker than the distilled water as high Ca^(2+)concentration suppressed the formation of free radicals.The presence of different waters also affects the impact of different minerals on nC_(16)H_(34) degradation,via its impact on mineral alterations.In the anhydrous nC_(16)H_(34)-mineral systems,calcite and clays catalyzed generation of low-molecular-weight(LMW)alkanes,particularly the clays.Quartz,feldspar,and calcite catalyzed generation of high-molecular-weight(HMW)alkanes and PAHs,whereas clays catalyzed the generation of LMW alkanes and mono-bicyclic aromatic hydrocarbons(M-BAHs).In the hydrous nC_(16)H_(34)-distilled water-mineral systems,all minerals but quartz promoted nC_(16)H_(34) degra-dation to generate more LMW alkanes,less HMW alkanes and PAHs.In the nC_(16)H_(34)-CaCl_(2) water-mineral systems,the promotion impact of minerals was weaker than that in the systems with distilled water.This study demonstrated the generation of different hydrocarbons with different fluorescence colors in the different nC_(16)H_(34)-water-mineral systems after heating for the same time,implying that fluorescence colors need to be interpreted carefully in investigation of hydrocarbon charging histories and oil origins in deeply buried reservoirs.Besides,the organic-inorganic interactions in different nC_(16)H_(34)-water-mineral systems proceeded in different pathways at different rates,which likely led to preservation of liquid hydrocarbons at different depth(temperature).Thus,quantitative investigations of the reaction kinetics in different hydrocarbon-water-rock systems are required to improve the prediction of hydrocar-bon evolution in deeply buried hydrocarbon reservoirs.