3D density inversion of gravity gradient data using the extrapolated Tikhonov regularization

We use the extrapolated Tikhonov regularization to deal with the ill-posed problem of 3D density inversion of gravity gradient data. The use of regularization parameters in the proposed method reduces the deviations between calculated and observed data. We also use the depth weighting function based on the eigenvector of gravity gradient tensor to eliminate undesired effects owing to the fast attenuation of the position function. Model data suggest that the extrapolated Tikhonov regularization in conjunction with the depth weighting function can effectively recover the 3D distribution of density anomalies. We conduct density inversion of gravity gradient data from the Australia Kauring test site and compare the inversion results with the published research results. The proposed inversion method can be used to obtain the 3D density distribution of underground anomalies.

利用点电流源和Tikhonov正则化的潜艇稳恒电场反演方法

为了评估潜艇水下腐蚀相关稳恒电场分布特性,从潜艇水下腐蚀相关稳恒电场产生机理出发,基于等效点电流源建立潜艇水下腐蚀相关稳恒电场正演模型,并利用Tikhonov正则化根据已知电场数据求解等效点电流源电流强度,对潜艇周围海水空间的腐蚀相关稳恒电场进行推算。将某型潜艇的COMSOL软件仿真结果作为模拟试验数据,对所提方法有效性进行验证。结果表明:由水深38 m电场值向水深42.5 m和水深33.5 m进行反演时,相对均方根误差、最大值相对误差、峰峰值相对误差均不超过0.06;由较近平面向较远平面进行反演时,即使推算深度达到45 m,相对均方根误差仍然在0.21以内;噪声标准差为实际电场最大值的0.1倍时,反演误差仍然小于0.1。该算法抗噪声能力强,精度较高,能较好地用于工程实践。

Application of Tikhonov regularization method to wind retrieval from scatterometer data Ⅱ: cyclone wind retrieval with consideration of rain

According to the conclusion of the simulation experiments in paper I, the Tikhonov regularization method is applied to cyclone wind retrieval with a rain-effect-considering geophysical model function （called CMF＋Rain）. The CMF＋Rain model which is based on the NASA scatterometer-2 （NSCAT2） GMF is presented to compensate for the effects of rain on cyclone wind retrieval. With the multiple solution scheme （MSS）, the noise of wind retrieval is effectively suppressed, but the influence of the background increases. It will cause a large wind direction error in ambiguity removal when the background error is large. However, this can be mitigated by the new ambiguity removal method of Tikhonov regularization as proved in the simulation experiments. A case study on an extratropical cyclone of hurricane observed with SeaWinds at 25-km resolution shows that the retrieved wind speed for areas with rain is in better agreement with that derived from the best track analysis for the GMF＋Rain model, but the wind direction obtained with the two-dimensional variational （2DVAR） ambiguity removal is incorrect. The new method of Tikhonov regularization effectively improves the performance of wind direction ambiguity removal through choosing appropriate regularization parameters and the retrieved wind speed is almost the same as that obtained from the 2DVAR.

Bathymetry inversion using the modifi ed gravitygeologic method:application of the rectangular prism model and Tikhonov regularization

Bathymetry data are usually obtained via single-beam or multibeam sounding;however,these methods exhibit low efficiency and coverage and are dependent on various parameters,including the condition of the vessel and sea state.To overcome these limitations,we propose a method for marine bathymetry inversion based on the satellite altimetry gravity anomaly data as a modification of the gravity-geologic method(GGM),which is a conventional terrain inversion method based on gravity data.In accordance with its principle,the modified method adopts a rectangular prism model for modeling the short-wavelength gravity anomaly and the Tikhonov regularization method to integrate the geophysical constraints,including the a priori water depth data and characteristics of the sea bottom relief.The a priori water depth data can be obtained based on the measurement data obtained from a ship,borehole information,etc.,and the existing bathymetry/terrain model can be considered as the initial model.Marquardt’s method is used during the inversion process,and the regularization parameter can be adaptively determined.The model test and application to the West Philippine Basin indicate the feasibility and eff ectiveness of the proposed method.The results indicate the capability of the proposed method to improve the overall accuracy of the water depth data.Then,the proposed method can be used to conduct a preliminary study of the ocean depths.Additionally,the results show that in the improved GGM,the density diff erence parameter has lost its original physical meaning,and it will not have a great impact on the inversion process.Based on the boundedness of the study area,the inversion result may exhibit a lower confi dence level near the margin than that near the center.Furthermore,the modifi ed GGM is time-and memory-intensive when compared with the conventional GGM.

