Constraints on Characteristics and Distribution of Gas Hydrate and Free Gas Using Broad-Band Processing of Three-Dimensional Seismic Data

Gas hydrate drilling expeditions in the Pearl River Mouth Basin,South China Sea,have identified concentrated gas hydrates with variable thickness.Moreover,free gas and the coexistence of gas hydrate and free gas have been confirmed by logging,coring,and production tests in the foraminifera-rich silty sediments with complex bottom-simulating reflectors(BSRs).The broad-band processing is conducted on conventional three-dimensional(3D)seismic data to improve the image and detection accuracy of gas hydratebearing layers and delineate the saturation and thickness of gas hydrate-and free gas-bearing sediments.Several geophysical attributes extracted along the base of the gas hydrate stability zone are used to demonstrate the variable distribution and the controlling factors for the differential enrichment of gas hydrate.The inverted gas hydrate saturation at the production zone is over 40% with a thickness of 90 m,showing the interbedded distribution with different boundaries between gas hydrate-and free gas-bearing layers.However,the gas hydrate saturation value at the adjacent canyon is 70%,with 30-m-thick patches and linear features.The lithological and fault controls on gas hydrate and free gas distributions are demonstrated by tracing each gas hydrate-bearing layer.Moreover,the BSR depths based on broad-band reprocessed 3D seismic data not only exhibit variations due to small-scale topographic changes caused by seafloor sedimentation and erosion but also show the upward shift of BSR and the blocky distribution of the coexistence of gas hydrate and free gas in the Pearl River Mouth Basin.

Extreme massive hydraulic fracturing in deep coalbed methane horizontal wells:A case study of the Linxing Block,eastern Ordos Basin,NW China

Deep coal seams show low permeability,low elastic modulus,high Poisson’s ratio,strong plasticity,high fracture initiation pressure,difficulty in fracture extension,and difficulty in proppants addition.We proposed the concept of large-scale stimulation by fracture network,balanced propagation and effective support of fracture network in fracturing design and developed the extreme massive hydraulic fracturing technique for deep coalbed methane(CBM)horizontal wells.This technique involves massive injection with high pumping rate+high-intensity proppant injection+perforation with equal apertures and limited flow+temporary plugging and diverting fractures+slick water with integrated variable viscosity+graded proppants with multiple sizes.The technique was applied in the pioneering test of a multi-stage fracturing horizontal well in deep CBM of Linxing Block,eastern margin of the Ordos Basin.The injection flow rate is 18 m^(3)/min,proppant intensity is 2.1 m^(3)/m,and fracturing fluid intensity is 16.5 m^(3)/m.After fracturing,a complex fracture network was formed,with an average fracture length of 205 m.The stimulated reservoir volume was 1987×10^(4)m^(3),and the peak gas production rate reached 6.0×10^(4)m^(3)/d,which achieved efficient development of deep CBM.

Probability Distribution Characteristics of Strong Nonlinear Waves Under Typhoon Conditions in the Northern South China Sea

The generation and propagation mechanism of strong nonlinear waves in the South China Sea is an essential research area. In this study, the third-generation wave model WAVEWATCH Ⅲ is employed to simulate wave fields under extreme sea states. The model, integrating the ST6 source term, is validated against observed data, demonstrating its credibility. The spatial distribution of the occurrence probability of strong nonlinear waves during typhoons is shown, and the waves in the straits and the northeastern part of the South China Sea show strong nonlinear characteristics. The high-order spectral model HOS-ocean is employed to simulate the random wave surface series beneath five different platform areas. The waves during the typhoon exhibit strong nonlinear characteristics, and freak waves exist. The space-varying probability model is established to describe the short-term probability distribution of nonlinear wave series. The exceedance probability distributions of the wave surface beneath different platform areas are compared and analyzed. The results show that with an increase in the platform area, the probability of a strong nonlinear wave beneath the platform increases.

Study on Free-Standing Riser Under Different Working Conditions Based on 3D Co-Rotational Beam Element

