In this paper, the Tikhonov regularization method was used to solve the nondegenerate compact hnear operator equation, which is a well-known ill-posed problem. Apart from the usual error level, the noise data were sup...In this paper, the Tikhonov regularization method was used to solve the nondegenerate compact hnear operator equation, which is a well-known ill-posed problem. Apart from the usual error level, the noise data were supposed to satisfy some additional monotonic condition. Moreover, with the assumption that the singular values of operator have power form, the improved convergence rates of the regularized solution were worked out.展开更多
Interaction between mesoscale perturbations of sea surface temperature(SSTmeso)and wind stress(WSmeso)has great influences on the ocean upwelling system and turbulent mixing in the atmospheric boundary layer.Using dai...Interaction between mesoscale perturbations of sea surface temperature(SSTmeso)and wind stress(WSmeso)has great influences on the ocean upwelling system and turbulent mixing in the atmospheric boundary layer.Using daily Quik-SCAT wind speed data and AMSR-E SST data,SSTmeso and WSmeso fields in the western coast of South America are extracted by using a locally weighted regression method(LOESS).The spatial patterns of SSTmeso and WSmeso indicate strong mesoscale SST-wind stress coupling in the region.The coupling coefficient between SSTmeso and WSmeso is about 0.0095 N/(m^2·℃)in winter and 0.0082 N/(m^2·℃)in summer.Based on mesoscale coupling relationships,the mesoscale perturbations of wind stress divergence(Div(WSmeso))and curl(Curl(WSmeso))can be obtained from the SST gradient perturbations,which can be further used to derive wind stress vector perturbations using the Tikhonov regularization method.The computational examples are presented in the western coast of South America and the patterns of the reconstructed WS meso are highly consistent with SSTmeso,but the amplitude can be underestimated significantly.By matching the spatially averaged maximum standard deviations of reconstructed WSmeso magnitude and observations,a reasonable magnitude of WSmeso can be obtained when a rescaling factor of 2.2 is used.As current ocean models forced by prescribed wind cannot adequately capture the mesoscale wind stress response,the empirical wind stress perturbation model developed in this study can be used to take into account the feedback effects of the mesoscale wind stress-SST coupling in ocean modeling.Further applications are discussed for taking into account the feedback effects of the mesoscale coupling in largescale climate models and the uncoupled ocean models.展开更多
In this paper,we consider a Cauchy problem of the time fractional diffusion equation(TFDE)in x∈[0,L].This problem is ubiquitous in science and engineering applications.The illposedness of the Cauchy problem is explai...In this paper,we consider a Cauchy problem of the time fractional diffusion equation(TFDE)in x∈[0,L].This problem is ubiquitous in science and engineering applications.The illposedness of the Cauchy problem is explained by its solution in frequency domain.Furthermore,the problem is formulated into a minimization problem with a modified Tikhonov regularization method.The gradient of the regularization functional based on an adjoint problem is deduced and the standard conjugate gradient method is presented for solving the minimization problem.The error estimates for the regularized solutions are obtained under Hp norm priori bound assumptions.Finally,numerical examples illustrate the effectiveness of the proposed method.展开更多
Downward continuation is a key step in processing airborne geomagnetic data. However,downward continuation is a typically ill-posed problem because its computation is unstable; thus, regularization methods are needed ...Downward continuation is a key step in processing airborne geomagnetic data. However,downward continuation is a typically ill-posed problem because its computation is unstable; thus, regularization methods are needed to realize effective continuation. According to the Poisson integral plane approximate relationship between observation and continuation data, the computation formulae combined with the fast Fourier transform（FFT）algorithm are transformed to a frequency domain for accelerating the computational speed. The iterative Tikhonov regularization method and the iterative Landweber regularization method are used in this paper to overcome instability and improve the precision of the results. The availability of these two iterative regularization methods in the frequency domain is validated by simulated geomagnetic data, and the continuation results show good precision.展开更多
The inverse problem to determine the vibrating velocity from known exteriorfield measurement pressure, involves the solution of a discrete ill-posed problem. To facilitate thecomputation of a meaningful approximate so...The inverse problem to determine the vibrating velocity from known exteriorfield measurement pressure, involves the solution of a discrete ill-posed problem. To facilitate thecomputation of a meaningful approximate solution possible, the indirect boundary element method(IBEM) code for investigating vibration velocity reconstruction and Tikhonov regularization methodby means of singular value decomposition (SVD) are used. The amount of regularization is determinedby a regularization parameter. Its optimal value is given by the L-curve approach. Numerical resultsindicate the reconstructed normal surface velocity is a good approximation to the real source.展开更多
In this paper,we consider the inverse problem for identifying the source term of the time-fractional equation with a hyper-Bessel operator.First,we prove that this inverse problem is ill-posed,and give the conditional...In this paper,we consider the inverse problem for identifying the source term of the time-fractional equation with a hyper-Bessel operator.First,we prove that this inverse problem is ill-posed,and give the conditional stability.Then,we give the optimal error bound for this inverse problem.Next,we use the fractional Tikhonov regularization method and the fractional Landweber iterative regularization method to restore the stability of the ill-posed problem,and give corresponding error estimates under different regularization parameter selection rules.Finally,we verify the effectiveness of the method through numerical examples.展开更多
Tikhonov regularization（TR） method has played a very important role in the gravity data and magnetic data process. In this paper, the Tikhonov regularization method with respect to the inversion of gravity data is d...Tikhonov regularization（TR） method has played a very important role in the gravity data and magnetic data process. In this paper, the Tikhonov regularization method with respect to the inversion of gravity data is discussed. and the extrapolated TR method（EXTR） is introduced to improve the fitting error. Furthermore, the effect of the parameters in the EXTR method on the fitting error, number of iterations, and inversion results are discussed in details. The computation results using a synthetic model with the same and different densities indicated that. compared with the TR method, the EXTR method not only achieves the a priori fitting error level set by the interpreter but also increases the fitting precision, although it increases the computation time and number of iterations. And the EXTR inversion results are more compact than the TR inversion results, which are more divergent. The range of the inversion data is closer to the default range of the model parameters, and the model features and default model density distribution agree well.展开更多
Nuclear Magnetic inversion is the basis of NMR Resonance （NMR） T2 logging interpretation. The regularization parameter selection of the penalty term directly influences the NMR T2 inversion result. We implemented b...Nuclear Magnetic inversion is the basis of NMR Resonance （NMR） T2 logging interpretation. The regularization parameter selection of the penalty term directly influences the NMR T2 inversion result. We implemented both norm smoothing and curvature smoothing methods for NMR T2 inversion, and compared the inversion results with respect to the optimal regular- ization parameters （（Xopt） which were selected by the dis- crepancy principle （DP）, generalized cross-validation （GCV）, S-curve, L-curve, and the slope of L-curve methods, respectively. The numerical results indicate that the DP method can lead to an oscillating or oversmoothed solution which is caused by an inaccurately estimated noise level. The （Xopt selected by the L-curve method is occa- sionally small or large which causes an undersmoothed or oversmoothed T2 distribution. The inversion results from GCV, S-curve and the slope of L-curve methods show satisfying inversion results. The slope of the L-curve method with less computation is more suitable for NMR T2 inversion. The inverted T2 distribution from norm smoothing is better than that from curvature smoothing when the noise level is high.展开更多
Two new regularization algorithms for solving the first-kind Volterra integral equation, which describes the pressure-rate deconvolution problem in well test data interpretation, are developed in this paper. The main ...Two new regularization algorithms for solving the first-kind Volterra integral equation, which describes the pressure-rate deconvolution problem in well test data interpretation, are developed in this paper. The main features of the problem are the strong nonuniform scale of the solution and large errors (up to 15%) in the input data. In both algorithms, the solution is represented as decomposition on special basic functions, which satisfy given a priori information on solution, and this idea allow us significantly to improve the quality of approximate solution and simplify solving the minimization problem. The theoretical details of the algorithms, as well as the results of numerical experiments for proving robustness of the algorithms, are presented.展开更多
A Cauchy problem for the elliptic equation with variable coefficients is considered. This problem is severely ill-posed. Then, we need use the regularization techniques to overcome its ill-posedness and get a stable n...A Cauchy problem for the elliptic equation with variable coefficients is considered. This problem is severely ill-posed. Then, we need use the regularization techniques to overcome its ill-posedness and get a stable numerical solution. In this paper, we use a modified Tikhonov regularization method to treat it. Under the a-priori bound assumptions for the exact solution, the convergence estimates of this method are established. Numerical results show that our method works well.展开更多
文摘In this paper, the Tikhonov regularization method was used to solve the nondegenerate compact hnear operator equation, which is a well-known ill-posed problem. Apart from the usual error level, the noise data were supposed to satisfy some additional monotonic condition. Moreover, with the assumption that the singular values of operator have power form, the improved convergence rates of the regularized solution were worked out.
基金Supported by the National Key Research and Development Program of China(No.2017YFC1404102(2017YFC1404100))the National Program on Global Change and Air-sea Interaction(No.GASI-IPOVAI-06)+3 种基金the National Natural Science Foundation of China(Nos.41490644(41490640),41690122(41690120))the Chinese Academy of Sciences Strategic Priority Project(No.XDA19060102)the NSFC Shandong Joint Fund for Marine Science Research Centers(No.U1406402)the Taishan Scholarship and the Recruitment Program of Global Experts。
文摘Interaction between mesoscale perturbations of sea surface temperature(SSTmeso)and wind stress(WSmeso)has great influences on the ocean upwelling system and turbulent mixing in the atmospheric boundary layer.Using daily Quik-SCAT wind speed data and AMSR-E SST data,SSTmeso and WSmeso fields in the western coast of South America are extracted by using a locally weighted regression method(LOESS).The spatial patterns of SSTmeso and WSmeso indicate strong mesoscale SST-wind stress coupling in the region.The coupling coefficient between SSTmeso and WSmeso is about 0.0095 N/(m^2·℃)in winter and 0.0082 N/(m^2·℃)in summer.Based on mesoscale coupling relationships,the mesoscale perturbations of wind stress divergence(Div(WSmeso))and curl(Curl(WSmeso))can be obtained from the SST gradient perturbations,which can be further used to derive wind stress vector perturbations using the Tikhonov regularization method.The computational examples are presented in the western coast of South America and the patterns of the reconstructed WS meso are highly consistent with SSTmeso,but the amplitude can be underestimated significantly.By matching the spatially averaged maximum standard deviations of reconstructed WSmeso magnitude and observations,a reasonable magnitude of WSmeso can be obtained when a rescaling factor of 2.2 is used.As current ocean models forced by prescribed wind cannot adequately capture the mesoscale wind stress response,the empirical wind stress perturbation model developed in this study can be used to take into account the feedback effects of the mesoscale wind stress-SST coupling in ocean modeling.Further applications are discussed for taking into account the feedback effects of the mesoscale coupling in largescale climate models and the uncoupled ocean models.
基金Supported by the National Natural Science Foundation of China(Grant No.11471253 and No.11571311)
文摘In this paper,we consider a Cauchy problem of the time fractional diffusion equation(TFDE)in x∈[0,L].This problem is ubiquitous in science and engineering applications.The illposedness of the Cauchy problem is explained by its solution in frequency domain.Furthermore,the problem is formulated into a minimization problem with a modified Tikhonov regularization method.The gradient of the regularization functional based on an adjoint problem is deduced and the standard conjugate gradient method is presented for solving the minimization problem.The error estimates for the regularized solutions are obtained under Hp norm priori bound assumptions.Finally,numerical examples illustrate the effectiveness of the proposed method.
基金supported by the National Natural Science Foundation of China(41304022,41174026,41104047)the National 973 Foundation(61322201,2013CB733303)+1 种基金the Key laboratory Foundation of Geo-space Environment and Geodesy of the Ministry of Education(13-01-08)the Youth Innovation Foundation of High Resolution Earth Observation(GFZX04060103-5-12)
文摘Downward continuation is a key step in processing airborne geomagnetic data. However,downward continuation is a typically ill-posed problem because its computation is unstable; thus, regularization methods are needed to realize effective continuation. According to the Poisson integral plane approximate relationship between observation and continuation data, the computation formulae combined with the fast Fourier transform（FFT）algorithm are transformed to a frequency domain for accelerating the computational speed. The iterative Tikhonov regularization method and the iterative Landweber regularization method are used in this paper to overcome instability and improve the precision of the results. The availability of these two iterative regularization methods in the frequency domain is validated by simulated geomagnetic data, and the continuation results show good precision.
文摘The inverse problem to determine the vibrating velocity from known exteriorfield measurement pressure, involves the solution of a discrete ill-posed problem. To facilitate thecomputation of a meaningful approximate solution possible, the indirect boundary element method(IBEM) code for investigating vibration velocity reconstruction and Tikhonov regularization methodby means of singular value decomposition (SVD) are used. The amount of regularization is determinedby a regularization parameter. Its optimal value is given by the L-curve approach. Numerical resultsindicate the reconstructed normal surface velocity is a good approximation to the real source.
基金supported by the National Natural Science Foundation of China(11961044)the Doctor Fund of Lan Zhou University of Technologythe Natural Science Foundation of Gansu Provice(21JR7RA214)。
文摘In this paper,we consider the inverse problem for identifying the source term of the time-fractional equation with a hyper-Bessel operator.First,we prove that this inverse problem is ill-posed,and give the conditional stability.Then,we give the optimal error bound for this inverse problem.Next,we use the fractional Tikhonov regularization method and the fractional Landweber iterative regularization method to restore the stability of the ill-posed problem,and give corresponding error estimates under different regularization parameter selection rules.Finally,we verify the effectiveness of the method through numerical examples.
基金supported by the National Scientific and Technological Plan(Nos.2009BAB43B00 and 2009BAB43B01)
文摘Tikhonov regularization（TR） method has played a very important role in the gravity data and magnetic data process. In this paper, the Tikhonov regularization method with respect to the inversion of gravity data is discussed. and the extrapolated TR method（EXTR） is introduced to improve the fitting error. Furthermore, the effect of the parameters in the EXTR method on the fitting error, number of iterations, and inversion results are discussed in details. The computation results using a synthetic model with the same and different densities indicated that. compared with the TR method, the EXTR method not only achieves the a priori fitting error level set by the interpreter but also increases the fitting precision, although it increases the computation time and number of iterations. And the EXTR inversion results are more compact than the TR inversion results, which are more divergent. The range of the inversion data is closer to the default range of the model parameters, and the model features and default model density distribution agree well.
基金funded by Shell International Exploration and Production Inc.(PT45371)the National Natural Science Foundation of China-China National Petroleum Corporation Petrochemical Engineering United Fund(U1262114)the National Natural Science Foundation of China(41272163)
文摘Nuclear Magnetic inversion is the basis of NMR Resonance （NMR） T2 logging interpretation. The regularization parameter selection of the penalty term directly influences the NMR T2 inversion result. We implemented both norm smoothing and curvature smoothing methods for NMR T2 inversion, and compared the inversion results with respect to the optimal regular- ization parameters （（Xopt） which were selected by the dis- crepancy principle （DP）, generalized cross-validation （GCV）, S-curve, L-curve, and the slope of L-curve methods, respectively. The numerical results indicate that the DP method can lead to an oscillating or oversmoothed solution which is caused by an inaccurately estimated noise level. The （Xopt selected by the L-curve method is occa- sionally small or large which causes an undersmoothed or oversmoothed T2 distribution. The inversion results from GCV, S-curve and the slope of L-curve methods show satisfying inversion results. The slope of the L-curve method with less computation is more suitable for NMR T2 inversion. The inverted T2 distribution from norm smoothing is better than that from curvature smoothing when the noise level is high.
文摘Two new regularization algorithms for solving the first-kind Volterra integral equation, which describes the pressure-rate deconvolution problem in well test data interpretation, are developed in this paper. The main features of the problem are the strong nonuniform scale of the solution and large errors (up to 15%) in the input data. In both algorithms, the solution is represented as decomposition on special basic functions, which satisfy given a priori information on solution, and this idea allow us significantly to improve the quality of approximate solution and simplify solving the minimization problem. The theoretical details of the algorithms, as well as the results of numerical experiments for proving robustness of the algorithms, are presented.
文摘A Cauchy problem for the elliptic equation with variable coefficients is considered. This problem is severely ill-posed. Then, we need use the regularization techniques to overcome its ill-posedness and get a stable numerical solution. In this paper, we use a modified Tikhonov regularization method to treat it. Under the a-priori bound assumptions for the exact solution, the convergence estimates of this method are established. Numerical results show that our method works well.