Some studies on the Tikhonov regularization method with additional assumptions for noise data

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.

Tikhonov正则化方法在航空γ测量数据处理中的应用

在航空γ测量中,当地形平坦且飞行高度变化不大时,采用传统的高度修正方法可得到良好的结果。实际测量中,地形崎岖不平、飞行高度突变等常见现象会导致高度修正后的结果仍存在偏差甚至错误。基于条带模型构设崎岖地形条件下的反演方程组,应用Tikhonov正则化方法求解可得到地面放射性核素的面活度浓度。Tikhonov正则化方法的应用结果表明:随机噪声较小时,采用广义交叉验证(generalized cross validation,GCV)方法选取正则化参数得到的反演结果较好;随机噪声较大时,采用L曲线方法选取正则化参数得到的反演结果较好。与传统高度修正方法的计算结果相比该方法好,且适用于飞行高度变化较大的情形。

A modified Tikhonov regularization method for a Cauchy problem of a time fractional diffusion equation

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.

Representing surface wind stress response to mesoscale SST perturbations in western coast of South America using Tikhonov regularization method

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.

用Tikhonov正则化方法同时反演对流扩散方程的对流速度和源函数

在给定两个附加观测数据的条件下,本文基于Tikhonov正则化方法研究了对流扩散方程的对流速度和源函数的同时反演问题.鉴于原问题是一个初始值非零的对流扩散方程,本文通过将初始值转化为源项得到了一个组合源项,首先将原问题转化为一个具有齐次条件的对流扩散问题.由于所得问题是不适定的,本文进而利用Tikhonov正则化方法构建了相应的极小化目标泛函,得到了问题最优解的存在性和应满足的必要条件.最后,对终端时刻较小的特殊情形,本文证明了最优解的唯一性和稳定性.

Application of Tikhonov Regularization Technique to the Kinetic Data of an Autocatalytic Reaction: Pyrolysis of N-Eicosane

A new technique based on Tikhonov regularization, for converting time-concentration data into concentration-reaction rate data, was applied to a novel pyrolysis investigation carried out by Susu and Kunugi [1]. The reaction which involves the thermal decomposition of n-eicosane using synthesis gas for K2CO3-catalyzed shift reaction was reported to be autocatalytic. This result was confirmed by applying Tikhonov regularization to the experimental data (conversion vs. time) presented by Susu and Kunugi [1]. Due to the ill-posed nature of the problem of obtaining reaction rates from experimental data, conventional methods will lead to noise amplification of the experimental data. Hence, Tikhonov regularization is preferably employed because it is entirely independent of reaction rate model and it also manages to keep noise amplification under control, thus, leading to more reliable results. This is shown by the agreement of the kinetic parameters obtained using the resulting conversion-reaction rate profile, with the Ostwald-type process for autocatalysis suggested by Susu and Kunugi [1].

VIBRATING VELOCITY RECONSTRUC-TION USING IBEM AND TIKHONOV REGULARIZATION

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.

分数阶椭圆方程反边值问题的分数Tikhonov正则化方法

该文研究了Tricomi-Gellerstedt-Keldysh型分数阶椭圆方程的反边值问题.对于该不适定问题,建立了条件稳定性结果.基于问题的不适定性,构造了分数Tikhonov正则化方法,以恢复解对测量数据的连续依赖性.在正则化参数的先验和后验选取规则下,分别给出并证明了相应的Hölder型收敛性结果.最后,通过两个数值例子验证了分数Tikhonov正则化方法的模拟效果.数值结果表明,该方法能稳定有效地处理文中反问题.

基于重复测量数据的Tikhonov型正则化方法

针对基于重复测量数据的非线性不适定问题提出了一种Tikhonov型正则化方法:利用“信号平均”重构带有凸罚项的Tikhonov型正则化解,其中正则化参数可利用先验选择规则或序列统计偏差原则确定,并在变分源条件下分析了该方法在这两种参数选择规则下的收敛性结果.

Convergence analysis on Browder-Tikhonov regularization for second-order evolution hemivariational inequality

This paper studies the Browder-Tikhonov regularization of a second-order evolution hemivariational inequality （SOEHVI） with non-coercive operators. With duality mapping, the regularized formulations and a derived first-order evolution hemivariational inequality （FOEHVI） for the problem considered are presented. By applying the Browder-Tikhonov regularization method to the derived FOEHVI, a sequence of regularized solutions to the regularized SOEHVI is constructed, and the strong convergence of the whole sequence of regularized solutions to a solution to the problem is proved.