When the free standing riser(FSR)is in service in the ocean,its mechanical properties are affected by various factors,including complex ocean current forces,buoyancy of the buoyancy can,and torque caused by the deflection of the upper floating body.These loads have a great influence on the deformation and internal force of the FSR.The static performance of FSR is investigated in this research under various working conditions.The finite element model of FSR is established based on the co-rotational method.The arc length approach is used to solve the model.The load is exerted in increments.The current load on the riser changes with the configuration of the riser.The accuracy of the numerical method is verified by Abaqus software.The calculation time is also compared.Then,the effects of uniform current,actual current and floating body yaw motion on FSR are studied by parameter analysis.Additionally,the influence of the FSR on the ocean current after the failure of part of the buoyancy can chamber is analyzed.The results show that the numerical model based on the co-rotational method can effectively simulate the large rotation and torsion behavior of FSR.This method has high computational efficiency and precision,and this method can quickly improve the efficiency of numerical calculation of static analysis of deep-water riser.The proposed technology may serve as an alternative to the existing proprietary commercial software,which uses a complex graphical user interface.

Bending Failure Mode and Prediction Method of the Compressive Strain Capacity of A Submarine Pipeline with Dent Defects

A dent is a common type of defects for submarine pipeline.For submarine pipelines,high hydrostatic pressure and internal pressure are the main loads.Once pipelines bend due to complex subsea conditions,the compression strain capacity may be exceeded.Research into the local buckling failure and accurate prediction of the compressive strain capacity are important.A finite element model of a pipeline with a dent is established.Local buckling failure under a bending moment is investigated,and the compressive strain capacity is calculated.The effects of different parameters on pipeline local buckling are analyzed.The results show that the dent depth,external pressure and internal pressure lead to different local buckling failure modes of the pipeline.A higher internal pressure indicates a larger compressive strain capacity,and the opposite is true for external pressure.When the ratio of external pressure to collapse pressure of intact pipeline is greater than 0.1,the deeper the dent,the greater the compressive strain capacity of the pipeline.And as the ratio is less than 0.1,the opposite is true.On the basis of these results,a regression equation for predicting the compressive strain capacity of a dented submarine pipeline is proposed,which can be referred to during the integrity assessment of a submarine pipeline.

Evaluation of the injection and plugging ability of a novel epoxy resin in cement cracks

Sustained casing pressure(SCP)is a crucial issue in the oil and gas production lifecycle.Epoxy resins,exhibiting exceptional compressive strength,ductility,and shear bonding strength,have the potential to form reliable barriers.The injectivity and sealing capacity of the epoxy resin is crucial parameters for the success of shallow remediation operations.This study aimed to develop and assess a novel solid-free resin sealant as an alternative to Portland cement for mitigating fluid leakage.The investigation evaluated the viscosity,compressive strength,and brittleness index of the epoxy resin sealant,as well as its tangential and normal shear strengths in conjunction with casing steel.The flow characteristics and sealing abilities of conventional cement and epoxy resin were comparatively analyzed in cracks.The results showed that the application of a viscosity reducer facilitated control over the curing time of the epoxy resin,ranging from 1.5 to 6 h,and reduced the initial viscosity from 865.53 to 118.71 m Pa,s.The mechanical properties of the epoxy resin initially increased with a rise in curing agent content before experiencing a minor decrease.The epoxy resin containing 30%curing agent exhibited optimal mechanical properties.After a 14-day curing period,the epoxy resin's compressive strength reached81.37 MPa,2.12 times higher than that of cement,whereas the elastic modulus of cement was 2.99 times greater than that of the epoxy resin.The brittleness index of epoxy resin is only 3.42,demonstrating high flexibility and toughness.The tangential and normal shear strengths of the epoxy resin exceeded those of cement by 3.17 and 2.82 times,respectively.In a 0.5 mm-wide crack,the injection pressure of the epoxy resin remained below 0.075 MPa,indicating superior injection and flow capabilities.Conversely,the injection pressure of cement surged dramatically to 2.61 MPa within 5 min.The breakthrough pressure of0.5 PV epoxy resin reached 7.53 MPa,decreasing the crack's permeability to 0.02 D,a mere 9.49%of the permeability observed following cement plugging.Upon sealing a 2 mm-wide crack using epoxy resin,the maximum breakthrough pressure attained 5.47 MPa,3.48 times of cement.These results suggest that epoxy resin sealant can be employed safely and effectively to seal cracks in the cement.

