International Journal of Geophysics

Volume 2018 (2018), Article ID 2913591, 8 pages

https://doi.org/10.1155/2018/2913591

## SEISGAMA: A Free C# Based Seismic Data Processing Software Platform

Correspondence should be addressed to Wiwit Suryanto; di.ca.mgu@sw

Received 4 July 2017; Accepted 24 December 2017; Published 30 January 2018

Academic Editor: Filippos Vallianatos

Copyright © 2018 Theodosius Marwan Irnaka et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Seismic reflection is one of the most popular methods in geophysical prospecting. Nevertheless, obtaining high resolution and accurate results requires a sophisticated processing stage. There are many open-source seismic reflection data processing software programs available; however, they often use a high-level programming language that decreases its overall performance, lacks intuitive user-interfaces, and is limited to a small set of tasks. These shortcomings reveal the need to develop new software using a programming language that is natively supported by Windows® operating systems, which uses a relatively medium-level programming language (such as C#) and can be enhanced by an intuitive user interface. SEISGAMA was designed to address this need and employs a modular concept, where each processing group is combined into one module to ensure continuous and easy development and documentation. SEISGAMA can perform basic seismic reflection processes. This ability is very useful, especially for educational purposes or during a quality control process (in the acquisition stage). Those processes can be easily carried out by users via specific menus on SEISGAMA’s main user interface. SEISGAMA has been tested, and its results have been verified using available theoretical frameworks and by comparison to similar commercial software.

#### 1. Introduction

Geophysics is a branch of physics that studies physical properties of the earth, such as mass density, seismic wave speed, electrical resistance, and magnetic properties [1–3] from surface measurements. These physical properties can be further used to image subsurface geological conditions of the earth. One of the most popular methods in geophysics is the seismic reflection method. This method is very sensitive to a rock’s acoustic impedance variation, which is characterized by seismic waves reflected in each interface of the subsurface layer or structure. Acoustic impedance is denoted by the product of rock density and seismic wave speed within a certain medium. Seismic reflection is very popular, especially in the oil and gas industry, because this method yields high resolution subsurface images of geological structures [4]. The seismic reflection method is also useful in engineering, environmental, geohazard, and ground-water exploration, as discussed by Steeples and Miller [5].

Recently, a large number of open-source software programs have emerged that can be used to process seismic reflection data; however, each program has its own disadvantages. For example, Seismic Unx [6] has a really good ability to process seismic reflection data, but it has no intuitive user interface and depends on a command-line interface. Templeton and Gough [7] developed a web-based seismic data processing tool based on the Seismic Unx platform. Through this web service, users can fill-out the form provided by the server, and the input data will be queued on a parallel processing scheduler application. In this sense, the online service is a good concept, but the performance is limited to user interaction via a form-based input with limited resources on the server side. DSISoft, which was developed by Beaty et al. [8], is a vertical seismic profiling data processor based on the MATLAB® programming language. MATLAB is a popular as programming language because it has a library of common mathematical and engineering functions. Unfortunately, it has a large computational cost as data size increases and is thus inefficient compared to a lower level programming language (e.g., C, Fortran, and C#). A comparison between MATLAB and C, made by Andrews [9], shows significant speed differences; He determines that a C program is up to 546 times faster than a MATLAB program. A more comprehensive seismic reflection data processing implementation in MATLAB was developed by Mousa and Al-Shuhail, [10]. Another example is SEISPRHO [11], software developed in Lazarus Pascal Language, which can be used to interpret seismic reflection and is enhanced with simple processing tools. This software is mainly developed for interpretation purposes.

The recent popularity of Python, another high-level programming language that is free and reliable, has many features to handle numerical problems and has high quality graphics. It is increasingly being used to develop geophysical data processing software [12, 13]. One of the advantages of Python is its ability to be installed on a server so that the application can be accessed remotely (e.g., [14]). Schwehr [15] has developed an application called* seismic-py*, which provides an infrastructure for creating and managing a Python library for each seismic data format. Rizvandi et al. [16] detailed the performance of an implementation of the Prestack Kirchhoff Time Migration in Python. There is currently no complete application that is built entirely with Python (i.e., ready to use without combining the Python library to another library), as this would require specialized Python programming skills, at least the basic ability to program in Python.

To the best of our knowledge, there is no other free seismic reflection data processing software program built with the C# programming language that has a complete seismic processing toolbox, can process raw data from field recordings, and includes an intuitive and appealing user interface. The C# programming language was considered due to its native compatibility with the most popular operating system, that is, Windows. Netmarketshare states that approximately 81% of desktop users have a Windows XP or newer operating system installed on their computers [17]. The C# programming language is officially supported and periodically updated by Windows and included on a .NET (dot NET) platform [18, 19]. Developing applications natively supported on Windows operating systems ensures long-term and continuous software development. Enhanced with an intuitive user interface, SEISGAMA can easily be used by beginner or advanced users, and it is particularly useful for undergraduate geophysics students as they are learning processes. The need for reliable and free seismic processing software has long been awaited, especially since the price of oil became variable. This software may ensure that oil exploration can be resumed with lower associated costs.

#### 2. Methods

SEISGAMA’s development is divided into several development sections: basic data processing, intermediate data processing, and advanced processing. This paper reports only the basic processing aspects of reflection seismic methods, and the advanced processing aspects will be discussed separately in another paper. The basic data processor that was developed in this research consists of amplitude correction, muting, - and - domain transform, velocity analysis, normal moveout (NMO) correction, and traces stacking. These processing steps combine into one trace with a high signal to noise ratio (SNR).

A more detailed explanation of the processing steps follows.

##### 2.1. Amplitude Correction

Amplitude correction was performed to compensate seismic signal weakening, which decreases exponentially as a function of time due to attenuation and geometrical spreading. Campillo et al. [20] explains that this phenomenon often holds in some regional seismic phases. The phenomenon is also described by geometrical spreading [21]. Amplitude attenuation may be exploited via an exponential or oscillatory equation; therefore, the amplitude correction can be formulated as an exponential equation (1). is seismic data after amplitude correction, is the data before amplitude correction, and is an exponential constant. This equation is an example of amplitude correction for a data-independent method. After amplitude correction, the signal has amplitude with no attenuation as the recording time is increased.

##### 2.2. Muting

The muting process removes unused parts of the seismic recording. This process is very important because in the seismic reflection method, only the reflection signal will be considered, and other signals, such as those from direct waves, ground-roll, or refraction, must be removed. The muting process is carried out by defining the area to be removed first, and then, SEISGAMA will automatically make the amplitude of that region equal to zero.

##### 2.3. - and - Transformation

Seismic data is typically represented in time-domain versus (range or) distance (-). The - transformation is carried out by using a 1D Fourier transform and is applied to the transformation for each seismic trace in the data. This will transform the time series data into frequency domain data. The results are then displayed in terms of distance (-axis) and frequency (-axis). The Fourier transformation is expressed byThe - transform utilizes the 2D Fourier transformation to transform a time-domain version of each seismic trace into the frequency domain in and transform the spatial domain of data into wave number (). In this transformation, seismic data in time versus distance will be transformed into frequency versus wave number (). The 2D Fourier transform is an extension of the 1D Fourier transformation and expressed by

##### 2.4. Velocity Analysis

The reflection from a horizontal reflector of a seismic wave on the record for a particular seismic source will follow a parabolic equationThese characteristics are used to estimate the wave propagation speed by assuming the velocity of the medium. Parabolic response patterns can be detected using a shot-gather-based seismic semblance calculation , as noted by Mousa and Al-Shuhail [10].

A seismic semblance is formulated as the quotient between stacked and prestacked energy subsequently multiplied by the number of seismic traces (7) for each time window (). Time window parameters, with width , are required to perform an calculation. Here, is the number of samples at the time window, ; is the number of the traces on seismic trace collections; is the amplitude at time and trace ; is stacked energy; and is prestacked energy.

##### 2.5. Normal Moveout Correction

A normal moveout (NMO) correction is carried out to restore the nonzero offset seismic reflectors to zero offset. The correction is done after we obtain and for each reflector during velocity analysis process. The time delay after NMO correction is defined in (9). A NMO correction is suitable for a reflector having a slope of less than 15°. When a reflector has higher slope than that threshold, another correction (Dip Move Out correction) is required.NMO corrections can be categorized by three classifications: if the NMO velocity () is too high, then it will have a shape of a downward convex reflector, if the is correct, then the reflector will be flat on , and if is too low, then the reflector will be an upward convex reflector. During this process, the NMO correction of seismic signals often causes a distortion in the frequency domain. To overcome this problem, a stretching factor (10) is a provision threshold that is required. If a signal is found to have a stretching factor that exceeds the threshold, then a muting correction will be performed.

##### 2.6. Stacking

Stacking is a process of merging many seismic traces that are already in a state of zero offset. The stacking process is done by quantifying the average amplitude of all signals (11). This stacking process is expected to improve SNR (signal to noise ratio) to factor [10].

#### 3. Implementation

SEISGAMA is written using C# programming language provided by Microsoft Visual Studio 2015. Some scientific computing functions and the display are implemented using libraries from Ilnumerics Computing Engine and Ilnumerics Visualization Engine [22]. SEISGAMA has a modular design that aims to facilitate future development of the application. Some modules that are available in the application include the following:* GamaseisProcessing* and* GamaseisView*.* GamaseisProcessing* is on a SEISGAMA module that contains the reading and writing data processing and some physics and mathematical calculations of the seismic reflection method.* GamaseisView* is a module that controls the graphical user interface, including windows, graphs, seismic sections, labels, and buttons.

*GamaseisProcessing* consists of scripts to read data in a SEG-Y format, performing amplitude correction, muting, velocity analysis, normal moveout correction, and stacking. The SEG-Y seismic data storage format is the most popular to date, and it is used by almost all commercial seismic software. The reading function of the SEG-Y format is already supported in any form; however, until now, this reading function has been limited to two-dimensional seismic reflection data. Other mathematical functions that have been included in this program are implemented in accordance with existing theory. Calculation of Fourier transforms on SEISGAMA uses libraries that have been prepared and optimized by Ilnumerics [22], which implements a special compilation of C language functions to implement the algorithm called the Fastest Fourier Transform from the West (FFTW) [23]. Utilizing a premade library, which has been tested and uses a lower level programming language, is expected to improve the computation time of the Fourier transform. Thus, the calculation of the transformation to - and - domain is mostly done using FFTW algorithm.

GamaseisView consists of a set of views, views that govern the interface between user interaction with the software. Figure 1 shows the main user interface of SEISGAMA where users can perform seismic data processing easily.