<?xml version="1.0" encoding="ISO-8859-1"?><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">
<front>
<journal-meta>
<journal-id>0122-5383</journal-id>
<journal-title><![CDATA[CT&F - Ciencia, Tecnología y Futuro]]></journal-title>
<abbrev-journal-title><![CDATA[C.T.F Cienc. Tecnol. Futuro]]></abbrev-journal-title>
<issn>0122-5383</issn>
<publisher>
<publisher-name><![CDATA[Instituto Colombiano del Petróleo (ICP) - ECOPETROL S.A.]]></publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id>S0122-53832018000200075</article-id>
<article-id pub-id-type="doi">10.29047/01225383.83</article-id>
<title-group>
<article-title xml:lang="en"><![CDATA[COMPARATIVE ANALYSIS OF 3D RTM IMPLEMENTATION STRATEGIES FOR AN EFFICIENT USE OF MEMORY IN A SINGLE GPU]]></article-title>
<article-title xml:lang="es"><![CDATA[ANÁLISIS COMPARATIVO DE LAS ESTRATEGIAS DE IMPLEMENTACIÓN RTM 3D PARA EL USO EFICIENTE DE LA MEMORIA EN UNA GPU]]></article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname><![CDATA[Salamanca]]></surname>
<given-names><![CDATA[William-A]]></given-names>
</name>
<xref ref-type="aff" rid="Aff"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname><![CDATA[Ramirez]]></surname>
<given-names><![CDATA[Ana-B]]></given-names>
</name>
<xref ref-type="aff" rid="Aff"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname><![CDATA[Vivas]]></surname>
<given-names><![CDATA[Flor-A]]></given-names>
</name>
<xref ref-type="aff" rid="Aff"/>
</contrib>
</contrib-group>
<aff id="Af1">
<institution><![CDATA[,Universidad Industrial de Santander  ]]></institution>
<addr-line><![CDATA[Bucaramanga ]]></addr-line>
<country>Colombia</country>
</aff>
<aff id="Af2">
<institution><![CDATA[,Instituto Colombiano del Petróleo  ]]></institution>
<addr-line><![CDATA[Bucaramanga ]]></addr-line>
<country>Colombia</country>
</aff>
<pub-date pub-type="pub">
<day>00</day>
<month>12</month>
<year>2018</year>
</pub-date>
<pub-date pub-type="epub">
<day>00</day>
<month>12</month>
<year>2018</year>
</pub-date>
<volume>8</volume>
<numero>2</numero>
<fpage>75</fpage>
<lpage>82</lpage>
<copyright-statement/>
<copyright-year/>
<self-uri xlink:href="http://www.scielo.org.co/scielo.php?script=sci_arttext&amp;pid=S0122-53832018000200075&amp;lng=en&amp;nrm=iso"></self-uri><self-uri xlink:href="http://www.scielo.org.co/scielo.php?script=sci_abstract&amp;pid=S0122-53832018000200075&amp;lng=en&amp;nrm=iso"></self-uri><self-uri xlink:href="http://www.scielo.org.co/scielo.php?script=sci_pdf&amp;pid=S0122-53832018000200075&amp;lng=en&amp;nrm=iso"></self-uri><abstract abstract-type="short" xml:lang="en"><p><![CDATA[ABSTRACT Reverse-Time Migration (RTM) is a two-way wave-equation based method used to generate images of the Earth's subsurface. RTM has been successfully used in seismic imaging as it allows defining complex structural areas. However, RTM is a highly computational expensive algorithm that requires the computation of both the source and the receiver wavefields for each shot. Fortunately, numerical methods that compute the wave propagation using the wave equation are highly parallelizable, so they can take leverage on GPU features. However, the main problem of a GPU-RTM implementation is memory management. To take advantage of the GPU computing capabilities, the transfers to host RAM memory storage, or more expensive hard disk storage must be avoided. We present the analysis of three different strategies to implement RTM using only the memory available on a single GPU: (1) Stored wavefield checkpointing (2) Backpropagation of source wavefield using stored boundaries, and (3) Backpropagation of source wavefield using the two last snapshots and random boundaries, showing that the large amount of memory required in the first two strategies becomes a restriction over the model size. The last method (using random boundary conditions) is shown as a suggested solution to the memory problem of using a single GPU.]]></p></abstract>
<abstract abstract-type="short" xml:lang="es"><p><![CDATA[RESUMEN La migración reversa en tiempo (RTM) es un método basado en la ecuación de onda bidireccional usado para generar imágenes del subsuelo. RTM ha sido empleado exitosamente en la exploración sísmica debido a que resuelve áreas de alta complejidad estructural. Sin embargo, RTM es un algoritmo con un alto costo computacional que requiere el cálculo del campo de la fuente y el campo de los receptores en cada disparo. Afortunadamente, los métodos numéricos que permiten la extrapolación del campo de onda son altamente paralelizables y se sacan provecho de la capacidad de cómputo de la GPU. Sin embargo, el principal problema de una implementación GPU-RTM es el manejo de memoria. Para sacar provecho de la capacidad de cómputo de la GPU, se evitaron las transferencias de memoria hacia la RAM del host u otras más costosas como las transferencias al disco duro. Se presenta el análisis de tres diferentes estrategias para implementar RTM usando únicamente la memoria disponible en una GPU: (1) Almacenar el campo en puntos de control, (2) Retropropagación del campo de la fuente almacenando las fronteras, y (3) Retropropagación del campo de la fuente usando los dos últimos snapshots y fronteras aleatorias, mostrando que la gran cantidad de memoria requerida por las dos primeras estrategias se convierte en una restricción sobre el tamaño del modelo. El último método (usando condiciones de frontera aleatoria) se presenta como la solución sugerida para el problema de memoria usando únicamente una GPU.]]></p></abstract>
<kwd-group>
<kwd lng="en"><![CDATA[Graphical Processing Units]]></kwd>
<kwd lng="en"><![CDATA[Reverse-Time Migration 3D]]></kwd>
<kwd lng="en"><![CDATA[Wavefield Computation Strategies]]></kwd>
<kwd lng="es"><![CDATA[Unidades de Procesamiento Gráfico (GPU)]]></kwd>
<kwd lng="es"><![CDATA[Migración Reversa en Tiempo 3D]]></kwd>
<kwd lng="es"><![CDATA[Estrategias de cálculo del campo de onda]]></kwd>
</kwd-group>
</article-meta>
</front><back>
<ref-list>
<ref id="B1">
<label>[1]</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Baysal]]></surname>
<given-names><![CDATA[E.]]></given-names>
</name>
<name>
<surname><![CDATA[Kosloff]]></surname>
<given-names><![CDATA[D. D.]]></given-names>
</name>
<name>
<surname><![CDATA[Sherwood]]></surname>
<given-names><![CDATA[J. W]]></given-names>
</name>
</person-group>
<article-title xml:lang=""><![CDATA[Reverse time migration]]></article-title>
<source><![CDATA[Geophysics]]></source>
<year>1983</year>
<volume>48</volume>
<numero>11</numero>
<issue>11</issue>
<page-range>1514-24</page-range></nlm-citation>
</ref>
<ref id="B2">
<label>[2]</label><nlm-citation citation-type="">
<collab>NVIDIA Corporation</collab>
<source><![CDATA[CUDA C best practices guide]]></source>
<year>2017</year>
</nlm-citation>
</ref>
<ref id="B3">
<label>[3]</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Clapp]]></surname>
<given-names><![CDATA[R. G]]></given-names>
</name>
<name>
<surname><![CDATA[Haohuan]]></surname>
<given-names><![CDATA[Fu]]></given-names>
</name>
<name>
<surname><![CDATA[Lind]]></surname>
<given-names><![CDATA[O]]></given-names>
</name>
</person-group>
<article-title xml:lang=""><![CDATA[Selecting the right hardware for reverse time migration]]></article-title>
<source><![CDATA[The Leading Edge]]></source>
<year>2010</year>
<volume>29</volume>
<numero>1</numero>
<issue>1</issue>
<page-range>48-58</page-range></nlm-citation>
</ref>
<ref id="B4">
<label>[4]</label><nlm-citation citation-type="confpro">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Dussaud]]></surname>
<given-names><![CDATA[E]]></given-names>
</name>
<name>
<surname><![CDATA[Total]]></surname>
<given-names><![CDATA[E]]></given-names>
</name>
<name>
<surname><![CDATA[Symes]]></surname>
<given-names><![CDATA[W. W]]></given-names>
</name>
<name>
<surname><![CDATA[Williamson]]></surname>
<given-names><![CDATA[P]]></given-names>
</name>
<name>
<surname><![CDATA[Lemaistre]]></surname>
<given-names><![CDATA[L]]></given-names>
</name>
<name>
<surname><![CDATA[Singer]]></surname>
<given-names><![CDATA[P]]></given-names>
</name>
<name>
<surname><![CDATA[Sa]]></surname>
<given-names><![CDATA[T]]></given-names>
</name>
</person-group>
<source><![CDATA[Computational strategies for reverse-time migration]]></source>
<year>2008</year>
<conf-name><![CDATA[ SEG Las Vegas 2008 Annual Meeting]]></conf-name>
<conf-loc> </conf-loc>
<page-range>2267-71</page-range></nlm-citation>
</ref>
<ref id="B5">
<label>[5]</label><nlm-citation citation-type="confpro">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Clapp]]></surname>
<given-names><![CDATA[R. G]]></given-names>
</name>
</person-group>
<source><![CDATA[Reverse time migration with random boundaries]]></source>
<year>2009</year>
<conf-name><![CDATA[ SEG Houston 2009 International Exposition and Annual Meeting]]></conf-name>
<conf-loc> </conf-loc>
<page-range>2809-13</page-range></nlm-citation>
</ref>
<ref id="B6">
<label>[6]</label><nlm-citation citation-type="confpro">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Pasalic]]></surname>
<given-names><![CDATA[D.]]></given-names>
</name>
<name>
<surname><![CDATA[McGarry]]></surname>
<given-names><![CDATA[R]]></given-names>
</name>
</person-group>
<source><![CDATA[Convolutional perfectly matched layer for isotropic and anisotropic acoustic wave equations]]></source>
<year>2010</year>
<conf-name><![CDATA[ SEG Annual Meeting]]></conf-name>
<conf-date>2010</conf-date>
<conf-loc> </conf-loc>
<page-range>2925-9</page-range></nlm-citation>
</ref>
<ref id="B7">
<label>[7]</label><nlm-citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Griewank]]></surname>
<given-names><![CDATA[A]]></given-names>
</name>
</person-group>
<source><![CDATA[Achieving logarithmic growth of temporal and spatial complexity in reverse automatic differentiation]]></source>
<year>2011</year>
<page-range>35-54</page-range><publisher-name><![CDATA[Optimization Methods and Software]]></publisher-name>
</nlm-citation>
</ref>
<ref id="B8">
<label>[8]</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Symes]]></surname>
<given-names><![CDATA[W. W]]></given-names>
</name>
</person-group>
<article-title xml:lang=""><![CDATA[Reverse time migration with optimal checkpointing]]></article-title>
<source><![CDATA[Geophysics]]></source>
<year>2007</year>
<volume>72</volume>
<numero>5</numero>
<issue>5</issue>
<page-range>SM213--SM221</page-range></nlm-citation>
</ref>
<ref id="B9">
<label>[9]</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Bo]]></surname>
<given-names><![CDATA[F]]></given-names>
</name>
<name>
<surname><![CDATA[Huazhong]]></surname>
<given-names><![CDATA[W]]></given-names>
</name>
</person-group>
<article-title xml:lang=""><![CDATA[Reverse time migration with source wavefield reconstruction strategy]]></article-title>
<source><![CDATA[Journal of Geophysics and Engineering]]></source>
<year>2012</year>
<volume>9</volume>
<numero>1</numero>
<issue>1</issue>
<page-range>69</page-range></nlm-citation>
</ref>
<ref id="B10">
<label>[10]</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Yang]]></surname>
<given-names><![CDATA[P.]]></given-names>
</name>
<name>
<surname><![CDATA[Brossier]]></surname>
<given-names><![CDATA[R.]]></given-names>
</name>
<name>
<surname><![CDATA[Métivier]]></surname>
<given-names><![CDATA[L.]]></given-names>
</name>
<name>
<surname><![CDATA[Virieux]]></surname>
<given-names><![CDATA[J]]></given-names>
</name>
</person-group>
<article-title xml:lang=""><![CDATA[Wavefield reconstruction in attenuating media: A checkpointing-assisted reverse-forward simulation method]]></article-title>
<source><![CDATA[Geophysics]]></source>
<year>2016</year>
<volume>81</volume>
<numero>6</numero>
<issue>6</issue>
<page-range>R349-62</page-range></nlm-citation>
</ref>
<ref id="B11">
<label>[11]</label><nlm-citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Gulunay]]></surname>
<given-names><![CDATA[N]]></given-names>
</name>
</person-group>
<source><![CDATA[FXDECON and complex wiener prediction filter]]></source>
<year>1986</year>
<page-range>279-81</page-range><publisher-name><![CDATA[SEG Technical Program Expanded Abstracts 1986]]></publisher-name>
</nlm-citation>
</ref>
<ref id="B12">
<label>[12]</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Naghizadeh]]></surname>
<given-names><![CDATA[M]]></given-names>
</name>
</person-group>
<article-title xml:lang=""><![CDATA[Seismic data interpolation and denoising in the frequency-wavenumber domain]]></article-title>
<source><![CDATA[Geophysics]]></source>
<year>2012</year>
<volume>77</volume>
<numero>2</numero>
<issue>2</issue>
<page-range>V71-80</page-range></nlm-citation>
</ref>
<ref id="B13">
<label>[13]</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Gholami]]></surname>
<given-names><![CDATA[A.]]></given-names>
</name>
<name>
<surname><![CDATA[Haghshenas]]></surname>
<given-names><![CDATA[H]]></given-names>
</name>
</person-group>
<article-title xml:lang=""><![CDATA[Curvelet-TV regularized Bregman iteration for seismic random noise attenuation]]></article-title>
<source><![CDATA[Journal of AppliedGeophysics]]></source>
<year>2014</year>
<volume>109</volume>
<page-range>233-41</page-range></nlm-citation>
</ref>
</ref-list>
</back>
</article>
