420 likes | 574 Views
COMPARISON OF WAVE EQUATION MIGRATION METHODS WITH PHASE ENCODING. Jianhua Yu. University of Utah. CONTENTS. . Introduction . Brief Description of Wave Equation Methods: SSF, PSPI, FFD (SSF+FD) . Phase Encoding Algorithm . Numerical Results . Conclusion and future work. CONTENTS.
E N D
COMPARISON OF WAVE EQUATION MIGRATION METHODS WITH PHASE ENCODING Jianhua Yu University of Utah
CONTENTS . Introduction . Brief Description of Wave Equation Methods: SSF, PSPI, FFD (SSF+FD) . Phase Encoding Algorithm . Numerical Results . Conclusion and future work
CONTENTS . Introduction . Brief Description of Wave Equation Methods: SSF, PSPI, FFD (SSF+FD) . Phase Encoding Algorithm . Numerical Results . Conclusion and future work
INTRODUCTION Migration methods includes two classes: . Ray-based migration method Features: a. Computational efficiency b. Capability of irregular acquisition geometry and target-processing
INTRODUCTION c. Depend on ray-tracing algorithm d. Less accurate image for complex area
INTRODUCTION Wave equation methods: Solution for two-way or one-way equations Features: a. Accurate wavefield extrapolation b. High quality image for complex area c. Expensive computational cost
OBJECTIVE Compare various wave-equation methods Compare phase encoding algorithms The final purpose is to develop an efficient wave equation method for 3D migration image in complex area
CONTENTS . Introduction . Brief Description of Wave Equation Methods: SSF, PSPI, FFD (SSF+FD) . Phase Encoding Algorithm . Numerical Results . Conclusion and future work
METHODOLOGY Phase-shift based wave equation methods: a. SSF b. PSPI c. FFD (SSF+FD)
Less memory requirement Compressed data in frequency domain Why use phase-shift based method?
Phase shift Method Basic equation for wavefield extrapolation is This equation can not handle velocity lateral variation !!
SSF METHOD Basic steps include: Split velocity field: Vi=V0+dVi Wavefiled extrapolation with V0 in frequency-wavenumber domain Phase correction with dVi in frequency-space domain
PSPI METHOD Basic steps of PSPI consist of: Introduce several reference velocities: Vi Wavefiled extrapolation with each Vi in frequency-wavenumber domain Interpolating reference wavefield in frequency- space domain
SSF+FD(FFD) METHOD Basic equation for FFD migration is: SSF FD FFD = SSF+FD
CONTENTS . Introduction . Brief Description of Wave Equation Methods: SSF, PSPI, FFD (SSF+FD) . Phase Encoding Algorithm . Numerical Results . Conclusion and future work
PHASE ENCODING ALGORITHM Where Sn is source term; Dn is nth shot gather; Ss and Ds are the encoded source term and super shot gather; n is the index of shot gather
can be obtained by following ways: Generate random coefficients with valuem of 1 or -1 Generate random coefficients with range [0,2 ] Generate random coefficients with Gaussian sampling HOW TO GET ?
Phase encoding Phase encoding HOW TO PHASE ENCODE ? Strategy I: G1 Gn Gm-n Gm K<M G’k G’1 Migrating Output
HOW TO PHASE ENCODE? Strategy II: G1 G2 Gi Gm Phase Encoding Gsuper Output Migrating
CONTENTS . Introduction . Brief Description of Wave Equation Methods: SSF, PSPI, FFD (SSF+FD) . Phase Encoding Algorithm . Numerical Results . Conclusion and future work
NUMERICAL RESULTS FOR SEG/EAGE SALT MODEL DATA X (km) 0 15 0 Depth (km) 3
TEST1: Comparison of Four Migration Methods X (km) X (km) 0 0 15 15 0 Depth (km) 3 Kirchhoff SSF Without Multiples
TEST1: Comparison of Four Migration Methods X (km) X (km) 0 0 15 15 0 Depth (km) 3 PSPI SSF+FD Without Multiples
TEST1: Comparison of Four Migration Methods X (km) X (km) 0 0 15 15 0 Depth (km) 3 Kirchhoff SSF With Multiples
TEST1: Comparison of Four Migration Methods X (km) X (km) 0 0 15 15 0 Depth (km) 3 PSPI SSF+FD With Multiples
TEST2: Comparison of Three Phase Encoding Methods X (km) X (km) 0 0 15 15 0 Depth (km) 3 (a) (c) 0 15 0 a. Random distribution b. Linear encoding c. Gaussian distribution Depth (km) 3 (b)
TEST3: Comparison of SSF Migration with Phase Encoding Strategy I X (km) X (km) 0 15 0 15 0 Depth (km) 3 (a). encode=2 (c). encode=10 0 Depth (km) 3 (b). encode=4 (d). no encoding(320)
TEST3: CPU Time of SSF Migration with Phase Encoding Strategy I
TEST3: Comparison of SSF Migration with Phase Encoding Strategy II X (km) X (km) 0 15 0 15 0 Depth (km) 3 (a). nmig=100 (c). nmig=60 0 Depth (km) 3 (b). nmig=80 (d). no encoding(320)
TEST3: CPU Time of SSF Migration with Phase Encoding Strategy II
TEST4: Comparison of SSF+FD Migration with Phase Encoding Strategy I X (km) X (km) 0 15 0 15 0 Depth (km) 3 (a). encode=2 (c). encode=10 0 Depth (km) 3 (b). encode=4 (d). no encoding(320)
TEST4: CPU Time of SSF+FD Migration with Phase Encoding Strategy I
TEST4: Comparison of SSF+FD Migration with Phase Encoding Strategy II X (km) X (km) 0 15 0 15 0 Depth (km) 3 (a). nmig=100 (c). nmig=80 0 Depth (km) 3 (b). nmig=80 (d). no encoding(320)
TEST4: CPU Time of SSF+FD Migration with Phase Encoding Strategy II
TEST5: CPU Time of SSF+FD Migration with Phase Encoding Strategy II 1 node 10 nodes
TEST6: Preliminary Result of 3D SSF Migration with Salt Model(zero offset data) X (km) X (km) 0 4 0 4 0 Depth (km) 10 Velocity slice (depth=1.6 km) Migration image (depth=1.6 km)
TEST4: Preliminary Result of 3D SSF Migration with Salt Model(zero offset data) X (km) X (km) 0 4 0 4 0 Depth (km) (c). nmig=80 10 Velocity slice (depth=2 km) Migration image (depth=2 km)
CONTENTS . Introduction . Brief Description of Wave Equation Methods: SSF, PSPI, FFD (SSF+FD) . Phase Encoding Algorithm . Numerical Results . Conclusion and future work
CONCLUSION The test results show: Phase encoding can reduce computational effort of wave equation migration by factor three or more SSF+FD yields best image result.
Future Work 3D prestack migration method is being tested on SALT model data We will mainly focus on developing an efficient wave equation method with phase encoding and other type compression technique for 3D migration
ACKNOWLEDGEMENTS I greatly appreciate the sponsors of UTAM Consortium for their financial support I also thank all the people who give me suggestions and help for this work