Coupling relationship and genetic mechanisms of shelf-edge delta and deep-water fan source-to-sink:A case study in Paleogene Zhuhai Formation in south subsag of Baiyun Sag,Pearl River Mouth Basin,China

The coupling relationship between shelf-edge deltas and deep-water fan sand bodies is a hot and cutting-edge field of international sedimentology and deep-water oil and gas exploration.Based on the newly acquired high-resolution 3D seismic,logging and core data of Pearl River Mouth Basin(PRMB),this paper dissected the shelf-edge delta to deep-water fan(SEDDF)depositional system in the Oligocene Zhuhai Formation of Paleogene in south subsag of Baiyun Sag,and revealed the complex coupling relationship from the continental shelf edge to deep-water fan sedimentation and its genetic mechanisms.The results show that during the deposition of the fourth to first members of the Zhuhai Formation,the scale of the SEDDF depositional system in the study area showed a pattern of first increasing and then decreasing,with deep-water fan developed in the third to first members and the largest plane distribution scale developed in the late stage of the second member.Based on the development of SEDDF depositional system along the source direction,three types of coupling relationships are divided,namely,deltas that are linked downdip to fans,deltas that lack downdip fans and fans that lack updip coeval deltas,with different depositional characteristics and genetic mechanisms.(1)Deltas that are linked downdip to fans:with the development of shelf-edge deltas in the shelf area and deep-water fans in the downdip slope area,and the strong source supply and relative sea level decline are the two key factors which control the development of this type of source-to-sink(S2S).The development of channels on the continental shelf edge is conducive to the formation of this type of S2S system even with weak source supply and high sea level.(2)Deltas that lack downdip fans:with the development of shelf edge deltas in shelf area,while deep water fans are not developed in the downdip slope area.The lack of“sources”and“channels”,and fluid transformation are the three main reasons for the formation of this type of S2S system.(3)Fans that lack updip coeval deltas:with the development of deep-water fans in continental slope area and the absence of updip coeval shelf edge deltas,which is jointly controlled by the coupling of fluid transformation at the shelf edge and the“channels”in the continental slope area.

Paleoenvironmental Characteristics of Paleogene Lacustrine Source Rocks in the Western Bozhong Sag,Bohai Bay Basin,China:Evidence from Biomarkers,Major and Trace Elements

The organic matter(OM)enrichment mechanisms and depositional environment characteristics of lacustrine source rocks in the western Bozhong Sag,Bohai Bay Basin in Northeast China remain controversial.To address these issues,based on Rock-Eval pyrolysis,kerogen macerals,H/C and O/C ratios,GC-MS,major and trace elements,the Dongying Formation Member(Mbr)3(E_(3)d_(3)),the Shahejie Formation mbrs 1 and 2(E_(2)s_(1+2)),and the Shahejie Mbr 3(E_(2)s_(3))source rocks in the western Bozhong Sag were studied.The above methods were used to reveal their geochemical properties,OM origins and depositional environments,all of which indicate that E_(2)s_(1+2)and E_(2)s_(3)are excellent source rocks,and that E_(3)d_(3)is of the second good quality.E_(3)d_(3)source rocks were formed under a warm and humid climate,mainly belong to fluvial/delta facies,the E_(3)d_(3)sediments formed under weakly oxidizing and freshwater conditions.Comparatively,the depositional environments of E_(2)s_(1+2)source rocks were arid and cold climate,representing saline or freshwater lacustrine facies,and the sediments of E_(2)s_(1+2)belong to anoxic or suboxic settings with large evaporation and salinity.During the period of E_(2)s_(3),the climate became warm and humid,indicating the freshwater lacustrine facies,and E_(2)s_(3)was characterized by freshwater and abundant algae.Moreover,compared with other intervals,the OM origin of E_(3)d_(3)source rocks has noticeable terrestrial input.The OM origin of the E_(2)s_(1+2)and E_(2)s_(3)are mainly plankton and bacteria.Tectonic subsidence and climate change have affected the changes of the depositional environment in the western Bozhong Sag,thus controlling the distribution of the source rocks,the geochemical characteristics in the three intervals of lacustrine source rocks have distinct differences.Overall,these factors are effective to evaluate the paleoenvironmental characteristics of source rocks by biomarkers,major and trace elements.The established models may have positive implications for research of lacustrine source rocks in offshore areas with few drillings.

