摘要 |
Method for processing seismic data includes: dividing an inversion into at least three frequency regions, region 1, region 2, and region 3; computing, via a computing processor, an inversion of source separation in region 2 using M=(GHG)-1GH D= VS-1 UHD, wherein, M is reflectivity of the data, G is the vibrator- sweep matrix, D is the input data matrix, U is a matrix of eigenvectors that span the data space, V is a matrix of eigenvectors that span the model space, S is a diagonal eigenvalue matrix whose diagonal elements are called singular values, and H is a conjugate transpose operator; selecting a stable solution within region 2 at a reference frequency (fRL), wherein fRL = f3 + δf, wherein δf is between 0 to 10 Hz, and wherein the solution at the frequency of fRL is chosen as a referenced solution; replacing singular values, S, within low-frequency range (f1 to f2) with S'=S+δS, where δs represents a fraction of singular value of the referenced solution and the fraction typically ranges between 1 to 30 percent; computing a generalized least-squares solution for the low-frequency range (f1 to f2) as M = V S'-1 UHD; selecting a stable solution within region 2 as a reference frequency (fRH), wherein fRH = f4 - δf, wherein δf is between 0 to 10 Hz and wherein the solution at the frequency of fRH is chosen as a referenced solution; replacing singular values, S, within high-frequency range (f5 to f6) with S'=S+δS, where δS represents a fraction of singular value of the referenced solution and the fraction typically ranges between 1 to 30 percent; and computing a generalized least-squares solution as, M = VS'-1 UHD. |