By conducting organic geochemical analysis of the samples taken from the drilled wells in Baiyun Sag of Pearl River Mouth Basin,China,the development characteristics of hydrocarbon source rocks in the sag are clarifie...By conducting organic geochemical analysis of the samples taken from the drilled wells in Baiyun Sag of Pearl River Mouth Basin,China,the development characteristics of hydrocarbon source rocks in the sag are clarified.Reconstruct the current geothermal field of the sag and restore the tectonic-thermal evolution process to predict the type,scale,and distribution of resources in Baiyun Sag through thermal pressure simulation experiments and numerical simulation.The Baiyun Sag is characterized by the development of Paleogene shallow lacustrine source rocks,which are deposited in a slightly oxidizing environment.The source rocks are mainly composed of terrestrial higher plants,with algae making a certain contribution,and are oil and gas source rocks.Current geothermal field of the sag was reconstructed,in which the range of geothermal gradients is(3.5–5.2)℃/100 m,showing an overall increasing trend from northwest to southeast,with significant differences in geothermal gradients across different sub-sags.Baiyun Sag has undergone two distinct periods of extensional process,the Eocene and Miocene,since the Cenozoic era.These two periods of heating and warming events have been identified,accelerating the maturation and evolution of source rocks.The main body of ancient basal heat flow value reached its highest at 13.82 Ma.The basin modelling results show that the maturity of source rocks is significantly higher in Baiyun main sub-sag than that in other sub-sags.The Eocene Wenchang Formation is currently in the stage of high maturity to over maturity,while the Eocene Enping Formation has reached the stage of maturity to high maturity.The rock thermal simulation experiment shows that the shallow lacustrine mudstone of the Wenchang Formation has a good potential of generating gas from kerogen cracking with high gas yield and long period of gas window.Shallow lacustrine mudstone of the Enping Formation has a good ability to generate light oil,and has ability to generate kerogen cracking gas in the late stage.The gas yield of shallow lacustrine mudstone of the Enping Formation is less than that of shallow lacustrine mudstone of the Wenchang Formation and the delta coal-bearing mudstone of the Enping Formation.The numerical simulation results indicate that the source rocks of Baiyun main sub-sag generate hydrocarbons earlier and have significantly higher hydrocarbon generation intensity than other sub-sags,with an average of about 1200×10^(4)t/km^(2).Oil and gas resources were mainly distributed in Baiyun main sub-sag and the main source rocks are distributed in the 3^(rd)and 4^(th)members of Wenchang Formation.Four favorable zones are selected for the division and evaluation of migration and aggregation units:No.(1)Panyu 30 nose-shaped structural belt,No.(3)Liuhua 29 nose-shaped uplift belt and Liwan 3 nose-shaped uplift belt,No.(2)gentle slope belt of Baiyun east sag,and No.(8)Baiyun 1 low-uplift.展开更多
Excessive carbon emissions have resulted in the greenhouse effect, causing considerable global climate change. Marine carbon storage has emerged as a crucial approach to addressing climate change. The Qiantang Sag(QTS...Excessive carbon emissions have resulted in the greenhouse effect, causing considerable global climate change. Marine carbon storage has emerged as a crucial approach to addressing climate change. The Qiantang Sag(QTS) in the East China Sea Shelf Basin, characterized by its extensive area, thick sedimentary strata, and optimal depth, presents distinct geological advantages for carbon dioxide(CO_(2)) storage. Focusing on the lower section of the Shimentan Formation in the Upper Cretaceous of the QTS, this study integrates seismic interpretation and drilling data with core and thin-section analysis. We reveal the vertical variation characteristics of the strata by providing a detailed stratigraphic description. We use petrophysical data to reveal the development characteristics of high-quality carbon-storage layers and favorable reservoircaprock combinations, thereby evaluating the geological conditions for CO_(2) storage in various stratigraphic sections. We identify Layer B of the lower Shimentan Formation as the most advantageous stratum for marine CO_(2) storage. Furthermore, we analyze the carbon emission trends in the adjacent Yangtze River Delta region. Considering the characteristics of the source and sink areas, we suggest a strong correlation between the carbon emission sources of the Yangtze River Delta and the CO_(2) storage area of the QTS, making the latter a priority area for conducting experiments on marine CO_(2) storage.展开更多
Depletion-induced stress change causes the redistribution of stress field in reservoirs,which can lead to the reorientation of principal stresses.Stress reorientation has a direct impact on fracture propagation of inf...Depletion-induced stress change causes the redistribution of stress field in reservoirs,which can lead to the reorientation of principal stresses.Stress reorientation has a direct impact on fracture propagation of infill wells.To understand the effect of stress reorientation on the propagation of infill well’s fractures,an integrated simulation workflow that combines the reservoir flow calculation and the infill well hydraulic fracturing modeling is adopted.The reservoir simulation is computed to examine the relationship between the extent of stress reversal region and reservoir properties.Then,the hydraulic fracturing model considering the altered stress field for production is built to characterize the stress evolution of secondary fracturing.Numerical simulations show that stress reorientation may occur due to the decreasing of the horizontal stresses in an elliptical region around the parent well.Also,the initial stress difference is the driving factor for stress reorientation.However,the bottom hole pressure,permeability and other properties connected with fluid flow control timing of the stress reorientation.The decrease of the horizontal stresses around the parent well lead to asymmetrical propagation of a hydraulic fracture of the infill well.The study provides insights on understanding the influence of stress reorientation to the infill well fracturing treatment and interference between parent and infill wells.展开更多
The molecular composition and biomarker distribution of various occurrences of organic matter in argillaceous source rocks developed in fresh and saline lacustrine environments were revealed by successive treatments o...The molecular composition and biomarker distribution of various occurrences of organic matter in argillaceous source rocks developed in fresh and saline lacustrine environments were revealed by successive treatments of solvent extraction followed by acid hydrolysis using gas chromatography–mass spectrometry. The free fraction obtained by solvent extraction provided abundant geochemical information concerning the sedimentary environment, thermal maturity and biogenic origin, and obvious differences existed between fresh and saline lacustrine source rock samples. Our research results indicate that the carbonate-mineral-bound(CM-bound) fraction released by successive acid hydrolysis could also serve as a significant biogenic indicator, as the bicyclic sesquiterpenoids, indicative of Botryococcus braunii origin, were specifically detected in quite high abundance in the acid-soluble fraction. In addition, the light end hydrocarbons were much better preserved in the acid-soluble fraction, and elemental sulfur was only detected in the CM-bound fraction, suggesting a relatively confined environment for the CM-bound fraction, which thus could preserve additional geochemical information compared to that of the free fraction. The CM-bound fraction also exhibited discernable differences between fresh and saline lacustrine samples. Therefore, it can be concluded that comprehensive analysis of free and CM-bound fractions in the argillaceous source rocks can provide a more authentic and objective interpretation of geologic conditions.展开更多
Organic matter(OM) is preserved as different occurrences in mudstones, which can affect the hydrocarbon generation process. However, little research has focused on hydrocarbon generation as a function of different occ...Organic matter(OM) is preserved as different occurrences in mudstones, which can affect the hydrocarbon generation process. However, little research has focused on hydrocarbon generation as a function of different occurrences of OM. This study collected a suite of mudstones in the Dongying Sag, Bohai Bay Basin, and conducted Rock-Eval Ⅵ pyrolysis after Soxhlet extraction and Na_(2)S_(2)O_(8) oxidation, aiming to quantify the OM with different occurrences and figure out the contributions of each occurrence of OM to the hydrocarbon generation. There are three types of occurrences of OM: soluble organic matter(SOM),mineral-bound organic matter(MOM), and particulate organic matter(POM). MOM is the most abundant among the three occurrence types of OM. SOM and MOM are the main hydrocarbon precursors, and their hydrocarbon contributions alternate with different kerogen types and layers. Additionally, MOMcontributed hydrocarbons are numerous at shallow depths;SOM-contributed hydrocarbons mainly occur at deep depths;and POM-contributed hydrocarbons change little with depth. These results demonstrate that MOM should be the main hydrocarbon precursor in shallow formations and that SOM is the main hydrocarbon contributor at deep depths.展开更多
According to the variable toe-to-heel well spacing, combined with the dislocation theory, discrete lattice method, and finite-element-method(FEM) based fluid-solid coupling, an integrated geological-engineering method...According to the variable toe-to-heel well spacing, combined with the dislocation theory, discrete lattice method, and finite-element-method(FEM) based fluid-solid coupling, an integrated geological-engineering method of volume fracturing for fan-shaped well pattern is proposed considering the geomechanical modeling, induced stress calculation, hydraulic fracturing simulation, and post-frac productivity evaluation. Besides, we propose the differential fracturing design for the conventional productivity-area and the potential production area for fan-shaped horizontal wells. After the fracturing of the conventional production area for H1 fan-shaped well platform, the research shows that the maximum reduction of the horizontal principal stress difference in the potential productivity-area is 0.2 MPa, which cannot cause the stress reversal, but this reduction is still conducive to the lateral propagation of hydraulic fractures. According to the optimized fracturing design, in zone-Ⅰ of the potential production area, only Well 2 is fractured, with a cluster spacing of 30 m and an injection rate of 12 m^(3)/min per stage;in zone-Ⅱ, Well 2 is fractured before Well 3, with a cluster spacing of 30 m and an injection rate of 12 m^(3)/min per stage. The swept area of the pore pressure drop in the potential production area is small, showing that the reservoir is not well developed. The hydraulic fracturing in the toe area can be improved by, for example, properly densifying the fractures and adjusting the fracture distribution, in order to enhance the swept volume and increase the reservoir utilization.展开更多
Reflection-based inversion that aims to reconstruct the low-to-intermediate wavenumbers of the subsurface model, can be a complementary to refraction-data-driven full-waveform inversion(FWI), especially for the deep t...Reflection-based inversion that aims to reconstruct the low-to-intermediate wavenumbers of the subsurface model, can be a complementary to refraction-data-driven full-waveform inversion(FWI), especially for the deep target area where diving waves cannot be acquired at the surface. Nevertheless, as a typical nonlinear inverse problem, reflection waveform inversion may easily suffer from the cycleskipping issue and have a slow convergence rate, if gradient-based first-order optimization methods are used. To improve the accuracy and convergence rate, we introduce the Hessian operator into reflection traveltime inversion(RTI) and reflection waveform inversion(RWI) in the framework of second-order optimization. A practical two-stage workflow is proposed to build the velocity model, in which Gauss-Newton RTI is first applied to mitigate the cycle-skipping problem and then Gauss-Newton RWI is employed to enhance the model resolution. To make the Gauss-Newton iterations more efficiently and robustly for large-scale applications, we introduce proper preconditioning for the Hessian matrix and design appropriate strategies to reduce the computational costs. The example of a real dataset from East China Sea demonstrates that the cascaded Hessian-based RTI and RWI have good potential to improve velocity model building and seismic imaging, especially for the deep targets.展开更多
Reliable prediction of the shale fracturing process is a challenging problem in exploiting deep shale oil and gas resources.Complex fracture networks need to be artificially created to employ deep shale oil and gas re...Reliable prediction of the shale fracturing process is a challenging problem in exploiting deep shale oil and gas resources.Complex fracture networks need to be artificially created to employ deep shale oil and gas reserves.Randomly distributed minerals and heterogeneities in shales significantly affect mechanical properties and fracturing behaviors in oil and gas exploitation.Describing the actual microstructure and associated heterogeneities in shales constitutes a significant challenge.The RFPA3D(rock failure process analysis parallel computing program)-based modeling approach is a promising numerical technique due to its unique capability to simulate the fracturing behavior of rocks.To improve traditional numerical technology and study crack propagation in shale on the microscopic scale,a combination of high-precision internal structure detection technology with the RFPA^(3D) numerical simulation method was developed to construct a real mineral structure-based modeling method.First,an improved digital image processing technique was developed to incorporate actual shale microstructures(focused ion beam scanning electron microscopy was used to capture shale microstructure images that reflect the distri-butions of different minerals)into the numerical model.Second,the effect of mineral inhomogeneity was considered by integrating the mineral statistical model obtained from the mineral nanoindentation experiments into the numerical model.By simulating a shale numerical model in which pyrite particles are wrapped by organic matter,the effects of shale microstructure and applied stress state on microcrack behavior and mechanical properties were investigated and analyzed.In this study,the effect of pyrite particles on fracture propagation was systematically analyzed and summarized for the first time.The results indicate that the distribution of minerals and initial defects dominated the fracture evolution and the failure mode.Cracks are generally initiated and propagated along the boundaries of hard mineral particles such as pyrite or in soft minerals such as organic matter.Locations with collections of hard minerals are more likely to produce complex fractures.This study provides a valuable method for un-derstanding the microfracture behavior of shales.展开更多
An improved preprocessed Yaroslavsky filter（IPYF）is proposed to avoid the nick effects and obtain a better denoising result when the noise variance is unknown.Different from its predecessors,the similarity between t...An improved preprocessed Yaroslavsky filter（IPYF）is proposed to avoid the nick effects and obtain a better denoising result when the noise variance is unknown.Different from its predecessors,the similarity between two pixels is calculated by shearlet features.The feature vector consists of initial denoised results by the non-subsampled shearlet transform hard thresholding（NSST-HT）and NSST coefficients,which can help allocate the averaging weights more reasonably.With the correct estimated noise variance,the NSST-HT can provide good denoised results as the initial estimation and high-frequency coefficients contribute large weights to preserve textures.In case of the incorrect estimated noise variance,the low-frequency coefficients will mitigate the nick effect in cartoon regions greatly,making the IPYF more robust than the original PYF.Detailed experimental results show that the IPYF is a very competitive method based on a comprehensive consideration involving peak signal to noise ratio（PSNR）,computing time,visual quality and method noise.展开更多
基金Supported by the National Oil and Gas Resource Evaluation Project for the 14th Five Year Plan of the Ministry of Natural Resources(QGYQZYPJ2022-3)China National Offshore Oil Corporation"14th Five Year Plan"Major Science and Technology Project(KJGG2022-0103-03)。
文摘By conducting organic geochemical analysis of the samples taken from the drilled wells in Baiyun Sag of Pearl River Mouth Basin,China,the development characteristics of hydrocarbon source rocks in the sag are clarified.Reconstruct the current geothermal field of the sag and restore the tectonic-thermal evolution process to predict the type,scale,and distribution of resources in Baiyun Sag through thermal pressure simulation experiments and numerical simulation.The Baiyun Sag is characterized by the development of Paleogene shallow lacustrine source rocks,which are deposited in a slightly oxidizing environment.The source rocks are mainly composed of terrestrial higher plants,with algae making a certain contribution,and are oil and gas source rocks.Current geothermal field of the sag was reconstructed,in which the range of geothermal gradients is(3.5–5.2)℃/100 m,showing an overall increasing trend from northwest to southeast,with significant differences in geothermal gradients across different sub-sags.Baiyun Sag has undergone two distinct periods of extensional process,the Eocene and Miocene,since the Cenozoic era.These two periods of heating and warming events have been identified,accelerating the maturation and evolution of source rocks.The main body of ancient basal heat flow value reached its highest at 13.82 Ma.The basin modelling results show that the maturity of source rocks is significantly higher in Baiyun main sub-sag than that in other sub-sags.The Eocene Wenchang Formation is currently in the stage of high maturity to over maturity,while the Eocene Enping Formation has reached the stage of maturity to high maturity.The rock thermal simulation experiment shows that the shallow lacustrine mudstone of the Wenchang Formation has a good potential of generating gas from kerogen cracking with high gas yield and long period of gas window.Shallow lacustrine mudstone of the Enping Formation has a good ability to generate light oil,and has ability to generate kerogen cracking gas in the late stage.The gas yield of shallow lacustrine mudstone of the Enping Formation is less than that of shallow lacustrine mudstone of the Wenchang Formation and the delta coal-bearing mudstone of the Enping Formation.The numerical simulation results indicate that the source rocks of Baiyun main sub-sag generate hydrocarbons earlier and have significantly higher hydrocarbon generation intensity than other sub-sags,with an average of about 1200×10^(4)t/km^(2).Oil and gas resources were mainly distributed in Baiyun main sub-sag and the main source rocks are distributed in the 3^(rd)and 4^(th)members of Wenchang Formation.Four favorable zones are selected for the division and evaluation of migration and aggregation units:No.(1)Panyu 30 nose-shaped structural belt,No.(3)Liuhua 29 nose-shaped uplift belt and Liwan 3 nose-shaped uplift belt,No.(2)gentle slope belt of Baiyun east sag,and No.(8)Baiyun 1 low-uplift.
基金Key Laboratory of Deep-time Geography and Environment Reconstruction and Applications of Ministry of Natural ResourcesChengdu University of Technology:DGERA20231110。
文摘Excessive carbon emissions have resulted in the greenhouse effect, causing considerable global climate change. Marine carbon storage has emerged as a crucial approach to addressing climate change. The Qiantang Sag(QTS) in the East China Sea Shelf Basin, characterized by its extensive area, thick sedimentary strata, and optimal depth, presents distinct geological advantages for carbon dioxide(CO_(2)) storage. Focusing on the lower section of the Shimentan Formation in the Upper Cretaceous of the QTS, this study integrates seismic interpretation and drilling data with core and thin-section analysis. We reveal the vertical variation characteristics of the strata by providing a detailed stratigraphic description. We use petrophysical data to reveal the development characteristics of high-quality carbon-storage layers and favorable reservoircaprock combinations, thereby evaluating the geological conditions for CO_(2) storage in various stratigraphic sections. We identify Layer B of the lower Shimentan Formation as the most advantageous stratum for marine CO_(2) storage. Furthermore, we analyze the carbon emission trends in the adjacent Yangtze River Delta region. Considering the characteristics of the source and sink areas, we suggest a strong correlation between the carbon emission sources of the Yangtze River Delta and the CO_(2) storage area of the QTS, making the latter a priority area for conducting experiments on marine CO_(2) storage.
基金the support provided by the Scientific Research and Technology Development Project of CNPC(Grant No.kt2017-19-01-1)the National Natural Science Foundation of China(Grant No.41772286,No.42077247 and No.42002271)+2 种基金Petro China Innovation Foundation(Grant No.2018D-5007-0202)Project funded by China Postdoctoral Science Foundation(Grant No.2021T140514)Fund of State Key Laboratory of Geomechanics and Geotechnical Engineering,Institute of Rock and Soil Mechanics,Chinese Academy of Sciences(Grant No.Z020009)。
文摘Depletion-induced stress change causes the redistribution of stress field in reservoirs,which can lead to the reorientation of principal stresses.Stress reorientation has a direct impact on fracture propagation of infill wells.To understand the effect of stress reorientation on the propagation of infill well’s fractures,an integrated simulation workflow that combines the reservoir flow calculation and the infill well hydraulic fracturing modeling is adopted.The reservoir simulation is computed to examine the relationship between the extent of stress reversal region and reservoir properties.Then,the hydraulic fracturing model considering the altered stress field for production is built to characterize the stress evolution of secondary fracturing.Numerical simulations show that stress reorientation may occur due to the decreasing of the horizontal stresses in an elliptical region around the parent well.Also,the initial stress difference is the driving factor for stress reorientation.However,the bottom hole pressure,permeability and other properties connected with fluid flow control timing of the stress reorientation.The decrease of the horizontal stresses around the parent well lead to asymmetrical propagation of a hydraulic fracture of the infill well.The study provides insights on understanding the influence of stress reorientation to the infill well fracturing treatment and interference between parent and infill wells.
基金funded by NSFC (National Science Foundation of China) No. 41372130
文摘The molecular composition and biomarker distribution of various occurrences of organic matter in argillaceous source rocks developed in fresh and saline lacustrine environments were revealed by successive treatments of solvent extraction followed by acid hydrolysis using gas chromatography–mass spectrometry. The free fraction obtained by solvent extraction provided abundant geochemical information concerning the sedimentary environment, thermal maturity and biogenic origin, and obvious differences existed between fresh and saline lacustrine source rock samples. Our research results indicate that the carbonate-mineral-bound(CM-bound) fraction released by successive acid hydrolysis could also serve as a significant biogenic indicator, as the bicyclic sesquiterpenoids, indicative of Botryococcus braunii origin, were specifically detected in quite high abundance in the acid-soluble fraction. In addition, the light end hydrocarbons were much better preserved in the acid-soluble fraction, and elemental sulfur was only detected in the CM-bound fraction, suggesting a relatively confined environment for the CM-bound fraction, which thus could preserve additional geochemical information compared to that of the free fraction. The CM-bound fraction also exhibited discernable differences between fresh and saline lacustrine samples. Therefore, it can be concluded that comprehensive analysis of free and CM-bound fractions in the argillaceous source rocks can provide a more authentic and objective interpretation of geologic conditions.
基金financially supported by the National Natural Science Foundation of China (Grant Nos. 41672115 and 41972126)the National Science and Technology Major Project of China(Grant No. 2016ZX05006001-003)。
文摘Organic matter(OM) is preserved as different occurrences in mudstones, which can affect the hydrocarbon generation process. However, little research has focused on hydrocarbon generation as a function of different occurrences of OM. This study collected a suite of mudstones in the Dongying Sag, Bohai Bay Basin, and conducted Rock-Eval Ⅵ pyrolysis after Soxhlet extraction and Na_(2)S_(2)O_(8) oxidation, aiming to quantify the OM with different occurrences and figure out the contributions of each occurrence of OM to the hydrocarbon generation. There are three types of occurrences of OM: soluble organic matter(SOM),mineral-bound organic matter(MOM), and particulate organic matter(POM). MOM is the most abundant among the three occurrence types of OM. SOM and MOM are the main hydrocarbon precursors, and their hydrocarbon contributions alternate with different kerogen types and layers. Additionally, MOMcontributed hydrocarbons are numerous at shallow depths;SOM-contributed hydrocarbons mainly occur at deep depths;and POM-contributed hydrocarbons change little with depth. These results demonstrate that MOM should be the main hydrocarbon precursor in shallow formations and that SOM is the main hydrocarbon contributor at deep depths.
基金Supported by National Natural Science Foundation of China (52104029,U2139204)PetroChina Science and Technology Innovation Foundation (2021 DQ02-0501)。
文摘According to the variable toe-to-heel well spacing, combined with the dislocation theory, discrete lattice method, and finite-element-method(FEM) based fluid-solid coupling, an integrated geological-engineering method of volume fracturing for fan-shaped well pattern is proposed considering the geomechanical modeling, induced stress calculation, hydraulic fracturing simulation, and post-frac productivity evaluation. Besides, we propose the differential fracturing design for the conventional productivity-area and the potential production area for fan-shaped horizontal wells. After the fracturing of the conventional production area for H1 fan-shaped well platform, the research shows that the maximum reduction of the horizontal principal stress difference in the potential productivity-area is 0.2 MPa, which cannot cause the stress reversal, but this reduction is still conducive to the lateral propagation of hydraulic fractures. According to the optimized fracturing design, in zone-Ⅰ of the potential production area, only Well 2 is fractured, with a cluster spacing of 30 m and an injection rate of 12 m^(3)/min per stage;in zone-Ⅱ, Well 2 is fractured before Well 3, with a cluster spacing of 30 m and an injection rate of 12 m^(3)/min per stage. The swept area of the pore pressure drop in the potential production area is small, showing that the reservoir is not well developed. The hydraulic fracturing in the toe area can be improved by, for example, properly densifying the fractures and adjusting the fracture distribution, in order to enhance the swept volume and increase the reservoir utilization.
基金supported by National Natural Science Foundation of China (42074157)the National Key Research and Development Program of China (2018YFC0310104)the Strategic Priority Research Program of the Chinese Academy of Science(XDA14010203)。
文摘Reflection-based inversion that aims to reconstruct the low-to-intermediate wavenumbers of the subsurface model, can be a complementary to refraction-data-driven full-waveform inversion(FWI), especially for the deep target area where diving waves cannot be acquired at the surface. Nevertheless, as a typical nonlinear inverse problem, reflection waveform inversion may easily suffer from the cycleskipping issue and have a slow convergence rate, if gradient-based first-order optimization methods are used. To improve the accuracy and convergence rate, we introduce the Hessian operator into reflection traveltime inversion(RTI) and reflection waveform inversion(RWI) in the framework of second-order optimization. A practical two-stage workflow is proposed to build the velocity model, in which Gauss-Newton RTI is first applied to mitigate the cycle-skipping problem and then Gauss-Newton RWI is employed to enhance the model resolution. To make the Gauss-Newton iterations more efficiently and robustly for large-scale applications, we introduce proper preconditioning for the Hessian matrix and design appropriate strategies to reduce the computational costs. The example of a real dataset from East China Sea demonstrates that the cascaded Hessian-based RTI and RWI have good potential to improve velocity model building and seismic imaging, especially for the deep targets.
基金supported by the Central Program of Basic Science of the National Natural Science Foundation of China(No.72088101)"The theory and application of resource and environment management in the digital economy era"+1 种基金The National Natural Science Foundation of China(No.41941018)Scientific research and technological development program of RIPED,"major research of basic geologic and synergy research of engineering practice on Gulong shale oil"(No.2021ycq01).
文摘Reliable prediction of the shale fracturing process is a challenging problem in exploiting deep shale oil and gas resources.Complex fracture networks need to be artificially created to employ deep shale oil and gas reserves.Randomly distributed minerals and heterogeneities in shales significantly affect mechanical properties and fracturing behaviors in oil and gas exploitation.Describing the actual microstructure and associated heterogeneities in shales constitutes a significant challenge.The RFPA3D(rock failure process analysis parallel computing program)-based modeling approach is a promising numerical technique due to its unique capability to simulate the fracturing behavior of rocks.To improve traditional numerical technology and study crack propagation in shale on the microscopic scale,a combination of high-precision internal structure detection technology with the RFPA^(3D) numerical simulation method was developed to construct a real mineral structure-based modeling method.First,an improved digital image processing technique was developed to incorporate actual shale microstructures(focused ion beam scanning electron microscopy was used to capture shale microstructure images that reflect the distri-butions of different minerals)into the numerical model.Second,the effect of mineral inhomogeneity was considered by integrating the mineral statistical model obtained from the mineral nanoindentation experiments into the numerical model.By simulating a shale numerical model in which pyrite particles are wrapped by organic matter,the effects of shale microstructure and applied stress state on microcrack behavior and mechanical properties were investigated and analyzed.In this study,the effect of pyrite particles on fracture propagation was systematically analyzed and summarized for the first time.The results indicate that the distribution of minerals and initial defects dominated the fracture evolution and the failure mode.Cracks are generally initiated and propagated along the boundaries of hard mineral particles such as pyrite or in soft minerals such as organic matter.Locations with collections of hard minerals are more likely to produce complex fractures.This study provides a valuable method for un-derstanding the microfracture behavior of shales.
基金Supported by Open Fund of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation(Southwest Petroleum University)(PL N1303)Open Fund of State Key Laboratory of Marine Geology(Tongji University)(MGK1412)+1 种基金Fundation of Graduate Innovation Center in NUAA(kfjj201430)the Fundamental Research Funds for the Central Universities
文摘An improved preprocessed Yaroslavsky filter（IPYF）is proposed to avoid the nick effects and obtain a better denoising result when the noise variance is unknown.Different from its predecessors,the similarity between two pixels is calculated by shearlet features.The feature vector consists of initial denoised results by the non-subsampled shearlet transform hard thresholding（NSST-HT）and NSST coefficients,which can help allocate the averaging weights more reasonably.With the correct estimated noise variance,the NSST-HT can provide good denoised results as the initial estimation and high-frequency coefficients contribute large weights to preserve textures.In case of the incorrect estimated noise variance,the low-frequency coefficients will mitigate the nick effect in cartoon regions greatly,making the IPYF more robust than the original PYF.Detailed experimental results show that the IPYF is a very competitive method based on a comprehensive consideration involving peak signal to noise ratio（PSNR）,computing time,visual quality and method noise.
