Hyperbolic estimation of sparse models from erratic data

We have developed a hyperbolic penalty function for image estimation. The center of a hyperbola is parabolic like that of an l(2) norm fitting. Its asymptotes are similar to l(1) norm fitting. A transition threshold must be chosen for regression equations of data fitting and another threshold for model regularization. We combined two methods: Newton's and a variant of conjugate gradient method to solve this problem in a manner we call the hyperbolic conjugate direction (HYCD) method. We tested examples of (1) velocity transform with strong noise (2) migration of aliased data, and (3) blocky interval velocity estimation. For the linear experiments we performed in this study, nonlinearity is introduced by the hyperbolic objective function, but the convexity of the sum of the hyperbolas assures the convergence of gradient methods. Because of the sufficiently reliable performance obtained on the three mainstream geophysical applications, we expect the HYCD solver method to become our default method.