The invention relates to the technical field of
geophysics,
seismic exploration, oil-gas exploration and tunnel and underground
engineering, and discloses a pre-stack generalized progressive-order Taylor expansion inversion method, which comprises the following steps of: establishing an initial model of
longitudinal wave velocity,
transverse wave velocity and density by using
logging information, and constructing a prior model through a damping least square method; constructing a first-order Taylor expansion form and a second-order Taylor expansion form of the synthetic seismic
record based on an accurate Zoeppritz equation, and introducing the first-order Taylor expansion form and the second-order Taylor expansion form as double targets into a multi-target arithmetic optimization
algorithm for joint solving; the balance of global search and local approximation is realized through an arithmetic operator, and a non-dominated solution set is maintained by means of an elite archive, so that an inversion result of a continuous progressive order from first-order approximation to second-order approximation is obtained; and finally, determining an elastic parameter inversion result according to a minimum oscillation and transverse consistency criterion. According to the method, the stability and prediction reliability of pre-stack AVA inversion can be improved.