Alkaline lacustrine shale is highly heterogeneous,and the complex relationship between the organicinorganic porosity network and hydrocarbon occurrence restricts the effectiveness of shale oil exploration and developm...Alkaline lacustrine shale is highly heterogeneous,and the complex relationship between the organicinorganic porosity network and hydrocarbon occurrence restricts the effectiveness of shale oil exploration and development.Herein,we investigated the Fengcheng Formation(P_(1)f)in Mahu Sag.This study integrated geochemistry,Soxhlet extraction,scanning electron microscopy,gas adsorption,and nuclear magnetic resonance T_(1)-T_(2)spectroscopy to elucidate the microscopic oil occurrence mechanisms in shales.Results indicate the presence of felsic shale,dolomitic shale,lime shale,and mixed shale within the P_(1)f.Matrix pores and microfractures associated with inorganic minerals are the predominant pore types in P_(1)f.Adsorbed oil primarily resides on the surfaces of organic matter and clay minerals,while free oil predominantly occupies inorganic pores and microfractures with larger pore sizes.Variations exist in the quantity and distribution of shale oil accumulation across different scales,where free oil and adsorbed oil are governed by dominant pores with diameters exceeding 10 nm and ineffective pores with diameters below 10 nm,respectively.Shale oil occurrence characteristics are influenced by organic matter,pore structure,and mineral composition.Felsic shale exhibits a high abundance of dominant pores,possesses the highest oil content,predominantly harbors free oil within these dominant pores,and demonstrates good mobility.Fluid occurrence in dolomitic shale and lime shale is intricate,with low oil content and a free oil to adsorbed oil ratio of 1:1.Mixed shale exhibits elevated clay mineral content and a scarcity of dominant pores.Moreover,ineffective pores contain increased bound water,resulting in medium oil content and limited mobility predominantly due to adsorption.Presently,shale oil mainly occurs in the dominant pores with a diameter larger than 10 nm in a free state.During the exploration and development of alkaline lacustrine shale oil resources,emphasis should be placed on identifying sweet spots within the felsic shale characterized by dominant pores.展开更多
The use of blended acquisition technology in marine seismic exploration has the advantages of high acquisition efficiency and low exploration costs.However,during acquisition,the primary source may be disturbed by adj...The use of blended acquisition technology in marine seismic exploration has the advantages of high acquisition efficiency and low exploration costs.However,during acquisition,the primary source may be disturbed by adjacent sources,resulting in blended noise that can adversely affect data processing and interpretation.Therefore,the de-blending method is needed to suppress blended noise and improve the quality of subsequent processing.Conventional de-blending methods,such as denoising and inversion methods,encounter challenges in parameter selection and entail high computational costs.In contrast,deep learning-based de-blending methods demonstrate reduced reliance on manual intervention and provide rapid calculation speeds post-training.In this study,we propose a Uformer network using a nonoverlapping window multihead attention mechanism designed for de-blending blended data in the common shot domain.We add the depthwise convolution to the feedforward network to improve Uformer’s ability to capture local context information.The loss function comprises SSIM and L1 loss.Our test results indicate that the Uformer outperforms convolutional neural networks and traditional denoising methods across various evaluation metrics,thus highlighting the effectiveness and advantages of Uformer in de-blending blended data.展开更多
In the leakage detection of reservoir dam bodies,traditional geophysical methods can only achieve one-time detection.Meanwhile,due to the non-uniqueness of geophysicalin version,how to improve the fidelity of geophysi...In the leakage detection of reservoir dam bodies,traditional geophysical methods can only achieve one-time detection.Meanwhile,due to the non-uniqueness of geophysicalin version,how to improve the fidelity of geophysical detection inversion profiles has become a key challenge in the industry.This study aims to construct a long-term real-time monitoring system using the time-lapse resistivity method,reveal the spatiotemporal evolution law of dam leakage,and provide technical support for accurate treatment.By integrating the Internet of Things(IoT),5G technology,and AI technology,real-time data acquisition,real-time transmission,and automatic inversionare realized.Through dynamic imaging analysis of the electrical anomaly characteristics of the leakage area and comparison between corresponding rainfall events,the leakage range and resistivity changes,the reliability and efficiency of dam leakagedete ction are significantly improved.This achieves long-termdynamic monitoring of dam leakage and provides a new perspective for the safe operation and maintenance of reservoirs.展开更多
Monitoring and delineating the spatial distribution of shale fracturing is fundamentally important to shale gas production. Standard monitoring methods, such as time-lapse seismic, cross-well seismic and micro-seismic...Monitoring and delineating the spatial distribution of shale fracturing is fundamentally important to shale gas production. Standard monitoring methods, such as time-lapse seismic, cross-well seismic and micro-seismic methods, are expensive, time- consuming, and do not show the changes in the formation with time. The resistivities of hydraulic fracturing fluid and reservoir rocks were measured. The results suggest that the injection fluid and consequently the injected reservoir are characterized by very low resistivity and high chargeability. This allows using of the controlled-source electromagnetic method （CSEM） to monitor shale gas hydraulic fracturing. Based on the geoelectrical model which was proposed according to the well-log and seismic data in the test area the change rule of the reacted electrical field was studied to account for the change of shale resistivity, and then the normalized residual resistivity method for time lapse processing was given. The time-domain electromagnetic method （TDEM） was used to continuously monitor the shale gas fracturing at the Fulin shale gas field in southern China. A high-power transmitter and multi-channel transient electromagnetic receiver array were adopted. 9 h time series of Ex component of 224 sites which were laid out on the surface and over three fracturing stages of a horizontal well at 2800 m depth was recorded. After data processing and calculation of the normalized resistivity residuals, the changes in the Ex signal were determined and a dynamic 3D image of the change in resistivity was constructed. This allows modeling the spatial distribution of the fracturing fluid. The model results suggest that TDEM is promising for monitoring hydraulic fracturing of shale.展开更多
Multistage deformation events have occurred in the northeastern Jiangshao Fault （Suture） Belt. The earliest two are ductile deformation events. The first is the ca. 820 Ma top-to-the-northwest ductile thrusting, whi...Multistage deformation events have occurred in the northeastern Jiangshao Fault （Suture） Belt. The earliest two are ductile deformation events. The first is the ca. 820 Ma top-to-the-northwest ductile thrusting, which directly resulted from the collision between the Cathaysia Old Land and the Chencai Arc （？） during the Late Neoproterozoic, and the Jiangnan Orogenic Belt that formed as the ocean closed between the Yangtze Plate and the jointed Cathaysia Old Land and the Chencai Arc due to continuous compression. The second is the ductile left-lateral strike-slipping that occurred in the latest Early Paleozoic. Since the Jinning period, all deformation events represent the reactivation or inversion of intraplate structures due to the collisions between the North China and Yangtze plates during the Triassic and between the Philippine Sea and Eurasian plates during the Cenozoic. In the Triassic, brittle right-lateral strike-slipping and subsequent top-to-the south thrusting occurred along the whole northeastern Jiangshao Fault Zone because of the collision between the North China and Yangtze plates. In the Late Mesozoic, regional extension took place across southeastern China. In the Cenozoic, the collision between the Philippine Sea and Eurasian plates resulted in brittle thrusts along the whole Jiangnan Old land in the Miocene. The Jiangshao Fault Belt is a weak zone in the crust with long history, and its reactivation is one of important characteristics of the deformation in South China; however, late-stage deformation events did not occur beyond the Jiangnan Old Land and most of them are parallel to the strike of the Old Land, which is similar to the Cenozoic deformation in Central Asia. In addition, the Jiangnan old Land is not a collisional boundary between the Yangtze Plate and Cathaysia Old Land in the Triassic.展开更多
Petroleum geophysicists recognize that many parameters related to oil and gas reservoirs are predicted using seismic attribute data. However, how best to optimize the seismic attributes, predict the character of thin ...Petroleum geophysicists recognize that many parameters related to oil and gas reservoirs are predicted using seismic attribute data. However, how best to optimize the seismic attributes, predict the character of thin sandstone reservoirs, and enhance the reservoir description accuracy is an important goal for geologists and geophysicists. Based on the theory of main component analysis, we present a new optimization method, called constrained main component analysis. Modeling estimates and real application in an oilfield show that it can enhance reservoir prediction accuracy and has better applicability.展开更多
3D inversion of borehole-surface electrical data for complex geo-electrical models is still a challenging problem in geophysical exploration. We have developed a program for 3D inversion to borehole-surface electrical...3D inversion of borehole-surface electrical data for complex geo-electrical models is still a challenging problem in geophysical exploration. We have developed a program for 3D inversion to borehole-surface electrical data based on the quasi-analytical approximation （QA） and re-weighted regularized conjugate gradient method （RRCG） algorithms using Visual Fortran 6.5. Application of the QA approximation to forward modeling and Frechet derivative computations speeds up the calculation dramatically. The trial calculation for synthetic data of theoretical model showed that the program is fast and highly precise.展开更多
The Triassic reservoir in the Jinan area of Tarim Oilfield consists largely of interbedded sand and shale. Because of the large overlap between sandstone and shale impedance, it is difficult to distinguish sandstone f...The Triassic reservoir in the Jinan area of Tarim Oilfield consists largely of interbedded sand and shale. Because of the large overlap between sandstone and shale impedance, it is difficult to distinguish sandstone from shale by acoustic impedance alone. Compared to acoustic impedance, elastic impedance contains more lithologic and physical information of the reservoir. Based on meticulous well-tie calibration, elastic impedance data volumes for 10°, 20°, and 30° emergence angles are obtained using pre-stack elastic impedance inversion. A non-linear statistical relationship between elastic impedance and shale content is set up by a PNN neural network. The non-linear mapping relationship is used to predict the reservoir shale content from elastic impedance, which will depict and predict the reservoir oil-bearing sands.展开更多
Secondary/minor structures occurring along the main fault surfaces are important indicators for judging the kinematic characteristics of faults.However,many factors can lead to the formation of these structures,which ...Secondary/minor structures occurring along the main fault surfaces are important indicators for judging the kinematic characteristics of faults.However,many factors can lead to the formation of these structures,which results in the difficulty for rapid judgment and application in the fields.A series of secondary faults/fractures developed due to the movement of main faults are the most important and widespread phenomena in the scope of brittle deformation.The morphology of the main fault surfaces is various,and former researchers mainly discussed the structures on the main even fault surfaces. However,the fluctuation of fault surfaces is the intrinsic character of the faults,and the intersection between the main fault and secondary faults/fractures can produce a series of kinematic indicators on the main fault surfaces.Based on previous studies and our observations,i.e.the structural traces of the P,R,R',T and X shears/faults along the main faults,some indicators which are rarely reported previously,are described in the paper.Furthermore,their reliabilities are also discussed,and more practical and reliable criteria are brought forward.We suggest that the simple application of congruous and incongruous steps without knowing their exact origins should be abandoned in the fields,and several types of indicators along one fault surface should be checked with each other as much as possible.Meanwhile,the origins of some other arcuate indicators on the fault surfaces are also discussed,and new models are brought forward.展开更多
Oil and gas seismic exploration have to adopt irregular seismic acquisition due to the increasingly complex exploration conditions to adapt to complex geological conditions and environments.However,the irregular seism...Oil and gas seismic exploration have to adopt irregular seismic acquisition due to the increasingly complex exploration conditions to adapt to complex geological conditions and environments.However,the irregular seismic acquisition is accompanied by the lack of acquisition data,which requires high-precision regularization.The sparse signal feature in the transform domain in compressed sensing theory is used in this paper to recover the missing signal,involving sparse transform base optimization and threshold modeling.First,this paper analyzes and compares the effects of six sparse transformation bases on the reconstruction accuracy and efficiency of irregular seismic data and establishes the quantitative relationship between sparse transformation and reconstruction accuracy and efficiency.Second,an adaptive threshold modeling method based on sparse coefficient is provided to improve the reconstruction accuracy.Test results show that the method has good adaptability to different seismic data and sparse transform bases.The f-x domain reconstruction method of effective frequency samples is studied to address the problem of low computational efficiency.The parallel computing strategy of curvelet transform combined with OpenMP is further proposed,which substantially improves the computational efficiency under the premise of ensuring the reconstruction accuracy.Finally,the actual acquisition data are used to verify the proposed method.The results indicate that the proposed method strategy can solve the regularization problem of irregular seismic data in production and improve the imaging quality of the target layer economically and efficiently.展开更多
The joint source scanning algorithm(SSA)scans locations and focal mechanisms of microseismic events simultaneously.Compared to the traditional source scanning algorithm,it yields much more events with extra informatio...The joint source scanning algorithm(SSA)scans locations and focal mechanisms of microseismic events simultaneously.Compared to the traditional source scanning algorithm,it yields much more events with extra information of focal mechanisms.The availability of more events and focal mechanisms make it possible to invert for a 2D gridded stress field.As a byproduct of hydrofracturing monitoring,the method offers a new way to extract stress field as a substitute to other more expensive technologies.This method is applied to a hydraulic fracturing dataset collected from one shale gas production field in the southeast of the Sichuan basin.A damped stress inversion is conducted to obtain a 2D stress field.five hydraulic-fracturing induced fractures can be determined from the result.The events associated with these fractures generally have relatively low stacked energy and are limited to the depth of horizontal well.One existing fault(possibly associated with the axis of the central Sichuan uplift)is also determined and the events associated with the existing fault generally have higher stacked energy and are more densely populated.The existing fault may also serve as a structural boundary where the rocks to the NW side are easier to be fractured while events on the other side are sparse with low stacked energy.The existing fault also divides the stress field into two regimes:the maximum compressional stress field to the NW and SE of the fault line are dominantly in NW-SE and N-S directions,respectively.展开更多
We present a 3D approach to numerical modeling of the borehole-surface electromagnetic （BSEM） method. The 3D electromagnetic response created by a vertical line current source in a layered medium is modeled using th...We present a 3D approach to numerical modeling of the borehole-surface electromagnetic （BSEM） method. The 3D electromagnetic response created by a vertical line current source in a layered medium is modeled using the 3D integral equation method. The modeling results are consistent with analytical solutions. 3D Born approximation inversion of BSEM data is also conducted for reservoir delineation. The inversion method is verified by a synthetic reservoir model.展开更多
Conventional land vertical seismic profiling （VSP） exploration usually uses P-wave sources and three-component geophones for receivers, emphasizing P- and converted S-waves. Previous studies show that both dynamite ...Conventional land vertical seismic profiling （VSP） exploration usually uses P-wave sources and three-component geophones for receivers, emphasizing P- and converted S-waves. Previous studies show that both dynamite borehole shots and vertical vibrations from controllable seismic sources at the surface will produce relatively strong pure P-waves and weaker pure S-waves. Interfaces with a large Poisson＇s ratio difference have a positive influence on the formation of strong transmitted converted S-waves. By a comparative analysis of pure S-waves from sources and converted downgoing S-waves, we believe that the main frequency of pure S-waves is usually lower than pure P-waves while the main frequency of downgoing converted S-waves is close to that of P-waves. We have studied zero-offset and offset VSP data from land P-wave sources. Results show that pure S-waves commonly exist in these data with differences in wave intensity. S-wave velocity can be obtained from the P-wave source zero-offset VSP data. Finally, we discuss the bright future of joint application of VSP P-and S-waves and the full use of S-waves in P-wave source VSP data.展开更多
Seismic velocity is important to migration of seismic data, interpretation of lithology and lithofacies as well as prediction of reservoir. The information of shear wave velocity is required to reduce the uncertainty ...Seismic velocity is important to migration of seismic data, interpretation of lithology and lithofacies as well as prediction of reservoir. The information of shear wave velocity is required to reduce the uncertainty for discriminating lithology, identifying fluid type in porous material and calculating gas saturation in reservoir prediction. Based on Zoeppritz equations, a numeral and scanning method was proposed in this paper. Shear wave velocity can be calculated with prestack converted wave data. The effects were demonstrated by inversion of theoretical and real seismic data.展开更多
Multi-wave exploration is an effective means for improving precision in the exploration and development of complex oil and gas reservoirs that are dense and have low permeability. However, convened wave data is charac...Multi-wave exploration is an effective means for improving precision in the exploration and development of complex oil and gas reservoirs that are dense and have low permeability. However, convened wave data is characterized by a low signal-to-noise ratio and low resolution, because the conventional deconvolution technology is easily affected by the frequency range limits, and there is limited scope for improving its resolution. The spectral inversion techniques is used to identify λ/8 thin layers and its breakthrough regarding band range limits has greatly improved the seismic resolution. The difficulty associated with this technology is how to use the stable inversion algorithm to obtain a high-precision reflection coefficient, and then to use this reflection coefficient to reconstruct broadband data for processing. In this paper, we focus on how to improve the vertical resolution of the converted PS-wave for multi-wave data processing. Based on previous research, we propose a least squares inversion algorithm with a total variation constraint, in which we uses the total variance as a priori information to solve under-determined problems, thereby improving the accuracy and stability of the inversion. Here, we simulate the Gaussian fitting amplitude spectrum to obtain broadband wavelet data, which we then process to obtain a higher resolution converted wave. We successfully apply the proposed inversion technology in the processing of high-resolution data from the Penglai region to obtain higher resolution convened wave data, which we then verify in a theoretical test. Improving the resolution of converted PS-wave data will provide more accurate data for subsequent velocity inversion and the extraction of reservoir reflection information.展开更多
Frequency attenuation occurs when seismic waves propagate through the porous reservoirs containing hydrocarbons. Current researches on the seismic frequency attenuation mainly focus on the post-stack domain instead of...Frequency attenuation occurs when seismic waves propagate through the porous reservoirs containing hydrocarbons. Current researches on the seismic frequency attenuation mainly focus on the post-stack domain instead of the prestack domain. Here we propose the frequency attenuation gradient vs. offset （FAGVO） based on the amplitude variation with offset and frequency attenuation integral equations. We derive the FAGVO equation that equals to zero in a full-elastic medium and is negative in a viscoelastic medium. FAGVO is affected by the viscosity of the medium, the coefficients of reflection, the frequency variation, and high-frequency attenuation. FAGVO uses the differences of partially stacked data to decrease the interference caused by subsurface strata affecting the frequency attenuation, highlights the frequency attenuation gradient anomalies in hydrocarbon-bearing reservoir pores, and finally realizes the hydrocarbon fluid identification. The method was verified using a two-dimensional wave equation forward model and was found to be cost effective. Furthermore, the method does not require well information, which can be applied in the stage of seismic exploration, especially, in the exploration of a none-well project.展开更多
3D seismic data recently acquired from the Ordos Basin shows three sets of regularly distributed fault systems,which overrides previous understanding that no faults were developed in this basin.Seismic interpretation ...3D seismic data recently acquired from the Ordos Basin shows three sets of regularly distributed fault systems,which overrides previous understanding that no faults were developed in this basin.Seismic interpretation suggests that the faults in the southwestern Ordos Basin have three basic characteristics,namely extreme micro-scale,distinct vertical stratification,and regularity of planar distribution.These NS-,NW-,and NE-trending fault systems developed in the Meso-Neoproterozoic e Lower Ordovician strata.Of these,the NS-trending fault system mainly consists of consequent and antithetic faults which show clear syndepositional deformation.The fault systems in the Carboniferous e Middle-Lower Triassic strata are not clear on seismic reflection profiles.The NW-and NE-trending fault systems are developed in the Upper Triassic e Middle Jurassic strata.Of these,the NW-trending fault system appears as a negative flower structure in sectional view and in an en echelon pattern in plan-view;they show transtensional deformation.A NE-trending fault system that developed in the Lower Cretaceous e Cenozoic strata shows a Y-shaped structural style and tension-shear properties.A comprehensive analysis of the regional stress fields at different geologic times is essential to determine the development,distribution direction,and intensity of the activity of fault systems in the Ordos Basin.Current exploration suggests three aspects in which the faults within the Ordos Basin are crucial to oil and gas accumulation.Firstly,these faults serve as vertical barriers that cause the formation of two sets of relatively independent petroleum systems in the Paleozoic and Mesozoic strata respectively;this is the basis for the‘upper oil and lower gas’distribution pattern.Secondly,the vertical communication of these faults is favorable for oil and gas migration,thus contributing to the typical characteristics of multiple oil and gas fields within the basin,i.e.oil and gas reservoirs with multiple superimposed strata.Finally,these faults and their associated fractures improve the permeability of Mesozoic tight reservoirs,providing favorable conditions for oil enrichment in areas around the fault systems.展开更多
The Longwangmiao Fm gas pool was found in 2012 in Moxi area in Leshan-Longnüsi Paleouplift,central Sichuan Basin.Previous geological studies showed that the Lower Cambrian Longwangmiao Fm reservoir is mainly pore...The Longwangmiao Fm gas pool was found in 2012 in Moxi area in Leshan-Longnüsi Paleouplift,central Sichuan Basin.Previous geological studies showed that the Lower Cambrian Longwangmiao Fm reservoir is mainly pore type reservoir of grain shoal facies,locally superimposed by the vuggy reservoir formed by Caledonian karstification,and the distribution of the reservoir has a critical control effect on the gas pool there.Therefore,the reservoir prediction in this area follows the approach of“looking for the overlap of shoal facies,karst,bright spots,and traps”.First,the favorable facies boundary and karstification effect range are defined based on sedimentary facies;on the basis of precise correlation of the top and bottom horizons of the reservoir,the favorable development area of the reservoir was then predicted according to the seismic response characteristics of“bright spots”of the Longwangmiao Fm reservoir;finally,favorable exploration areas were selected according to the effective configuration of reservoir and trap.Due to the horizontal change in reservoir and lithology,Longwangmiao Fm top has no uniform interface features on the seismic profile,making it difficult to trace.In the actual interpretation process,the underlying reference horizon is sought out first,then the bottom boundary of Longwangmiao Fm is correlated,and finally its top boundary is correlated.The prediction results indicated that the Longwangmiao reservoir distributes in band shape in a wide range around the Paleouplift;and the exploration prospects mainly involve three domains:lithologic-stratigraphic traps near the pinchout line of Longwangmiao Fm;lithologic traps formed by isolated shoal body;and structural and lithologic traps at the north flank of the Paleouplift.展开更多
The high resolution 3D nonlinear integrated inversion method is based on nonlinear theory. Under layer control, the log data from several wells （or all wells） in the study area and seismic trace data adjacent to the...The high resolution 3D nonlinear integrated inversion method is based on nonlinear theory. Under layer control, the log data from several wells （or all wells） in the study area and seismic trace data adjacent to the wells are input to a network with multiple inputs and outputs and are integratedly trained to obtain an adaptive weight function of the entire study area. Integrated nonlinear mapping relationships are built and updated by the lateral and vertical geologic variations of the reservoirs. Therefore, the inversion process and its inversion results can be constrained and controlled and a stable seismic inversion section with high resolution with velocity inversion, impedance inversion, and density inversion sections, can be gained. Good geologic effects have been obtained in model computation tests and real data processing, which verified that this method has high precision, good practicality, and can be used for quantitative reservoir analysis.展开更多
基金financially supported by the State Key Laboratory of Shale Oil and Gas Enrichment Mechanisms and Efficient Development(33550000-22-ZC0613-0006)National Natural Science Foundation of China(42202133)+2 种基金CNPC Innovation Fund(2022DQ02-0106)Strategic Cooperation Technology Projects of the CNPC and CUPB(ZLZX2020-01-05)Key Laboratory of Tectonics and Petroleum Resources(China University of Geosciences),Ministry of Education,China(TPR-2023-05)。
文摘Alkaline lacustrine shale is highly heterogeneous,and the complex relationship between the organicinorganic porosity network and hydrocarbon occurrence restricts the effectiveness of shale oil exploration and development.Herein,we investigated the Fengcheng Formation(P_(1)f)in Mahu Sag.This study integrated geochemistry,Soxhlet extraction,scanning electron microscopy,gas adsorption,and nuclear magnetic resonance T_(1)-T_(2)spectroscopy to elucidate the microscopic oil occurrence mechanisms in shales.Results indicate the presence of felsic shale,dolomitic shale,lime shale,and mixed shale within the P_(1)f.Matrix pores and microfractures associated with inorganic minerals are the predominant pore types in P_(1)f.Adsorbed oil primarily resides on the surfaces of organic matter and clay minerals,while free oil predominantly occupies inorganic pores and microfractures with larger pore sizes.Variations exist in the quantity and distribution of shale oil accumulation across different scales,where free oil and adsorbed oil are governed by dominant pores with diameters exceeding 10 nm and ineffective pores with diameters below 10 nm,respectively.Shale oil occurrence characteristics are influenced by organic matter,pore structure,and mineral composition.Felsic shale exhibits a high abundance of dominant pores,possesses the highest oil content,predominantly harbors free oil within these dominant pores,and demonstrates good mobility.Fluid occurrence in dolomitic shale and lime shale is intricate,with low oil content and a free oil to adsorbed oil ratio of 1:1.Mixed shale exhibits elevated clay mineral content and a scarcity of dominant pores.Moreover,ineffective pores contain increased bound water,resulting in medium oil content and limited mobility predominantly due to adsorption.Presently,shale oil mainly occurs in the dominant pores with a diameter larger than 10 nm in a free state.During the exploration and development of alkaline lacustrine shale oil resources,emphasis should be placed on identifying sweet spots within the felsic shale characterized by dominant pores.
基金supported by the National Natural Science Foundation of China(Research on Dynamic Location of Receiving Points and Wave Field Separation Technology Based on Deep Learning in OBN Seismic Exploration,No.42074140)the Sinopec Geophysical Corporation,Project of OBC/OBN Seismic Data Wave Field Characteristics Analysis and Ghost Wave Suppression(No.SGC-202206)。
文摘The use of blended acquisition technology in marine seismic exploration has the advantages of high acquisition efficiency and low exploration costs.However,during acquisition,the primary source may be disturbed by adjacent sources,resulting in blended noise that can adversely affect data processing and interpretation.Therefore,the de-blending method is needed to suppress blended noise and improve the quality of subsequent processing.Conventional de-blending methods,such as denoising and inversion methods,encounter challenges in parameter selection and entail high computational costs.In contrast,deep learning-based de-blending methods demonstrate reduced reliance on manual intervention and provide rapid calculation speeds post-training.In this study,we propose a Uformer network using a nonoverlapping window multihead attention mechanism designed for de-blending blended data in the common shot domain.We add the depthwise convolution to the feedforward network to improve Uformer’s ability to capture local context information.The loss function comprises SSIM and L1 loss.Our test results indicate that the Uformer outperforms convolutional neural networks and traditional denoising methods across various evaluation metrics,thus highlighting the effectiveness and advantages of Uformer in de-blending blended data.
文摘In the leakage detection of reservoir dam bodies,traditional geophysical methods can only achieve one-time detection.Meanwhile,due to the non-uniqueness of geophysicalin version,how to improve the fidelity of geophysical detection inversion profiles has become a key challenge in the industry.This study aims to construct a long-term real-time monitoring system using the time-lapse resistivity method,reveal the spatiotemporal evolution law of dam leakage,and provide technical support for accurate treatment.By integrating the Internet of Things(IoT),5G technology,and AI technology,real-time data acquisition,real-time transmission,and automatic inversionare realized.Through dynamic imaging analysis of the electrical anomaly characteristics of the leakage area and comparison between corresponding rainfall events,the leakage range and resistivity changes,the reliability and efficiency of dam leakagedete ction are significantly improved.This achieves long-termdynamic monitoring of dam leakage and provides a new perspective for the safe operation and maintenance of reservoirs.
基金supported by NSFC(Grant No.U1562109 and 41774082)the National Major Research Plan(Grant No.2016YFC0601100and 2016ZX05004)the Project of Scientific Research and Technological Development,CNPC(Grant No.2017D-5006-16)
文摘Monitoring and delineating the spatial distribution of shale fracturing is fundamentally important to shale gas production. Standard monitoring methods, such as time-lapse seismic, cross-well seismic and micro-seismic methods, are expensive, time- consuming, and do not show the changes in the formation with time. The resistivities of hydraulic fracturing fluid and reservoir rocks were measured. The results suggest that the injection fluid and consequently the injected reservoir are characterized by very low resistivity and high chargeability. This allows using of the controlled-source electromagnetic method （CSEM） to monitor shale gas hydraulic fracturing. Based on the geoelectrical model which was proposed according to the well-log and seismic data in the test area the change rule of the reacted electrical field was studied to account for the change of shale resistivity, and then the normalized residual resistivity method for time lapse processing was given. The time-domain electromagnetic method （TDEM） was used to continuously monitor the shale gas fracturing at the Fulin shale gas field in southern China. A high-power transmitter and multi-channel transient electromagnetic receiver array were adopted. 9 h time series of Ex component of 224 sites which were laid out on the surface and over three fracturing stages of a horizontal well at 2800 m depth was recorded. After data processing and calculation of the normalized resistivity residuals, the changes in the Ex signal were determined and a dynamic 3D image of the change in resistivity was constructed. This allows modeling the spatial distribution of the fracturing fluid. The model results suggest that TDEM is promising for monitoring hydraulic fracturing of shale.
基金funded by the Nonprofit Special Research Program"The formation and destruction of northeastern segment of Cathaysia-the Yangtze Plate Suture Zone and their mineralization"(No.200811015)from the Ministry of Land and Resourcethe Land Resource Survey Project of the Ministry of Land and Natural Resources,China"The convergence and breakup process of main blocks of China and their geological background for mineralization"(Nos.1212011121064,1212011121068)from the China Geological Survey
文摘Multistage deformation events have occurred in the northeastern Jiangshao Fault （Suture） Belt. The earliest two are ductile deformation events. The first is the ca. 820 Ma top-to-the-northwest ductile thrusting, which directly resulted from the collision between the Cathaysia Old Land and the Chencai Arc （？） during the Late Neoproterozoic, and the Jiangnan Orogenic Belt that formed as the ocean closed between the Yangtze Plate and the jointed Cathaysia Old Land and the Chencai Arc due to continuous compression. The second is the ductile left-lateral strike-slipping that occurred in the latest Early Paleozoic. Since the Jinning period, all deformation events represent the reactivation or inversion of intraplate structures due to the collisions between the North China and Yangtze plates during the Triassic and between the Philippine Sea and Eurasian plates during the Cenozoic. In the Triassic, brittle right-lateral strike-slipping and subsequent top-to-the south thrusting occurred along the whole northeastern Jiangshao Fault Zone because of the collision between the North China and Yangtze plates. In the Late Mesozoic, regional extension took place across southeastern China. In the Cenozoic, the collision between the Philippine Sea and Eurasian plates resulted in brittle thrusts along the whole Jiangnan Old land in the Miocene. The Jiangshao Fault Belt is a weak zone in the crust with long history, and its reactivation is one of important characteristics of the deformation in South China; however, late-stage deformation events did not occur beyond the Jiangnan Old Land and most of them are parallel to the strike of the Old Land, which is similar to the Cenozoic deformation in Central Asia. In addition, the Jiangnan old Land is not a collisional boundary between the Yangtze Plate and Cathaysia Old Land in the Triassic.
文摘Petroleum geophysicists recognize that many parameters related to oil and gas reservoirs are predicted using seismic attribute data. However, how best to optimize the seismic attributes, predict the character of thin sandstone reservoirs, and enhance the reservoir description accuracy is an important goal for geologists and geophysicists. Based on the theory of main component analysis, we present a new optimization method, called constrained main component analysis. Modeling estimates and real application in an oilfield show that it can enhance reservoir prediction accuracy and has better applicability.
文摘3D inversion of borehole-surface electrical data for complex geo-electrical models is still a challenging problem in geophysical exploration. We have developed a program for 3D inversion to borehole-surface electrical data based on the quasi-analytical approximation （QA） and re-weighted regularized conjugate gradient method （RRCG） algorithms using Visual Fortran 6.5. Application of the QA approximation to forward modeling and Frechet derivative computations speeds up the calculation dramatically. The trial calculation for synthetic data of theoretical model showed that the program is fast and highly precise.
文摘The Triassic reservoir in the Jinan area of Tarim Oilfield consists largely of interbedded sand and shale. Because of the large overlap between sandstone and shale impedance, it is difficult to distinguish sandstone from shale by acoustic impedance alone. Compared to acoustic impedance, elastic impedance contains more lithologic and physical information of the reservoir. Based on meticulous well-tie calibration, elastic impedance data volumes for 10°, 20°, and 30° emergence angles are obtained using pre-stack elastic impedance inversion. A non-linear statistical relationship between elastic impedance and shale content is set up by a PNN neural network. The non-linear mapping relationship is used to predict the reservoir shale content from elastic impedance, which will depict and predict the reservoir oil-bearing sands.
基金funded by the National Basic Research Progam of China(Nos. 2007CB411306 and 2001CB409810)China Geological Survey(No.1212010611806)the National Natural Science Foundation of China(No.40702032)
文摘Secondary/minor structures occurring along the main fault surfaces are important indicators for judging the kinematic characteristics of faults.However,many factors can lead to the formation of these structures,which results in the difficulty for rapid judgment and application in the fields.A series of secondary faults/fractures developed due to the movement of main faults are the most important and widespread phenomena in the scope of brittle deformation.The morphology of the main fault surfaces is various,and former researchers mainly discussed the structures on the main even fault surfaces. However,the fluctuation of fault surfaces is the intrinsic character of the faults,and the intersection between the main fault and secondary faults/fractures can produce a series of kinematic indicators on the main fault surfaces.Based on previous studies and our observations,i.e.the structural traces of the P,R,R',T and X shears/faults along the main faults,some indicators which are rarely reported previously,are described in the paper.Furthermore,their reliabilities are also discussed,and more practical and reliable criteria are brought forward.We suggest that the simple application of congruous and incongruous steps without knowing their exact origins should be abandoned in the fields,and several types of indicators along one fault surface should be checked with each other as much as possible.Meanwhile,the origins of some other arcuate indicators on the fault surfaces are also discussed,and new models are brought forward.
基金supported by the National Science and Technology Major project(No.2016ZX05024001003)the Innovation Consortium Project of China Petroleum,and the Southwest Petroleum University(No.2020CX010201).
文摘Oil and gas seismic exploration have to adopt irregular seismic acquisition due to the increasingly complex exploration conditions to adapt to complex geological conditions and environments.However,the irregular seismic acquisition is accompanied by the lack of acquisition data,which requires high-precision regularization.The sparse signal feature in the transform domain in compressed sensing theory is used in this paper to recover the missing signal,involving sparse transform base optimization and threshold modeling.First,this paper analyzes and compares the effects of six sparse transformation bases on the reconstruction accuracy and efficiency of irregular seismic data and establishes the quantitative relationship between sparse transformation and reconstruction accuracy and efficiency.Second,an adaptive threshold modeling method based on sparse coefficient is provided to improve the reconstruction accuracy.Test results show that the method has good adaptability to different seismic data and sparse transform bases.The f-x domain reconstruction method of effective frequency samples is studied to address the problem of low computational efficiency.The parallel computing strategy of curvelet transform combined with OpenMP is further proposed,which substantially improves the computational efficiency under the premise of ensuring the reconstruction accuracy.Finally,the actual acquisition data are used to verify the proposed method.The results indicate that the proposed method strategy can solve the regularization problem of irregular seismic data in production and improve the imaging quality of the target layer economically and efficiently.
基金the National Key R&D Program of China(Nos.2018YFC1503401,and 2016ZX05023004)National Natural Science Foundation of China(Nos.4217040570,and 41674059).
文摘The joint source scanning algorithm(SSA)scans locations and focal mechanisms of microseismic events simultaneously.Compared to the traditional source scanning algorithm,it yields much more events with extra information of focal mechanisms.The availability of more events and focal mechanisms make it possible to invert for a 2D gridded stress field.As a byproduct of hydrofracturing monitoring,the method offers a new way to extract stress field as a substitute to other more expensive technologies.This method is applied to a hydraulic fracturing dataset collected from one shale gas production field in the southeast of the Sichuan basin.A damped stress inversion is conducted to obtain a 2D stress field.five hydraulic-fracturing induced fractures can be determined from the result.The events associated with these fractures generally have relatively low stacked energy and are limited to the depth of horizontal well.One existing fault(possibly associated with the axis of the central Sichuan uplift)is also determined and the events associated with the existing fault generally have higher stacked energy and are more densely populated.The existing fault may also serve as a structural boundary where the rocks to the NW side are easier to be fractured while events on the other side are sparse with low stacked energy.The existing fault also divides the stress field into two regimes:the maximum compressional stress field to the NW and SE of the fault line are dominantly in NW-SE and N-S directions,respectively.
文摘We present a 3D approach to numerical modeling of the borehole-surface electromagnetic （BSEM） method. The 3D electromagnetic response created by a vertical line current source in a layered medium is modeled using the 3D integral equation method. The modeling results are consistent with analytical solutions. 3D Born approximation inversion of BSEM data is also conducted for reservoir delineation. The inversion method is verified by a synthetic reservoir model.
文摘Conventional land vertical seismic profiling （VSP） exploration usually uses P-wave sources and three-component geophones for receivers, emphasizing P- and converted S-waves. Previous studies show that both dynamite borehole shots and vertical vibrations from controllable seismic sources at the surface will produce relatively strong pure P-waves and weaker pure S-waves. Interfaces with a large Poisson＇s ratio difference have a positive influence on the formation of strong transmitted converted S-waves. By a comparative analysis of pure S-waves from sources and converted downgoing S-waves, we believe that the main frequency of pure S-waves is usually lower than pure P-waves while the main frequency of downgoing converted S-waves is close to that of P-waves. We have studied zero-offset and offset VSP data from land P-wave sources. Results show that pure S-waves commonly exist in these data with differences in wave intensity. S-wave velocity can be obtained from the P-wave source zero-offset VSP data. Finally, we discuss the bright future of joint application of VSP P-and S-waves and the full use of S-waves in P-wave source VSP data.
基金National 973 Key Basic Research Development Program (2005CB422104)SINOPEC's Scientific and Techno-logical Development Program(P05063)
文摘Seismic velocity is important to migration of seismic data, interpretation of lithology and lithofacies as well as prediction of reservoir. The information of shear wave velocity is required to reduce the uncertainty for discriminating lithology, identifying fluid type in porous material and calculating gas saturation in reservoir prediction. Based on Zoeppritz equations, a numeral and scanning method was proposed in this paper. Shear wave velocity can be calculated with prestack converted wave data. The effects were demonstrated by inversion of theoretical and real seismic data.
基金supported by the China National Petroleum Corporation Scientific research and technology development project(Nos.2013E-38-08)
文摘Multi-wave exploration is an effective means for improving precision in the exploration and development of complex oil and gas reservoirs that are dense and have low permeability. However, convened wave data is characterized by a low signal-to-noise ratio and low resolution, because the conventional deconvolution technology is easily affected by the frequency range limits, and there is limited scope for improving its resolution. The spectral inversion techniques is used to identify λ/8 thin layers and its breakthrough regarding band range limits has greatly improved the seismic resolution. The difficulty associated with this technology is how to use the stable inversion algorithm to obtain a high-precision reflection coefficient, and then to use this reflection coefficient to reconstruct broadband data for processing. In this paper, we focus on how to improve the vertical resolution of the converted PS-wave for multi-wave data processing. Based on previous research, we propose a least squares inversion algorithm with a total variation constraint, in which we uses the total variance as a priori information to solve under-determined problems, thereby improving the accuracy and stability of the inversion. Here, we simulate the Gaussian fitting amplitude spectrum to obtain broadband wavelet data, which we then process to obtain a higher resolution converted wave. We successfully apply the proposed inversion technology in the processing of high-resolution data from the Penglai region to obtain higher resolution convened wave data, which we then verify in a theoretical test. Improving the resolution of converted PS-wave data will provide more accurate data for subsequent velocity inversion and the extraction of reservoir reflection information.
基金supported by the Major National Science&Technology Project(No.2016ZX05003-003)the National High Technology Research and Development Program of China(863 Program)(No.2013AA064201)
文摘Frequency attenuation occurs when seismic waves propagate through the porous reservoirs containing hydrocarbons. Current researches on the seismic frequency attenuation mainly focus on the post-stack domain instead of the prestack domain. Here we propose the frequency attenuation gradient vs. offset （FAGVO） based on the amplitude variation with offset and frequency attenuation integral equations. We derive the FAGVO equation that equals to zero in a full-elastic medium and is negative in a viscoelastic medium. FAGVO is affected by the viscosity of the medium, the coefficients of reflection, the frequency variation, and high-frequency attenuation. FAGVO uses the differences of partially stacked data to decrease the interference caused by subsurface strata affecting the frequency attenuation, highlights the frequency attenuation gradient anomalies in hydrocarbon-bearing reservoir pores, and finally realizes the hydrocarbon fluid identification. The method was verified using a two-dimensional wave equation forward model and was found to be cost effective. Furthermore, the method does not require well information, which can be applied in the stage of seismic exploration, especially, in the exploration of a none-well project.
基金project entitled Seismic Identification and Accumulation Control of Strike-Slip Faults in Superimposed Basins inWest-central Part of China initiated by the Bureau of Geophysical Prospecting,CNPC(No.:03-02-2022).
文摘3D seismic data recently acquired from the Ordos Basin shows three sets of regularly distributed fault systems,which overrides previous understanding that no faults were developed in this basin.Seismic interpretation suggests that the faults in the southwestern Ordos Basin have three basic characteristics,namely extreme micro-scale,distinct vertical stratification,and regularity of planar distribution.These NS-,NW-,and NE-trending fault systems developed in the Meso-Neoproterozoic e Lower Ordovician strata.Of these,the NS-trending fault system mainly consists of consequent and antithetic faults which show clear syndepositional deformation.The fault systems in the Carboniferous e Middle-Lower Triassic strata are not clear on seismic reflection profiles.The NW-and NE-trending fault systems are developed in the Upper Triassic e Middle Jurassic strata.Of these,the NW-trending fault system appears as a negative flower structure in sectional view and in an en echelon pattern in plan-view;they show transtensional deformation.A NE-trending fault system that developed in the Lower Cretaceous e Cenozoic strata shows a Y-shaped structural style and tension-shear properties.A comprehensive analysis of the regional stress fields at different geologic times is essential to determine the development,distribution direction,and intensity of the activity of fault systems in the Ordos Basin.Current exploration suggests three aspects in which the faults within the Ordos Basin are crucial to oil and gas accumulation.Firstly,these faults serve as vertical barriers that cause the formation of two sets of relatively independent petroleum systems in the Paleozoic and Mesozoic strata respectively;this is the basis for the‘upper oil and lower gas’distribution pattern.Secondly,the vertical communication of these faults is favorable for oil and gas migration,thus contributing to the typical characteristics of multiple oil and gas fields within the basin,i.e.oil and gas reservoirs with multiple superimposed strata.Finally,these faults and their associated fractures improve the permeability of Mesozoic tight reservoirs,providing favorable conditions for oil enrichment in areas around the fault systems.
基金Major Project of National Science and Technology(No.2011ZX05004-005)Major Project of Exploration of PetroChina“Hydrocarbon Evaluation of Leshan-Longnüsi Palaeohigh in Sichuan Basin and Research on Exploration Support Technologies”(No.2012ZD01-03-01)。
文摘The Longwangmiao Fm gas pool was found in 2012 in Moxi area in Leshan-Longnüsi Paleouplift,central Sichuan Basin.Previous geological studies showed that the Lower Cambrian Longwangmiao Fm reservoir is mainly pore type reservoir of grain shoal facies,locally superimposed by the vuggy reservoir formed by Caledonian karstification,and the distribution of the reservoir has a critical control effect on the gas pool there.Therefore,the reservoir prediction in this area follows the approach of“looking for the overlap of shoal facies,karst,bright spots,and traps”.First,the favorable facies boundary and karstification effect range are defined based on sedimentary facies;on the basis of precise correlation of the top and bottom horizons of the reservoir,the favorable development area of the reservoir was then predicted according to the seismic response characteristics of“bright spots”of the Longwangmiao Fm reservoir;finally,favorable exploration areas were selected according to the effective configuration of reservoir and trap.Due to the horizontal change in reservoir and lithology,Longwangmiao Fm top has no uniform interface features on the seismic profile,making it difficult to trace.In the actual interpretation process,the underlying reference horizon is sought out first,then the bottom boundary of Longwangmiao Fm is correlated,and finally its top boundary is correlated.The prediction results indicated that the Longwangmiao reservoir distributes in band shape in a wide range around the Paleouplift;and the exploration prospects mainly involve three domains:lithologic-stratigraphic traps near the pinchout line of Longwangmiao Fm;lithologic traps formed by isolated shoal body;and structural and lithologic traps at the north flank of the Paleouplift.
基金supported by the Key Project of the National Natural Scientific Foundation(Grant No.40839909)
文摘The high resolution 3D nonlinear integrated inversion method is based on nonlinear theory. Under layer control, the log data from several wells （or all wells） in the study area and seismic trace data adjacent to the wells are input to a network with multiple inputs and outputs and are integratedly trained to obtain an adaptive weight function of the entire study area. Integrated nonlinear mapping relationships are built and updated by the lateral and vertical geologic variations of the reservoirs. Therefore, the inversion process and its inversion results can be constrained and controlled and a stable seismic inversion section with high resolution with velocity inversion, impedance inversion, and density inversion sections, can be gained. Good geologic effects have been obtained in model computation tests and real data processing, which verified that this method has high precision, good practicality, and can be used for quantitative reservoir analysis.
