The Kaleidoscope Project: Breaking the Sound Barrier in ... - prace
The Kaleidoscope Project: Breaking the Sound Barrier in ... - prace
The Kaleidoscope Project: Breaking the Sound Barrier in ... - prace
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
<strong>Break<strong>in</strong>g</strong> <strong>the</strong> <strong>Sound</strong> <strong>Barrier</strong> <strong>in</strong><br />
Seismic Imag<strong>in</strong>g:<br />
<strong>The</strong> <strong>Kaleidoscope</strong><br />
<strong>Project</strong><br />
PRACE Amsterdam September 2008
ContentsContents<br />
Oil Exploration <strong>in</strong> <strong>the</strong> Gulf of Mexico<br />
Seismic Methods<br />
<strong>The</strong> <strong>Kaleidoscope</strong> <strong>Project</strong><br />
Conclud<strong>in</strong>g Remarks
Oil Exploration <strong>in</strong> <strong>the</strong> Gulf of Mexico.<br />
Cartoon from Greg Guidri, Shell<br />
From Greg Guidri, Shell
<strong>The</strong> Future of Oil Exploration <strong>in</strong> GoM<br />
Cartoon from Greg Guidri, Shell
Boundary conditions <strong>in</strong> GOM<br />
•A competitive environment like GOM where<br />
Significant rema<strong>in</strong><strong>in</strong>g reserves<br />
Well cost over 150 MMUS$<br />
Data is multiclient<br />
requires for an Operator to develop a Strong Competitive<br />
Advantage based <strong>in</strong> Propietary State-of-<strong>the</strong>-Art Technology.<br />
•That technology must focus on Imag<strong>in</strong>g and Drill<strong>in</strong>g.<br />
Stena Drillmax <strong>Project</strong> a new concept <strong>in</strong> deepwater<br />
drill<strong>in</strong>g<br />
<strong>The</strong> <strong>Kaleidoscope</strong> <strong>Project</strong> seek<strong>in</strong>g Imag<strong>in</strong>g Quality and<br />
Fast Turnaround.
<strong>The</strong> Seismic Method: Introduction<br />
• We can f<strong>in</strong>d directly some natural resources (i.e.<br />
wood) but Hydrocarbons (HC) are produced from<br />
reservoirs under <strong>the</strong> earth<br />
• Opportunities for direct observation of subsurface<br />
rocks are limited to outcrops, m<strong>in</strong>es, tunnels and<br />
caves, which toge<strong>the</strong>r offer only a glimpse of <strong>the</strong><br />
<strong>in</strong>ner Earth<br />
• After geologists identify a potential zone for HC<br />
presence, we need <strong>in</strong>direct methods to explore<br />
subsurface. Geophysical groups of methods used<br />
by <strong>the</strong> <strong>in</strong>dustry:<br />
Gravimetry<br />
Electromagnetic<br />
Seismic
Basic Oil Geology<br />
Some geologic structures that form traps:<br />
SALT DOME<br />
Anticl<strong>in</strong>e<br />
Fault trap<br />
GOM
<strong>The</strong> Seismic Method: Introduction
<strong>The</strong> Seismic Method: Introduction<br />
Mar<strong>in</strong>e adquisition
<strong>The</strong> Seismic Method: Initial data
<strong>The</strong> Seismic Method: F<strong>in</strong>al results
<strong>The</strong> Seismic Method <strong>in</strong> HPC Context<br />
• Seismic Imag<strong>in</strong>g is a Grand Challenge for HPC<br />
Solves an Inverse Problem: Observed Data Creates a Model:<br />
no unique solution and iterative process and depends on<br />
human <strong>in</strong>teraction for <strong>in</strong>termediate evaluations/iterations<br />
Algorithms Based <strong>in</strong> partial differential equations. Algorithms<br />
have evolutioned <strong>in</strong> time accord<strong>in</strong>g to availability of HW/SF<br />
We can asume some historical algorithms as shortcuts<br />
accord<strong>in</strong>g to technology avaible at that time<br />
Applications designed for run <strong>in</strong> cluster LINUX and tailored for a<br />
particular processor<br />
Academic world <strong>in</strong>terested <strong>in</strong> evolution and new solutions<br />
• Bad news: Different E&P projects require different algorithms.<br />
No magical recipes.
Subsalt Plays <strong>in</strong> GOM<br />
Cartoon from Greg Guidry, Shell
Alternatives for Imag<strong>in</strong>g <strong>in</strong> GOM<br />
• Alternatives for Subsalt Imag<strong>in</strong>g<br />
– One Way Wave Equation: dip limitation<br />
– Kirchhoff PSDM: problems with complex<br />
structures where multi-pah<strong>in</strong>g occurs<br />
– RTM: Not known <strong>the</strong>oretical limitation<br />
• <strong>The</strong> limitation of RTM: Operational, it is a<br />
PETAFLOPS application<br />
– Need for at least <strong>in</strong>crease <strong>in</strong> one Order of<br />
Magnitude <strong>in</strong> Comput<strong>in</strong>g Power
Comparative Generic Algorithm<br />
2000<br />
2003<br />
2006<br />
2008-9<br />
Kirchhoff<br />
Migration<br />
Common<br />
Azimuth<br />
Migration<br />
Shot<br />
Profile<br />
Migration<br />
Reverse Time<br />
Migration<br />
Kernel<br />
Rays<br />
1-way<br />
1-way<br />
2-way<br />
Dip limitation<br />
None<br />
< 90 o<br />
< 90 o<br />
None<br />
Disk per shot<br />
(relative)<br />
1<br />
0.1<br />
1<br />
> 100<br />
Computation time<br />
per shot (relative)<br />
0.1<br />
0.1<br />
1<br />
10<br />
Better Image
Subsalt Plays and RTM <strong>in</strong> GoM<br />
Industry idea: <strong>The</strong> future of Subsalt Imag<strong>in</strong>g is RTM
<strong>Kaleidoscope</strong> <strong>Project</strong>: <strong>The</strong> vision<br />
• Start to work with <strong>the</strong> idea of no computer<br />
limitation: To have a complete RTM when <strong>the</strong><br />
petascale threshold were really commercially<br />
avaible<br />
• To use <strong>the</strong> PSDM sources of value schema for<br />
search colaborators and def<strong>in</strong>e <strong>the</strong> project<br />
methodology
<strong>Kaleidoscope</strong>: Collaboration Model<br />
Stanford Exploration <strong>Project</strong><br />
3DGeo<br />
Collaboration<br />
to promote<br />
R+D<br />
IBM BSC CSIC<br />
Repsol
<strong>Kaleidoscope</strong>: <strong>Project</strong> Methodology<br />
Model<br />
Algorithm<br />
Capacity
<strong>Kaleidoscope</strong>: Capacity<br />
Model<br />
Algorithm<br />
Capacity
Peta-scale capacity: Our vision<br />
Widely available peta-scale capacity will require:<br />
• Multi-element (multicore) processors, given that CPU<br />
frequency has reached a limit. Any multi-element processor<br />
will be difficult to program.<br />
• This multi-element processor must have very low power<br />
consumption to make peta-scale capacity economically<br />
viable and technically feasible.<br />
• Widely available peta-scale capacity requires a very cheap<br />
processor. That implies mass production <strong>in</strong> figures larger<br />
than <strong>the</strong> current production for processors <strong>in</strong> <strong>the</strong> PC<br />
<strong>in</strong>dustry.
<strong>The</strong> Cell/B.E. Processor<br />
• Co-developed by Sony, IBM and<br />
Toshiba, it was presented <strong>in</strong><br />
February 2005.<br />
• Every processor is <strong>in</strong>itially composed<br />
by 9 subunits that work<br />
<strong>in</strong>dependently.<br />
• It differs completely from <strong>the</strong><br />
philosophy of multicore processors.<br />
• Currently marketed <strong>in</strong> <strong>the</strong> Sony<br />
PlayStation 3 and <strong>in</strong> <strong>the</strong> IBM QS2x<br />
Blades.<br />
Courtesy IBM
<strong>The</strong> Mare Nostrum Supercomputer<br />
Not <strong>the</strong> most powerful world supercomputer but <strong>the</strong> most beaufitul <strong>in</strong>deed
<strong>Kaleidoscope</strong>: Algorithm<br />
Model<br />
Algorithm<br />
Capacity
Next Generation of Imag<strong>in</strong>g Tools<br />
Anisotropic One-pass and Two-pass Shot Profile<br />
Migration<br />
Plane-wave Migration & Migration <strong>in</strong> Tilted<br />
Coord<strong>in</strong>ates (PWSPM)<br />
Anisotropic Reverse Time Migration (RTM)<br />
Automated Wavepath Velocity Model Updat<strong>in</strong>g<br />
Waveform Inversion of Seismic Data<br />
IBM Power Optimized software<br />
Cell/B.E. Optimized software
<strong>Kaleidoscope</strong>: Model<br />
Model<br />
Algorithm<br />
Capacity
Depth and Time Migration<br />
Real Case: IMPORTANCE OF VELOCITY MODEL<br />
Time Process<strong>in</strong>g<br />
Depth Process<strong>in</strong>g
Model<br />
• Repsol YPF has <strong>the</strong> practical experience needed<br />
by <strong>the</strong> project: RESULTS<br />
• Repsol YPF provides experience <strong>in</strong> Model<br />
Build<strong>in</strong>g: “trench’s process<strong>in</strong>g”
Conclud<strong>in</strong>gs remarks<br />
• <strong>The</strong> <strong>Kaleidoscope</strong> <strong>Project</strong> encompasses a simultaneous<br />
research <strong>in</strong> Hardware and Software to achieve a Petascale<br />
solution to Seismic Imag<strong>in</strong>g.<br />
• Software Research focuses on <strong>the</strong> Quality of Algorithms<br />
avoid<strong>in</strong>g shorcuts or tradeoffs to ensure <strong>the</strong> maximum<br />
possible Image quality regardless of <strong>the</strong> comput<strong>in</strong>g power<br />
required.<br />
• In parallel, Hardware research focuses on <strong>the</strong><br />
development of peta-scale supercomputer from Commodity<br />
Off-<strong>the</strong>-Shelf technology aimed to massive markets.<br />
• Speed and Power are guaranteed while ensur<strong>in</strong>g low cost,<br />
future improvement, and cont<strong>in</strong>uity on <strong>the</strong> solution.<br />
• <strong>The</strong> <strong>Kaleidoscope</strong> <strong>Project</strong> has anticipated this trend and<br />
is creat<strong>in</strong>g <strong>the</strong> new standards for this technological leap<br />
and new methodology <strong>in</strong> R+D <strong>in</strong>dustry/academia projects
Thank You !<br />
Visit us at www.repsolypf.com