利用GPS垂直位移反演区域陆地水储量变化的TSVD-Tikhonov正则化方法

利用GPS垂直位移反演区域陆地水储量变化(TWSC)属于典型的病态问题,其关键是如何进行稳定求解并提高反演结果的精度和可靠性.本文引入TSVD-Tikhonov组合正则化方法对利用GPS垂直位移反演区域TWSC的病态问题进行求解,并以四川省TWSC反演为例进行分析与验证.首先,通过数值模拟对TSVD、Tikhonov和TSVD-Tikhonov正则化方法采用不同正则化参数选取策略(RMSE最小准则、GCV法和L-curve法)进行反演,结果显示基于TSVD-Tikhonov正则化反演的TWSC比单独使用TSVD或Tikhonov正则化反演结果的精度和可靠性更高,这三种正则化方法反演2005年1月至12月的TWSC差值的平均STD分别为14.97 mm、7.03 mm和5.04 mm.其次,利用中国地壳运动观测网络(CMONOC)的72个GPS测站的垂直位移数据,基于TSVD-Tikhonov正则化反演了四川省2010年12月至2021年2月的TWSC时间序列,结果表明GPS反演的TWSC与GRACE/GFO Mascon模型(JPL、CSR和GSFC)的空间分布特征及季节性变化符合较好,但其TWSC信号的振幅比GRACE/GFO Mascon模型更强.最后,采用广义三角帽方法(GTCH)融合不同类型的降水、蒸散发和径流数据,并根据水量平衡方程计算的dTWSC/dt序列(PER-dS/dt)对GPS反演的dTWSC/dt序列(GPS-dS/dt)和GRACE/GFO Mascon模型融合的dTWSC/dt序列(GRACE/GFO-dS/dt)进行验证,结果表明这三类dTWSC/dt序列的季节性变化符合较好,平滑后GPS-dS/dt和GRACE/GFO-dS/dt序列与PER-dS/dt序列的相关系数分别为0.78和0.87,但GPS相比GRACE/GFO对降水变化的响应更为敏感.本文研究证明了TSVD-Tikhonov组合正则化方法能够提高GPS垂直位移反演区域TWSC的精度和可靠性,同时也表明GPS观测数据对局部水质量负荷变化更为敏感,可作为GRACE/GFO反演区域TWSC的有益补充.

Numerical estimation of choice of the regularization parameter for NMR T2 inversion

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.

Downward continuation of airborne geomagnetic data based on two iterative regularization methods in the frequency domain

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.

基于分数阶Tikhonov正则化的激光吸收光谱燃烧场二维重建光路优化研究

为探究基于激光吸收光谱技术的燃烧场二维测量光路布置方式,实现有限投影下更精确的燃烧场二维重建,根据分数阶微积分理论,提出一种基于分数阶Tikhonov正则化的光路优化方法.将经典的整数阶Tikhonov正则化推广到分数阶模式,建立了基于分数阶Tikhonov正则化的光路设计目标函数.利用遗传算法分析(0,1)范围内不同阶数的计算结果,得到最佳光路布置方式.采用近红外波段7185.6 cm^(-1)的H_(2)O特征吸收谱线结合20条测试光路对10×10离散化网格区域进行计算,对比分析五种光路布置方式对多种分布模型的重建结果,结果表明,基于分数阶Tikhonov正则化的光路布置方式具有最佳重建效果.研究结果对有限投影条件下激光吸收光谱二维测量光路的优化设计理论研究具有重要意义,可以促进激光吸收光谱技术在复杂发动机燃烧场二维重建及燃烧效率提升方面的应用.

New Regularization Algorithms for Solving the Deconvolution Problem in Well Test Data Interpretation

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.

Efficient multiuser detector based on box-constrained dichotomous coordinate descent and regularization

The presented iterative multiuser detection technique was based on joint deregularized and box-constrained solution to quadratic optimization with iterations similar to that used in the nonstationary Tikhonov iterated algorithm.The deregularization maximized the energy of the solution,which was opposite to the Tikhonov regularization where the energy was minimized.However,combined with box-constraints,the deregularization forced the solution to be close to the binary set.It further exploited the box-constrained dichotomous coordinate descent algorithm and adapted it to the nonstationary iterative Tikhonov regularization to present an efficient detector.As a result,the worst-case and average complexity are reduced down as K2.8 and K2.5 floating point operation per second,respectively.The development improves the "efficient frontier" in multiuser detection,which is illustrated by simulation results.In addition,most operations in the detector are additions and bit-shifts.This makes the proposed technique attractive for fixed-point hardware implementation.