Model-driven and data-driven inversions are two prominent methods for obtaining P-wave impedance,which is significant in reservoir description and identification.Based on proper initial models,most model-driven method...Model-driven and data-driven inversions are two prominent methods for obtaining P-wave impedance,which is significant in reservoir description and identification.Based on proper initial models,most model-driven methods primarily use the limited frequency bandwidth information of seismic data and can invert P-wave impedance with high accuracy,but not high resolution.Conventional data-driven methods mainly employ the information from well-log data and can provide high-accuracy and highresolution P-wave impedance owing to the superior nonlinear curve fitting capacity of neural networks.However,these methods require a significant number of training samples,which are frequently insufficient.To obtain P-wave impedance with both high accuracy and high resolution,we propose a model-data-driven inversion method using Res Nets and the normalized zero-lag cross-correlation objective function which is effective for avoiding local minima and suppressing random noise.By using initial models and training samples,the proposed model-data-driven method can invert P-wave impedance with satisfactory accuracy and resolution.Tests on synthetic and field data demonstrate the proposed method’s efficacy and practicability.展开更多
The present research work attempted to delineate and characterize the reservoir facies from the Dawson Canyon Formation in the Penobscot field,Scotian Basin.An integrated study of instantaneous frequency,P-impedance,v...The present research work attempted to delineate and characterize the reservoir facies from the Dawson Canyon Formation in the Penobscot field,Scotian Basin.An integrated study of instantaneous frequency,P-impedance,volume of clay and neutron-porosity attributes,and structural framework was done to unravel the Late Cretaceous depositional system and reservoir facies distribution patterns within the study area.Fault strikes were found in the EW and NEE-SWW directions indicating the dominant course of tectonic activities during the Late Cretaceous period in the region.P-impedance was estimated using model-based seismic inversion.Petrophysical properties such as the neutron porosity(NPHI)and volume of clay(VCL)were estimated using the multilayer perceptron neural network with high accuracy.Comparatively,a combination of low instantaneous frequency(15-30 Hz),moderate to high impedance(7000-9500 gm/cc*m/s),low neutron porosity(27%-40%)and low volume of clay(40%-60%),suggests fair-to-good sandstone development in the Dawson Canyon Formation.After calibration with the welllog data,it is found that further lowering in these attribute responses signifies the clean sandstone facies possibly containing hydrocarbons.The present study suggests that the shale lithofacies dominates the Late Cretaceous deposition(Dawson Canyon Formation)in the Penobscot field,Scotian Basin.Major faults and overlying shale facies provide structural and stratigraphic seals and act as a suitable hydrocarbon entrapment mechanism in the Dawson Canyon Formation's reservoirs.The present research advocates the integrated analysis of multi-attributes estimated using different methods to minimize the risk involved in hydrocarbon exploration.展开更多
Electrical impedance tomography (EIT) aims to reconstruct the conductivity distribution using the boundary measured voltage potential. Traditional regularization based method would suffer from error propagation due to...Electrical impedance tomography (EIT) aims to reconstruct the conductivity distribution using the boundary measured voltage potential. Traditional regularization based method would suffer from error propagation due to the iteration process. The statistical inverse problem method uses statistical inference to estimate unknown parameters. In this article, we develop a nonlinear weighted anisotropic total variation (NWATV) prior density function based on the recently proposed NWATV regularization method. We calculate the corresponding posterior density function, i.e., the solution of the EIT inverse problem in the statistical sense, via a modified Markov chain Monte Carlo (MCMC) sampling. We do numerical experiment to validate the proposed approach.展开更多
Seismic impedance inversion is an important technique for structure identification and reservoir prediction.Model-based and data-driven impedance inversion are the commonly used inversion methods.In practice,the geoph...Seismic impedance inversion is an important technique for structure identification and reservoir prediction.Model-based and data-driven impedance inversion are the commonly used inversion methods.In practice,the geophysical inversion problem is essentially an ill-posedness problem,which means that there are many solutions corresponding to the same seismic data.Therefore,regularization schemes,which can provide stable and unique inversion results to some extent,have been introduced into the objective function as constrain terms.Among them,given a low-frequency initial impedance model is the most commonly used regularization method,which can provide a smooth and stable solution.However,this model-based inversion method relies heavily on the initial model and the inversion result is band limited to the effective frequency bandwidth of seismic data,which cannot effectively improve the seismic vertical resolution and is difficult to be applied to complex structural regions.Therefore,we propose a data-driven approach for high-resolution impedance inversion based on the bidirectional long short-term memory recurrent neural network,which regards seismic data as time-series rather than image-like patches.Compared with the model-based inversion method,the data-driven approach provides higher resolution inversion results,which demonstrates the effectiveness of the data-driven method for recovering the high-frequency components.However,judging from the inversion results for characterization the spatial distribution of thin-layer sands,the accuracy of high-frequency components is difficult to guarantee.Therefore,we add the model constraint to the objective function to overcome the shortages of relying only on the data-driven schemes.First,constructing the supervisor1 based on the bidirectional long short-term memory recurrent neural network,which provides the predicted impedance with higher resolution.Then,convolution constraint as supervisor2 is introduced into the objective function to guarantee the reliability and accuracy of the inversion results,which makes the synthetic seismic data obtained from the inversion result consistent with the input data.Finally,we test the proposed scheme based on the synthetic and field seismic data.Compared to model-based and purely data-driven impedance inversion methods,the proposed approach provides more accurate and reliable inversion results while with higher vertical resolution and better spatial continuity.The inversion results accurately characterize the spatial distribution relationship of thin sands.The model tests demonstrate that the model-constrained and data-driven impedance inversion scheme can effectively improve the thin-layer structure characterization based on the seismic data.Moreover,tests on the oil field data indicate the practicality and adaptability of the proposed method.展开更多
Although the ambiguity of seismic inversion is widely recognized in both theory and practice, so far as a concrete inversion example is concerned, there is not any objective, controllable method or any standard for ho...Although the ambiguity of seismic inversion is widely recognized in both theory and practice, so far as a concrete inversion example is concerned, there is not any objective, controllable method or any standard for how to evaluate and determine its ambiguity and reliability, especially for the high frequency components beyond the effective seismic frequency band. Taking log-constrained impedance inversion as an example, a new appraisal method is proposed on the basis of analyzing a simple geological model. Firstly, the inverted impedance model is transformed to a reflection coefficient series. Secondly, the maximum effective frequency of the real seismic data is chosen as a cutoff point and the reflection coefficient series is decomposed into two components by low-pass and high-pass filters. Thirdly, the geometrical reflection characteristics of the high-frequency components and that of the real seismic data are compared and analyzed. Then, the reliability of the inverted impedance model is appraised according to the similarity of geometrical characteristics between the high-frequency components and the real seismic data. The new method avoids some subjectivity in appraising the inverted result, and helps to enhance the reliability of reservoir prediction by impedance inversion technology.展开更多
Considering the uncertainty of the electrical axis for two-dimensional audo-magnetotelluric(AMT) data processing, an AMT inversion method with the Central impedance tensor was presented. First, we present a calculatio...Considering the uncertainty of the electrical axis for two-dimensional audo-magnetotelluric(AMT) data processing, an AMT inversion method with the Central impedance tensor was presented. First, we present a calculation expression of the Central impedance tensor in AMT, which can be considered as the arithmetic mean of TE-polarization mode and TM-polarization mode in the twodimensional geo-electrical model. Second, a least-squares iterative inversion algorithm is established, based on a smoothnessconstrained model, and an improved L-curve method is adopted to determine the best regularization parameters. We then test the above inversion method with synthetic data and field data. The test results show that this two-dimensional AMT inversion scheme for the responses of Central impedance is effective and can reconstruct reasonable two-dimensional subsurface resistivity structures. We conclude that the Central impedance tensor is a useful tool for two-dimensional inversion of AMT data.展开更多
Deep learning is widely used for seismic impedance inversion,but few work provides in-depth research and analysis on designing the architectures of deep neural networks and choosing the network hyperparameters.This pa...Deep learning is widely used for seismic impedance inversion,but few work provides in-depth research and analysis on designing the architectures of deep neural networks and choosing the network hyperparameters.This paper is dedicated to comprehensively studying on the significant aspects of deep neural networks that affect the inversion results.We experimentally reveal how network hyperparameters and architectures affect the inversion performance,and develop a series of methods which are proven to be effective in reconstructing high-frequency information in the estimated impedance model.Experiments demonstrate that the proposed multi-scale architecture is helpful to reconstruct more high-frequency details than a conventional network.Besides,the reconstruction of high-frequency information can be further promoted by introducing a perceptual loss and a generative adversarial network from the computer vision perspective.More importantly,the experimental results provide valuable references for designing proper network architectures in the seismic inversion problem.展开更多
Connolly (1999 ) 有弹性的阻抗(EI ) 方程是 P 波浪速度, S 波浪速度,密度,和发生角度的功能。常规倒置方法能仅仅直接基于这个方程提取 P 速度, S 速度,和密度数据,在不同发生角度的有弹性的阻抗不在一样的规模,它使比较困难...Connolly (1999 ) 有弹性的阻抗(EI ) 方程是 P 波浪速度, S 波浪速度,密度,和发生角度的功能。常规倒置方法能仅仅直接基于这个方程提取 P 速度, S 速度,和密度数据,在不同发生角度的有弹性的阻抗不在一样的规模,它使比较困难。我们基于灰色的等建议一个新有弹性的阻抗方程。(1999 ) 地址的 Zoeppritz 近似使用 Lam é参数常规倒置方法的缺乏。这个方程被使正常化在不同角度统一有弹性的阻抗尺寸并且为倒置使用了。Lam é参数能从用在 Lam é参数和有弹性的阻抗之间的线性关系从倒置获得的有弹性的阻抗数据直接被提取。有弹性的参数用这个新方法提取了的申请例子表演更稳定、正确并且能恢复水库信息很好。新方法是常规方法上的改进基于 Connolly 的方程。展开更多
Elastic impedance inversion with high efficiency and high stability has become one of the main directions of seismic pre-stack inversion. The nonlinear elastic impedance inversion method based on a fast Markov chain M...Elastic impedance inversion with high efficiency and high stability has become one of the main directions of seismic pre-stack inversion. The nonlinear elastic impedance inversion method based on a fast Markov chain Monte Carlo (MCMC) method is proposed in this paper, combining conventional MCMC method based on global optimization with a preconditioned conjugate gradient (PCG) algorithm based on local optimization, so this method does not depend strongly on the initial model. It converges to the global optimum quickly and efficiently on the condition that effi- ciency and stability of inversion are both taken into consid- eration at the same time. The test data verify the feasibility and robustness of the method, and based on this method, we extract the effective pore-fluid bulk modulus, which is applied to reservoir fluid identification and detection, and consequently, a better result has been achieved.展开更多
The classical elastic impedance(EI) inversion method,however,is based on the L2-norm misfit function and considerably sensitive to outliers,assuming the noise of the seismic data to be the Guassian-distribution.So we ...The classical elastic impedance(EI) inversion method,however,is based on the L2-norm misfit function and considerably sensitive to outliers,assuming the noise of the seismic data to be the Guassian-distribution.So we have developed a more robust elastic impedance inversion based on the L1-norm misfit function,and the noise is assumed to be non-Gaussian.Meanwhile,some regularization methods including the sparse constraint regularization and elastic impedance point constraint regularization are incorporated to improve the ill-posed characteristics of the seismic inversion problem.Firstly,we create the Ll-norm misfit objective function of pre-stack inversion problem based on the Bayesian scheme within the sparse constraint regularization and elastic impedance point constraint regularization.And then,we obtain more robust elastic impedances of different angles which are less sensitive to outliers in seismic data by using the IRLS strategy.Finally,we extract the P-wave and S-wave velocity and density by using the more stable parameter extraction method.Tests on synthetic data show that the P-wave and S-wave velocity and density parameters are still estimated reasonable with moderate noise.A test on the real data set shows that compared to the results of the classical elastic impedance inversion method,the estimated results using the proposed method can get better lateral continuity and more distinct show of the gas,verifying the feasibility and stability of the method.展开更多
This paper presents an electrical impedance tomography(EIT)method using a partial-differential-equationconstrained optimization approach.The forward problem in the inversion framework is described by a complete electr...This paper presents an electrical impedance tomography(EIT)method using a partial-differential-equationconstrained optimization approach.The forward problem in the inversion framework is described by a complete electrodemodel(CEM),which seeks the electric potential within the domain and at surface electrodes considering the contact impedance between them.The finite element solution of the electric potential has been validated using a commercial code.The inverse medium problem for reconstructing the unknown electrical conductivity profile is formulated as an optimization problem constrained by the CEM.The method seeks the optimal solution of the domain’s electrical conductivity to minimize a Lagrangian functional consisting of a least-squares objective functional and a regularization term.Enforcing the stationarity of the Lagrangian leads to state,adjoint,and control problems,which constitute the Karush-Kuhn-Tucker(KKT)first-order optimality conditions.Subsequently,the electrical conductivity profile of the domain is iteratively updated by solving the KKT conditions in the reduced space of the control variable.Numerical results show that the relative error of the measured and calculated electric potentials after the inversion is less than 1%,demonstrating the successful reconstruction of heterogeneous electrical conductivity profiles using the proposed EIT method.This method thus represents an application framework for nondestructive evaluation of structures and geotechnical site characterization.展开更多
Simultaneous waveform inversion was used to predict lithofacies and fluid type across the field. Very often, characterizing reservoirs in terms of lithology and fluid type using conventional methods is replete with un...Simultaneous waveform inversion was used to predict lithofacies and fluid type across the field. Very often, characterizing reservoirs in terms of lithology and fluid type using conventional methods is replete with uncertainties, especially in marginal fields. An approach is employed in this study that integrated rock physics and waveform inverse modelling for lithology and fluid-type characterization to appropriately identify potential hydrocarbon saturated zones and their corresponding lithology. Seismic and well-log data were analyzed using Hampson Russel software. The method adopted includes lithofacies and fluid content analysis using rock physics parameters and seismic simultaneous inverse modelling. Rock physics analysis identified 2 broad reservoirs namely: HDZ1 and HDZ2 reservoirs. Results from the inverse modelling showed that low values of acoustic impedance from 19,743 to 20,487 (ft/s)(g/cc) reflect hydrocarbon-bearing reservoirs while medium to high values shows brine and shale respectively, with brine zone ranging from 20,487 to 22,531 (ft/s)(g/cc) and shale above 22,531 (ft/s)(g/cc). Two lithofacies were identified from inversion analysis of Vp/Vs and Mu-Rho, namely: sand and shale with VpVs 1.95 values respectively. Mu-Rho > 12.29 (GPa)(g/cc) and <12.29 (GPa) (g/cc) represent sand and shale respectively. From 3D volume, it was observed that a high accumulation of hydrocarbon was observed to be saturated at the north to the eastern part of the field forming a meandering channel. Sands were mainly distributed around the northeastern to the southwestern part of the field, that tends to be away from Well 029. This was also validated by the volume of rigidity modulus (Mu-Rho) showing high values indicating sands fall within the northeastern part of the field.展开更多
For a typical marine shale reservoir in the Jiaoshiba area, Sichuan Basin of China, P-impedance is sensitive for identifying lithology but not suitable for indicating good shale reservoirs. In comparison, density is a...For a typical marine shale reservoir in the Jiaoshiba area, Sichuan Basin of China, P-impedance is sensitive for identifying lithology but not suitable for indicating good shale reservoirs. In comparison, density is an important quantity, which is sensitive for identifying the organic-rich mud shale from non-organic-rich mud shale. Due to the poor data quality and incidence angle range, density cannot be easily inverted by directly solving the ill-posed pre-stack seismic inversion in this area. Meanwhile, the traditional density regularizations implemented by directly using the more robust P-impedance inversion tend to be inaccurate for recovering density for this shale reservoir. In this paper, we combine the P-impedance and the minus uranium to construct the pseudo-P-impedance(PIp) at well locations. The PIp is observed to be sensitive for identifying organic-rich mud shale and has a good correlation with density in this area. We employ the PIp–density relation into the pre-stack inversion framework to estimate density. Three types of regularization are tested on both numerical and field data: These are no regularization, traditional regularization and the proposed approach. It is observed that the proposed method is better for recovering the density of organic-rich mud shale in the Jiaoshiba area.展开更多
Comprehensive inversion of logging and seismic data presented in this paper is a method to improve seismic data resolution. It involves using ample high-frequency information and complete low-frequency information of ...Comprehensive inversion of logging and seismic data presented in this paper is a method to improve seismic data resolution. It involves using ample high-frequency information and complete low-frequency information of known logging to make up for the lack of limited bandwidth of practical seismic recording, obtaining an approximate reflection coefficient sequence (or wave impedance) of high resolution by iterative inversion and providing more reliable seismic evidence for further lithologic interpretation and lateral tracking, correlation and prediction of thin reservoir. The comprehensive inversion can be realized in the following steps: (1) to establish an initial model of higher resolution; (2) to obtain wavelets, and (3) to constrain iterative inversion. The key to this inversion lies in building an initial model. It is assumed from our experience that when the initial model is properly given, iterative inversion can be quickly converged to the ideal result.展开更多
Accurate estimation of fracture density and orientation is of great significance for seismic characterization of fractured reservoirs.Here,we propose a novel methodology to estimate fracture density and orientation fr...Accurate estimation of fracture density and orientation is of great significance for seismic characterization of fractured reservoirs.Here,we propose a novel methodology to estimate fracture density and orientation from azimuthal elastic impedance(AEI)difference using singular value decomposition(SVD).Based on Hudson's model,we first derive the AEI equation containing fracture density in HTI media,and then obtain basis functions and singular values from the normalized AEI difference utilizing SVD.Analysis shows that the basis function changing with azimuth is related to fracture orientation,fracture density is the linearly weighted sum of singular values,and the first singular value contributes the most to fracture density.Thus,we develop an SVD-based fracture density and orientation inversion approach constrained by smooth prior elastic parameters.Synthetic example shows that fracture density and orientation can be stably estimated,and the correlation coefficient between the true value and the estimated fracture density is above 0.85 even when an S/N ratio of 2.Field data example shows that the estimated fracture orientation is consistent with the interpretation of image log data,and the estimated fracture density reliably indicates fractured gas-bearing reservoir,which could help to guide the exploration and development of fractured reservoirs.展开更多
Fluid and effective fracture identification in reservoirs is a crucial part of reservoir prediction.The frequency-dependent AVO inversion algorithms have proven to be effective for identifying fluid through its disper...Fluid and effective fracture identification in reservoirs is a crucial part of reservoir prediction.The frequency-dependent AVO inversion algorithms have proven to be effective for identifying fluid through its dispersion property.However,the conventional frequency-dependent AVO inversion algorithms based on Smith&Gidlow and Aki&Richards approximations do not consider the acquisition azimuth of seismic data and neglect the effect of seismic anisotropic dispersion in the actual medium.The aligned fractures in the subsurface medium induce anisotropy.The seismic anisotropy should be considered while accounting for the seismic dispersion properties through fluid-saturated fractured reservoirs.Anisotropy in such reservoirs is frequency-related due to wave-induced fluid-flow(WIFF)between interconnected fractures and pores.It can be used to identify fluid and effective fractures(fluid-saturated)by using azimuthal seismic data via anisotropic dispersion properties.In this paper,based on Rüger’s equation,we derived an analytical expression in the frequency domain for the frequencydependent AVOAz inversion in terms of fracture orientation,dispersion gradient of isotropic background rock,anisotropic dispersion gradient,and the dispersion at a normal incident angle.The frequency-dependent AVOAz equation utilizes azimuthal seismic data and considers the effect of both isotropic and anisotropic dispersion.Reassigned Gabor Transform(RGT)is used to achieve highresolution frequency division data.We then propose the frequency-dependent AVOAz inversion method to identify fluid and characterize effective fractures in fractured porous reservoirs.Through application to high-qualified seismic data of dolomite and carbonate reservoirs,the results show that the method is useful for identifying fluid and effective fractures in fluid-saturated fractured rocks.展开更多
Using a bottom simulating reflector(BSR)on a seismic profile to identify marine gas hydrate is a traditional seismic exploration method.However,owing to the abundance differences between the gas hydrate and free gas i...Using a bottom simulating reflector(BSR)on a seismic profile to identify marine gas hydrate is a traditional seismic exploration method.However,owing to the abundance differences between the gas hydrate and free gas in different regions,the BSR may be unremarkable on the seismic profile and invisible in certain cases.With the improvement of exploration precision,difficulty arises in meeting the requirements of distinguishing the abundance differences in the gas hydrate based on BSR.Hence,we studied other sensitive attributes to ascertain the existence of gas hydrate and its abundance variations,eventually improving the success rate of drilling and productivity.In this paper,we analyzed the contradiction between the seismic profile data and drilling sampling data from the Blake Ridge.We extracted different attributes and performed multi-parameter constraint analysis based on the prestack elastic wave impedance inversion.Then,we compared the analysis results with the drilling sampling data.Eventually,we determined five sensitive attributes that can better indicate the existence of gas hydrate and its abundance variations.This method overcomes the limitations of recognizing the gas hydrate methods based on BSR or single inversion attribute.Moreover,the conclusions can notably improve the identification accuracy of marine gas hydrate and provide excellent reference significance for the recognition of marine gas hydrate.Notably,the different geological features of reservoirs feature different sensitivities to the prestacking attributes when using the prestack elastic inversion in different areas.展开更多
Seismic inversion and basic theory are briefly presented and the main idea of this method is introduced. Both non-linear wave equation inversion technique and Complete Utilization of Samples Information (CUSI) neural ...Seismic inversion and basic theory are briefly presented and the main idea of this method is introduced. Both non-linear wave equation inversion technique and Complete Utilization of Samples Information (CUSI) neural network analysis are used in lithological interpretation in Jibei coal field. The prediction results indicate that this method can provide reliable data for thin coal exploitation and promising area evaluation.展开更多
基金financially supported by the Important National Science&Technology Specific Project of China(Grant No.2017ZX05018-005)
文摘Model-driven and data-driven inversions are two prominent methods for obtaining P-wave impedance,which is significant in reservoir description and identification.Based on proper initial models,most model-driven methods primarily use the limited frequency bandwidth information of seismic data and can invert P-wave impedance with high accuracy,but not high resolution.Conventional data-driven methods mainly employ the information from well-log data and can provide high-accuracy and highresolution P-wave impedance owing to the superior nonlinear curve fitting capacity of neural networks.However,these methods require a significant number of training samples,which are frequently insufficient.To obtain P-wave impedance with both high accuracy and high resolution,we propose a model-data-driven inversion method using Res Nets and the normalized zero-lag cross-correlation objective function which is effective for avoiding local minima and suppressing random noise.By using initial models and training samples,the proposed model-data-driven method can invert P-wave impedance with satisfactory accuracy and resolution.Tests on synthetic and field data demonstrate the proposed method’s efficacy and practicability.
文摘The present research work attempted to delineate and characterize the reservoir facies from the Dawson Canyon Formation in the Penobscot field,Scotian Basin.An integrated study of instantaneous frequency,P-impedance,volume of clay and neutron-porosity attributes,and structural framework was done to unravel the Late Cretaceous depositional system and reservoir facies distribution patterns within the study area.Fault strikes were found in the EW and NEE-SWW directions indicating the dominant course of tectonic activities during the Late Cretaceous period in the region.P-impedance was estimated using model-based seismic inversion.Petrophysical properties such as the neutron porosity(NPHI)and volume of clay(VCL)were estimated using the multilayer perceptron neural network with high accuracy.Comparatively,a combination of low instantaneous frequency(15-30 Hz),moderate to high impedance(7000-9500 gm/cc*m/s),low neutron porosity(27%-40%)and low volume of clay(40%-60%),suggests fair-to-good sandstone development in the Dawson Canyon Formation.After calibration with the welllog data,it is found that further lowering in these attribute responses signifies the clean sandstone facies possibly containing hydrocarbons.The present study suggests that the shale lithofacies dominates the Late Cretaceous deposition(Dawson Canyon Formation)in the Penobscot field,Scotian Basin.Major faults and overlying shale facies provide structural and stratigraphic seals and act as a suitable hydrocarbon entrapment mechanism in the Dawson Canyon Formation's reservoirs.The present research advocates the integrated analysis of multi-attributes estimated using different methods to minimize the risk involved in hydrocarbon exploration.
文摘Electrical impedance tomography (EIT) aims to reconstruct the conductivity distribution using the boundary measured voltage potential. Traditional regularization based method would suffer from error propagation due to the iteration process. The statistical inverse problem method uses statistical inference to estimate unknown parameters. In this article, we develop a nonlinear weighted anisotropic total variation (NWATV) prior density function based on the recently proposed NWATV regularization method. We calculate the corresponding posterior density function, i.e., the solution of the EIT inverse problem in the statistical sense, via a modified Markov chain Monte Carlo (MCMC) sampling. We do numerical experiment to validate the proposed approach.
基金funded by R&D Department of China National Petroleum Corporation(2022DQ0604-04)the Strategic Cooperation Technology Projects of CNPC and CUPB(ZLZX2020-03)the Science Research and Technology Development of PetroChina(2021DJ1206).
文摘Seismic impedance inversion is an important technique for structure identification and reservoir prediction.Model-based and data-driven impedance inversion are the commonly used inversion methods.In practice,the geophysical inversion problem is essentially an ill-posedness problem,which means that there are many solutions corresponding to the same seismic data.Therefore,regularization schemes,which can provide stable and unique inversion results to some extent,have been introduced into the objective function as constrain terms.Among them,given a low-frequency initial impedance model is the most commonly used regularization method,which can provide a smooth and stable solution.However,this model-based inversion method relies heavily on the initial model and the inversion result is band limited to the effective frequency bandwidth of seismic data,which cannot effectively improve the seismic vertical resolution and is difficult to be applied to complex structural regions.Therefore,we propose a data-driven approach for high-resolution impedance inversion based on the bidirectional long short-term memory recurrent neural network,which regards seismic data as time-series rather than image-like patches.Compared with the model-based inversion method,the data-driven approach provides higher resolution inversion results,which demonstrates the effectiveness of the data-driven method for recovering the high-frequency components.However,judging from the inversion results for characterization the spatial distribution of thin-layer sands,the accuracy of high-frequency components is difficult to guarantee.Therefore,we add the model constraint to the objective function to overcome the shortages of relying only on the data-driven schemes.First,constructing the supervisor1 based on the bidirectional long short-term memory recurrent neural network,which provides the predicted impedance with higher resolution.Then,convolution constraint as supervisor2 is introduced into the objective function to guarantee the reliability and accuracy of the inversion results,which makes the synthetic seismic data obtained from the inversion result consistent with the input data.Finally,we test the proposed scheme based on the synthetic and field seismic data.Compared to model-based and purely data-driven impedance inversion methods,the proposed approach provides more accurate and reliable inversion results while with higher vertical resolution and better spatial continuity.The inversion results accurately characterize the spatial distribution relationship of thin sands.The model tests demonstrate that the model-constrained and data-driven impedance inversion scheme can effectively improve the thin-layer structure characterization based on the seismic data.Moreover,tests on the oil field data indicate the practicality and adaptability of the proposed method.
基金supported by the Major Basic Research Development Program of China’s 973 Project(grant No.2007CB209608)the Science and Technology Innovation Foundation of CNPC(grant No.2010D-5006-0301)
文摘Although the ambiguity of seismic inversion is widely recognized in both theory and practice, so far as a concrete inversion example is concerned, there is not any objective, controllable method or any standard for how to evaluate and determine its ambiguity and reliability, especially for the high frequency components beyond the effective seismic frequency band. Taking log-constrained impedance inversion as an example, a new appraisal method is proposed on the basis of analyzing a simple geological model. Firstly, the inverted impedance model is transformed to a reflection coefficient series. Secondly, the maximum effective frequency of the real seismic data is chosen as a cutoff point and the reflection coefficient series is decomposed into two components by low-pass and high-pass filters. Thirdly, the geometrical reflection characteristics of the high-frequency components and that of the real seismic data are compared and analyzed. Then, the reliability of the inverted impedance model is appraised according to the similarity of geometrical characteristics between the high-frequency components and the real seismic data. The new method avoids some subjectivity in appraising the inverted result, and helps to enhance the reliability of reservoir prediction by impedance inversion technology.
基金supported by National Natural Science Foundation of China (grant 41674080)Higher School Doctor Subject Special Scientific Research Foundation (grant 20110162120064)
文摘Considering the uncertainty of the electrical axis for two-dimensional audo-magnetotelluric(AMT) data processing, an AMT inversion method with the Central impedance tensor was presented. First, we present a calculation expression of the Central impedance tensor in AMT, which can be considered as the arithmetic mean of TE-polarization mode and TM-polarization mode in the twodimensional geo-electrical model. Second, a least-squares iterative inversion algorithm is established, based on a smoothnessconstrained model, and an improved L-curve method is adopted to determine the best regularization parameters. We then test the above inversion method with synthetic data and field data. The test results show that this two-dimensional AMT inversion scheme for the responses of Central impedance is effective and can reconstruct reasonable two-dimensional subsurface resistivity structures. We conclude that the Central impedance tensor is a useful tool for two-dimensional inversion of AMT data.
基金supported by the National Natural Science Foundation of China under Grant No.42050104
文摘Deep learning is widely used for seismic impedance inversion,but few work provides in-depth research and analysis on designing the architectures of deep neural networks and choosing the network hyperparameters.This paper is dedicated to comprehensively studying on the significant aspects of deep neural networks that affect the inversion results.We experimentally reveal how network hyperparameters and architectures affect the inversion performance,and develop a series of methods which are proven to be effective in reconstructing high-frequency information in the estimated impedance model.Experiments demonstrate that the proposed multi-scale architecture is helpful to reconstruct more high-frequency details than a conventional network.Besides,the reconstruction of high-frequency information can be further promoted by introducing a perceptual loss and a generative adversarial network from the computer vision perspective.More importantly,the experimental results provide valuable references for designing proper network architectures in the seismic inversion problem.
基金the sponsorship of the National Basic Research Program of China (973 Program,2013CB228604,2014CB239201)the National Oil and Gas Major Projects of China (2011ZX05014-001-010HZ,2011ZX05014-001-006-XY570) for their funding of this research
文摘Elastic impedance inversion with high efficiency and high stability has become one of the main directions of seismic pre-stack inversion. The nonlinear elastic impedance inversion method based on a fast Markov chain Monte Carlo (MCMC) method is proposed in this paper, combining conventional MCMC method based on global optimization with a preconditioned conjugate gradient (PCG) algorithm based on local optimization, so this method does not depend strongly on the initial model. It converges to the global optimum quickly and efficiently on the condition that effi- ciency and stability of inversion are both taken into consid- eration at the same time. The test data verify the feasibility and robustness of the method, and based on this method, we extract the effective pore-fluid bulk modulus, which is applied to reservoir fluid identification and detection, and consequently, a better result has been achieved.
基金Projects(U1562215,41674130,41404088)supported by the National Natural Science Foundation of ChinaProjects(2013CB228604,2014CB239201)supported by the National Basic Research Program of China+1 种基金Projects(2016ZX05027004-001,2016ZX05002006-009)supported by the National Oil and Gas Major Projects of ChinaProject(15CX08002A)supported by the Fundamental Research Funds for the Central Universities,China
文摘The classical elastic impedance(EI) inversion method,however,is based on the L2-norm misfit function and considerably sensitive to outliers,assuming the noise of the seismic data to be the Guassian-distribution.So we have developed a more robust elastic impedance inversion based on the L1-norm misfit function,and the noise is assumed to be non-Gaussian.Meanwhile,some regularization methods including the sparse constraint regularization and elastic impedance point constraint regularization are incorporated to improve the ill-posed characteristics of the seismic inversion problem.Firstly,we create the Ll-norm misfit objective function of pre-stack inversion problem based on the Bayesian scheme within the sparse constraint regularization and elastic impedance point constraint regularization.And then,we obtain more robust elastic impedances of different angles which are less sensitive to outliers in seismic data by using the IRLS strategy.Finally,we extract the P-wave and S-wave velocity and density by using the more stable parameter extraction method.Tests on synthetic data show that the P-wave and S-wave velocity and density parameters are still estimated reasonable with moderate noise.A test on the real data set shows that compared to the results of the classical elastic impedance inversion method,the estimated results using the proposed method can get better lateral continuity and more distinct show of the gas,verifying the feasibility and stability of the method.
基金funded by the National Research Foundation of Korea,the Grant from a Basic Science and Engineering Research Project(NRF-2017R1C1B200497515)and the Grant from Basic Laboratory Support Project(NRF-2020R1A4A101882611).
文摘This paper presents an electrical impedance tomography(EIT)method using a partial-differential-equationconstrained optimization approach.The forward problem in the inversion framework is described by a complete electrodemodel(CEM),which seeks the electric potential within the domain and at surface electrodes considering the contact impedance between them.The finite element solution of the electric potential has been validated using a commercial code.The inverse medium problem for reconstructing the unknown electrical conductivity profile is formulated as an optimization problem constrained by the CEM.The method seeks the optimal solution of the domain’s electrical conductivity to minimize a Lagrangian functional consisting of a least-squares objective functional and a regularization term.Enforcing the stationarity of the Lagrangian leads to state,adjoint,and control problems,which constitute the Karush-Kuhn-Tucker(KKT)first-order optimality conditions.Subsequently,the electrical conductivity profile of the domain is iteratively updated by solving the KKT conditions in the reduced space of the control variable.Numerical results show that the relative error of the measured and calculated electric potentials after the inversion is less than 1%,demonstrating the successful reconstruction of heterogeneous electrical conductivity profiles using the proposed EIT method.This method thus represents an application framework for nondestructive evaluation of structures and geotechnical site characterization.
文摘Simultaneous waveform inversion was used to predict lithofacies and fluid type across the field. Very often, characterizing reservoirs in terms of lithology and fluid type using conventional methods is replete with uncertainties, especially in marginal fields. An approach is employed in this study that integrated rock physics and waveform inverse modelling for lithology and fluid-type characterization to appropriately identify potential hydrocarbon saturated zones and their corresponding lithology. Seismic and well-log data were analyzed using Hampson Russel software. The method adopted includes lithofacies and fluid content analysis using rock physics parameters and seismic simultaneous inverse modelling. Rock physics analysis identified 2 broad reservoirs namely: HDZ1 and HDZ2 reservoirs. Results from the inverse modelling showed that low values of acoustic impedance from 19,743 to 20,487 (ft/s)(g/cc) reflect hydrocarbon-bearing reservoirs while medium to high values shows brine and shale respectively, with brine zone ranging from 20,487 to 22,531 (ft/s)(g/cc) and shale above 22,531 (ft/s)(g/cc). Two lithofacies were identified from inversion analysis of Vp/Vs and Mu-Rho, namely: sand and shale with VpVs 1.95 values respectively. Mu-Rho > 12.29 (GPa)(g/cc) and <12.29 (GPa) (g/cc) represent sand and shale respectively. From 3D volume, it was observed that a high accumulation of hydrocarbon was observed to be saturated at the north to the eastern part of the field forming a meandering channel. Sands were mainly distributed around the northeastern to the southwestern part of the field, that tends to be away from Well 029. This was also validated by the volume of rigidity modulus (Mu-Rho) showing high values indicating sands fall within the northeastern part of the field.
基金NSFC and Sinopec Joint Key Project (U1663207)the China Geology Survey Project (DD20160195)+2 种基金973 Program (2014CB239104)National Key S&T Projects (2017ZX05049002)China Postdoctoral Science Foundation for the financial support
文摘For a typical marine shale reservoir in the Jiaoshiba area, Sichuan Basin of China, P-impedance is sensitive for identifying lithology but not suitable for indicating good shale reservoirs. In comparison, density is an important quantity, which is sensitive for identifying the organic-rich mud shale from non-organic-rich mud shale. Due to the poor data quality and incidence angle range, density cannot be easily inverted by directly solving the ill-posed pre-stack seismic inversion in this area. Meanwhile, the traditional density regularizations implemented by directly using the more robust P-impedance inversion tend to be inaccurate for recovering density for this shale reservoir. In this paper, we combine the P-impedance and the minus uranium to construct the pseudo-P-impedance(PIp) at well locations. The PIp is observed to be sensitive for identifying organic-rich mud shale and has a good correlation with density in this area. We employ the PIp–density relation into the pre-stack inversion framework to estimate density. Three types of regularization are tested on both numerical and field data: These are no regularization, traditional regularization and the proposed approach. It is observed that the proposed method is better for recovering the density of organic-rich mud shale in the Jiaoshiba area.
文摘Comprehensive inversion of logging and seismic data presented in this paper is a method to improve seismic data resolution. It involves using ample high-frequency information and complete low-frequency information of known logging to make up for the lack of limited bandwidth of practical seismic recording, obtaining an approximate reflection coefficient sequence (or wave impedance) of high resolution by iterative inversion and providing more reliable seismic evidence for further lithologic interpretation and lateral tracking, correlation and prediction of thin reservoir. The comprehensive inversion can be realized in the following steps: (1) to establish an initial model of higher resolution; (2) to obtain wavelets, and (3) to constrain iterative inversion. The key to this inversion lies in building an initial model. It is assumed from our experience that when the initial model is properly given, iterative inversion can be quickly converged to the ideal result.
基金sponsorship of the National Natural Science Foundation of China(41674130,U19B2008)the Postgraduate Innovation Project in China University of Petroleum(East China)(YCX2021016)for their funding this research。
文摘Accurate estimation of fracture density and orientation is of great significance for seismic characterization of fractured reservoirs.Here,we propose a novel methodology to estimate fracture density and orientation from azimuthal elastic impedance(AEI)difference using singular value decomposition(SVD).Based on Hudson's model,we first derive the AEI equation containing fracture density in HTI media,and then obtain basis functions and singular values from the normalized AEI difference utilizing SVD.Analysis shows that the basis function changing with azimuth is related to fracture orientation,fracture density is the linearly weighted sum of singular values,and the first singular value contributes the most to fracture density.Thus,we develop an SVD-based fracture density and orientation inversion approach constrained by smooth prior elastic parameters.Synthetic example shows that fracture density and orientation can be stably estimated,and the correlation coefficient between the true value and the estimated fracture density is above 0.85 even when an S/N ratio of 2.Field data example shows that the estimated fracture orientation is consistent with the interpretation of image log data,and the estimated fracture density reliably indicates fractured gas-bearing reservoir,which could help to guide the exploration and development of fractured reservoirs.
基金supported by the National Major Science and Technology Project of China(2016ZX05004003)the National Natural Science Foundation of China(41574103,41974120,U20B2015)Open Fund of State Key Laboratory of Coal Resources and Safe Mining(Grant No.SKLCRSM19KFA08)。
文摘Fluid and effective fracture identification in reservoirs is a crucial part of reservoir prediction.The frequency-dependent AVO inversion algorithms have proven to be effective for identifying fluid through its dispersion property.However,the conventional frequency-dependent AVO inversion algorithms based on Smith&Gidlow and Aki&Richards approximations do not consider the acquisition azimuth of seismic data and neglect the effect of seismic anisotropic dispersion in the actual medium.The aligned fractures in the subsurface medium induce anisotropy.The seismic anisotropy should be considered while accounting for the seismic dispersion properties through fluid-saturated fractured reservoirs.Anisotropy in such reservoirs is frequency-related due to wave-induced fluid-flow(WIFF)between interconnected fractures and pores.It can be used to identify fluid and effective fractures(fluid-saturated)by using azimuthal seismic data via anisotropic dispersion properties.In this paper,based on Rüger’s equation,we derived an analytical expression in the frequency domain for the frequencydependent AVOAz inversion in terms of fracture orientation,dispersion gradient of isotropic background rock,anisotropic dispersion gradient,and the dispersion at a normal incident angle.The frequency-dependent AVOAz equation utilizes azimuthal seismic data and considers the effect of both isotropic and anisotropic dispersion.Reassigned Gabor Transform(RGT)is used to achieve highresolution frequency division data.We then propose the frequency-dependent AVOAz inversion method to identify fluid and characterize effective fractures in fractured porous reservoirs.Through application to high-qualified seismic data of dolomite and carbonate reservoirs,the results show that the method is useful for identifying fluid and effective fractures in fluid-saturated fractured rocks.
基金supported by the National Natural Science Foundation of China (No. 41230318)
文摘Using a bottom simulating reflector(BSR)on a seismic profile to identify marine gas hydrate is a traditional seismic exploration method.However,owing to the abundance differences between the gas hydrate and free gas in different regions,the BSR may be unremarkable on the seismic profile and invisible in certain cases.With the improvement of exploration precision,difficulty arises in meeting the requirements of distinguishing the abundance differences in the gas hydrate based on BSR.Hence,we studied other sensitive attributes to ascertain the existence of gas hydrate and its abundance variations,eventually improving the success rate of drilling and productivity.In this paper,we analyzed the contradiction between the seismic profile data and drilling sampling data from the Blake Ridge.We extracted different attributes and performed multi-parameter constraint analysis based on the prestack elastic wave impedance inversion.Then,we compared the analysis results with the drilling sampling data.Eventually,we determined five sensitive attributes that can better indicate the existence of gas hydrate and its abundance variations.This method overcomes the limitations of recognizing the gas hydrate methods based on BSR or single inversion attribute.Moreover,the conclusions can notably improve the identification accuracy of marine gas hydrate and provide excellent reference significance for the recognition of marine gas hydrate.Notably,the different geological features of reservoirs feature different sensitivities to the prestacking attributes when using the prestack elastic inversion in different areas.
文摘Seismic inversion and basic theory are briefly presented and the main idea of this method is introduced. Both non-linear wave equation inversion technique and Complete Utilization of Samples Information (CUSI) neural network analysis are used in lithological interpretation in Jibei coal field. The prediction results indicate that this method can provide reliable data for thin coal exploitation and promising area evaluation.