Study on the clogging mechanism of punching screen in sand control by the punching structure parameters

As an independent sand control unit or a common protective shell of a high-quality screen,the punching screen is the outermost sand retaining unit of the sand control pipe which is used in geothermal well or oil and gas well.However,most screens only consider the influence of the internal sand retaining medium parameters in the sand control performance design while ignoring the influence of the plugging of the punching screen on the overall sand retaining performance of the screen.To explore the clogging mechanism of the punching screen,this paper established the clogging mechanism calculation model of a single punching screen sand control unit by using the computational fluid mechanics-discrete element method(CFD-DEM)combined method.According to the combined motion of particles and fluids,the influence of the internal flow state on particle motion and accumulation was analyzed.The results showed that(1)the clogging process of the punching sand control unit is divided into three stages:initial clogging,aggravation of clogging and stability of clogging.In the initial stage of blockage,coarse particles form a loose bridge structure,and blockage often occurs preferentially at the streamline gathering place below chamfering inside the sand control unit.In the stage of blockage intensification,the particle mass develops into a relatively complete sand bridge,which develops from both ends of the opening to the center of the opening.In the stable plugging stage,the sand deposits show a“fan shape”and form a“V-shaped”gully inside the punching slot element.(2)Under a certain reservoir particle-size distribution,The slit length and opening height have a large influence on the permeability and blockage rate,while the slit width size has little influence on the permeability and blockage rate.The microscopic clogging mechanism and its law of the punching screen prevention unit are proposed in this study,which has some field guidance significance for the design of punching screen and sand prevention selection.

Application of C_(30)tetracyclic polyprenoids as effective biomarker in oil-to source rock correlation in the ZhuⅢdepression,Zhujiangkou Basin,northern South China

The northern South China Sea,including the Zhujiangkou Basin and the Beibuwan Basin,developed high-quality lacustrine source rocks during the Eocene rifting period.These source rocks are vital for hydrocarbon generation in the northern South China Sea.The ZhuⅠdepression in the Zhujiangkou Basin and the Beibuwan Basin typically exhibit high abundance of C_(30)4-methyl steranes.However,shales in the Eocene Wenchang Formation in the ZhuⅢdepression of the Zhujiangkou Basin contains lower quantities of high-quality lacustrine source rocks with 4-methyl steranes,which often co-elute with some pentacyclic triterpanes in gas chromatography-mass spectrometry(GC-MS).Therefore,the single 4-methylsterane parameter based on GC-MS cannot accurately distinguish organic source in the deep to semi-deep water lacustrine source rocks of the Wenchang Formation from other source rocks,thus impeding the recognition of their contributions to petroleum reservoirs.In this study,GC-MS of aliphatic hydrocarbons,palynofacies and algal identification,as well as stable carbon isotope data of organic matter were used to identify the algal species and construct the paleoclimate during deposition of the Wenchang Formation source rocks in the ZhuⅢdepression of the Zhujiangkou Basin.It is suggested that during the Wenchang Formation period,freshwater green algae prevailed in the lake,which is likely account for the relatively low content of 4-methyl steranes in the high-quality lacustrine source rocks.Controlled by the algal species,it is proposed that the content of C_(30)tetracyclic polyprenoids(TPP)can better indicate the quality of the Wenchang source rocks than C_(30)4-methyl steranes.Consequently,a relationship between the TPP index and the quality of the lacustrine source rocks in the Wenchang Formation of the ZhuⅢdepression was established.A higher TPP index indicates higher organic matter abundance and hydrogen index of the lacustrine source rocks.When applied to the origin analysis of oils in the ZhuⅢdepression,it is believed that the organic-rich deep lacustrine source rocks in the Wenchang Formation made great contribution to the transitional zone crude oils in the Wenchang A and Wenchang B depressions.

