The soybean aphid, Aphis glycines Matsumura(Hemiptera: Aphididae), is one of the greatest threats to soybean production, and both trend analysis and periodic analysis of its population dynamics are important for integ...The soybean aphid, Aphis glycines Matsumura(Hemiptera: Aphididae), is one of the greatest threats to soybean production, and both trend analysis and periodic analysis of its population dynamics are important for integrated pest management(IPM). Based on systematically investigating soybean aphid populations in the field from 2018 to 2020, this study adopted the inverse logistic model for the first time, and combined it with the classical logistic model to describe the changes in seasonal population abundance from colonization to extinction in the field. Then, the increasing and decreasing phases of the population fluctuation were divided by calculating the inflection points of the models, which exhibited distinct seasonal trends of the soybean aphid populations in each year. In addition, multifactor logistic models were then established for the first time, in which the abundance of soybean aphids in the field changed with time and relevant environmental conditions. This model enabled the prediction of instantaneous aphid abundance at a given time based on relevant meteorological data. Taken as a whole, the successful approaches implemented in this study could be used to build a theoretical framework for practical IPM strategies for controlling soybean aphids.展开更多
In order to explore the influence of welding parameters and to investigate the Al alloy (AA) nugget formation process, a comprehensive model involving electrical-thermal-mechanical and metallurgical analysis was estab...In order to explore the influence of welding parameters and to investigate the Al alloy (AA) nugget formation process, a comprehensive model involving electrical-thermal-mechanical and metallurgical analysis was established to numerically display the resistance spot welding (RSW) process within multiple fields and understand the AA-RSW physics. A multi-disciplinary finite element method (FEM) framework and a empirical sub-model were built to analyze the affecting factors on weld nugget and the underlying nature of welding physics with dynamic simulation procedure. Specifically, a counter-intuitive phenomenon of the resistance time-variation caused by the transient inverse virtual variation (TIVV) effect was highlighted and analyzed on the basis of welding current and temperature distribution simulation. The empirical model describing the TIVV phenomenon was used for modifying the dynamic resistance simulation during the AA spot welding process. The numerical and experimental results show that the proposed multi-field FEM model agrees with the measured AA welding feature, and the modified dynamic resistance model captures the physics of nugget growth and the electrical-thermal behavior under varying welding current and fluctuating heat input.展开更多
This paper presents trajectory tracking control works concerning quadrotor aerial robot with rigid cross structure. The quadrotor consists of four propellers which are two paired clockwise rotate and anticlockwise rot...This paper presents trajectory tracking control works concerning quadrotor aerial robot with rigid cross structure. The quadrotor consists of four propellers which are two paired clockwise rotate and anticlockwise rotate. A nonlinear dynamic model of the quadrotor is provided, and a controller based on the improved dynamic inverse is synthesized for the purpose of stabilization and trajectory tracking. The proposed control strategy has been tested in simulation that can balance the deviation of model inaccuracy well.展开更多
The Bayesian method of statistical analysis has been applied to the parameter identification problem. A method is presented to identify parameters of dynamic models with the Bayes estimators of measurement frequencies...The Bayesian method of statistical analysis has been applied to the parameter identification problem. A method is presented to identify parameters of dynamic models with the Bayes estimators of measurement frequencies. This is based on the solution of an inverse generalized evaluate problem. The stochastic nature of test data is considered and a normal distribution is used for the measurement frequencies. An additional feature is that the engineer's confidence in the measurement frequencies is quantified and incorporated into the identification procedure. A numerical example demonstrates the efficiency of the method.展开更多
Models of marine ecosystem dynamics play an important role in revealing the evolution mechanisms of marine ecosystems and in forecasting their future changes. Most traditional ecological dynamics models are establishe...Models of marine ecosystem dynamics play an important role in revealing the evolution mechanisms of marine ecosystems and in forecasting their future changes. Most traditional ecological dynamics models are established based on basic physical and biological laws, and have obvious dynamic characteristics and ecological significance. However, they are not flexible enough for the variability of environment conditions and ecological processes found in offshore marine areas, where it is often difficult to obtain parameters for the model, and the precision of the model is often low. In this paper, a new modeling method is introduced, which aims to establish an evolution model of marine ecosystems by coupling statistics with differential dynamics. Firstly, we outline the basic concept and method of inverse modeling of marine ecosystems. Then we set up a statistical dynamics model of marine ecosystems evolution according to annual ecological observation data from Jiaozhou Bay. This was done under the forcing conditions of sea surface temperature and surface irradiance and considering the state variables of phytoplankton, zooplankton and nutrients. This model is dynamic, makes the best of field observation data, and the average predicted precision can reach 90% or higher. A simpler model can be easily obtained through eliminating the terms with smaller contributions according to the weight coefficients of model differential items. The method proposed in this paper avoids the difficulties of obtaining and optimizing parameters, which exist in traditional research, and it provides a new path for research of marine ecological dynamics.展开更多
With the high-tech industrialization of earth observation satellite remote sensing and the implementation of digital earth strategy,the energy and natural resources have been decided to be the key research fields in C...With the high-tech industrialization of earth observation satellite remote sensing and the implementation of digital earth strategy,the energy and natural resources have been decided to be the key research fields in China.In these fields,from the model based on topology data,through simple feature data model to rule-based data model,the basic spatial analysis algorithms have been developed展开更多
The expressions of matrix construction by using the singular value decomposition (SVD) are applied to the physics parameter identification of dynamic model. Then, based upon to the characteristics of a kind of matrix ...The expressions of matrix construction by using the singular value decomposition (SVD) are applied to the physics parameter identification of dynamic model. Then, based upon to the characteristics of a kind of matrix construction method, the orders of the parameter identification model can be reduced. After reducing, the mathematics and physics correspondence relations between the subsystem and the original system are distinct. the condensation errors can be avoided. The numerical example shows the benefit of the presented methodology.展开更多
This work investigated the applicability of heterogeneous and pseudo-homogeneous models to predict the dynamic behavior of a fixed-bed catalytic reactor. Some issues concerning the dynamic behavior of the system were ...This work investigated the applicability of heterogeneous and pseudo-homogeneous models to predict the dynamic behavior of a fixed-bed catalytic reactor. Some issues concerning the dynamic behavior of the system were discussed, such as the prediction of the inverse response phenomenon. The proposed models (Het- erogeneous I and II and Pseudo-homogeneous) were able to predict with qualitative similarity the main characteristics of the dynamic behavior of a fixed-bed catalytic reactor, including the inverse response. The computational time demanded for the solution of the heterogeneous models was 10 to 50% longer than in the case of the pseudo-homogeneous model, making the use of the former suitable for applications where computational time is not the major restriction (off-line applications). On the other hand, when on-line applications are required, the simplified model (Pseudo-homogeneous model) showed to be a good alternative because this model was able to predict (qualitatively) the dynamics of the reactor using a faster and easier numerical solution.展开更多
This study addresses the parameter identification problem in a system of time-dependent quasi-linear partial differential equations(PDEs).Using the integral equation method,we prove the uniqueness of the inverse probl...This study addresses the parameter identification problem in a system of time-dependent quasi-linear partial differential equations(PDEs).Using the integral equation method,we prove the uniqueness of the inverse problem in nonlinear PDEs.Moreover,using the method of successive approximations,we develop a novel iterative algorithm to estimate sorption isotherms.The stability results of the algorithm are proven under both a priori and a posteriori stopping rules.A numerical example is given to show the efficiency and robustness of the proposed new approach.展开更多
SiO_(2)-particle reinforced silicon rubber composite(SP-RSRC)is a widely utilized material that offers shock absorption protection to various engineering structures in impact environments.This paper presents a compreh...SiO_(2)-particle reinforced silicon rubber composite(SP-RSRC)is a widely utilized material that offers shock absorption protection to various engineering structures in impact environments.This paper presents a comprehensive investigation of the mechanical behavior of SP-RSRC under various strain rates,employing a combination of experimental,theoretical,and numerical analyses.Firstly,quasi-static and dynamic compression tests were performed on SP-RSRC utilizing a universal testing machine and split Hopkinson pressure bar(SHPB)apparatus.Nonlinear stress-strain relationships of SP-RSRC were obtained for strain rates ranging from 1×10^(−3) to 3065 s^(−1).The results indicated that the composite showed evident strain rate sensitivity,along with nonlinearity.Then,a nonlinear visco-hyperelastic constitutive model was developed,consisting of a hyperelastic component utilizing the 3rd-order Ogden energy function and a viscous component employing a rate-dependent relaxation time scheme.The model accurately characterized the dynamic mechanical response of SP-RSRC,effectively mitigating the challenge of calibrating an excessive number of material parameters inherent in conventional viscoelastic models.Furthermore,the simplified rubber material(SRM)model,integrated within the LS-DYNA software,was chosen to depict the mechanical properties of SP-RSRC in numerical simulations.The parameters of the SRM model were further calibrated based on the strain-stress relationships of SP-RSRC,as predicted by the developed nonlinear visco-hyperelastic constitutive model.Finally,an inverse ballistic experiment using a single-stage air gun was conducted for SP-RSRC.Numerical simulations of SHPB experiments and the inverse ballistic experiment were then performed,and the reliability of the calibrated SRM model was verified by comparing the results of experiments and numerical simulations.This study offers a valuable reference for the utilization of SP-RSRC in the realm of impact protection.展开更多
Offshore carbon capture, utilization, and storage(OCCUS) is regarded as a crucial technology for mitigating greenhouse gas emissions.Quantitative monitoring maps of sealed carbon dioxide are necessary in a comprehensi...Offshore carbon capture, utilization, and storage(OCCUS) is regarded as a crucial technology for mitigating greenhouse gas emissions.Quantitative monitoring maps of sealed carbon dioxide are necessary in a comprehensive OCCUS project. A potential high-resolution method for the aforementioned purpose lies in the full-waveform inversion(FWI) of time-lapse seismic data. However, practical applications of FWI are severely restricted by the well-known cycle-skipping problem. A new time-lapse FWI method using cross-correlation-based dynamic time warping(CDTW) is proposed to detect changes in the subsurface property due to carbon dioxide(CO_(2)) injection and address the aforementioned issue. The proposed method, namely CDTW, which combines the advantages of cross-correlation and dynamic time warping, is employed in the automatic estimation of the discrepancy between the seismic signals simulated using the baseline/initial model and those acquired. The proposed FWI method can then back-project the estimated discrepancy to the subsurface space domain, thereby facilitating retrieval of the induced subsurface property change by taking the difference between the inverted baseline and monitor models. Numerical results on pairs of signals prove that CDTW can obtain reliable shifts under amplitude modulation and noise contamination conditions. The performance of CDTW substantially outperforms that of the conventional dynamic time warping method. The proposed time-lapse fullwaveform inversion(FWI) method is applied to the Frio-2 CO_(2) storage model. The baseline and monitor models are inverted from the corresponding time-lapse seismic data. The changes in velocity due to CO_(2) injection are reconstructed by the difference between the baseline and the monitor models.展开更多
In this paper we first present a CG-type method for inverse eigenvalue problem of constructing real and symmetric matrices M,D and K for the quadratic pencil Q(λ)=λ^(2)M+λD+K,so that Q(λ)has a prescribed subset of...In this paper we first present a CG-type method for inverse eigenvalue problem of constructing real and symmetric matrices M,D and K for the quadratic pencil Q(λ)=λ^(2)M+λD+K,so that Q(λ)has a prescribed subset of eigenvalues and eigenvectors.This method can determine the solvability of the inverse eigenvalue problem automatically.We then consider the least squares model for updating a quadratic pencil Q(λ).More precisely,we update the model coefficient matrices M,C and K so that(i)the updated model reproduces the measured data,(ii)the symmetry of the original model is preserved,and(iii)the difference between the analytical triplet(M,D,K)and the updated triplet(M_(new),D_(new),K_(new))is minimized.In this paper a computationally efficient method is provided for such model updating and numerical examples are given to illustrate the effectiveness of the proposed method.展开更多
为了在地震资料和裂缝储层特征之间建立联系,对裂缝储层采用了等效介质模型.而传统的等效介质模型未充分考虑非完全弹性介质理论和基于频变各向异性理论的双相或多相流体假设,也不能对实际裂缝储层中的地震波频散和衰减现象提供准确合...为了在地震资料和裂缝储层特征之间建立联系,对裂缝储层采用了等效介质模型.而传统的等效介质模型未充分考虑非完全弹性介质理论和基于频变各向异性理论的双相或多相流体假设,也不能对实际裂缝储层中的地震波频散和衰减现象提供准确合理的解释,并且储层参数的反演研究对裂缝储层的定性预测和定量描述举足轻重.为此,本文首先根据所提出的黏弹性Chapman-Kelvin动态等效介质模型,该模型考虑了耦合的双相流体假设、黏弹性理论、喷射流以及斑块效应,并在此基础上分析了裂缝储层参数(主要为裂缝密度、裂缝长度、孔隙度和含水饱和度)对地震波频变特征的影响.然后基于黏弹性Chapman-Kelvin模型与Schoenberg和Protazio概括的Zoeppritz方程所计算出的频变反射系数,分析了反射PP波和PS波的频变AVAZ(Amplitude Versus Angle and Azimuth)特性和PP波频变反射系数与裂缝储层参数的关系.同时考虑到发生地震频散时,反射系数和频率产生关系,构建了在角度、方位和时间域内的新型正演方程.最后,基于PP波频变反射系数对裂缝密度、裂缝长度、孔隙度和含水饱和度的变化有较好的敏感性特点,进行了两种反演方法研究:其一是基于贝叶斯理论直接反演方法,其中以对数绝对范数作为似然函数和高斯分布,即L2范数度量作为先验约束;其二是基于频变反射系数的MCMC(Markov Chain Monte Carlo)随机反演方法.通过频变AVAZ进行裂缝储层参数的反演研究,结果表明,MCMC随机反演方法在缺失先验的储层参数信息时,反演结果的不确定性较强.当存在有效且足够的先验信息时,反演结果的可靠性进一步提升.而基于贝叶斯理论直接反演方法,则显示出基于频变AVAZ来区别大尺度裂缝和微尺度裂隙的潜力.展开更多
Approximate Dynamic Inversion (ADI) is basically an approximation of exact dynamic inversionor feedback linearisation, which converts a nonlinear system to an equivalent linear structure.This method can be widely appl...Approximate Dynamic Inversion (ADI) is basically an approximation of exact dynamic inversionor feedback linearisation, which converts a nonlinear system to an equivalent linear structure.This method can be widely applied for controlling minimum phase, nonaffine-in-control systems.For applying the ADI method, a fast dynamic subsystem for deriving explicit inversion of thenonaffine equation is required. With full state feedback, ADI may be expressed in the same way asa Proportional Integral (PI) controller with only knowledge of the sign of control effectiveness andalso without any approximation. The Model Reference Adaptive Controller (MRAC) augmentedwith the PI method is an adaptive control technique where the PI parameters are updated/tunedas per the control methodology based on the MRAC-Massachusetts Institute of Technology (MIT)rule so that the plant is capable to follow the reference model. The main objective of this paperis to find the relationship between ADI and MRAC augmented with a PI controller.展开更多
A technique for inversion of a nonlinear dynamical model from the observations has been developed, and examined by using the Lorenz system. The results show that the Lorenz system can be accurately inversed by using t...A technique for inversion of a nonlinear dynamical model from the observations has been developed, and examined by using the Lorenz system. The results show that the Lorenz system can be accurately inversed by using the observation data set. This technique will be broadly used in establishing the nonlinear dynamical model from an actual data set.展开更多
基金supported by the Chinese National Special Fund for Agro-scientific Research in the Public Interest (201003025 and 201103022)the National Key Research and Development Program of China (2018YFD0201004)the Discipline Construction Project of Liaoning Academy of Agricultural Sciences, China (2019DD082612)。
文摘The soybean aphid, Aphis glycines Matsumura(Hemiptera: Aphididae), is one of the greatest threats to soybean production, and both trend analysis and periodic analysis of its population dynamics are important for integrated pest management(IPM). Based on systematically investigating soybean aphid populations in the field from 2018 to 2020, this study adopted the inverse logistic model for the first time, and combined it with the classical logistic model to describe the changes in seasonal population abundance from colonization to extinction in the field. Then, the increasing and decreasing phases of the population fluctuation were divided by calculating the inflection points of the models, which exhibited distinct seasonal trends of the soybean aphid populations in each year. In addition, multifactor logistic models were then established for the first time, in which the abundance of soybean aphids in the field changed with time and relevant environmental conditions. This model enabled the prediction of instantaneous aphid abundance at a given time based on relevant meteorological data. Taken as a whole, the successful approaches implemented in this study could be used to build a theoretical framework for practical IPM strategies for controlling soybean aphids.
基金Projects (11202125, 61175038) supported by the National Natural Science Foundation of China
文摘In order to explore the influence of welding parameters and to investigate the Al alloy (AA) nugget formation process, a comprehensive model involving electrical-thermal-mechanical and metallurgical analysis was established to numerically display the resistance spot welding (RSW) process within multiple fields and understand the AA-RSW physics. A multi-disciplinary finite element method (FEM) framework and a empirical sub-model were built to analyze the affecting factors on weld nugget and the underlying nature of welding physics with dynamic simulation procedure. Specifically, a counter-intuitive phenomenon of the resistance time-variation caused by the transient inverse virtual variation (TIVV) effect was highlighted and analyzed on the basis of welding current and temperature distribution simulation. The empirical model describing the TIVV phenomenon was used for modifying the dynamic resistance simulation during the AA spot welding process. The numerical and experimental results show that the proposed multi-field FEM model agrees with the measured AA welding feature, and the modified dynamic resistance model captures the physics of nugget growth and the electrical-thermal behavior under varying welding current and fluctuating heat input.
文摘This paper presents trajectory tracking control works concerning quadrotor aerial robot with rigid cross structure. The quadrotor consists of four propellers which are two paired clockwise rotate and anticlockwise rotate. A nonlinear dynamic model of the quadrotor is provided, and a controller based on the improved dynamic inverse is synthesized for the purpose of stabilization and trajectory tracking. The proposed control strategy has been tested in simulation that can balance the deviation of model inaccuracy well.
文摘The Bayesian method of statistical analysis has been applied to the parameter identification problem. A method is presented to identify parameters of dynamic models with the Bayes estimators of measurement frequencies. This is based on the solution of an inverse generalized evaluate problem. The stochastic nature of test data is considered and a normal distribution is used for the measurement frequencies. An additional feature is that the engineer's confidence in the measurement frequencies is quantified and incorporated into the identification procedure. A numerical example demonstrates the efficiency of the method.
基金Supported by the National Basic Research Program of China (973 Program) (No. 2010CB428703)Oceanic Science Fund for Young Scholar of SOA (Nos. 2010225, 2010118)+1 种基金Public Science and Technology Research Funds Projects of Ocean of China (Nos. 201005008, 201005009)Open Fund of MOIDAT (No. 201011)
文摘Models of marine ecosystem dynamics play an important role in revealing the evolution mechanisms of marine ecosystems and in forecasting their future changes. Most traditional ecological dynamics models are established based on basic physical and biological laws, and have obvious dynamic characteristics and ecological significance. However, they are not flexible enough for the variability of environment conditions and ecological processes found in offshore marine areas, where it is often difficult to obtain parameters for the model, and the precision of the model is often low. In this paper, a new modeling method is introduced, which aims to establish an evolution model of marine ecosystems by coupling statistics with differential dynamics. Firstly, we outline the basic concept and method of inverse modeling of marine ecosystems. Then we set up a statistical dynamics model of marine ecosystems evolution according to annual ecological observation data from Jiaozhou Bay. This was done under the forcing conditions of sea surface temperature and surface irradiance and considering the state variables of phytoplankton, zooplankton and nutrients. This model is dynamic, makes the best of field observation data, and the average predicted precision can reach 90% or higher. A simpler model can be easily obtained through eliminating the terms with smaller contributions according to the weight coefficients of model differential items. The method proposed in this paper avoids the difficulties of obtaining and optimizing parameters, which exist in traditional research, and it provides a new path for research of marine ecological dynamics.
文摘With the high-tech industrialization of earth observation satellite remote sensing and the implementation of digital earth strategy,the energy and natural resources have been decided to be the key research fields in China.In these fields,from the model based on topology data,through simple feature data model to rule-based data model,the basic spatial analysis algorithms have been developed
文摘The expressions of matrix construction by using the singular value decomposition (SVD) are applied to the physics parameter identification of dynamic model. Then, based upon to the characteristics of a kind of matrix construction method, the orders of the parameter identification model can be reduced. After reducing, the mathematics and physics correspondence relations between the subsystem and the original system are distinct. the condensation errors can be avoided. The numerical example shows the benefit of the presented methodology.
文摘This work investigated the applicability of heterogeneous and pseudo-homogeneous models to predict the dynamic behavior of a fixed-bed catalytic reactor. Some issues concerning the dynamic behavior of the system were discussed, such as the prediction of the inverse response phenomenon. The proposed models (Het- erogeneous I and II and Pseudo-homogeneous) were able to predict with qualitative similarity the main characteristics of the dynamic behavior of a fixed-bed catalytic reactor, including the inverse response. The computational time demanded for the solution of the heterogeneous models was 10 to 50% longer than in the case of the pseudo-homogeneous model, making the use of the former suitable for applications where computational time is not the major restriction (off-line applications). On the other hand, when on-line applications are required, the simplified model (Pseudo-homogeneous model) showed to be a good alternative because this model was able to predict (qualitatively) the dynamics of the reactor using a faster and easier numerical solution.
基金supported by the National Natural Science Foundation of China(No.12171036)Beijing Natural Science Foundation(No.Z210001)the NSF of China No.11971221,Guangdong NSF Major Fund No.2021ZDZX1001,the Shenzhen Sci-Tech Fund Nos.RCJC20200714114556020,JCYJ20200109115422828 and JCYJ20190809150413261,National Center for Applied Mathematics Shenzhen,and SUSTech International Center for Mathematics.
文摘This study addresses the parameter identification problem in a system of time-dependent quasi-linear partial differential equations(PDEs).Using the integral equation method,we prove the uniqueness of the inverse problem in nonlinear PDEs.Moreover,using the method of successive approximations,we develop a novel iterative algorithm to estimate sorption isotherms.The stability results of the algorithm are proven under both a priori and a posteriori stopping rules.A numerical example is given to show the efficiency and robustness of the proposed new approach.
文摘SiO_(2)-particle reinforced silicon rubber composite(SP-RSRC)is a widely utilized material that offers shock absorption protection to various engineering structures in impact environments.This paper presents a comprehensive investigation of the mechanical behavior of SP-RSRC under various strain rates,employing a combination of experimental,theoretical,and numerical analyses.Firstly,quasi-static and dynamic compression tests were performed on SP-RSRC utilizing a universal testing machine and split Hopkinson pressure bar(SHPB)apparatus.Nonlinear stress-strain relationships of SP-RSRC were obtained for strain rates ranging from 1×10^(−3) to 3065 s^(−1).The results indicated that the composite showed evident strain rate sensitivity,along with nonlinearity.Then,a nonlinear visco-hyperelastic constitutive model was developed,consisting of a hyperelastic component utilizing the 3rd-order Ogden energy function and a viscous component employing a rate-dependent relaxation time scheme.The model accurately characterized the dynamic mechanical response of SP-RSRC,effectively mitigating the challenge of calibrating an excessive number of material parameters inherent in conventional viscoelastic models.Furthermore,the simplified rubber material(SRM)model,integrated within the LS-DYNA software,was chosen to depict the mechanical properties of SP-RSRC in numerical simulations.The parameters of the SRM model were further calibrated based on the strain-stress relationships of SP-RSRC,as predicted by the developed nonlinear visco-hyperelastic constitutive model.Finally,an inverse ballistic experiment using a single-stage air gun was conducted for SP-RSRC.Numerical simulations of SHPB experiments and the inverse ballistic experiment were then performed,and the reliability of the calibrated SRM model was verified by comparing the results of experiments and numerical simulations.This study offers a valuable reference for the utilization of SP-RSRC in the realm of impact protection.
文摘Offshore carbon capture, utilization, and storage(OCCUS) is regarded as a crucial technology for mitigating greenhouse gas emissions.Quantitative monitoring maps of sealed carbon dioxide are necessary in a comprehensive OCCUS project. A potential high-resolution method for the aforementioned purpose lies in the full-waveform inversion(FWI) of time-lapse seismic data. However, practical applications of FWI are severely restricted by the well-known cycle-skipping problem. A new time-lapse FWI method using cross-correlation-based dynamic time warping(CDTW) is proposed to detect changes in the subsurface property due to carbon dioxide(CO_(2)) injection and address the aforementioned issue. The proposed method, namely CDTW, which combines the advantages of cross-correlation and dynamic time warping, is employed in the automatic estimation of the discrepancy between the seismic signals simulated using the baseline/initial model and those acquired. The proposed FWI method can then back-project the estimated discrepancy to the subsurface space domain, thereby facilitating retrieval of the induced subsurface property change by taking the difference between the inverted baseline and monitor models. Numerical results on pairs of signals prove that CDTW can obtain reliable shifts under amplitude modulation and noise contamination conditions. The performance of CDTW substantially outperforms that of the conventional dynamic time warping method. The proposed time-lapse fullwaveform inversion(FWI) method is applied to the Frio-2 CO_(2) storage model. The baseline and monitor models are inverted from the corresponding time-lapse seismic data. The changes in velocity due to CO_(2) injection are reconstructed by the difference between the baseline and the monitor models.
基金Research supported by National Natural Science Foundation of China(10571047 and 10861005)Provincial Natural Science Foundation of Guangxi(0991238)。
文摘In this paper we first present a CG-type method for inverse eigenvalue problem of constructing real and symmetric matrices M,D and K for the quadratic pencil Q(λ)=λ^(2)M+λD+K,so that Q(λ)has a prescribed subset of eigenvalues and eigenvectors.This method can determine the solvability of the inverse eigenvalue problem automatically.We then consider the least squares model for updating a quadratic pencil Q(λ).More precisely,we update the model coefficient matrices M,C and K so that(i)the updated model reproduces the measured data,(ii)the symmetry of the original model is preserved,and(iii)the difference between the analytical triplet(M,D,K)and the updated triplet(M_(new),D_(new),K_(new))is minimized.In this paper a computationally efficient method is provided for such model updating and numerical examples are given to illustrate the effectiveness of the proposed method.
文摘为了在地震资料和裂缝储层特征之间建立联系,对裂缝储层采用了等效介质模型.而传统的等效介质模型未充分考虑非完全弹性介质理论和基于频变各向异性理论的双相或多相流体假设,也不能对实际裂缝储层中的地震波频散和衰减现象提供准确合理的解释,并且储层参数的反演研究对裂缝储层的定性预测和定量描述举足轻重.为此,本文首先根据所提出的黏弹性Chapman-Kelvin动态等效介质模型,该模型考虑了耦合的双相流体假设、黏弹性理论、喷射流以及斑块效应,并在此基础上分析了裂缝储层参数(主要为裂缝密度、裂缝长度、孔隙度和含水饱和度)对地震波频变特征的影响.然后基于黏弹性Chapman-Kelvin模型与Schoenberg和Protazio概括的Zoeppritz方程所计算出的频变反射系数,分析了反射PP波和PS波的频变AVAZ(Amplitude Versus Angle and Azimuth)特性和PP波频变反射系数与裂缝储层参数的关系.同时考虑到发生地震频散时,反射系数和频率产生关系,构建了在角度、方位和时间域内的新型正演方程.最后,基于PP波频变反射系数对裂缝密度、裂缝长度、孔隙度和含水饱和度的变化有较好的敏感性特点,进行了两种反演方法研究:其一是基于贝叶斯理论直接反演方法,其中以对数绝对范数作为似然函数和高斯分布,即L2范数度量作为先验约束;其二是基于频变反射系数的MCMC(Markov Chain Monte Carlo)随机反演方法.通过频变AVAZ进行裂缝储层参数的反演研究,结果表明,MCMC随机反演方法在缺失先验的储层参数信息时,反演结果的不确定性较强.当存在有效且足够的先验信息时,反演结果的可靠性进一步提升.而基于贝叶斯理论直接反演方法,则显示出基于频变AVAZ来区别大尺度裂缝和微尺度裂隙的潜力.
文摘Approximate Dynamic Inversion (ADI) is basically an approximation of exact dynamic inversionor feedback linearisation, which converts a nonlinear system to an equivalent linear structure.This method can be widely applied for controlling minimum phase, nonaffine-in-control systems.For applying the ADI method, a fast dynamic subsystem for deriving explicit inversion of thenonaffine equation is required. With full state feedback, ADI may be expressed in the same way asa Proportional Integral (PI) controller with only knowledge of the sign of control effectiveness andalso without any approximation. The Model Reference Adaptive Controller (MRAC) augmentedwith the PI method is an adaptive control technique where the PI parameters are updated/tunedas per the control methodology based on the MRAC-Massachusetts Institute of Technology (MIT)rule so that the plant is capable to follow the reference model. The main objective of this paperis to find the relationship between ADI and MRAC augmented with a PI controller.
基金Project supported by the National Natural Science Foundation of China.
文摘A technique for inversion of a nonlinear dynamical model from the observations has been developed, and examined by using the Lorenz system. The results show that the Lorenz system can be accurately inversed by using the observation data set. This technique will be broadly used in establishing the nonlinear dynamical model from an actual data set.