In the process of solving inverse gaphysical problems, the conjugate gradient (CG) method is one of the classical and valuable main algorithms. Based on classical conjugate gradient standard algorithm, the present paper made further development and deduced the standard algorithm forms for figuring out damping(IRLS)solution under the meaning of least absolute deviation(LAD), which resulted in more generalization of the application of the CG method. This method is called generalized conjugate gradient (GCG) method.In order to render the functions of GCG method more understandable. the paper has given two examples and compared the calculations with the calculations obtained by the generally recognized singular value decompeition(SVD)algorithm for figuring out morbid problems. The results show that the GCG method is characterized by strong capacity for solving morbid problems, high precision and fast operation speed. In addition, the GCG method obviously has two more features:(1) simple algorithm and flexible programming; (2) the capacity for maintaining spare characteristics of the coefficient matrix.
[1] 何宝侃、周熙襄、钟本善:地球物理反演问题中的.优化方法,地质出版社,1980年。
[2] Fulton Koehler and M. Turban Tuner: The use of the conjugate-gdiem gorithm in the computation of predictive deovolution operators, Geophysics, Vol. 50, No.12, P.2752-2758, 1985.
[3] Enders A. Robinson: Seismic time一invariant oonwlutional model, Geophysics, Vol. 50, No.12, P.2742-2752, 1985.
[4] Jeffrey R. Thoraon and Jon F. Clearbout: Vekxity - stack end slant-stack stochastic inversion, Geophysics, Vol. 50,No.12, P.2727-2741, 1986.
[5] Cluck S. and Fabre N et al: Robust multichannel stratigraphic inversion of seismic tracts, Paper prevented at the 51st EAEG meeting-West Berlin, June, 1989.
[6] Ruben D. Martinet end Bruce E. Cornish, et al: Complex~characterization by multiperameter constrained. Prexnted at the SEG/EAEG R}rch Workshop on Reservoir Geophysics, Dallas, Texas, July 31-August 3,1988.