Dynamic Analysis of A Deepwater Drilling Riser with A New Hang-off System

The safety of risers in hang-off states is a vital challenge in offshore oil and gas engineering.A new hang-off system installed on top of risers is proposed for improving the security of risers.This approach leads to a challenging problem:coupling the dynamics of risers with a new hang-off system combined with multiple structures and complex constraints.To accurately analyze the dynamic responses of the coupled system,a coupled dynamic model is established based on the Euler-Bernoulli beam-column theory and penalty function method.A comprehensive analysis method is proposed for coupled dynamic analysis by combining the finite element method and the Newmarkβmethod.An analysis program is also developed in MATLAB for dynamic simulation.The simulation results show that the dynamic performances of the risers at the top part are significantly improved by the new hang-off system,especially the novel design,which includes the centralizer and articulation joint.The bending moment and lateral deformation of the risers at the top part decrease,while the hang-off joint experiences a great bending moment at the bottom of the lateral restraint area which requires particular attention in design and application.The platform navigation speed range under the safety limits of risers expands with the new hang-off system in use.

Geochemical characteristics and origins of natural gases in the eastern Cote d’Ivoire Basin,West Africa

The gas sources in the eastern Cote d’Ivoire Basin(Tano Basin)are seldom reported and remain controversial due to multiple sets of potential source rocks and poorly documented geochemical characteristics of natural gases.The marine source rock potential from the Upper Albian to Turonian as well as the molecular composition and the stable carbon isotope composition of natural gases in the eastern Cote d’Ivoire Basin were studied in detail to investigate the origins of natural gases.The total organic carbon(TOC),hydrogen index(HI),and generation potential(S_(1)+S_(2))of source rocks indicate that both sapropelic source rocks and humic source rocks developed during the late Albian,whereas sapropelic source rocks developed during the Cenomanian and the Turonian.The normal order ofδ^(13)CH_(4)<δ^(13)C_(2)H6<δ^(13)C_(3)H_(8)(δ^(13)C_(1)<δ^(13)C_(2)<δ^(13)C_(3)),the relationship between C_(2)/C_(3)molar ratio andδ^(13)C_(2)-δ^(13)C_(3),and the plot ofδ^(13)C_(1)versus C_(1)/(C_(2)+C_(3))collectively show that the natural gases are thermogenic due to the primary cracking of kerogen,including the typical oil-associated gases from Well D-1,the mixed oil-associated gases and coal-derived gases from Well G-1 and Well L-1.Based on the plot ofδ^(13)C_(1)versusδ^(13)C_(2)and the established relationship betweenδ^(13)C_(1)and equivalent vitrinite reflectance(Ro),we proposed that the natural gases are in a mature stage(Ro generally varies from 1.0%to 1.3%).Combined with results of basin modelling and oil-to-source correlation,the transitional to marine source rocks during the late Albian were thought to have made a great contribution to the natural gases.Our study will make a better understanding on petroleum system in the eastern Cote d’Ivoire Basin.

Origin of condensate oil and natural gas in Bozhong 19-6 gas field,Bohai Bay Basin

