Method of tracking discontinuous seismic horizons in seismic images,申请号US201213604964-传众专利搜索

Method of tracking discontinuous seismic horizons in seismic images

申请公布号：US8811677(B2)

申请号：US201213604964

申请日期：2012.09.06

申请公布日期：2014.08.19

申请人：

Total SA

发明人：Guillon Sebastien;Zinck Guillaume;Donias Marc;Lavialle Olivier

分类号：G06K9/00

主分类号：G06K9/00

代理机构：

Patterson Thuente Pedersen, P.A.

代理人：Patterson Thuente Pedersen, P.A.

地址：Paris FR

摘要：The invention targets in particular a method for searching for a seismic horizon in a seismic image of the subsoil. The method notably comprises the designation of two points belonging to the horizon sought, the iterative search for the best solution for the equation of the seismic horizon notably by virtue of a breakdown of this solution into two components, a pseudo-continuous component and a jump component.

主权项：1. A method for searching for a seismic horizon in a seismic image of the subsoil, the method comprising: designating two points of respective coordinates x1, y1 and xN, yN in the seismic image as being two points belonging to the horizon sought, an integer number N of discrete abscissae x1, x2, . . . , xN being defined between the abscissae x1, xN of the two designated points; considering a fault position at a discrete abscissa xnα where nα is an integer between 1and N−1 and a fault amplitude Cd in a discrete domain; by taking two functions f and g0 which can be derived over the interval [x1, xN], initializing a fault amplitude in a continuous domain with the value C0=Cd−[f(xnα+1)−f(xnα)], a pseudo-continuous component {tilde over (τ)}0 (xn) and a jump component {circumflex over (τ)}0(xn) defined over the N discrete abscissae according to {tilde over (τ)}0(xn)=f(xn) and {circumflex over (τ)}0(xn)=g0(xn)+C0·H(n−nα), where H(·) designates the Heaviside function; performing a number of iterations of a sequence of computation steps for an iteration index k initialized at k=0 then incremented by units, the sequence comprising, for an index k: computing a residue rk(xn)=∇{tilde over (τ)}k(xn)+∇gk(xn)−p(xn, {tilde over (τ)}k(xn)+{circumflex over (τ)}k(xn)), where ∇ designates the gradient operator and p(xn, y) designates the tangent of an inclination estimated for a position of discrete abscissa xn and of ordinate y in the seismic image;solving a Poisson equation Δ(δ{tilde over (τ)}k)=−div(rk) to determine an update term δ{tilde over (τ)}k(xn), with conditions at the Dirichlet limits: {⁢δ⁢⁢τ~0⁡(x1)=y1-f⁡(x1)-g0⁡(x1)⁢δ⁢⁢τ~0⁡(xN)=yN-f⁡(xN)-g0⁡(xN)-C0 on the first iteration and {⁢δ⁢⁢τ~k⁡(x1)=gk-1⁡(x1)-gk⁡(x1)⁢δ⁢⁢τ~k⁡(xN)=gk-1⁡(xN)-gk⁡(xN)+Ck-1-Ck on each iteration of index k≧1, where gk is a function that can be derived over the interval [x1, xN];updating the pseudo-continuous component according to {tilde over (τ)}k+1(xn)={tilde over (τ)}k(xn)+δ{tilde over (τ)}k(xn);updating the fault amplitude in the continuous domain according to Ck+1=Cd−└{tilde over (τ)}k+1(xnα+1)−{tilde over (τ)}k+1(xnα)┘;updating the jump component according to {tilde over (τ)}k+1(xn)=gk+1(xn)+Ck+1·H(n−nα);if a final value of the iteration index k is reached, computing a function τnα,Cd(xn) representing the ordinate in the seismic image of a horizon estimated as a function of the discrete abscissa xn, from a sum of the pseudo-continuous component and of the jump component; andif the final value of the iteration index k is not reached, executing the next iteration of the sequence.