The two-component cold atom systems with anisotropic hopping amplitudes can be phenomenologically described by a two-dimensional Ising-XY coupled model with spatial anisotropy.At low temperatures,theoretical predictio...The two-component cold atom systems with anisotropic hopping amplitudes can be phenomenologically described by a two-dimensional Ising-XY coupled model with spatial anisotropy.At low temperatures,theoretical predictions[Phys.Rev.A 72053604(2005)]and[arXiv:0706.1609]indicate the existence of a topological ordered phase characterized by Ising and XY disorder but with 2XY ordering.However,due to ergodic difficulties faced by Monte Carlo methods at low temperatures,this topological phase has not been numerically explored.We propose a linear cluster updating Monte Carlo method,which flips spins without rejection in the anisotropy limit but does not change the energy.Using this scheme and conventional Monte Carlo methods,we succeed in revealing the nature of topological phases with half-vortices and domain walls.In the constructed global phase diagram,Ising and XY-type transitions are very close to each other and differ significantly from the schematic phase diagram reported earlier.We also propose and explore a wide range of quantities,including magnetism,superfluidity,specific heat,susceptibility,and even percolation susceptibility,and obtain consistent and reliable results.Furthermore,we observed first-order transitions characterized by common intersection points in magnetizations for different system sizes,as opposed to the conventional phase transition where Binder cumulants of various sizes share common intersections.The critical exponents of different types of phase transitions are reasonably fitted.The results are useful to help cold atom experiments explore the half-vortex topological phase.展开更多
Straightness error is an important parameter in measuring high-precision shafts. New generation geometrical product speeifieation（GPS） requires the measurement uncertainty characterizing the reliability of the resul...Straightness error is an important parameter in measuring high-precision shafts. New generation geometrical product speeifieation（GPS） requires the measurement uncertainty characterizing the reliability of the results should be given together when the measurement result is given. Nowadays most researches on straightness focus on error calculation and only several research projects evaluate the measurement uncertainty based on ＂The Guide to the Expression of Uncertainty in Measurement（GUM）＂. In order to compute spatial straightness error（SSE） accurately and rapidly and overcome the limitations of GUM, a quasi particle swarm optimization（QPSO） is proposed to solve the minimum zone SSE and Monte Carlo Method（MCM） is developed to estimate the measurement uncertainty. The mathematical model of minimum zone SSE is formulated. In QPSO quasi-random sequences are applied to the generation of the initial position and velocity of particles and their velocities are modified by the constriction factor approach. The flow of measurement uncertainty evaluation based on MCM is proposed, where the heart is repeatedly sampling from the probability density function（PDF） for every input quantity and evaluating the model in each case. The minimum zone SSE of a shaft measured on a Coordinate Measuring Machine（CMM） is calculated by QPSO and the measurement uncertainty is evaluated by MCM on the basis of analyzing the uncertainty contributors. The results show that the uncertainty directly influences the product judgment result. Therefore it is scientific and reasonable to consider the influence of the uncertainty in judging whether the parts are accepted or rejected, especially for those located in the uncertainty zone. The proposed method is especially suitable when the PDF of the measurand cannot adequately be approximated by a Gaussian distribution or a scaled and shifted t-distribution and the measurement model is non-linear.展开更多
A two-dimensional （2D） full band self-consistent ensemble Monte Carlo （MC） method for solving the quantum Boltzmann equation, including collision broadening and quantum potential corrections, is developed to exten...A two-dimensional （2D） full band self-consistent ensemble Monte Carlo （MC） method for solving the quantum Boltzmann equation, including collision broadening and quantum potential corrections, is developed to extend the MC method to the study of nano-scale semiconductor devices with obvious quantum mechanical （QM） effects. The quantum effects both in real space and momentum space in nano-scale semiconductor devices can be simulated. The effective mobility in the inversion layer of n and p channel MOSFET is simulated and compared with experimental data to verify this method. With this method 50nm ultra thin body silicon on insulator MOSFET are simulated. Results indicate that this method can be used to simulate the 2D QM effects in semiconductor devices including tunnelling effect.展开更多
In loosely coupled or large-scale problems with high dominance ratios,slow fission source convergence can take extremely long time,reducing Monte Carlo(MC)criticality calculation efficiency.Although various accelerati...In loosely coupled or large-scale problems with high dominance ratios,slow fission source convergence can take extremely long time,reducing Monte Carlo(MC)criticality calculation efficiency.Although various acceleration methods have been developed,some methods cannot reduce convergence times,whereas others have been limited to specific problem geometries.In this study,a new fission source convergence acceleration(FSCA)method,the forced propagation(FP)method,has been proposed,which forces the fission source to propagate and accelerate fission source convergence.Additionally,some stabilization techniques have been designed to render the method more practical.The resulting stabilized method was then successfully implemented in the MC transport code,and its feasibility and effectiveness were tested using the modified OECD/NEA,one-dimensional slab benchmark,and the Hoogenboom full-core problem.The comparison results showed that the FP method was able to achieve efficient FSCA.展开更多
One of the most critical and complicated steps in mine design is a selection of suitable mining method based upon geological,geotechnical,geographical,safety and economical parameters.The aim of this study is developi...One of the most critical and complicated steps in mine design is a selection of suitable mining method based upon geological,geotechnical,geographical,safety and economical parameters.The aim of this study is developing a Monte Carlo simulation to selection the optimum mining method by using effective and major criteria and at the same time,taking subjective judgments of decision makers into consideration.Proposed approach is based on the combination of Monte Carlo simulation with conventional Analytic Hierarchy Process(AHP).Monte Carlo simulation is used to determine the confdence level of each alternative’s score,is calculated by AHP,with the respect to the variance of decision makers’opinion.The proposed method is applied for Jajarm Bauxite Mine in Iran and eventually the most appropriate mining methods for this mine are ranked.展开更多
A stratified sampling Monte Carlo method to analyze the reliability of structural systems is presented. Introducing a small exploratory simulation, this method overcomes the difficulties for getting the systematic sam...A stratified sampling Monte Carlo method to analyze the reliability of structural systems is presented. Introducing a small exploratory simulation, this method overcomes the difficulties for getting the systematic sampling probability of all the strata. Several useful and efficient stratification methods are given and the strategies of stratification and simulation are studied. A general conclusion has been presented corresponding to actual engineering structures. The strict theoretical proof has been given,and it is especially effective to solve probabilistic integration. Statistic error of evaluating failure probability is reduced obviously. Especially in highly non-linear and nonreonvex problems, it is more accurate than other methods. Compared with other variance reduction techniques, this method can obtain a more obvious variance reduction and an increased sampling efficiency. Moreover, without strict limiting condition, it is convenient to use. This method is especially suitable to solve the reliability problem of structural systems with multiple failure modes and highly non-linear safety margin equations.展开更多
A new hybrid method, Monte-Carlo-Heat-Flux （MCHF） method, was presented to analyze the radiative heat transfer of participating medium in a three-dimensional rectangular enclosure using combined the Monte-Carlo meth...A new hybrid method, Monte-Carlo-Heat-Flux （MCHF） method, was presented to analyze the radiative heat transfer of participating medium in a three-dimensional rectangular enclosure using combined the Monte-Carlo method with the heat flux method. Its accuracy and reliability was proved by comparing the computational results with exact results from classical ＂Zone Method＂.展开更多
Monte-Carlo method is used for estimating coalbed methane （CBM） resources in key coal mining areas of China. Monte-Carlo method is shown to be superior to the traditional volumetric method with constant parameters i...Monte-Carlo method is used for estimating coalbed methane （CBM） resources in key coal mining areas of China. Monte-Carlo method is shown to be superior to the traditional volumetric method with constant parameters in the calculation of CBM resources. The focus of the article is to introduce the main algorithm and the realization of functions estimated by Monte-Carlo method, including selection of parameters, determination of distribution function, generation of pseudo-random numbers, and evaluation of the parameters corresponding to pseudo-random numbers. A specified software on the basis of Monte-Carlo method is developed using Visual C＋＋ for the assessment of the CBM resources. A case study shows that calculation results using Monte-Carlo method have smaller error range in comparison with those using volumetric method.展开更多
Tooth modification technique is widely used in gear industry to improve the meshing performance of gearings. However, few of the present studies on tooth modification considers the influence of inevitable random error...Tooth modification technique is widely used in gear industry to improve the meshing performance of gearings. However, few of the present studies on tooth modification considers the influence of inevitable random errors on gear modification effects. In order to investigate the uncertainties of tooth modification amount variations on system＇s dynamic behaviors of a helical planetary gears, an analytical dynamic model including tooth modification parameters is proposed to carry out a deterministic analysis on the dynamics of a helical planetary gear. The dynamic meshing forces as well as the dynamic transmission errors of the sun-planet 1 gear pair with and without tooth modifications are computed and compared to show the effectiveness of tooth modifications on gear dynamics enhancement. By using response surface method, a fitted regression model for the dynamic transmission error（DTE） fluctuations is established to quantify the relationship between modification amounts and DTE fluctuations. By shifting the inevitable random errors arousing from manufacturing and installing process to tooth modification amount variations, a statistical tooth modification model is developed and a methodology combining Monte Carlo simulation and response surface method is presented for uncertainty analysis of tooth modifications. The uncertainly analysis reveals that the system＇s dynamic behaviors do not obey the normal distribution rule even though the design variables are normally distributed. In addition, a deterministic modification amount will not definitely achieve an optimal result for both static and dynamic transmission error fluctuation reduction simultaneously.展开更多
Measurement uncertainty plays an important role in laser tracking measurement analyses. In the present work, the guides to the expression of uncertainty in measurement(GUM) uncertainty framework(GUF) and its supplemen...Measurement uncertainty plays an important role in laser tracking measurement analyses. In the present work, the guides to the expression of uncertainty in measurement(GUM) uncertainty framework(GUF) and its supplement, the Monte Carlo method, were used to estimate the uncertainty of task-specific laser tracker measurements. First, the sources of error in laser tracker measurement were analyzed in detail, including instruments, measuring network fusion, measurement strategies, measurement process factors(such as the operator), measurement environment, and task-specific data processing. Second, the GUM and Monte Carlo methods and their application to laser tracker measurement were presented. Finally, a case study involving the uncertainty estimation of a cylindricity measurement process using the GUF and Monte Carlo methods was illustrated. The expanded uncertainty results(at 95% confidence levels) obtained with the Monte Carlo method are 0.069 mm(least-squares criterion) and 0.062 mm(minimum zone criterion), respectively, while with the GUM uncertainty framework, none but the result of least-squares criterion can be got, which is 0.071 mm. Thus, the GUM uncertainty framework slightly underestimates the overall uncertainty by 10%. The results demonstrate that the two methods have different characteristics in task-specific uncertainty evaluations of laser tracker measurements. The results indicate that the Monte Carlo method is a practical tool for applying the principle of propagation of distributions and does not depend on the assumptions and limitations required by the law of propagation of uncertainties(GUF). These features of the Monte Carlo method reduce the risk of an unreliable measurement of uncertainty estimation, particularly in cases of complicated measurement models, without the need to evaluate partial derivatives. In addition, the impact of sampling strategy and evaluation method on the uncertainty of the measurement results can also be taken into account with Monte Carlo method, which plays a guiding role in measurement planning.展开更多
Considering the ocean water＇s optical attenuation and the roughness of the sea surface, we analyze the security of continuous-variable （CV） quantum key distribution （QKD） based Mr-to-water channel. The effects of...Considering the ocean water＇s optical attenuation and the roughness of the sea surface, we analyze the security of continuous-variable （CV） quantum key distribution （QKD） based Mr-to-water channel. The effects of the absorp- tion and scattering on the transmittance of underwater quantum channel and the maximum secure transmission distance are studied. Considering the roughness of the sea surface, we simulate the performance bounds of CV QKD with different wind speeds using the Monte Carlo method. The results show that even if the secret key rate gradually reduces as the wind speed increases, the maximum transmission distance will not be affected obviously. Compared to the works regarding short-distance underwater optical communication, our research represents a significant step towards establishing secure communication between air platform and submarine vehicle.展开更多
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.展开更多
基金Project supported by the Hefei National Research Center for Physical Sciences at the Microscale (Grant No.KF2021002)the Natural Science Foundation of Shanxi Province,China (Grant Nos.202303021221029 and 202103021224051)+2 种基金the National Natural Science Foundation of China (Grant Nos.11975024,12047503,and 12275263)the Anhui Provincial Supporting Program for Excellent Young Talents in Colleges and Universities (Grant No.gxyq ZD2019023)the National Key Research and Development Program of China (Grant No.2018YFA0306501)。
文摘The two-component cold atom systems with anisotropic hopping amplitudes can be phenomenologically described by a two-dimensional Ising-XY coupled model with spatial anisotropy.At low temperatures,theoretical predictions[Phys.Rev.A 72053604(2005)]and[arXiv:0706.1609]indicate the existence of a topological ordered phase characterized by Ising and XY disorder but with 2XY ordering.However,due to ergodic difficulties faced by Monte Carlo methods at low temperatures,this topological phase has not been numerically explored.We propose a linear cluster updating Monte Carlo method,which flips spins without rejection in the anisotropy limit but does not change the energy.Using this scheme and conventional Monte Carlo methods,we succeed in revealing the nature of topological phases with half-vortices and domain walls.In the constructed global phase diagram,Ising and XY-type transitions are very close to each other and differ significantly from the schematic phase diagram reported earlier.We also propose and explore a wide range of quantities,including magnetism,superfluidity,specific heat,susceptibility,and even percolation susceptibility,and obtain consistent and reliable results.Furthermore,we observed first-order transitions characterized by common intersection points in magnetizations for different system sizes,as opposed to the conventional phase transition where Binder cumulants of various sizes share common intersections.The critical exponents of different types of phase transitions are reasonably fitted.The results are useful to help cold atom experiments explore the half-vortex topological phase.
基金supported by National Natural Science Foundation of China (Grant No. 51075198)Jiangsu Provincial Natural Science Foundation of China (Grant No. BK2010479)+2 种基金Innovation Research of Nanjing Institute of Technology, China (Grant No. CKJ20100008)Jiangsu Provincial Foundation of 333 Talents Engineering of ChinaJiangsu Provincial Foundation of Six Talented Peak of China
文摘Straightness error is an important parameter in measuring high-precision shafts. New generation geometrical product speeifieation（GPS） requires the measurement uncertainty characterizing the reliability of the results should be given together when the measurement result is given. Nowadays most researches on straightness focus on error calculation and only several research projects evaluate the measurement uncertainty based on ＂The Guide to the Expression of Uncertainty in Measurement（GUM）＂. In order to compute spatial straightness error（SSE） accurately and rapidly and overcome the limitations of GUM, a quasi particle swarm optimization（QPSO） is proposed to solve the minimum zone SSE and Monte Carlo Method（MCM） is developed to estimate the measurement uncertainty. The mathematical model of minimum zone SSE is formulated. In QPSO quasi-random sequences are applied to the generation of the initial position and velocity of particles and their velocities are modified by the constriction factor approach. The flow of measurement uncertainty evaluation based on MCM is proposed, where the heart is repeatedly sampling from the probability density function（PDF） for every input quantity and evaluating the model in each case. The minimum zone SSE of a shaft measured on a Coordinate Measuring Machine（CMM） is calculated by QPSO and the measurement uncertainty is evaluated by MCM on the basis of analyzing the uncertainty contributors. The results show that the uncertainty directly influences the product judgment result. Therefore it is scientific and reasonable to consider the influence of the uncertainty in judging whether the parts are accepted or rejected, especially for those located in the uncertainty zone. The proposed method is especially suitable when the PDF of the measurand cannot adequately be approximated by a Gaussian distribution or a scaled and shifted t-distribution and the measurement model is non-linear.
基金Project supported by the Special Foundation for State Major Basic Research Program of China （Grant No G2000035602） and the National Natural Science Foundation of China （Grant No 90307006）.
文摘A two-dimensional （2D） full band self-consistent ensemble Monte Carlo （MC） method for solving the quantum Boltzmann equation, including collision broadening and quantum potential corrections, is developed to extend the MC method to the study of nano-scale semiconductor devices with obvious quantum mechanical （QM） effects. The quantum effects both in real space and momentum space in nano-scale semiconductor devices can be simulated. The effective mobility in the inversion layer of n and p channel MOSFET is simulated and compared with experimental data to verify this method. With this method 50nm ultra thin body silicon on insulator MOSFET are simulated. Results indicate that this method can be used to simulate the 2D QM effects in semiconductor devices including tunnelling effect.
基金supported by the National Natural Science Foundation of China(Nos.11775126,11545013,11605101)the Young Elite Scientists Sponsorship Program by CAST(No.2016QNRC001)+1 种基金Science Challenge Project by MIIT of China(No.TZ2018001)Tsinghua University,Initiative Scientific Research Program。
文摘In loosely coupled or large-scale problems with high dominance ratios,slow fission source convergence can take extremely long time,reducing Monte Carlo(MC)criticality calculation efficiency.Although various acceleration methods have been developed,some methods cannot reduce convergence times,whereas others have been limited to specific problem geometries.In this study,a new fission source convergence acceleration(FSCA)method,the forced propagation(FP)method,has been proposed,which forces the fission source to propagate and accelerate fission source convergence.Additionally,some stabilization techniques have been designed to render the method more practical.The resulting stabilized method was then successfully implemented in the MC transport code,and its feasibility and effectiveness were tested using the modified OECD/NEA,one-dimensional slab benchmark,and the Hoogenboom full-core problem.The comparison results showed that the FP method was able to achieve efficient FSCA.
文摘One of the most critical and complicated steps in mine design is a selection of suitable mining method based upon geological,geotechnical,geographical,safety and economical parameters.The aim of this study is developing a Monte Carlo simulation to selection the optimum mining method by using effective and major criteria and at the same time,taking subjective judgments of decision makers into consideration.Proposed approach is based on the combination of Monte Carlo simulation with conventional Analytic Hierarchy Process(AHP).Monte Carlo simulation is used to determine the confdence level of each alternative’s score,is calculated by AHP,with the respect to the variance of decision makers’opinion.The proposed method is applied for Jajarm Bauxite Mine in Iran and eventually the most appropriate mining methods for this mine are ranked.
文摘A stratified sampling Monte Carlo method to analyze the reliability of structural systems is presented. Introducing a small exploratory simulation, this method overcomes the difficulties for getting the systematic sampling probability of all the strata. Several useful and efficient stratification methods are given and the strategies of stratification and simulation are studied. A general conclusion has been presented corresponding to actual engineering structures. The strict theoretical proof has been given,and it is especially effective to solve probabilistic integration. Statistic error of evaluating failure probability is reduced obviously. Especially in highly non-linear and nonreonvex problems, it is more accurate than other methods. Compared with other variance reduction techniques, this method can obtain a more obvious variance reduction and an increased sampling efficiency. Moreover, without strict limiting condition, it is convenient to use. This method is especially suitable to solve the reliability problem of structural systems with multiple failure modes and highly non-linear safety margin equations.
基金financially supported by the National Natural Science Foundation of China (No.50464004)
文摘A new hybrid method, Monte-Carlo-Heat-Flux （MCHF） method, was presented to analyze the radiative heat transfer of participating medium in a three-dimensional rectangular enclosure using combined the Monte-Carlo method with the heat flux method. Its accuracy and reliability was proved by comparing the computational results with exact results from classical ＂Zone Method＂.
基金supported by the National Natural Science Foundation of China (No. 40730422)the Foundation of China University of Mining & Technology (No. OF061003)State Key Laboratory of Geological Processes and Mineral Resources (No. GPMR2007-11)
文摘Monte-Carlo method is used for estimating coalbed methane （CBM） resources in key coal mining areas of China. Monte-Carlo method is shown to be superior to the traditional volumetric method with constant parameters in the calculation of CBM resources. The focus of the article is to introduce the main algorithm and the realization of functions estimated by Monte-Carlo method, including selection of parameters, determination of distribution function, generation of pseudo-random numbers, and evaluation of the parameters corresponding to pseudo-random numbers. A specified software on the basis of Monte-Carlo method is developed using Visual C＋＋ for the assessment of the CBM resources. A case study shows that calculation results using Monte-Carlo method have smaller error range in comparison with those using volumetric method.
基金Supported by National Natural Science Foundation of China(Grant No.51375013)Anhui Provincial Natural Science Foundation of China(Grant No.1208085ME64)
文摘Tooth modification technique is widely used in gear industry to improve the meshing performance of gearings. However, few of the present studies on tooth modification considers the influence of inevitable random errors on gear modification effects. In order to investigate the uncertainties of tooth modification amount variations on system＇s dynamic behaviors of a helical planetary gears, an analytical dynamic model including tooth modification parameters is proposed to carry out a deterministic analysis on the dynamics of a helical planetary gear. The dynamic meshing forces as well as the dynamic transmission errors of the sun-planet 1 gear pair with and without tooth modifications are computed and compared to show the effectiveness of tooth modifications on gear dynamics enhancement. By using response surface method, a fitted regression model for the dynamic transmission error（DTE） fluctuations is established to quantify the relationship between modification amounts and DTE fluctuations. By shifting the inevitable random errors arousing from manufacturing and installing process to tooth modification amount variations, a statistical tooth modification model is developed and a methodology combining Monte Carlo simulation and response surface method is presented for uncertainty analysis of tooth modifications. The uncertainly analysis reveals that the system＇s dynamic behaviors do not obey the normal distribution rule even though the design variables are normally distributed. In addition, a deterministic modification amount will not definitely achieve an optimal result for both static and dynamic transmission error fluctuation reduction simultaneously.
基金Project(51318010402)supported by General Armament Department Pre-Research Program of China
文摘Measurement uncertainty plays an important role in laser tracking measurement analyses. In the present work, the guides to the expression of uncertainty in measurement(GUM) uncertainty framework(GUF) and its supplement, the Monte Carlo method, were used to estimate the uncertainty of task-specific laser tracker measurements. First, the sources of error in laser tracker measurement were analyzed in detail, including instruments, measuring network fusion, measurement strategies, measurement process factors(such as the operator), measurement environment, and task-specific data processing. Second, the GUM and Monte Carlo methods and their application to laser tracker measurement were presented. Finally, a case study involving the uncertainty estimation of a cylindricity measurement process using the GUF and Monte Carlo methods was illustrated. The expanded uncertainty results(at 95% confidence levels) obtained with the Monte Carlo method are 0.069 mm(least-squares criterion) and 0.062 mm(minimum zone criterion), respectively, while with the GUM uncertainty framework, none but the result of least-squares criterion can be got, which is 0.071 mm. Thus, the GUM uncertainty framework slightly underestimates the overall uncertainty by 10%. The results demonstrate that the two methods have different characteristics in task-specific uncertainty evaluations of laser tracker measurements. The results indicate that the Monte Carlo method is a practical tool for applying the principle of propagation of distributions and does not depend on the assumptions and limitations required by the law of propagation of uncertainties(GUF). These features of the Monte Carlo method reduce the risk of an unreliable measurement of uncertainty estimation, particularly in cases of complicated measurement models, without the need to evaluate partial derivatives. In addition, the impact of sampling strategy and evaluation method on the uncertainty of the measurement results can also be taken into account with Monte Carlo method, which plays a guiding role in measurement planning.
基金Supported by the National Natural Science Foundation of China under Grant No 61572529
文摘Considering the ocean water＇s optical attenuation and the roughness of the sea surface, we analyze the security of continuous-variable （CV） quantum key distribution （QKD） based Mr-to-water channel. The effects of the absorp- tion and scattering on the transmittance of underwater quantum channel and the maximum secure transmission distance are studied. Considering the roughness of the sea surface, we simulate the performance bounds of CV QKD with different wind speeds using the Monte Carlo method. The results show that even if the secret key rate gradually reduces as the wind speed increases, the maximum transmission distance will not be affected obviously. Compared to the works regarding short-distance underwater optical communication, our research represents a significant step towards establishing secure communication between air platform and submarine vehicle.
基金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.