The discovery of the Bozhong 19-6 gas field has opened a new frontier for deep gas exploration in the Bohai Bay Basin,with a great potential for further gas exploration.However,poor understanding of oil and gas origin has been limiting the exploration progress in this area.To clarify the origin of condensate oil and gas in Bozhong 19-6 gas field,this study adequately utilized the organic geochemical analysis data to investigate the composition and geochemical characteristics of condensate oil and natural gas,and analyzed the relationship between condensate oil and the three sets of source rocks in the nearby subsags.Results show that the lighter components dominate the condensate oil,with a forward type predominance.The parent material of crude oil was primarily deposited in a shallow,clay-rich,low-salinity,weakly reducing aquatic environment.The condensate and natural gas have similar parent source characteristics and maturity,with Ro ranging from 1.4%to 1.6%.Both are products of high maturity stage,indicating that they are hydrocarbon compounds produced by the same group of source rocks in the same stage.Oil-sources correlation shows that condensate oil and gas mainly originate from the source rocks of the third member of Shahejie Formation in the nearby subsags of the Bozhong 19-6 structural belt.

Time-Lapse Full-Waveform Inversion Using Cross-Correlation-Based Dynamic Time Warping

Offshore carbon capture, utilization, and storage(OCCUS) is regarded as a crucial technology for mitigating greenhouse gas emissions.Quantitative monitoring maps of sealed carbon dioxide are necessary in a comprehensive OCCUS project. A potential high-resolution method for the aforementioned purpose lies in the full-waveform inversion(FWI) of time-lapse seismic data. However, practical applications of FWI are severely restricted by the well-known cycle-skipping problem. A new time-lapse FWI method using cross-correlation-based dynamic time warping(CDTW) is proposed to detect changes in the subsurface property due to carbon dioxide(CO_(2)) injection and address the aforementioned issue. The proposed method, namely CDTW, which combines the advantages of cross-correlation and dynamic time warping, is employed in the automatic estimation of the discrepancy between the seismic signals simulated using the baseline/initial model and those acquired. The proposed FWI method can then back-project the estimated discrepancy to the subsurface space domain, thereby facilitating retrieval of the induced subsurface property change by taking the difference between the inverted baseline and monitor models. Numerical results on pairs of signals prove that CDTW can obtain reliable shifts under amplitude modulation and noise contamination conditions. The performance of CDTW substantially outperforms that of the conventional dynamic time warping method. The proposed time-lapse fullwaveform inversion(FWI) method is applied to the Frio-2 CO_(2) storage model. The baseline and monitor models are inverted from the corresponding time-lapse seismic data. The changes in velocity due to CO_(2) injection are reconstructed by the difference between the baseline and the monitor models.

Application of Feature, Event, and Process Methods to Leakage Scenario Development for Offshore CO_(2) Geological Storage

Offshore carbon dioxide(CO_(2)) geological storage(OCGS) represents a significant strategy for addressing climate change by curtailing greenhouse gas emissions. Nonetheless, the risk of CO_(2) leakage poses a substantial concern associated with this technology. This study introduces an innovative approach for establishing OCGS leakage scenarios, involving four pivotal stages, namely, interactive matrix establishment, risk matrix evaluation, cause–effect analysis, and scenario development, which has been implemented in the Pearl River Estuary Basin in China. The initial phase encompassed the establishment of an interaction matrix for OCGS systems based on features, events, and processes. Subsequent risk matrix evaluation and cause–effect analysis identified key system components, specifically CO_(2) injection and faults/features. Building upon this analysis, two leakage risk scenarios were successfully developed, accompanied by the corresponding mitigation measures. In addition, this study introduces the application of scenario development to risk assessment, including scenario numerical simulation and quantitative assessment. Overall, this research positively contributes to the sustainable development and safe operation of OCGS projects and holds potential for further refinement and broader application to diverse geographical environments and project requirements. This comprehensive study provides valuable insights into the establishment of OCGS leakage scenarios and demonstrates their practical application to risk assessment, laying the foundation for promoting the sustainable development and safe operation of ocean CO_(2) geological storage projects while proposing possibilities for future improvements and broader applications to different contexts.

