We apply stochastic seismic inversion and Bayesian facies classification for porosity modeling and igneous rock identification in the presalt interval of the Santos Basin. This integration of seismic and well-derived ...We apply stochastic seismic inversion and Bayesian facies classification for porosity modeling and igneous rock identification in the presalt interval of the Santos Basin. This integration of seismic and well-derived information enhances reservoir characterization. Stochastic inversion and Bayesian classification are powerful tools because they permit addressing the uncertainties in the model. We used the ES-MDA algorithm to achieve the realizations equivalent to the percentiles P10, P50, and P90 of acoustic impedance, a novel method for acoustic inversion in presalt. The facies were divided into five: reservoir 1,reservoir 2, tight carbonates, clayey rocks, and igneous rocks. To deal with the overlaps in acoustic impedance values of facies, we included geological information using a priori probability, indicating that structural highs are reservoir-dominated. To illustrate our approach, we conducted porosity modeling using facies-related rock-physics models for rock-physics inversion in an area with a well drilled in a coquina bank and evaluated the thickness and extension of an igneous intrusion near the carbonate-salt interface. The modeled porosity and the classified seismic facies are in good agreement with the ones observed in the wells. Notably, the coquinas bank presents an improvement in the porosity towards the top. The a priori probability model was crucial for limiting the clayey rocks to the structural lows. In Well B, the hit rate of the igneous rock in the three scenarios is higher than 60%, showing an excellent thickness-prediction capability.展开更多
Stochastic seismic inversion is the combination of geostatistics and seismic inversion technology which integrates information from seismic records, well logs, and geostatistics into a posterior probability density fu...Stochastic seismic inversion is the combination of geostatistics and seismic inversion technology which integrates information from seismic records, well logs, and geostatistics into a posterior probability density function （PDF） of subsurface models. The Markov chain Monte Carlo （MCMC） method is used to sample the posterior PDF and the subsurface model characteristics can be inferred by analyzing a set of the posterior PDF samples. In this paper, we first introduce the stochastic seismic inversion theory, discuss and analyze the four key parameters： seismic data signal-to-noise ratio （S/N）, variogram, the posterior PDF sample number, and well density, and propose the optimum selection of these parameters. The analysis results show that seismic data S/N adjusts the compromise between the influence of the seismic data and geostatistics on the inversion results, the variogram controls the smoothness of the inversion results, the posterior PDF sample number determines the reliability of the statistical characteristics derived from the samples, and well density influences the inversion uncertainty. Finally, the comparison between the stochastic seismic inversion and the deterministic model based seismic inversion indicates that the stochastic seismic inversion can provide more reliable information of the subsurface character.展开更多
Young's modulus and Poisson's ratio are crucial parameters for reservoir characterization and rock brittleness evaluation.Conventional methods often rely on indirect computation or approximations of the Zoeppr...Young's modulus and Poisson's ratio are crucial parameters for reservoir characterization and rock brittleness evaluation.Conventional methods often rely on indirect computation or approximations of the Zoeppritz equations to estimate Young's modulus,which can introduce cumulative errors and reduce the accuracy of inversion results.To address these issues,this paper introduces the analytical solution of the Zoeppritz equation into the inversion process.The equation is re-derived and expressed in terms of Young's modulus,Poisson's ratio,and density.Within the Bayesian framework,we construct an objective function for the joint inversion of PP and PS waves.Traditional gradient-based algorithms often suffer from low precision and the computational complexity.In this study,we address limitations of conventional approaches related to low precision and complicated code by using Circle chaotic mapping,Levy flights,and Gaussian mutation to optimize the quantum particle swarm optimization(QPSO),named improved quantum particle swarm optimization(IQPSO).The IQPSO demonstrates superior global optimization capabilities.We test the proposed inversion method with both synthetic and field data.The test results demonstrate the proposed method's feasibility and effectiveness,indicating an improvement in inversion accuracy over traditional methods.展开更多
The reliable estimation of the wavenumber space(k-space)of the plates remains a longterm concern for acoustic modeling and structural dynamic behavior characterization.Most current analyses of wavenumber identificatio...The reliable estimation of the wavenumber space(k-space)of the plates remains a longterm concern for acoustic modeling and structural dynamic behavior characterization.Most current analyses of wavenumber identification methods are based on the deterministic hypothesis.To this end,an inverse method is proposed for identifying wave propagation characteristics of twodimensional structures under stochastic conditions,such as wavenumber space,dispersion curves,and band gaps.The proposed method is developed based on an algebraic identification scheme in the polar coordinate system framework,thus named Algebraic K-Space Identification(AKSI)technique.Additionally,a model order estimation strategy and a wavenumber filter are proposed to ensure that AKSI is successfully applied.The main benefit of AKSI is that it is a reliable and fast method under four stochastic conditions:(A)High level of signal noise;(B)Small perturbation caused by uncertainties in measurement points’coordinates;(C)Non-periodic sampling;(D)Unknown structural periodicity.To validate the proposed method,we numerically benchmark AKSI and three other inverse methods to extract dispersion curves on three plates under stochastic conditions.One experiment is then performed on an isotropic steel plate.These investigations demonstrate that AKSI is a good in-situ k-space estimator under stochastic conditions.展开更多
With the development of computational power, there has been an increased focus on data-fitting related seismic inversion techniques for high fidelity seismic velocity model and image, such as full-waveform inversion a...With the development of computational power, there has been an increased focus on data-fitting related seismic inversion techniques for high fidelity seismic velocity model and image, such as full-waveform inversion and least squares migration. However, though more advanced than conventional methods, these data fitting methods can be very expensive in terms of computational cost. Recently, various techniques to optimize these data-fitting seismic inversion problems have been implemented to cater for the industrial need for much improved efficiency. In this study, we propose a general stochastic conjugate gradient method for these data-fitting related inverse problems. We first prescribe the basic theory of our method and then give synthetic examples. Our numerical experiments illustrate the potential of this method for large-size seismic inversion application.展开更多
In the conventional stochastic inversion method,the spatial structure information of underground strata is usually characterized by variograms.However,effectively characterizing the heterogeneity of complex strata is ...In the conventional stochastic inversion method,the spatial structure information of underground strata is usually characterized by variograms.However,effectively characterizing the heterogeneity of complex strata is difficult.In this paper,multiple parameters are used to fully explore the underground formation information in the known seismic reflection and well log data.The spatial structure characteristics of complex underground reservoirs are described more comprehensively using multiple statistical characteristic parameters.We propose a prestack seismic stochastic inversion method based on prior information on statistical characteristic parameters.According to the random medium theory,this method obtains several statistical characteristic parameters from known seismic and logging data,constructs a prior information model that meets the spatial structure characteristics of the underground strata,and integrates multiparameter constraints into the likelihood function to construct the objective function.The very fast quantum annealing algorithm is used to optimize and update the objective function to obtain the fi nal inversion result.The model test shows that compared with the traditional prior information model construction method,the prior information model based on multiple parameters in this paper contains more detailed stratigraphic information,which can better describe complex underground reservoirs.A real data analysis shows that the stochastic inversion method proposed in this paper can effectively predict the geophysical characteristics of complex underground reservoirs and has a high resolution.展开更多
The Bayesian inversion method is a stochastic approach based on the Bayesian theory.With the development of sampling algorithms and computer technologies,the Bayesian inversion method has been widely used in geophysic...The Bayesian inversion method is a stochastic approach based on the Bayesian theory.With the development of sampling algorithms and computer technologies,the Bayesian inversion method has been widely used in geophysical inversion problems.In this study,we conduct inversion experiments using crosshole seismic travel-time data to examine the characteristics and performance of the stochastic Bayesian inversion based on the Markov chain Monte Carlo sampling scheme and the traditional deterministic inversion with Tikhonov regularization.Velocity structures with two different spatial variations are considered,one with a chessboard pattern and the other with an interface mimicking the Mohorovicicdiscontinuity(Moho).Inversions are carried out with different scenarios of model discretization and source–receiver configurations.Results show that the Bayesian method yields more robust single-model estimations than the deterministic method,with smaller model errors.In addition,the Bayesian method provides the posterior probabilistic distribution function of the model space,which can help us evaluate the quality of the inversion result.展开更多
The major storage space types in the carbonate reservoir in the Ordovician in the TZ45 area are secondary dissolution caves.For the prediction of caved carbonate reservoir,post-stack methods are commonly used in the o...The major storage space types in the carbonate reservoir in the Ordovician in the TZ45 area are secondary dissolution caves.For the prediction of caved carbonate reservoir,post-stack methods are commonly used in the oilfield at present since pre-stack inversion is always limited by poor seismic data quality and insufficient logging data.In this paper,based on amplitude preserved seismic data processing and rock-physics analysis,pre-stack inversion is employed to predict the caved carbonate reservoir in TZ45 area by seriously controlling the quality of inversion procedures.These procedures mainly include angle-gather conversion,partial stack,wavelet estimation,low-frequency model building and inversion residual analysis.The amplitude-preserved data processing method can achieve high quality data based on the principle that they are very consistent with the synthetics.Besides,the foundation of pre-stack inversion and reservoir prediction criterion can be established by the connection between reservoir property and seismic reflection through rock-physics analysis.Finally,the inversion result is consistent with drilling wells in most cases.It is concluded that integrated with amplitude-preserved processing and rock-physics,pre-stack inversion can be effectively applied in the caved carbonate reservoir prediction.展开更多
The extensive application of pre-stack depth migration has produced huge volumes of seismic data,which allows for the possibility of developing seismic inversions of reservoir properties from seismic data in the depth...The extensive application of pre-stack depth migration has produced huge volumes of seismic data,which allows for the possibility of developing seismic inversions of reservoir properties from seismic data in the depth domain.It is difficult to estimate seismic wavelets directly from seismic data due to the nonstationarity of the data in the depth domain.We conduct a velocity transformation of seismic data to make the seismic data stationary and then apply the ridge regression method to estimate a constant seismic wavelet.The estimated constant seismic wavelet is constructed as a set of space-variant seismic wavelets dominated by velocities at different spatial locations.Incorporating the weighted superposition principle,a synthetic seismogram is generated by directly employing the space-variant seismic wavelets in the depth domain.An inversion workflow based on the model-driven method is developed in the depth domain by incorporating the nonlinear conjugate gradient algorithm,which avoids additional data conversions between the time and depth domains.The impedance inversions of the synthetic and field seismic data in the depth domain show good results,which demonstrates that seismic inversion in the depth domain is feasible.The approach provides an alternative for forward numerical analyses and elastic property inversions of depth-domain seismic data.It is advantageous for further studies concerning the stability,accuracy,and efficiency of seismic inversions in the depth domain.展开更多
Joint PP–PS inversion offers better accuracy and resolution than conventional P-wave inversion. P-and S-wave elastic moduli determined through data inversions are key parameters for reservoir evaluation and fluid cha...Joint PP–PS inversion offers better accuracy and resolution than conventional P-wave inversion. P-and S-wave elastic moduli determined through data inversions are key parameters for reservoir evaluation and fluid characterization. In this paper, starting with the exact Zoeppritz equation that relates P-and S-wave moduli, a coefficient that describes the reflections of P-and converted waves is established. This method effectively avoids error introduced by approximations or indirect calculations, thus improving the accuracy of the inversion results. Considering that the inversion problem is ill-posed and that the forward operator is nonlinear, prior constraints on the model parameters and modified low-frequency constraints are also introduced to the objective function to make the problem more tractable. This modified objective function is solved over many iterations to continuously optimize the background values of the velocity ratio, which increases the stability of the inversion process. Tests of various models show that the method effectively improves the accuracy and stability of extracting P and S-wave moduli from underdetermined data. This method can be applied to provide inferences for reservoir exploration and fluid extraction.展开更多
This paper considers a concrete stochastic nonlinear system with stochastic unmeasurable inverse dynamics. Motivated by the concept of integral input-to-state stability （iISS） in deterministic systems and stochastic...This paper considers a concrete stochastic nonlinear system with stochastic unmeasurable inverse dynamics. Motivated by the concept of integral input-to-state stability （iISS） in deterministic systems and stochastic input-to-state stability （SISS） in stochastic systems, a concept of stochastic integral input-to-state stability （SiISS） using Lyapunov functions is first introduced. A constructive strategy is proposed to design a dynamic output feedback control law, which drives the state to the origin almost surely while keeping all other closed-loop signals almost surely bounded. At last, a simulation is given to verify the effectiveness of the control law.展开更多
Classical multi-channel technology can significantly reduce the pre-stack seismic inversion uncertainty, especially for complex geology such as high dipping structures. However, due to the consideration of complex str...Classical multi-channel technology can significantly reduce the pre-stack seismic inversion uncertainty, especially for complex geology such as high dipping structures. However, due to the consideration of complex structure or reflection features, the existing multi-channel inversion methods have to adopt the highly time-consuming strategy of arranging seismic data trace-by-trace, limiting its wide application in pre-stack inversion. A fast pre-stack multi-channel inversion constrained by seismic reflection features has been proposed to address this issue. The key to our method is to re-characterize the reflection features to directly constrain the pre-stack inversion through a Hadamard product operator without rearranging the seismic data. The seismic reflection features can reflect the distribution of the stratum reflection interface, and we obtained them from the post-stack profile by searching the shortest local Euclidean distance between adjacent seismic traces. Instead of directly constructing a large-size reflection features constraint operator advocated by the conventional methods, through decomposing the reflection features along the vertical and horizontal direction at a particular sampling point, we have constructed a computationally well-behaved constraint operator represented by the vertical and horizontal partial derivatives. Based on the Alternating Direction Method of Multipliers (ADMM) optimization, we have derived a fast algorithm for solving the objective function, including Hadamard product operators. Compared with the conventional reflection features constrained inversion, the proposed method is more efficient and accurate, proved on the Overthrust model and a field data set.展开更多
With the development of exploration of oil and gas resources,the requirements for seismic inversion results are getting more accurate.In particular,it is hoped that the distribution patterns of oil and gas reservoirs ...With the development of exploration of oil and gas resources,the requirements for seismic inversion results are getting more accurate.In particular,it is hoped that the distribution patterns of oil and gas reservoirs can be finely characterized,and the seismic inversion results can clearly characterize the location of stratigraphic boundaries and meet the needs of accurate geological description.Specifically,for pre-stack AVO inversion,it is required to be able to distinguish smaller geological targets in the depth or time domain,and clearly depict the vertical boundaries of the geological objects.In response to the above requirements,we introduce the preprocessing regularization of the adaptive edge-preserving smooth filter into the pre-stack AVO elastic parameter inversion to clearly invert the position of layer boundary and improve the accuracy of the inversion results.展开更多
This work explores the inverse stochastic resonance(ISR) induced by bounded noise and the multiple inverse stochastic resonance induced by time delay by constructing a modular neural network, where the modified Oja’s...This work explores the inverse stochastic resonance(ISR) induced by bounded noise and the multiple inverse stochastic resonance induced by time delay by constructing a modular neural network, where the modified Oja’s synaptic learning rule is employed to characterize synaptic plasticity in this network. Meanwhile, the effects of synaptic plasticity on the ISR dynamics are investigated. Through numerical simulations, it is found that the mean firing rate curve under the influence of bounded noise has an inverted bell-like shape, which implies the appearance of ISR. Moreover, synaptic plasticity with smaller learning rate strengthens this ISR phenomenon, while synaptic plasticity with larger learning rate weakens or even destroys it. On the other hand, the mean firing rate curve under the influence of time delay is found to exhibit a decaying oscillatory process, which represents the emergence of multiple ISR. However, the multiple ISR phenomenon gradually weakens until it disappears with increasing noise amplitude. On the same time, synaptic plasticity with smaller learning rate also weakens this multiple ISR phenomenon, while synaptic plasticity with larger learning rate strengthens it. Furthermore, we find that changes of synaptic learning rate can induce the emergence of ISR phenomenon. We hope these obtained results would provide new insights into the study of ISR in neuroscience.展开更多
Inverse stochastic resonance(ISR)is a phenomenon in which the firing activity of a neuron is inhibited at a certain noise level.In this paper,the effects of potassium channel blockage on ISR in single Hodgkin-Huxley n...Inverse stochastic resonance(ISR)is a phenomenon in which the firing activity of a neuron is inhibited at a certain noise level.In this paper,the effects of potassium channel blockage on ISR in single Hodgkin-Huxley neurons and in small-world networks were investigated.For the single neuron,the ion channel noise-induced ISR phenomenon can occur only in a certain small range of potassium channel blockage ratio.Bifurcation analysis showed that this small range is the bistable region regulated by the external bias current.For small-world networks,the effect of non-homogeneous network blockage on ISR was investigated.The network blockage ratio was used to represent the proportion of potassium-channel-blocked neurons to total network neurons.It is found that an increase in network blockage ratio at small coupling strengths results in shorter ISR duration.When the coupling strength is increased,the ISR is more significant in the case of a large network blockage ratio.The ISR phenomenon is determined by the network blockage ratio,the coupling strength,and the ion channel noise.Our results will provide new perspectives on the observation of ISR in neuroscience experiments.展开更多
Variation of reservoir physical properties can cause changes in its elastic parameters. However, this is not a simple linear relation. Furthermore, the lack of observations, data overlap, noise interference, and ideal...Variation of reservoir physical properties can cause changes in its elastic parameters. However, this is not a simple linear relation. Furthermore, the lack of observations, data overlap, noise interference, and idealized models increases the uncertainties of the inversion result. Thus, we propose an inversion method that is different from traditional statistical rock physics modeling. First, we use deterministic and stochastic rock physics models considering the uncertainties of elastic parameters obtained by prestack seismic inversion and introduce weighting coefficients to establish a weighted statistical relation between reservoir and elastic parameters. Second, based on the weighted statistical relation, we use Markov chain Monte Carlo simulations to generate the random joint distribution space of reservoir and elastic parameters that serves as a sample solution space of an objective function. Finally, we propose a fast solution criterion to maximize the posterior probability density and obtain reservoir parameters. The method has high efficiency and application potential.展开更多
Inversion of Young’s modulus,Poisson’s ratio and density from pre-stack seismic data has been proved to be feasible and effective.However,the existing methods do not take full advantage of the prior information.With...Inversion of Young’s modulus,Poisson’s ratio and density from pre-stack seismic data has been proved to be feasible and effective.However,the existing methods do not take full advantage of the prior information.Without considering the lateral continuity of the inversion results,these methods need to invert the reflectivity first.In this paper,we propose multi-gather simultaneous inversion for pre-stack seismic data.Meanwhile,the total variation（TV）regularization,L1 norm regularization and initial model constraint are used.In order to solve the objective function contains L1norm,TV norm and L2 norm,we develop an algorithm based on split Bregman iteration.The main advantages of our method are as follows：（1）The elastic parameters are calculated directly from objective function rather than from their reflectivity,therefore the stability and accuracy of the inversion process can be ensured.（2）The inversion results are more in accordance with the prior geological information.（3）The lateral continuity of the inversion results are improved.The proposed method is illustrated by theoretical model data and experimented with a 2-D field data.展开更多
Numerical simulation of concrete-faced rockfill dams(CFRDs)considering the spatial variability of rockfill has become a popular research topic in recent years.In order to determine uncertain rockfill properties effici...Numerical simulation of concrete-faced rockfill dams(CFRDs)considering the spatial variability of rockfill has become a popular research topic in recent years.In order to determine uncertain rockfill properties efficiently and reliably,this study developed an uncertainty inversion analysis method for rockfill material parameters using the stacking ensemble strategy and Jaya optimizer.The comprehensive implementation process of the proposed model was described with an illustrative CFRD example.First,the surrogate model method using the stacking ensemble algorithm was used to conduct the Monte Carlo stochastic finite element calculations with reduced computational cost and improved accuracy.Afterwards,the Jaya algorithm was used to inversely calculate the combination of the coefficient of variation of rockfill material parameters.This optimizer obtained higher accuracy and more significant uncertainty reduction than traditional optimizers.Overall,the developed model effectively identified the random parameters of rockfill materials.This study provided scientific references for uncertainty analysis of CFRDs.In addition,the proposed method can be applied to other similar engineering structures.展开更多
Heavy oil has high density and viscosity, and exhibits viscoelasticity. Gassmann's theory is not suitable for materials saturated with viscoelastic fluids. Directly applying such model leads to unreliable results ...Heavy oil has high density and viscosity, and exhibits viscoelasticity. Gassmann's theory is not suitable for materials saturated with viscoelastic fluids. Directly applying such model leads to unreliable results for seismic inversion of heavy oil reservoir. To describe the viscoelastic behavior of heavy oil, we modeled the elastic properties of heavy oil with varying viscosity and frequency using the Cole-Cole-Maxwell (CCM) model. Then, we used a CCoherent Potential Approximation (CPA) instead of the Gassmann equations to account for the fluid effect, by extending the single-phase fluid condition to two-phase fluid (heavy oil and water) condition, so that partial saturation of heavy oil can be considered. This rock physics model establishes the relationship between the elastic modulus of reservoir rock and viscosity, frequency and saturation. The viscosity of the heavy oil and the elastic moduli and porosity of typical reservoir rock samples were measured in laboratory, which were used for calibration of the rock physics model. The well-calibrated frequency-variant CPA model was applied to the prediction of the P- and S-wave velocities in the seismic frequency range (1–100 Hz) and the inversion of petrophysical parameters for a heavy oil reservoir. The pre-stack inversion results of elastic parameters are improved compared with those results using the CPA model in the sonic logging frequency (∼10 kHz), or conventional rock physics model such as the Xu-Payne model. In addition, the inversion of the porosity of the reservoir was conducted with the simulated annealing method, and the result fits reasonably well with the logging curve and depicts the location of the heavy oil reservoir on the time slice. The application of the laboratory-calibrated CPA model provides better results with the velocity dispersion correction, suggesting the important role of accurate frequency dependent rock physics models in the seismic prediction of heavy oil reservoirs.展开更多
基金Equinor for financing the R&D projectthe Institute of Science and Technology of Petroleum Geophysics of Brazil for supporting this research。
文摘We apply stochastic seismic inversion and Bayesian facies classification for porosity modeling and igneous rock identification in the presalt interval of the Santos Basin. This integration of seismic and well-derived information enhances reservoir characterization. Stochastic inversion and Bayesian classification are powerful tools because they permit addressing the uncertainties in the model. We used the ES-MDA algorithm to achieve the realizations equivalent to the percentiles P10, P50, and P90 of acoustic impedance, a novel method for acoustic inversion in presalt. The facies were divided into five: reservoir 1,reservoir 2, tight carbonates, clayey rocks, and igneous rocks. To deal with the overlaps in acoustic impedance values of facies, we included geological information using a priori probability, indicating that structural highs are reservoir-dominated. To illustrate our approach, we conducted porosity modeling using facies-related rock-physics models for rock-physics inversion in an area with a well drilled in a coquina bank and evaluated the thickness and extension of an igneous intrusion near the carbonate-salt interface. The modeled porosity and the classified seismic facies are in good agreement with the ones observed in the wells. Notably, the coquinas bank presents an improvement in the porosity towards the top. The a priori probability model was crucial for limiting the clayey rocks to the structural lows. In Well B, the hit rate of the igneous rock in the three scenarios is higher than 60%, showing an excellent thickness-prediction capability.
基金supported by the National Major Science and Technology Project of China on Development of Big Oil-Gas Fields and Coalbed Methane (No. 2008ZX05010-002)
文摘Stochastic seismic inversion is the combination of geostatistics and seismic inversion technology which integrates information from seismic records, well logs, and geostatistics into a posterior probability density function （PDF） of subsurface models. The Markov chain Monte Carlo （MCMC） method is used to sample the posterior PDF and the subsurface model characteristics can be inferred by analyzing a set of the posterior PDF samples. In this paper, we first introduce the stochastic seismic inversion theory, discuss and analyze the four key parameters： seismic data signal-to-noise ratio （S/N）, variogram, the posterior PDF sample number, and well density, and propose the optimum selection of these parameters. The analysis results show that seismic data S/N adjusts the compromise between the influence of the seismic data and geostatistics on the inversion results, the variogram controls the smoothness of the inversion results, the posterior PDF sample number determines the reliability of the statistical characteristics derived from the samples, and well density influences the inversion uncertainty. Finally, the comparison between the stochastic seismic inversion and the deterministic model based seismic inversion indicates that the stochastic seismic inversion can provide more reliable information of the subsurface character.
基金supported by Fundamental Research Funds for the Central Universities,CHD300102264715National Key Research and Development Program of China under Grant 2021YFA0716902Natural Science Basic Research Program of Shaanxi 2024JCYBMS-199。
文摘Young's modulus and Poisson's ratio are crucial parameters for reservoir characterization and rock brittleness evaluation.Conventional methods often rely on indirect computation or approximations of the Zoeppritz equations to estimate Young's modulus,which can introduce cumulative errors and reduce the accuracy of inversion results.To address these issues,this paper introduces the analytical solution of the Zoeppritz equation into the inversion process.The equation is re-derived and expressed in terms of Young's modulus,Poisson's ratio,and density.Within the Bayesian framework,we construct an objective function for the joint inversion of PP and PS waves.Traditional gradient-based algorithms often suffer from low precision and the computational complexity.In this study,we address limitations of conventional approaches related to low precision and complicated code by using Circle chaotic mapping,Levy flights,and Gaussian mutation to optimize the quantum particle swarm optimization(QPSO),named improved quantum particle swarm optimization(IQPSO).The IQPSO demonstrates superior global optimization capabilities.We test the proposed inversion method with both synthetic and field data.The test results demonstrate the proposed method's feasibility and effectiveness,indicating an improvement in inversion accuracy over traditional methods.
基金supported by the Lyon Acoustics Center of Lyon University,Francefunded by the China Scholarship Council(CSC)。
文摘The reliable estimation of the wavenumber space(k-space)of the plates remains a longterm concern for acoustic modeling and structural dynamic behavior characterization.Most current analyses of wavenumber identification methods are based on the deterministic hypothesis.To this end,an inverse method is proposed for identifying wave propagation characteristics of twodimensional structures under stochastic conditions,such as wavenumber space,dispersion curves,and band gaps.The proposed method is developed based on an algebraic identification scheme in the polar coordinate system framework,thus named Algebraic K-Space Identification(AKSI)technique.Additionally,a model order estimation strategy and a wavenumber filter are proposed to ensure that AKSI is successfully applied.The main benefit of AKSI is that it is a reliable and fast method under four stochastic conditions:(A)High level of signal noise;(B)Small perturbation caused by uncertainties in measurement points’coordinates;(C)Non-periodic sampling;(D)Unknown structural periodicity.To validate the proposed method,we numerically benchmark AKSI and three other inverse methods to extract dispersion curves on three plates under stochastic conditions.One experiment is then performed on an isotropic steel plate.These investigations demonstrate that AKSI is a good in-situ k-space estimator under stochastic conditions.
基金partially supported by the National Natural Science Foundation of China (No.41230318)
文摘With the development of computational power, there has been an increased focus on data-fitting related seismic inversion techniques for high fidelity seismic velocity model and image, such as full-waveform inversion and least squares migration. However, though more advanced than conventional methods, these data fitting methods can be very expensive in terms of computational cost. Recently, various techniques to optimize these data-fitting seismic inversion problems have been implemented to cater for the industrial need for much improved efficiency. In this study, we propose a general stochastic conjugate gradient method for these data-fitting related inverse problems. We first prescribe the basic theory of our method and then give synthetic examples. Our numerical experiments illustrate the potential of this method for large-size seismic inversion application.
基金the National Science Foundation of China(No.42074136 and U19B2008)the Major National Science and Technology Projects(No.2016ZX05027004-001 and 2016ZX05002-005-009)+1 种基金the Fundamental Research Funds for the Central Universities(No.19CX02007A)China Postdoctoral Science Foundation(No.2020M672170).
文摘In the conventional stochastic inversion method,the spatial structure information of underground strata is usually characterized by variograms.However,effectively characterizing the heterogeneity of complex strata is difficult.In this paper,multiple parameters are used to fully explore the underground formation information in the known seismic reflection and well log data.The spatial structure characteristics of complex underground reservoirs are described more comprehensively using multiple statistical characteristic parameters.We propose a prestack seismic stochastic inversion method based on prior information on statistical characteristic parameters.According to the random medium theory,this method obtains several statistical characteristic parameters from known seismic and logging data,constructs a prior information model that meets the spatial structure characteristics of the underground strata,and integrates multiparameter constraints into the likelihood function to construct the objective function.The very fast quantum annealing algorithm is used to optimize and update the objective function to obtain the fi nal inversion result.The model test shows that compared with the traditional prior information model construction method,the prior information model based on multiple parameters in this paper contains more detailed stratigraphic information,which can better describe complex underground reservoirs.A real data analysis shows that the stochastic inversion method proposed in this paper can effectively predict the geophysical characteristics of complex underground reservoirs and has a high resolution.
基金supported by the National Natural Science Foundation of China (grant nos. 41930103 and 41674052)
文摘The Bayesian inversion method is a stochastic approach based on the Bayesian theory.With the development of sampling algorithms and computer technologies,the Bayesian inversion method has been widely used in geophysical inversion problems.In this study,we conduct inversion experiments using crosshole seismic travel-time data to examine the characteristics and performance of the stochastic Bayesian inversion based on the Markov chain Monte Carlo sampling scheme and the traditional deterministic inversion with Tikhonov regularization.Velocity structures with two different spatial variations are considered,one with a chessboard pattern and the other with an interface mimicking the Mohorovicicdiscontinuity(Moho).Inversions are carried out with different scenarios of model discretization and source–receiver configurations.Results show that the Bayesian method yields more robust single-model estimations than the deterministic method,with smaller model errors.In addition,the Bayesian method provides the posterior probabilistic distribution function of the model space,which can help us evaluate the quality of the inversion result.
基金supported by National Basic Research Program(2006CB202304)of Chinaco-supported by the National Basic Research Program of China(Grant No.2011CB201103)the National Science and Technology Major Project of China(Grant No.2011ZX05004003)
文摘The major storage space types in the carbonate reservoir in the Ordovician in the TZ45 area are secondary dissolution caves.For the prediction of caved carbonate reservoir,post-stack methods are commonly used in the oilfield at present since pre-stack inversion is always limited by poor seismic data quality and insufficient logging data.In this paper,based on amplitude preserved seismic data processing and rock-physics analysis,pre-stack inversion is employed to predict the caved carbonate reservoir in TZ45 area by seriously controlling the quality of inversion procedures.These procedures mainly include angle-gather conversion,partial stack,wavelet estimation,low-frequency model building and inversion residual analysis.The amplitude-preserved data processing method can achieve high quality data based on the principle that they are very consistent with the synthetics.Besides,the foundation of pre-stack inversion and reservoir prediction criterion can be established by the connection between reservoir property and seismic reflection through rock-physics analysis.Finally,the inversion result is consistent with drilling wells in most cases.It is concluded that integrated with amplitude-preserved processing and rock-physics,pre-stack inversion can be effectively applied in the caved carbonate reservoir prediction.
基金supported by the National Natural Science Foundation of China(No.41574130,41874143 and 41374134)the National Science and Technology Major Project of China(No.2016ZX05014-001-009)the Sichuan Provincial Youth Science&Technology Innovative Research Group Fund(No.2016TD0023)
文摘The extensive application of pre-stack depth migration has produced huge volumes of seismic data,which allows for the possibility of developing seismic inversions of reservoir properties from seismic data in the depth domain.It is difficult to estimate seismic wavelets directly from seismic data due to the nonstationarity of the data in the depth domain.We conduct a velocity transformation of seismic data to make the seismic data stationary and then apply the ridge regression method to estimate a constant seismic wavelet.The estimated constant seismic wavelet is constructed as a set of space-variant seismic wavelets dominated by velocities at different spatial locations.Incorporating the weighted superposition principle,a synthetic seismogram is generated by directly employing the space-variant seismic wavelets in the depth domain.An inversion workflow based on the model-driven method is developed in the depth domain by incorporating the nonlinear conjugate gradient algorithm,which avoids additional data conversions between the time and depth domains.The impedance inversions of the synthetic and field seismic data in the depth domain show good results,which demonstrates that seismic inversion in the depth domain is feasible.The approach provides an alternative for forward numerical analyses and elastic property inversions of depth-domain seismic data.It is advantageous for further studies concerning the stability,accuracy,and efficiency of seismic inversions in the depth domain.
基金supported by the National Science and Technology Major Project(No.2016ZX05047-002-001)
文摘Joint PP–PS inversion offers better accuracy and resolution than conventional P-wave inversion. P-and S-wave elastic moduli determined through data inversions are key parameters for reservoir evaluation and fluid characterization. In this paper, starting with the exact Zoeppritz equation that relates P-and S-wave moduli, a coefficient that describes the reflections of P-and converted waves is established. This method effectively avoids error introduced by approximations or indirect calculations, thus improving the accuracy of the inversion results. Considering that the inversion problem is ill-posed and that the forward operator is nonlinear, prior constraints on the model parameters and modified low-frequency constraints are also introduced to the objective function to make the problem more tractable. This modified objective function is solved over many iterations to continuously optimize the background values of the velocity ratio, which increases the stability of the inversion process. Tests of various models show that the method effectively improves the accuracy and stability of extracting P and S-wave moduli from underdetermined data. This method can be applied to provide inferences for reservoir exploration and fluid extraction.
基金supported by National Natural Science Foundation of China (No. 60774010, 10971256, and 60974028)Jiangsu"Six Top Talents" (No. 07-A-020)+2 种基金Natural Science Foundation of Jiangsu Province (No. BK2009083)Program for Fundamental Research of Natural Sciences in Universities of Jiangsu Province(No.07KJB510114)Natural Science Foundation of Xuzhou Normal University (No. 08XLB20)
文摘This paper considers a concrete stochastic nonlinear system with stochastic unmeasurable inverse dynamics. Motivated by the concept of integral input-to-state stability （iISS） in deterministic systems and stochastic input-to-state stability （SISS） in stochastic systems, a concept of stochastic integral input-to-state stability （SiISS） using Lyapunov functions is first introduced. A constructive strategy is proposed to design a dynamic output feedback control law, which drives the state to the origin almost surely while keeping all other closed-loop signals almost surely bounded. At last, a simulation is given to verify the effectiveness of the control law.
基金We would like to acknowledge the sponsorship of the National Natural Science Foundation of China(42004092,42030103,41974119)Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology(Qingdao)(Grant No.2021QNLM020001-6)Young Elite Scientists Sponsorship Program by CAST(2021QNRC001).
文摘Classical multi-channel technology can significantly reduce the pre-stack seismic inversion uncertainty, especially for complex geology such as high dipping structures. However, due to the consideration of complex structure or reflection features, the existing multi-channel inversion methods have to adopt the highly time-consuming strategy of arranging seismic data trace-by-trace, limiting its wide application in pre-stack inversion. A fast pre-stack multi-channel inversion constrained by seismic reflection features has been proposed to address this issue. The key to our method is to re-characterize the reflection features to directly constrain the pre-stack inversion through a Hadamard product operator without rearranging the seismic data. The seismic reflection features can reflect the distribution of the stratum reflection interface, and we obtained them from the post-stack profile by searching the shortest local Euclidean distance between adjacent seismic traces. Instead of directly constructing a large-size reflection features constraint operator advocated by the conventional methods, through decomposing the reflection features along the vertical and horizontal direction at a particular sampling point, we have constructed a computationally well-behaved constraint operator represented by the vertical and horizontal partial derivatives. Based on the Alternating Direction Method of Multipliers (ADMM) optimization, we have derived a fast algorithm for solving the objective function, including Hadamard product operators. Compared with the conventional reflection features constrained inversion, the proposed method is more efficient and accurate, proved on the Overthrust model and a field data set.
基金support of China national key project 41904130 and key research project 041020080060.
文摘With the development of exploration of oil and gas resources,the requirements for seismic inversion results are getting more accurate.In particular,it is hoped that the distribution patterns of oil and gas reservoirs can be finely characterized,and the seismic inversion results can clearly characterize the location of stratigraphic boundaries and meet the needs of accurate geological description.Specifically,for pre-stack AVO inversion,it is required to be able to distinguish smaller geological targets in the depth or time domain,and clearly depict the vertical boundaries of the geological objects.In response to the above requirements,we introduce the preprocessing regularization of the adaptive edge-preserving smooth filter into the pre-stack AVO elastic parameter inversion to clearly invert the position of layer boundary and improve the accuracy of the inversion results.
基金the National Natural Science Foundation of China(Grant No.11972217).
文摘This work explores the inverse stochastic resonance(ISR) induced by bounded noise and the multiple inverse stochastic resonance induced by time delay by constructing a modular neural network, where the modified Oja’s synaptic learning rule is employed to characterize synaptic plasticity in this network. Meanwhile, the effects of synaptic plasticity on the ISR dynamics are investigated. Through numerical simulations, it is found that the mean firing rate curve under the influence of bounded noise has an inverted bell-like shape, which implies the appearance of ISR. Moreover, synaptic plasticity with smaller learning rate strengthens this ISR phenomenon, while synaptic plasticity with larger learning rate weakens or even destroys it. On the other hand, the mean firing rate curve under the influence of time delay is found to exhibit a decaying oscillatory process, which represents the emergence of multiple ISR. However, the multiple ISR phenomenon gradually weakens until it disappears with increasing noise amplitude. On the same time, synaptic plasticity with smaller learning rate also weakens this multiple ISR phenomenon, while synaptic plasticity with larger learning rate strengthens it. Furthermore, we find that changes of synaptic learning rate can induce the emergence of ISR phenomenon. We hope these obtained results would provide new insights into the study of ISR in neuroscience.
基金the National Natural Science Foundation of China(No.12175080)the Fundamental Research Funds for the Central Universities(No.CCNU22JC009),China.
文摘Inverse stochastic resonance(ISR)is a phenomenon in which the firing activity of a neuron is inhibited at a certain noise level.In this paper,the effects of potassium channel blockage on ISR in single Hodgkin-Huxley neurons and in small-world networks were investigated.For the single neuron,the ion channel noise-induced ISR phenomenon can occur only in a certain small range of potassium channel blockage ratio.Bifurcation analysis showed that this small range is the bistable region regulated by the external bias current.For small-world networks,the effect of non-homogeneous network blockage on ISR was investigated.The network blockage ratio was used to represent the proportion of potassium-channel-blocked neurons to total network neurons.It is found that an increase in network blockage ratio at small coupling strengths results in shorter ISR duration.When the coupling strength is increased,the ISR is more significant in the case of a large network blockage ratio.The ISR phenomenon is determined by the network blockage ratio,the coupling strength,and the ion channel noise.Our results will provide new perspectives on the observation of ISR in neuroscience experiments.
基金supported by the National Science and Technology Major Project(No.2011 ZX05007-006)the 973 Program of China(No.2013CB228604)the Major Project of Petrochina(No.2014B-0610)
文摘Variation of reservoir physical properties can cause changes in its elastic parameters. However, this is not a simple linear relation. Furthermore, the lack of observations, data overlap, noise interference, and idealized models increases the uncertainties of the inversion result. Thus, we propose an inversion method that is different from traditional statistical rock physics modeling. First, we use deterministic and stochastic rock physics models considering the uncertainties of elastic parameters obtained by prestack seismic inversion and introduce weighting coefficients to establish a weighted statistical relation between reservoir and elastic parameters. Second, based on the weighted statistical relation, we use Markov chain Monte Carlo simulations to generate the random joint distribution space of reservoir and elastic parameters that serves as a sample solution space of an objective function. Finally, we propose a fast solution criterion to maximize the posterior probability density and obtain reservoir parameters. The method has high efficiency and application potential.
基金supported by the National Natural Science Foundation of China (Nos.61775030,61571096,41301460,61362018,and 41274127)the key projects of Hunan Provincial Department of Education (No.16A174)
文摘Inversion of Young’s modulus,Poisson’s ratio and density from pre-stack seismic data has been proved to be feasible and effective.However,the existing methods do not take full advantage of the prior information.Without considering the lateral continuity of the inversion results,these methods need to invert the reflectivity first.In this paper,we propose multi-gather simultaneous inversion for pre-stack seismic data.Meanwhile,the total variation（TV）regularization,L1 norm regularization and initial model constraint are used.In order to solve the objective function contains L1norm,TV norm and L2 norm,we develop an algorithm based on split Bregman iteration.The main advantages of our method are as follows：（1）The elastic parameters are calculated directly from objective function rather than from their reflectivity,therefore the stability and accuracy of the inversion process can be ensured.（2）The inversion results are more in accordance with the prior geological information.（3）The lateral continuity of the inversion results are improved.The proposed method is illustrated by theoretical model data and experimented with a 2-D field data.
基金supported by the National Natural Science Foundation of China(Grants No.51879185 and 52179139)the Open Fund of the Hubei Key Laboratory of Construction and Management in Hydropower Engineering(Grant No.2020KSD06).
文摘Numerical simulation of concrete-faced rockfill dams(CFRDs)considering the spatial variability of rockfill has become a popular research topic in recent years.In order to determine uncertain rockfill properties efficiently and reliably,this study developed an uncertainty inversion analysis method for rockfill material parameters using the stacking ensemble strategy and Jaya optimizer.The comprehensive implementation process of the proposed model was described with an illustrative CFRD example.First,the surrogate model method using the stacking ensemble algorithm was used to conduct the Monte Carlo stochastic finite element calculations with reduced computational cost and improved accuracy.Afterwards,the Jaya algorithm was used to inversely calculate the combination of the coefficient of variation of rockfill material parameters.This optimizer obtained higher accuracy and more significant uncertainty reduction than traditional optimizers.Overall,the developed model effectively identified the random parameters of rockfill materials.This study provided scientific references for uncertainty analysis of CFRDs.In addition,the proposed method can be applied to other similar engineering structures.
基金supported by NSFC(41930425)Science Foundation of China University of Petroleum,Beijing(No.2462020YXZZ008)+1 种基金R&D Department of China National Petroleum Corporation(Investigations on fundamental experiments and advanced theoretical methods in geophysical prospecting applications(2022DQ0604-01)the Strategic Cooperation Technology Projects of CNPC and CUPB(ZLZX2020-03)and NSFC(42274142).
文摘Heavy oil has high density and viscosity, and exhibits viscoelasticity. Gassmann's theory is not suitable for materials saturated with viscoelastic fluids. Directly applying such model leads to unreliable results for seismic inversion of heavy oil reservoir. To describe the viscoelastic behavior of heavy oil, we modeled the elastic properties of heavy oil with varying viscosity and frequency using the Cole-Cole-Maxwell (CCM) model. Then, we used a CCoherent Potential Approximation (CPA) instead of the Gassmann equations to account for the fluid effect, by extending the single-phase fluid condition to two-phase fluid (heavy oil and water) condition, so that partial saturation of heavy oil can be considered. This rock physics model establishes the relationship between the elastic modulus of reservoir rock and viscosity, frequency and saturation. The viscosity of the heavy oil and the elastic moduli and porosity of typical reservoir rock samples were measured in laboratory, which were used for calibration of the rock physics model. The well-calibrated frequency-variant CPA model was applied to the prediction of the P- and S-wave velocities in the seismic frequency range (1–100 Hz) and the inversion of petrophysical parameters for a heavy oil reservoir. The pre-stack inversion results of elastic parameters are improved compared with those results using the CPA model in the sonic logging frequency (∼10 kHz), or conventional rock physics model such as the Xu-Payne model. In addition, the inversion of the porosity of the reservoir was conducted with the simulated annealing method, and the result fits reasonably well with the logging curve and depicts the location of the heavy oil reservoir on the time slice. The application of the laboratory-calibrated CPA model provides better results with the velocity dispersion correction, suggesting the important role of accurate frequency dependent rock physics models in the seismic prediction of heavy oil reservoirs.
