Optical multilayer thin film structures have been widely used in numerous photonic applications.However,existing inverse design methods have many drawbacks because they either fail to quickly adapt to different design...Optical multilayer thin film structures have been widely used in numerous photonic applications.However,existing inverse design methods have many drawbacks because they either fail to quickly adapt to different design targets,or are difficult to suit for different types of structures,e.g.,designing for different materials at each layer.These methods also cannot accommodate versatile design situations under different angles and polarizations.In addition,how to benefit practical fabrications and manufacturing has not been extensively considered yet.In this work,we introduce OptoGPT(Opto Generative Pretrained Transformer),a decoder-only transformer,to solve all these drawbacks and issues simultaneously.展开更多
In this paper,we investigate the reverse order law for Drazin inverse of three bound-ed linear operators under some commutation relations.Moreover,the Drazin invertibility of sum is also obtained for two bounded linea...In this paper,we investigate the reverse order law for Drazin inverse of three bound-ed linear operators under some commutation relations.Moreover,the Drazin invertibility of sum is also obtained for two bounded linear operators and its expression is presented.展开更多
Flow units(FU)rock typing is a common technique for characterizing reservoir flow behavior,producing reliable porosity and permeability estimation even in complex geological settings.However,the lateral extrapolation ...Flow units(FU)rock typing is a common technique for characterizing reservoir flow behavior,producing reliable porosity and permeability estimation even in complex geological settings.However,the lateral extrapolation of FU away from the well into the whole reservoir grid is commonly a difficult task and using the seismic data as constraints is rarely a subject of study.This paper proposes a workflow to generate numerous possible 3D volumes of flow units,porosity and permeability below the seismic resolution limit,respecting the available seismic data at larger scales.The methodology is used in the Mero Field,a Brazilian presalt carbonate reservoir located in the Santos Basin,who presents a complex and heterogenic geological setting with different sedimentological processes and diagenetic history.We generated metric flow units using the conventional core analysis and transposed to the well log data.Then,given a Markov chain Monte Carlo algorithm,the seismic data and the well log statistics,we simulated acoustic impedance,decametric flow units(DFU),metric flow units(MFU),porosity and permeability volumes in the metric scale.The aim is to estimate a minimum amount of MFU able to calculate realistic scenarios porosity and permeability scenarios,without losing the seismic lateral control.In other words,every porosity and permeability volume simulated produces a synthetic seismic that match the real seismic of the area,even in the metric scale.The achieved 3D results represent a high-resolution fluid flow reservoir modelling considering the lateral control of the seismic during the process and can be directly incorporated in the dynamic characterization workflow.展开更多
Surface-wave inversion is a powerful tool for revealing the Earth's internal structure.However,aside from shear-wave velocity(v_(S)),other parameters can influence the inversion outcomes,yet these have not been sy...Surface-wave inversion is a powerful tool for revealing the Earth's internal structure.However,aside from shear-wave velocity(v_(S)),other parameters can influence the inversion outcomes,yet these have not been systematically discussed.This study investigates the influence of various parameter assumptions on the results of surface-wave inversion,including the compressional and shear velocity ratio(v_(P)/v_(S)),shear-wave attenuation(Q_(S)),density(ρ),Moho interface,and sedimentary layer.We constructed synthetic models to generate dispersion data and compared the obtained results with different parameter assumptions with those of the true model.The results indicate that the v_(P)/v_(S) ratio,Q_(S),and density(ρ) have minimal effects on absolute velocity values and perturbation patterns in the inversion.Conversely,assumptions about the Moho interface and sedimentary layer significantly influenced absolute velocity values and perturbation patterns.Introducing an erroneous Mohointerface depth in the initial model of the inversion significantly affected the v_(S) model near that depth,while using a smooth initial model results in relatively minor deviations.The assumption on the sedimentary layer not only affects shallow structure results but also impacts the result at greater depths.Non-linear inversion methods outperform linear inversion methods,particularly for the assumptions of the Moho interface and sedimentary layer.Joint inversion with other data types,such as receiver functions or Rayleigh wave ellipticity,and using data from a broader period range or higher-mode surface waves,can mitigate these deviations.Furthermore,incorporating more accurate prior information can improve inversion results.展开更多
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.展开更多
The reaction order plays a crucial role in evaluating the response rate of acid-rock.However,the conventional two-scale model typically assumes that the reaction order is constant as one,which can lead to significant ...The reaction order plays a crucial role in evaluating the response rate of acid-rock.However,the conventional two-scale model typically assumes that the reaction order is constant as one,which can lead to significant deviations from reality.To address this issue,this study proposes a novel multi-order dynamic model for acid-rock reaction by combining rotating disk experimental data with theoretical derivation.Through numerical simulations,this model allows for the investigation of the impact of acidification conditions on different orders of reaction,thereby providing valuable insights for on-site construction.The analysis reveals that higher response orders require higher optimal acid liquid flow rates,and lower optimal H+diffusion coefficients,and demonstrate no significant correlation with acid concentration.Consequently,it is recommended to increase the displacement and use high-viscosity acid for reservoirs with high calcite content,while reducing the displacement and using low-viscosity acid for reservoirs with high dolomite content.展开更多
Many important problems in science and engineering require solving the so-called parametric partial differential equations(PDEs),i.e.,PDEs with different physical parameters,boundary conditions,shapes of computational...Many important problems in science and engineering require solving the so-called parametric partial differential equations(PDEs),i.e.,PDEs with different physical parameters,boundary conditions,shapes of computational domains,etc.Typical reduced order modeling techniques accelerate the solution of the parametric PDEs by projecting them onto a linear trial manifold constructed in the ofline stage.These methods often need a predefined mesh as well as a series of precomputed solution snapshots,and may struggle to balance between the efficiency and accuracy due to the limitation of the linear ansatz.Utilizing the nonlinear representation of neural networks(NNs),we propose the Meta-Auto-Decoder(MAD)to construct a nonlinear trial manifold,whose best possible performance is measured theoretically by the decoder width.Based on the meta-learning concept,the trial manifold can be learned in a mesh-free and unsupervised way during the pre-training stage.Fast adaptation to new(possibly heterogeneous)PDE parameters is enabled by searching on this trial manifold,and optionally fine-tuning the trial manifold at the same time.Extensive numerical experiments show that the MAD method exhibits a faster convergence speed without losing the accuracy than other deep learning-based methods.展开更多
Unmanned aerial vehicle transient electromagnetic(UAV-TEM)is a novel airborne exploration method that offers advantages such as low cost,simple operation,high exploration efficiency and suitability for near-surface ex...Unmanned aerial vehicle transient electromagnetic(UAV-TEM)is a novel airborne exploration method that offers advantages such as low cost,simple operation,high exploration efficiency and suitability for near-surface exploration in complex terrain areas.To improve the accuracy of data interpretation in this method,the authors conducted a systematic three-dimensional(3D)forward modeling and inversion of the UAV-TEM.This study utilized the finite element method based on unstructured tetrahedral elements and employed the second-order backward Euler method for time discretization.This allowed for accurate 3D modeling and accounted for the effects of complex terrain.Based on these,the influence characteristics of flight altitudes and the sizes,burial depths,and resistivities of anomalies are compared and analyzed to explore the UAV-TEM systems’exploration capability.Lastly,four typical geoelectrical models of landslides are designed,and the inversion method based on the Gauss-Newton optimization method is used to image the landslide models and analyze the imaging effect of the UAV-TEM method on landslide geohazards.Numerical results showed that UAV-TEM could have better exploration resolution and fine imaging of nearsurface structures,providing important technical support for monitoring,early warning,and preventing landslides and other geological hazards.展开更多
Methane generation in landfills and its inadequate management represent the major avoidable source of anthropogenic methane today. This paper models methane production and the potential resources expected (electrical ...Methane generation in landfills and its inadequate management represent the major avoidable source of anthropogenic methane today. This paper models methane production and the potential resources expected (electrical energy production and potential carbon credits from avoided CH4 emissions) from its proper management in a municipal solid waste landfill located in Ouagadougou, Burkina Faso. The modeling was carried out using two first-order decay (FOD) models (LandGEM V3.02 and SWANA) using parameters evaluated on the basis of the characteristics of the waste admitted to the landfill and weather data for the site. At the same time, production data have been collected since 2016 in order to compare them with the model results. The results obtained from these models were compared to experimental one. For the simulation of methane production, the SWANA model showed better consistency with experimental data, with a coefficient of determination (R²) of 0.59 compared with the LandGEM model, which obtained a coefficient of 0.006. Thus, despite the low correlation values linked to the poor consistency of experimental data, the SWANA model models methane production much better than the LandGEM model. Thus, despite the low correlation values linked to the poor consistency of the experimental data, the SWANA model models methane production much better than the LandGEM V3.02 model. It was noted that the poor consistency of the experimental data justifies these low coefficients, and that they can be improved in the future thanks to ongoing in situ measurements. According to the SWANA model prediction, in 27 years of operation a biogas plant with 33% electrical efficiency using biogas from the Polesgo landfill would avoid 1,340 GgCO2e. Also, the evaluation of revenues due to electricity and carbon credit gave a total revenue derived from methane production of US$27.38 million at a cost of US$10.5/tonne CO2e.展开更多
The medium for next-generation communication is considered as fiber for fast,secure communication and switching capability.Mode division and space division multiplexing provide an excellent switching capability with h...The medium for next-generation communication is considered as fiber for fast,secure communication and switching capability.Mode division and space division multiplexing provide an excellent switching capability with high data transmission rate.In this work,the authors have approached an inverse modeling technique using regression-based machine learning to design a weakly coupled few-mode fiber for facilitating mode division multiplexing.The technique is adapted to predict the accurate profile parameters for the proposed few-mode fiber to obtain the maximum number of modes.It is for a three-ring-core few-mode fiber for guiding five,ten,fifteen,and twenty modes.Three types of regression models namely ordinary least-square linear multi-output regression,k-nearest neighbors of multi-output regression,and ID3 algorithm-based decision trees for multi-output regression are used for predicting the multiple profile parameters.It is observed that the ID3-based decision tree for multioutput regression is the robust,highly-accurate machine learning model for fast modeling of FMFs.The proposed fiber claims to be an efficient candidate for the next-generation 5G and 6G backhaul networks using mode division multiplexing.展开更多
The outbreak of COVID-19 in 2019 resulted in numerous infections and deaths. In order to better study the transmission of COVID-19, this article adopts an improved fractional-order SIR model. Firstly, the properties o...The outbreak of COVID-19 in 2019 resulted in numerous infections and deaths. In order to better study the transmission of COVID-19, this article adopts an improved fractional-order SIR model. Firstly, the properties of the model are studied, including the feasible domain and bounded solutions of the system. Secondly, the stability of the system is discussed, among other things. Then, the GMMP method is introduced to obtain numerical solutions for the COVID-19 system and combined with the improved MH-NMSS-PSO parameter estimation method to fit the real data of Delhi, India from April 1, 2020 to June 30, 2020. The results show that the fitting effect is quite ideal. Finally, long-term predictions were made on the number of infections. We accurately estimate that the peak number of infections in Delhi, India, can reach around 2.1 million. This paper also compares the fitting performance of the integer-order COVID-19 model and the fractional-order COVID-19 model using the real data from Delhi. The results indicate that the fractional-order model with different orders, as we proposed, performs the best.展开更多
This paper investigates an improved SIR model for COVID-19 based on the Caputo fractional derivative. Firstly, the properties of the model are studied, including the feasible domain and bounded solutions of the system...This paper investigates an improved SIR model for COVID-19 based on the Caputo fractional derivative. Firstly, the properties of the model are studied, including the feasible domain and bounded solutions of the system. Secondly, the stability of the system is discussed, among other things. Then, the GMMP method is introduced to obtain numerical solutions for the COVID-19 system. Numerical simulations were conducted using MATLAB, and the results indicate that our model is valuable for studying virus transmission.展开更多
The concept of generalized order statistics has been introduced as a unified approach to a variety of models of ordered random variables with different interpretations. In this paper, we develop methodology for constr...The concept of generalized order statistics has been introduced as a unified approach to a variety of models of ordered random variables with different interpretations. In this paper, we develop methodology for constructing inference based on n selected generalized order statistics (GOS) from inverse Weibull distribution (IWD), Bayesian and non-Bayesian approaches have been used to obtain the estimators of the parameters and reliability function. We have examined Bayes estimates under various losses such as the balanced squared error (balanced SEL) and balanced LINEX loss functions are considered. We show that Bayes estimate under balanced SEL and balanced LINEX loss functions are more general, which include the symmetric and asymmetric losses as special cases. This was done under assumption of discrete-continuous mixture prior for the unknown model parameters. The parametric bootstrap method has been used to construct confidence interval for the parameters and reliability function. Progressively type-II censored and k-record values as a special case of GOS are considered. Finally a practical example using real data set was used for illustration.展开更多
This paper presents a unique and formal method of quantifying the similarity or distance between sedimentary facies successions from measured sections in outcrop or drilled wells and demonstrates its first application...This paper presents a unique and formal method of quantifying the similarity or distance between sedimentary facies successions from measured sections in outcrop or drilled wells and demonstrates its first application in inverse stratigraphic modeling. A sedimentary facies succession is represented with a string of symbols, or facies codes in its natural vertical order, in which each symbol brings with it one attribute such as thickness for the facies. These strings are called attributed strings. A similarity measure is defined between the attributed strings based on a syntactic pattern-recognition technique. A dynamic programming algorithm is used to calculate the similarity. Inverse stratigraphic modeling aims to generate quantitative 3D facies models based on forward stratigraphic modeling that honors observed datasets. One of the key techniques in inverse stratigraphic modeling is how to quantify the similarity or distance between simulated and observed sedimentary facies successions at data locations in order for the forward model to condition the simulation results to the observed dataset such as measured sections or drilled wells. This quantification technique comparing sedimentary successions is demonstrated in the form of a cost function based on the defined distance in our inverse stratigraphic modeling implemented with forward modeling optimization.展开更多
Identifying the unknown geometric and material information of a multi-shield object by analyzing the radiation signature measurements is always an important problem in national and global security. In order to identif...Identifying the unknown geometric and material information of a multi-shield object by analyzing the radiation signature measurements is always an important problem in national and global security. In order to identify the unknown shielding layer thicknesses of a source/shield system with gamma-ray spectra, we have developed a derivative-free inverse radiation transport model based on a differential evolution algorithm with global and local neighbourhoods(IRT-DEGL). In the present paper, the IRT-DEGL model is further extended for estimating the unknown thicknesses with random initial guesses and material mass densities of multi-shielding layers as well as their combinations. Using the detected gamma-ray spectra,the illustration of inverse studies is implemented and the main influence factors for inverse results are also analyzed.展开更多
By using aluminum alloys,the properties of the material in sheet hydroforming were obtained based on the identification of parameters for constitutive models by inverse modeling in which the friction coefficients were...By using aluminum alloys,the properties of the material in sheet hydroforming were obtained based on the identification of parameters for constitutive models by inverse modeling in which the friction coefficients were also considered in 2D and 3D simulations.With consideration of identified simulation parameters by inverse modeling,some key process parameters including tool dimensions and pre-bulging on the forming processes in sheet hydroforming were investigated and optimized.Based on the optimized parameters,the sheet hydroforming process can be analyzed more accurately to improve the robust design.It proves that the results from simulation based on the identified parameters are in good agreement with those from experiments.展开更多
Errors inevitably exist in numerical weather prediction (NWP) due to imperfect numeric and physical parameterizations. To eliminate these errors, by considering NWP as an inverse problem, an unknown term in the pred...Errors inevitably exist in numerical weather prediction (NWP) due to imperfect numeric and physical parameterizations. To eliminate these errors, by considering NWP as an inverse problem, an unknown term in the prediction equations can be estimated inversely by using the past data, which are presumed to represent the imperfection of the NWP model (model error, denoted as ME). In this first paper of a two-part series, an iteration method for obtaining the MEs in past intervals is presented, and the results from testing its convergence in idealized experiments are reported. Moreover, two batches of iteration tests were applied in the global forecast system of the Global and Regional Assimilation and Prediction System (GRAPES-GFS) for July-August 2009 and January-February 2010. The datasets associated with the initial conditions and sea surface temperature (SST) were both based on NCEP (National Centers for Environmental Prediction) FNL (final) data. The results showed that 6th h forecast errors were reduced to 10% of their original value after a 20-step iteration. Then, off-line forecast error corrections were estimated linearly based on the 2-month mean MEs and compared with forecast errors. The estimated error corrections agreed well with the forecast errors, but the linear growth rate of the estimation was steeper than the forecast error. The advantage of this iteration method is that the MEs can provide the foundation for online correction. A larger proportion of the forecast errors can be expected to be canceled out by properly introducing the model error correction into GRAPES-GFS.展开更多
Model error is one of the key factors restricting the accuracy of numerical weather prediction (NWP). Considering the continuous evolution of the atmosphere, the observed data (ignoring the measurement error) can ...Model error is one of the key factors restricting the accuracy of numerical weather prediction (NWP). Considering the continuous evolution of the atmosphere, the observed data (ignoring the measurement error) can be viewed as a series of solutions of an accurate model governing the actual atmosphere. Model error is represented as an unknown term in the accurate model, thus NWP can be considered as an inverse problem to uncover the unknown error term. The inverse problem models can absorb long periods of observed data to generate model error correction procedures. They thus resolve the deficiency and faultiness of the NWP schemes employing only the initial-time data. In this study we construct two inverse problem models to estimate and extrapolate the time-varying and spatial-varying model errors in both the historical and forecast periods by using recent observations and analogue phenomena of the atmosphere. Numerical experiment on Burgers' equation has illustrated the substantial forecast improvement using inverse problem algorithms. The proposed inverse problem methods of suppressing NWP errors will be useful in future high accuracy applications of NWP.展开更多
基金the National Science Foundation(PFI-008513 and FET-2309403)for the support of this work.
文摘Optical multilayer thin film structures have been widely used in numerous photonic applications.However,existing inverse design methods have many drawbacks because they either fail to quickly adapt to different design targets,or are difficult to suit for different types of structures,e.g.,designing for different materials at each layer.These methods also cannot accommodate versatile design situations under different angles and polarizations.In addition,how to benefit practical fabrications and manufacturing has not been extensively considered yet.In this work,we introduce OptoGPT(Opto Generative Pretrained Transformer),a decoder-only transformer,to solve all these drawbacks and issues simultaneously.
基金supported by the NNSF of China(12261065)the NSF of Inner Mongolia(2022MS01005)+1 种基金the Basic Science Research Fund of the Universities Directly under the Inner Mongolia Autonomous Re-gion(JY20220084)the Program for Innovative Research Team in Universities of Inner Mongolia Autonomous Region(NMGIRT2317).
文摘In this paper,we investigate the reverse order law for Drazin inverse of three bound-ed linear operators under some commutation relations.Moreover,the Drazin invertibility of sum is also obtained for two bounded linear operators and its expression is presented.
文摘Flow units(FU)rock typing is a common technique for characterizing reservoir flow behavior,producing reliable porosity and permeability estimation even in complex geological settings.However,the lateral extrapolation of FU away from the well into the whole reservoir grid is commonly a difficult task and using the seismic data as constraints is rarely a subject of study.This paper proposes a workflow to generate numerous possible 3D volumes of flow units,porosity and permeability below the seismic resolution limit,respecting the available seismic data at larger scales.The methodology is used in the Mero Field,a Brazilian presalt carbonate reservoir located in the Santos Basin,who presents a complex and heterogenic geological setting with different sedimentological processes and diagenetic history.We generated metric flow units using the conventional core analysis and transposed to the well log data.Then,given a Markov chain Monte Carlo algorithm,the seismic data and the well log statistics,we simulated acoustic impedance,decametric flow units(DFU),metric flow units(MFU),porosity and permeability volumes in the metric scale.The aim is to estimate a minimum amount of MFU able to calculate realistic scenarios porosity and permeability scenarios,without losing the seismic lateral control.In other words,every porosity and permeability volume simulated produces a synthetic seismic that match the real seismic of the area,even in the metric scale.The achieved 3D results represent a high-resolution fluid flow reservoir modelling considering the lateral control of the seismic during the process and can be directly incorporated in the dynamic characterization workflow.
基金supported by the Special Fund of the Institute of Geophysics, China Earthquake Administration (No. DQJB21B32)the National Key R&D Program of China (No. 2022YFF0800601)。
文摘Surface-wave inversion is a powerful tool for revealing the Earth's internal structure.However,aside from shear-wave velocity(v_(S)),other parameters can influence the inversion outcomes,yet these have not been systematically discussed.This study investigates the influence of various parameter assumptions on the results of surface-wave inversion,including the compressional and shear velocity ratio(v_(P)/v_(S)),shear-wave attenuation(Q_(S)),density(ρ),Moho interface,and sedimentary layer.We constructed synthetic models to generate dispersion data and compared the obtained results with different parameter assumptions with those of the true model.The results indicate that the v_(P)/v_(S) ratio,Q_(S),and density(ρ) have minimal effects on absolute velocity values and perturbation patterns in the inversion.Conversely,assumptions about the Moho interface and sedimentary layer significantly influenced absolute velocity values and perturbation patterns.Introducing an erroneous Mohointerface depth in the initial model of the inversion significantly affected the v_(S) model near that depth,while using a smooth initial model results in relatively minor deviations.The assumption on the sedimentary layer not only affects shallow structure results but also impacts the result at greater depths.Non-linear inversion methods outperform linear inversion methods,particularly for the assumptions of the Moho interface and sedimentary layer.Joint inversion with other data types,such as receiver functions or Rayleigh wave ellipticity,and using data from a broader period range or higher-mode surface waves,can mitigate these deviations.Furthermore,incorporating more accurate prior information can improve inversion results.
基金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.
基金financially supported by the National Natural Science Foundation of China(Project No.51874336)the National Key Technologies Research and Development Program of China during the 13th Five-Year Plan Period(Project No.2017ZX005030005)。
文摘The reaction order plays a crucial role in evaluating the response rate of acid-rock.However,the conventional two-scale model typically assumes that the reaction order is constant as one,which can lead to significant deviations from reality.To address this issue,this study proposes a novel multi-order dynamic model for acid-rock reaction by combining rotating disk experimental data with theoretical derivation.Through numerical simulations,this model allows for the investigation of the impact of acidification conditions on different orders of reaction,thereby providing valuable insights for on-site construction.The analysis reveals that higher response orders require higher optimal acid liquid flow rates,and lower optimal H+diffusion coefficients,and demonstrate no significant correlation with acid concentration.Consequently,it is recommended to increase the displacement and use high-viscosity acid for reservoirs with high calcite content,while reducing the displacement and using low-viscosity acid for reservoirs with high dolomite content.
基金supported by the National Key R&D Program of China under Grant No.2021ZD0110400.
文摘Many important problems in science and engineering require solving the so-called parametric partial differential equations(PDEs),i.e.,PDEs with different physical parameters,boundary conditions,shapes of computational domains,etc.Typical reduced order modeling techniques accelerate the solution of the parametric PDEs by projecting them onto a linear trial manifold constructed in the ofline stage.These methods often need a predefined mesh as well as a series of precomputed solution snapshots,and may struggle to balance between the efficiency and accuracy due to the limitation of the linear ansatz.Utilizing the nonlinear representation of neural networks(NNs),we propose the Meta-Auto-Decoder(MAD)to construct a nonlinear trial manifold,whose best possible performance is measured theoretically by the decoder width.Based on the meta-learning concept,the trial manifold can be learned in a mesh-free and unsupervised way during the pre-training stage.Fast adaptation to new(possibly heterogeneous)PDE parameters is enabled by searching on this trial manifold,and optionally fine-tuning the trial manifold at the same time.Extensive numerical experiments show that the MAD method exhibits a faster convergence speed without losing the accuracy than other deep learning-based methods.
基金Supported by Key Research and Development Project of Guangxi Pr ovince(No.AB21196028).
文摘Unmanned aerial vehicle transient electromagnetic(UAV-TEM)is a novel airborne exploration method that offers advantages such as low cost,simple operation,high exploration efficiency and suitability for near-surface exploration in complex terrain areas.To improve the accuracy of data interpretation in this method,the authors conducted a systematic three-dimensional(3D)forward modeling and inversion of the UAV-TEM.This study utilized the finite element method based on unstructured tetrahedral elements and employed the second-order backward Euler method for time discretization.This allowed for accurate 3D modeling and accounted for the effects of complex terrain.Based on these,the influence characteristics of flight altitudes and the sizes,burial depths,and resistivities of anomalies are compared and analyzed to explore the UAV-TEM systems’exploration capability.Lastly,four typical geoelectrical models of landslides are designed,and the inversion method based on the Gauss-Newton optimization method is used to image the landslide models and analyze the imaging effect of the UAV-TEM method on landslide geohazards.Numerical results showed that UAV-TEM could have better exploration resolution and fine imaging of nearsurface structures,providing important technical support for monitoring,early warning,and preventing landslides and other geological hazards.
文摘Methane generation in landfills and its inadequate management represent the major avoidable source of anthropogenic methane today. This paper models methane production and the potential resources expected (electrical energy production and potential carbon credits from avoided CH4 emissions) from its proper management in a municipal solid waste landfill located in Ouagadougou, Burkina Faso. The modeling was carried out using two first-order decay (FOD) models (LandGEM V3.02 and SWANA) using parameters evaluated on the basis of the characteristics of the waste admitted to the landfill and weather data for the site. At the same time, production data have been collected since 2016 in order to compare them with the model results. The results obtained from these models were compared to experimental one. For the simulation of methane production, the SWANA model showed better consistency with experimental data, with a coefficient of determination (R²) of 0.59 compared with the LandGEM model, which obtained a coefficient of 0.006. Thus, despite the low correlation values linked to the poor consistency of experimental data, the SWANA model models methane production much better than the LandGEM model. Thus, despite the low correlation values linked to the poor consistency of the experimental data, the SWANA model models methane production much better than the LandGEM V3.02 model. It was noted that the poor consistency of the experimental data justifies these low coefficients, and that they can be improved in the future thanks to ongoing in situ measurements. According to the SWANA model prediction, in 27 years of operation a biogas plant with 33% electrical efficiency using biogas from the Polesgo landfill would avoid 1,340 GgCO2e. Also, the evaluation of revenues due to electricity and carbon credit gave a total revenue derived from methane production of US$27.38 million at a cost of US$10.5/tonne CO2e.
文摘The medium for next-generation communication is considered as fiber for fast,secure communication and switching capability.Mode division and space division multiplexing provide an excellent switching capability with high data transmission rate.In this work,the authors have approached an inverse modeling technique using regression-based machine learning to design a weakly coupled few-mode fiber for facilitating mode division multiplexing.The technique is adapted to predict the accurate profile parameters for the proposed few-mode fiber to obtain the maximum number of modes.It is for a three-ring-core few-mode fiber for guiding five,ten,fifteen,and twenty modes.Three types of regression models namely ordinary least-square linear multi-output regression,k-nearest neighbors of multi-output regression,and ID3 algorithm-based decision trees for multi-output regression are used for predicting the multiple profile parameters.It is observed that the ID3-based decision tree for multioutput regression is the robust,highly-accurate machine learning model for fast modeling of FMFs.The proposed fiber claims to be an efficient candidate for the next-generation 5G and 6G backhaul networks using mode division multiplexing.
文摘The outbreak of COVID-19 in 2019 resulted in numerous infections and deaths. In order to better study the transmission of COVID-19, this article adopts an improved fractional-order SIR model. Firstly, the properties of the model are studied, including the feasible domain and bounded solutions of the system. Secondly, the stability of the system is discussed, among other things. Then, the GMMP method is introduced to obtain numerical solutions for the COVID-19 system and combined with the improved MH-NMSS-PSO parameter estimation method to fit the real data of Delhi, India from April 1, 2020 to June 30, 2020. The results show that the fitting effect is quite ideal. Finally, long-term predictions were made on the number of infections. We accurately estimate that the peak number of infections in Delhi, India, can reach around 2.1 million. This paper also compares the fitting performance of the integer-order COVID-19 model and the fractional-order COVID-19 model using the real data from Delhi. The results indicate that the fractional-order model with different orders, as we proposed, performs the best.
文摘This paper investigates an improved SIR model for COVID-19 based on the Caputo fractional derivative. Firstly, the properties of the model are studied, including the feasible domain and bounded solutions of the system. Secondly, the stability of the system is discussed, among other things. Then, the GMMP method is introduced to obtain numerical solutions for the COVID-19 system. Numerical simulations were conducted using MATLAB, and the results indicate that our model is valuable for studying virus transmission.
文摘The concept of generalized order statistics has been introduced as a unified approach to a variety of models of ordered random variables with different interpretations. In this paper, we develop methodology for constructing inference based on n selected generalized order statistics (GOS) from inverse Weibull distribution (IWD), Bayesian and non-Bayesian approaches have been used to obtain the estimators of the parameters and reliability function. We have examined Bayes estimates under various losses such as the balanced squared error (balanced SEL) and balanced LINEX loss functions are considered. We show that Bayes estimate under balanced SEL and balanced LINEX loss functions are more general, which include the symmetric and asymmetric losses as special cases. This was done under assumption of discrete-continuous mixture prior for the unknown model parameters. The parametric bootstrap method has been used to construct confidence interval for the parameters and reliability function. Progressively type-II censored and k-record values as a special case of GOS are considered. Finally a practical example using real data set was used for illustration.
基金financially was supported by Colorado School of Minessupported by the Science and Technology Ministry of China (2016ZX05033003)+1 种基金China Academy of Sciences (XDA14010204)Sinopec (G5800-15-ZS-KJB016)
文摘This paper presents a unique and formal method of quantifying the similarity or distance between sedimentary facies successions from measured sections in outcrop or drilled wells and demonstrates its first application in inverse stratigraphic modeling. A sedimentary facies succession is represented with a string of symbols, or facies codes in its natural vertical order, in which each symbol brings with it one attribute such as thickness for the facies. These strings are called attributed strings. A similarity measure is defined between the attributed strings based on a syntactic pattern-recognition technique. A dynamic programming algorithm is used to calculate the similarity. Inverse stratigraphic modeling aims to generate quantitative 3D facies models based on forward stratigraphic modeling that honors observed datasets. One of the key techniques in inverse stratigraphic modeling is how to quantify the similarity or distance between simulated and observed sedimentary facies successions at data locations in order for the forward model to condition the simulation results to the observed dataset such as measured sections or drilled wells. This quantification technique comparing sedimentary successions is demonstrated in the form of a cost function based on the defined distance in our inverse stratigraphic modeling implemented with forward modeling optimization.
基金supported by the National Natural Science Foundation of China(Nos.11605163 and 21504085)the China Academy of Engineering Physics Foundation for Development of Science and Technology(No.201580103014 and No.2015B0301063)+1 种基金the Foundation for Special Talents in China Academy of Engineering Physics(No.TP201502-3)the Sichuan Science and Technology Development Foundation for Young Scientists(No.2017Q0050)
文摘Identifying the unknown geometric and material information of a multi-shield object by analyzing the radiation signature measurements is always an important problem in national and global security. In order to identify the unknown shielding layer thicknesses of a source/shield system with gamma-ray spectra, we have developed a derivative-free inverse radiation transport model based on a differential evolution algorithm with global and local neighbourhoods(IRT-DEGL). In the present paper, the IRT-DEGL model is further extended for estimating the unknown thicknesses with random initial guesses and material mass densities of multi-shielding layers as well as their combinations. Using the detected gamma-ray spectra,the illustration of inverse studies is implemented and the main influence factors for inverse results are also analyzed.
基金Project(9901351) supported by the Hydromechanical Deep Drawing without a Draw DieProject(1057001) supported by the National Natural Science Foundation of China
文摘By using aluminum alloys,the properties of the material in sheet hydroforming were obtained based on the identification of parameters for constitutive models by inverse modeling in which the friction coefficients were also considered in 2D and 3D simulations.With consideration of identified simulation parameters by inverse modeling,some key process parameters including tool dimensions and pre-bulging on the forming processes in sheet hydroforming were investigated and optimized.Based on the optimized parameters,the sheet hydroforming process can be analyzed more accurately to improve the robust design.It proves that the results from simulation based on the identified parameters are in good agreement with those from experiments.
基金funded by the National Natural Science Foundation Science Fund for Youth (Grant No.41405095)the Key Projects in the National Science and Technology Pillar Program during the Twelfth Fiveyear Plan Period (Grant No.2012BAC22B02)the National Natural Science Foundation Science Fund for Creative Research Groups (Grant No.41221064)
文摘Errors inevitably exist in numerical weather prediction (NWP) due to imperfect numeric and physical parameterizations. To eliminate these errors, by considering NWP as an inverse problem, an unknown term in the prediction equations can be estimated inversely by using the past data, which are presumed to represent the imperfection of the NWP model (model error, denoted as ME). In this first paper of a two-part series, an iteration method for obtaining the MEs in past intervals is presented, and the results from testing its convergence in idealized experiments are reported. Moreover, two batches of iteration tests were applied in the global forecast system of the Global and Regional Assimilation and Prediction System (GRAPES-GFS) for July-August 2009 and January-February 2010. The datasets associated with the initial conditions and sea surface temperature (SST) were both based on NCEP (National Centers for Environmental Prediction) FNL (final) data. The results showed that 6th h forecast errors were reduced to 10% of their original value after a 20-step iteration. Then, off-line forecast error corrections were estimated linearly based on the 2-month mean MEs and compared with forecast errors. The estimated error corrections agreed well with the forecast errors, but the linear growth rate of the estimation was steeper than the forecast error. The advantage of this iteration method is that the MEs can provide the foundation for online correction. A larger proportion of the forecast errors can be expected to be canceled out by properly introducing the model error correction into GRAPES-GFS.
基金Project supported by the Special Scientific Research Project for Public Interest(Grant No.GYHY201206009)the Fundamental Research Funds for the Central Universities,China(Grant Nos.lzujbky-2012-13 and lzujbky-2013-11)the National Basic Research Program of China(Grant Nos.2012CB955902 and 2013CB430204)
文摘Model error is one of the key factors restricting the accuracy of numerical weather prediction (NWP). Considering the continuous evolution of the atmosphere, the observed data (ignoring the measurement error) can be viewed as a series of solutions of an accurate model governing the actual atmosphere. Model error is represented as an unknown term in the accurate model, thus NWP can be considered as an inverse problem to uncover the unknown error term. The inverse problem models can absorb long periods of observed data to generate model error correction procedures. They thus resolve the deficiency and faultiness of the NWP schemes employing only the initial-time data. In this study we construct two inverse problem models to estimate and extrapolate the time-varying and spatial-varying model errors in both the historical and forecast periods by using recent observations and analogue phenomena of the atmosphere. Numerical experiment on Burgers' equation has illustrated the substantial forecast improvement using inverse problem algorithms. The proposed inverse problem methods of suppressing NWP errors will be useful in future high accuracy applications of NWP.