TBAB hydrate formation and growth in a microdevice under static and dynamic conditions

The natural gas hydrate has become one of the most promising future green energy sources on the earth.The natural gas hydrates mostly exist in the sediments with porous structure, so a solid understanding of the hydrate formation and growth processes in the porous medium is of significance for the exploitation of natural gas hydrate. The micro-packed bed device is one of the efficient microfluidic devices in the engineering field, but it has been rarely used for the hydrate-based research. In this study, a transparent micro-packed bed device filled with glass beads was developed to mimic the porous condition of sediments, and used to in-situ visualize the hydrate formation and growth habits in the pore spaces under both static and dynamic conditions. For the static experiment, two types of hydrate growth patterns in porous medium were observed and identified in the micro-packed bed device, which were the graincoating growth and pore-filling growth. For the dynamic condition, the hydrate formation, growth,distribution habits and hydrate blockage phenomena in the pore spaces were in-situ visually captured.The impacts of flowrate and subcooling on the pressure variation of the micro-packed bed and the duration of the hydrate growth under dynamic flow condition in pores were in-situ monitored and analyzed. The higher flowrate could result in the faster hydrate growth and more severe blockage in pores, but the effect of subcooling condition might be less significant at the high flowrate.

A semi-analytical pressure and rate transient analysis model for inner boundary and propped fractures exhibiting dynamic behavior under long-term production conditions

The loss of hydrocarbon production caused by the dynamic behavior of the inner boundary and propped fractures under long-term production conditions has been widely reported in recent studies.However,the quantitative relationships for the variations of the inner boundary and propped fractures have not been determined and incorporated in the semi-analytical models for the pressure and rate transient analysis.This work focuses on describing the variations of the inner boundary and propped fractures and capturing the typical characteristics from the pressure transient curves.A generalized semi-analytical model was developed to characterize the dynamic behavior of the inner boundary and propped fractures under long-term production conditions.The pressure-dependent length shrinkage coefficients,which quantify the length changes of the inner zone and propped fractures,are modified and incorporated into this multi-zone semi-analytical model.With simultaneous numerical iterations and numerical inversions in Laplace and real-time space,the transient solutions to pressure and rate behavior are determined in just a few seconds.The dynamic behavior of the inner boundary and propped fractures on transient pressure curves is divided into five periods:fracture bilinear flow(FR1),dynamic PFs flow(FR2),inner-area linear flow(FR3),dynamic inner boundary flow(FR4),and outer-area dominated linear flow(FR5).The early hump during FR2 period and a positive upward shift during FR4period are captured on the log-log pressure transient curves,reflecting the dynamic behavior of the inner boundary and propped fractures during the long-term production period.The transient pressure behavior will exhibit greater positive upward trend and the flow rate will be lower with the shrinkage of the inner boundary.The pressure derivative curve will be upward earlier as the inner boundary shrinks more rapidly.The lower permeability caused by the closure of un-propped fractures in the inner zone results in greater upward in pressure derivative curves.If the permeability loss for the dynamic behavior of the inner boundary caused by the closure of un-propped fractures is neglected,the flow rate will be overestimated in the later production period.

Production induced fracture closure of deep shale gas well under thermo-hydro-mechanical conditions

