Patent ID: 11860326
Assignee: INSTITUTE OF GEOLOGY AND GEOPHYSICS, CHINESE ACADEMY OF SCIENCES
Field: Measurement (Instruments)
Classification: CPC G | IPC G

Claim 3:
4. The fault characterization method for precise navigation of deep oil and gas based on image fusion as described in claim 1 is characterized in that the dip filtering, specifically:
based on the reshaped seismic waveform data;
represent a local plane wave as:, ∂
     
      u
      lpw
     
    
    
     ∂
     x
    
   
   +
   
    
     θ
     lpw
    
    ⁢
    
     
      ∂
      
       u
       lpw
      
     
     
      ∂
      
       t
       lpw
      
     
    
   
  
  =
  0
 

wherein, ulpw(x, tlpw) represents a wave field function, x represents an offset distance or a horizontal distance, θlpw represents a local dip angle of a continuous spatiotemporal domain, and tlpw represents a time of the local plane wave;
the local plane wave solution operator is:

dres=D*ssc

wherein, dres represents a residual, ssc represents a seismic channel, and D represents a plane-wave destruction filter:, D
  =
  
   [
   
    
     
      I
     
     
      0
     
     
      0
     
     
      …
     
     
      0
     
    
    
     
      
       -
       
        P
        
         1
         ,
         2
        
       
      
     
     
      I
     
     
      0
     
     
      …
     
     
      0
     
    
    
     
      0
     
     
      
       -
       
        P
        
         2
         ,
         3
        
       
      
     
     
      I
     
     
      …
     
     
      0
     
    
    
     
      …
     
     
      …
     
     
      …
     
     
      ⋱
     
     
      …
     
    
    
     
      0
     
     
      0
     
     
      …
     
     
      
       -
       
        P
        
         
          N
          -
          1
         
         ,
         N
        
       
      
     
     
      I
     
    
   
   ]
  
 

wherein, Pi, j represents a prediction factor for number of seismic channels from the i-th to the jth channel, which is a function of the local dip angle, θlpw and I and Pi, j are both M×1-dimensional matrix;
the local dip angle θlpw is estimated by minimizing an objective function of prediction error d:

θlpw=argmin||D(θlpw)||22+ξreg2||sreg[θlpw]||22

wherein, ξreg represents a regularization proportion coefficient and sreg represents a reshaping regularization operator.