In the existing landslide susceptibility prediction(LSP)models,the influences of random errors in landslide conditioning factors on LSP are not considered,instead the original conditioning factors are directly taken a...In the existing landslide susceptibility prediction(LSP)models,the influences of random errors in landslide conditioning factors on LSP are not considered,instead the original conditioning factors are directly taken as the model inputs,which brings uncertainties to LSP results.This study aims to reveal the influence rules of the different proportional random errors in conditioning factors on the LSP un-certainties,and further explore a method which can effectively reduce the random errors in conditioning factors.The original conditioning factors are firstly used to construct original factors-based LSP models,and then different random errors of 5%,10%,15% and 20%are added to these original factors for con-structing relevant errors-based LSP models.Secondly,low-pass filter-based LSP models are constructed by eliminating the random errors using low-pass filter method.Thirdly,the Ruijin County of China with 370 landslides and 16 conditioning factors are used as study case.Three typical machine learning models,i.e.multilayer perceptron(MLP),support vector machine(SVM)and random forest(RF),are selected as LSP models.Finally,the LSP uncertainties are discussed and results show that:(1)The low-pass filter can effectively reduce the random errors in conditioning factors to decrease the LSP uncertainties.(2)With the proportions of random errors increasing from 5%to 20%,the LSP uncertainty increases continuously.(3)The original factors-based models are feasible for LSP in the absence of more accurate conditioning factors.(4)The influence degrees of two uncertainty issues,machine learning models and different proportions of random errors,on the LSP modeling are large and basically the same.(5)The Shapley values effectively explain the internal mechanism of machine learning model predicting landslide sus-ceptibility.In conclusion,greater proportion of random errors in conditioning factors results in higher LSP uncertainty,and low-pass filter can effectively reduce these random errors.展开更多
Real-time capabilities and computational efficiency are provided by parallel image processing utilizing OpenMP. However, race conditions can affect the accuracy and reliability of the outcomes. This paper highlights t...Real-time capabilities and computational efficiency are provided by parallel image processing utilizing OpenMP. However, race conditions can affect the accuracy and reliability of the outcomes. This paper highlights the importance of addressing race conditions in parallel image processing, specifically focusing on color inverse filtering using OpenMP. We considered three solutions to solve race conditions, each with distinct characteristics: #pragma omp atomic: Protects individual memory operations for fine-grained control. #pragma omp critical: Protects entire code blocks for exclusive access. #pragma omp parallel sections reduction: Employs a reduction clause for safe aggregation of values across threads. Our findings show that the produced images were unaffected by race condition. However, it becomes evident that solving the race conditions in the code makes it significantly faster, especially when it is executed on multiple cores.展开更多
Reservoir inversion by production history matching is an important way to decrease the uncertainty of the reservoir description. Ensemble Kalman filter (EnKF) is a new data assimilation method. There are two problem...Reservoir inversion by production history matching is an important way to decrease the uncertainty of the reservoir description. Ensemble Kalman filter (EnKF) is a new data assimilation method. There are two problems have to be solved for the standard EnKF. One is the inconsistency between the updated model and the updated dynamical variables for nonlinear problems, another is the filter divergence caused by the small ensemble size. We improved the EnKF to overcome these two problems. We use the half iterative EnKF (HIEnKF) for reservoir inversion by doing history matching. During the H1EnKF process, the prediction data are obtained by rerunning the reservoir simulator using the updated model. This can guarantee that the updated dynamical variables are consistent with the updated model. The updated model can nonlinearly affect the prediction data. It is proved that HIEnKF is similar to the first iteration of the EnRML method. Covariance localization is introduced to alleviate filter divergence and spurious correlations caused by the small ensemble size. By defining the shape and size of the correlation area, spurious correlation between the gridblocks far apart is alleviated. More freedom of the model ensemble is preserved. The results of history matching and inverse problem obtained from the HIEnKF with covariance localization are improved. The results show that the model freedom increases with a decrease in the correlation length. Therefore the production data can be matched better. But too small a correlation length can lose some reservoir information and this would cause big errors in the reservoir model estimation.展开更多
According to the effective sampling of particles and the particles impoverishment caused by re-sampling in particle filter,an improved particle filtering algorithm based on observation inversion optimal sampling was p...According to the effective sampling of particles and the particles impoverishment caused by re-sampling in particle filter,an improved particle filtering algorithm based on observation inversion optimal sampling was proposed. Firstly,virtual observations were generated from the latest observation,and two sampling strategies were presented. Then,the previous time particles were sampled by utilizing the function inversion relationship between observation and system state. Finally,the current time particles were generated on the basis of the previous time particles and the system one-step state transition model. By the above method,sampling particles can make full use of the latest observation information and the priori modeling information,so that they further approximate the true state. The theoretical analysis and experimental results show that the new algorithm filtering accuracy and real-time outperform obviously the standard particle filter,the extended Kalman particle filter and the unscented particle filter.展开更多
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.展开更多
As an important indicator parameter of fluid identification,fluid factor has always been a concern for scholars.However,when predicting Russell fluid factor or effective pore-fluid bulk modulus,it is necessary to intr...As an important indicator parameter of fluid identification,fluid factor has always been a concern for scholars.However,when predicting Russell fluid factor or effective pore-fluid bulk modulus,it is necessary to introduce a new rock skeleton parameter which is the dry-rock VP/VS ratio squared(DVRS).In the process of fluid factor calculation or inversion,the existing methods take this parameter as a static constant,which has been estimated in advance,and then apply it to the fluid factor calculation and inversion.The fluid identification analysis based on a portion of the Marmousi 2 model and numerical forward modeling test show that,taking the DVRS as a static constant will limit the identification ability of fluid factor and reduce the inversion accuracy.To solve the above problems,we proposed a new method to regard the DVRS as a dynamic variable varying with depth and lithology for the first time,then apply it to fluid factor calculation and inversion.Firstly,the exact Zoeppritz equations are rewritten into a new form containing the fluid factor and DVRS of upper and lower layers.Next,the new equations are applied to the four parameters simultaneous inversion based on the generalized nonlinear inversion(GNI)method.The testing results on a portion of the Marmousi 2 model and field data show that dynamic DVRS can significantly improve the fluid factor identification ability,effectively suppress illusion.Both synthetic and filed data tests also demonstrate that the GNI method based on Bayesian deterministic inversion(BDI)theory can successfully solve the above four parameter simultaneous inversion problem,and taking the dynamic DVRS as a target inversion parameter can effectively improve the inversion accuracy of fluid factor.All these results completely verified the feasibility and effectiveness of the proposed method.展开更多
Aiming at the Four-Dimensional Variation source term inversion algorithm proposed earlier,the observation error regularization factor is introduced to improve the prediction accuracy of the diffusion model,and an impr...Aiming at the Four-Dimensional Variation source term inversion algorithm proposed earlier,the observation error regularization factor is introduced to improve the prediction accuracy of the diffusion model,and an improved Four-Dimensional Variation source term inversion algorithm with observation error regularization(OER-4DVAR STI model)is formed.Firstly,by constructing the inversion process and basic model of OER-4DVAR STI model,its basic principle and logical structure are studied.Secondly,the observation error regularization factor estimation method based on Bayesian optimization is proposed,and the error factor is separated and optimized by two parameters:error statistical time and deviation degree.Finally,the scientific,feasible and advanced nature of the OER-4DVAR STI model are verified by numerical simulation and tracer test data.The experimental results show that OER-4DVAR STI model can better reverse calculate the hazard source term information under the conditions of high atmospheric stability and flat underlying surface.Compared with the previous inversion algorithm,the source intensity estimation accuracy of OER-4DVAR STI model is improved by about 46.97%,and the source location estimation accuracy is improved by about 26.72%.展开更多
The global bathymetry models are usually of low accuracy over the coastline of polar areas due to the harsh climatic environment and the complex topography.Satellite altimetric gravity data can be a supplement and pla...The global bathymetry models are usually of low accuracy over the coastline of polar areas due to the harsh climatic environment and the complex topography.Satellite altimetric gravity data can be a supplement and plays a key role in bathymetry modeling over these regions.The Synthetic Aperture Radar(SAR)altimeters in the missions like CryoSat-2 and Sentinel-3A/3B can relieve waveform contamination that existed in conventional altimeters and provide data with improved accuracy and spatial resolution.In this study,we investigate the potential application of SAR altimetric gravity data in enhancing coastal bathymetry,where the effects on local bathymetry modeling introduced from SAR altimetry data are quantified and evaluated.Furthermore,we study the effects on bathymetry modeling by using different scale factor calculation approaches,where a partition-wise scheme is implemented.The numerical experiment over the South Sandwich Islands near Antarctica suggests that using SARbased altimetric gravity data improves local coastal bathymetry modeling,compared with the model calculated without SAR altimetry data by a magnitude of 3:55 m within 10 km of offshore areas.Moreover,by using the partition-wise scheme for scale factor calculation,the quality of the coastal bathymetry model is improved by 7.34 m compared with the result derived from the traditional method.These results indicate the superiority of using SAR altimetry data in coastal bathymetry inversion.展开更多
基金This work is funded by the National Natural Science Foundation of China(Grant Nos.42377164 and 52079062)the National Science Fund for Distinguished Young Scholars of China(Grant No.52222905).
文摘In the existing landslide susceptibility prediction(LSP)models,the influences of random errors in landslide conditioning factors on LSP are not considered,instead the original conditioning factors are directly taken as the model inputs,which brings uncertainties to LSP results.This study aims to reveal the influence rules of the different proportional random errors in conditioning factors on the LSP un-certainties,and further explore a method which can effectively reduce the random errors in conditioning factors.The original conditioning factors are firstly used to construct original factors-based LSP models,and then different random errors of 5%,10%,15% and 20%are added to these original factors for con-structing relevant errors-based LSP models.Secondly,low-pass filter-based LSP models are constructed by eliminating the random errors using low-pass filter method.Thirdly,the Ruijin County of China with 370 landslides and 16 conditioning factors are used as study case.Three typical machine learning models,i.e.multilayer perceptron(MLP),support vector machine(SVM)and random forest(RF),are selected as LSP models.Finally,the LSP uncertainties are discussed and results show that:(1)The low-pass filter can effectively reduce the random errors in conditioning factors to decrease the LSP uncertainties.(2)With the proportions of random errors increasing from 5%to 20%,the LSP uncertainty increases continuously.(3)The original factors-based models are feasible for LSP in the absence of more accurate conditioning factors.(4)The influence degrees of two uncertainty issues,machine learning models and different proportions of random errors,on the LSP modeling are large and basically the same.(5)The Shapley values effectively explain the internal mechanism of machine learning model predicting landslide sus-ceptibility.In conclusion,greater proportion of random errors in conditioning factors results in higher LSP uncertainty,and low-pass filter can effectively reduce these random errors.
文摘Real-time capabilities and computational efficiency are provided by parallel image processing utilizing OpenMP. However, race conditions can affect the accuracy and reliability of the outcomes. This paper highlights the importance of addressing race conditions in parallel image processing, specifically focusing on color inverse filtering using OpenMP. We considered three solutions to solve race conditions, each with distinct characteristics: #pragma omp atomic: Protects individual memory operations for fine-grained control. #pragma omp critical: Protects entire code blocks for exclusive access. #pragma omp parallel sections reduction: Employs a reduction clause for safe aggregation of values across threads. Our findings show that the produced images were unaffected by race condition. However, it becomes evident that solving the race conditions in the code makes it significantly faster, especially when it is executed on multiple cores.
基金support from the Shandong Natural Science Foundation(Grant No.ZR2010EM053)the Fundamental Research Funds for the Central Universities(Grant No.10CX04042A)
文摘Reservoir inversion by production history matching is an important way to decrease the uncertainty of the reservoir description. Ensemble Kalman filter (EnKF) is a new data assimilation method. There are two problems have to be solved for the standard EnKF. One is the inconsistency between the updated model and the updated dynamical variables for nonlinear problems, another is the filter divergence caused by the small ensemble size. We improved the EnKF to overcome these two problems. We use the half iterative EnKF (HIEnKF) for reservoir inversion by doing history matching. During the H1EnKF process, the prediction data are obtained by rerunning the reservoir simulator using the updated model. This can guarantee that the updated dynamical variables are consistent with the updated model. The updated model can nonlinearly affect the prediction data. It is proved that HIEnKF is similar to the first iteration of the EnRML method. Covariance localization is introduced to alleviate filter divergence and spurious correlations caused by the small ensemble size. By defining the shape and size of the correlation area, spurious correlation between the gridblocks far apart is alleviated. More freedom of the model ensemble is preserved. The results of history matching and inverse problem obtained from the HIEnKF with covariance localization are improved. The results show that the model freedom increases with a decrease in the correlation length. Therefore the production data can be matched better. But too small a correlation length can lose some reservoir information and this would cause big errors in the reservoir model estimation.
基金Project(60634030) supported by the Key Project of the National Natural Science Foundation of ChinaProject(60702066) supported by the National Natural Science Foundation of China+1 种基金Project (2007ZC53037) supported by Aviation Science Foundation of ChinaProject(CASC0214) supported by the Space-Flight Innovation Foundation of China
文摘According to the effective sampling of particles and the particles impoverishment caused by re-sampling in particle filter,an improved particle filtering algorithm based on observation inversion optimal sampling was proposed. Firstly,virtual observations were generated from the latest observation,and two sampling strategies were presented. Then,the previous time particles were sampled by utilizing the function inversion relationship between observation and system state. Finally,the current time particles were generated on the basis of the previous time particles and the system one-step state transition model. By the above method,sampling particles can make full use of the latest observation information and the priori modeling information,so that they further approximate the true state. The theoretical analysis and experimental results show that the new algorithm filtering accuracy and real-time outperform obviously the standard particle filter,the extended Kalman particle filter and the unscented particle filter.
基金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(41904116,41874156,42074167 and 42204135)the Natural Science Foundation of Hunan Province(2020JJ5168)the China Postdoctoral Science Foundation(2021M703629)for their funding of this research.
文摘As an important indicator parameter of fluid identification,fluid factor has always been a concern for scholars.However,when predicting Russell fluid factor or effective pore-fluid bulk modulus,it is necessary to introduce a new rock skeleton parameter which is the dry-rock VP/VS ratio squared(DVRS).In the process of fluid factor calculation or inversion,the existing methods take this parameter as a static constant,which has been estimated in advance,and then apply it to the fluid factor calculation and inversion.The fluid identification analysis based on a portion of the Marmousi 2 model and numerical forward modeling test show that,taking the DVRS as a static constant will limit the identification ability of fluid factor and reduce the inversion accuracy.To solve the above problems,we proposed a new method to regard the DVRS as a dynamic variable varying with depth and lithology for the first time,then apply it to fluid factor calculation and inversion.Firstly,the exact Zoeppritz equations are rewritten into a new form containing the fluid factor and DVRS of upper and lower layers.Next,the new equations are applied to the four parameters simultaneous inversion based on the generalized nonlinear inversion(GNI)method.The testing results on a portion of the Marmousi 2 model and field data show that dynamic DVRS can significantly improve the fluid factor identification ability,effectively suppress illusion.Both synthetic and filed data tests also demonstrate that the GNI method based on Bayesian deterministic inversion(BDI)theory can successfully solve the above four parameter simultaneous inversion problem,and taking the dynamic DVRS as a target inversion parameter can effectively improve the inversion accuracy of fluid factor.All these results completely verified the feasibility and effectiveness of the proposed method.
基金Ministry of Science and Technology of the People’s Republic of China for its support and guidance(Grant No.2018YFC0214100)。
文摘Aiming at the Four-Dimensional Variation source term inversion algorithm proposed earlier,the observation error regularization factor is introduced to improve the prediction accuracy of the diffusion model,and an improved Four-Dimensional Variation source term inversion algorithm with observation error regularization(OER-4DVAR STI model)is formed.Firstly,by constructing the inversion process and basic model of OER-4DVAR STI model,its basic principle and logical structure are studied.Secondly,the observation error regularization factor estimation method based on Bayesian optimization is proposed,and the error factor is separated and optimized by two parameters:error statistical time and deviation degree.Finally,the scientific,feasible and advanced nature of the OER-4DVAR STI model are verified by numerical simulation and tracer test data.The experimental results show that OER-4DVAR STI model can better reverse calculate the hazard source term information under the conditions of high atmospheric stability and flat underlying surface.Compared with the previous inversion algorithm,the source intensity estimation accuracy of OER-4DVAR STI model is improved by about 46.97%,and the source location estimation accuracy is improved by about 26.72%.
基金supported by the National Natural Science Foundation of China(No.42004008)the Natural Science Foundation of Jiangsu Province,China(No.BK20190498)+1 种基金the Fundamental Research Funds for the Central Universities(No.B220202055)the State Scholarship Fund from Chinese Scholarship Council(No.201306270014).
文摘The global bathymetry models are usually of low accuracy over the coastline of polar areas due to the harsh climatic environment and the complex topography.Satellite altimetric gravity data can be a supplement and plays a key role in bathymetry modeling over these regions.The Synthetic Aperture Radar(SAR)altimeters in the missions like CryoSat-2 and Sentinel-3A/3B can relieve waveform contamination that existed in conventional altimeters and provide data with improved accuracy and spatial resolution.In this study,we investigate the potential application of SAR altimetric gravity data in enhancing coastal bathymetry,where the effects on local bathymetry modeling introduced from SAR altimetry data are quantified and evaluated.Furthermore,we study the effects on bathymetry modeling by using different scale factor calculation approaches,where a partition-wise scheme is implemented.The numerical experiment over the South Sandwich Islands near Antarctica suggests that using SARbased altimetric gravity data improves local coastal bathymetry modeling,compared with the model calculated without SAR altimetry data by a magnitude of 3:55 m within 10 km of offshore areas.Moreover,by using the partition-wise scheme for scale factor calculation,the quality of the coastal bathymetry model is improved by 7.34 m compared with the result derived from the traditional method.These results indicate the superiority of using SAR altimetry data in coastal bathymetry inversion.