Deep shale gas reservoirs have geological characteristics of high temperature,high pressure,high stress,and inferior ability to pass through fluids.The multi-stage fractured horizontal well is the key to exploiting the deep shale gas reservoir.However,during the production process,the effectiveness of the hydraulic fracture network decreases with the closure of fractures,which accelerates the decline of shale gas production.In this paper,we addressed the problems of unclear fracture closure mechanisms and low accuracy of shale gas production prediction during deep shale gas production.Then we established the fluid—solid—heat coupled model coupling the deformation and fluid flow among the fracture surface,proppant and the shale matrix.When the fluid—solid—heat coupled model was applied to the fracture network,it was well solved by our numerical method named discontinuous discrete fracture method.Compared with the conventional discrete fracture method,the discontinuous discrete fracture method can describe the three-dimensional morphology of the fracture while considering the effect of the change of fracture surface permeation coefficient on the coupled fracture—matrix flow and describing the displacement discontinuity across the fracture.Numerical simulations revealed that the degree of fracture closure increases as the production time proceeds,and the degree of closure of the secondary fractures is higher than that of the primary fractures.Shale creep and proppant embedment both increase the degree of fracture closure.The reduction in fracture surface permeability due to proppant embedment reduces the rate of fluid transfer between matrix and fracture,which has often been overlooked in the past.However,it significantly impacts shale gas production,with calculations showing a 24.7%cumulative three-year yield reduction.This study is helpful to understand the mechanism of hydraulic fracture closure.Therefore,it provides the theoretical guidance for maintaining the long-term effectiveness of hydraulic fractures.

Influence of pressure disturbance wave on dynamic response characteristics of liquid film seal for multiphase pump

Slug flow or high GVF(Gas Volume Fraction)conditions can cause pressure disturbance waves and alternating loads at the boundary of mechanical seals for multiphase pumps,endangering the safety of multiphase pump units.The mechanical seal model is simplified by using periodic boundary conditions and numerical calculations are carried out based on the Zwart-Gerber-Belamri cavitation model.UDF(User Define Function)programs such as structural dynamics equations,alternating load equations,and pressure disturbance equations are embedded in numerical calculations,and the dynamic response characteristics of mechanical seal are studied using layered dynamic mesh technology.The results show that when the pressure disturbance occurs at the inlet,as the amplitude and period of the disturbance increase,the film thickness gradually decreases.And the fundamental reason for the hysteresis of the film thickness change is that the pressure in the high-pressure area cannot be restored in a timely manner.The maximum value of leakage and the minimum value of axial velocity are independent of the disturbance period and determined by the disturbance amplitude.The mutual interference between enhanced waves does not have a significant impact on the film thickness,while the front wave in the attenuated wave has a promoting effect on the subsequent film thickness changes,and the fluctuation of the liquid film cavitation rate and axial velocity under the attenuated wave condition deviates from the initial values.Compared with pressure disturbance conditions,alternating load conditions have a more significant impact on film thickness and leakage.During actual operation,it is necessary to avoid alternating load conditions in multiphase pump mechanical seals.

Probabilistic seismic inversion based on physics-guided deep mixture density network

Deterministic inversion based on deep learning has been widely utilized in model parameters estimation.Constrained by logging data,seismic data,wavelet and modeling operator,deterministic inversion based on deep learning can establish nonlinear relationships between seismic data and model parameters.However,seismic data lacks low-frequency and contains noise,which increases the non-uniqueness of the solutions.The conventional inversion method based on deep learning can only establish the deterministic relationship between seismic data and parameters,and cannot quantify the uncertainty of inversion.In order to quickly quantify the uncertainty,a physics-guided deep mixture density network(PG-DMDN)is established by combining the mixture density network(MDN)with the deep neural network(DNN).Compared with Bayesian neural network(BNN)and network dropout,PG-DMDN has lower computing cost and shorter training time.A low-frequency model is introduced in the training process of the network to help the network learn the nonlinear relationship between narrowband seismic data and low-frequency impedance.In addition,the block constraints are added to the PG-DMDN framework to improve the horizontal continuity of the inversion results.To illustrate the benefits of proposed method,the PG-DMDN is compared with existing semi-supervised inversion method.Four synthetic data examples of Marmousi II model are utilized to quantify the influence of forward modeling part,low-frequency model,noise and the pseudo-wells number on inversion results,and prove the feasibility and stability of the proposed method.In addition,the robustness and generality of the proposed method are verified by the field seismic data.