Awali Priyono

Follower

Following

Public Views

Josiah Heyman

University of Texas at El Paso (UTEP)

Nina Glick Schiller

The University of Manchester

Armando Marques-Guedes

UNL - New University of Lisbon

Nando Sigona

University of Birmingham

Anna Horolets

University of Warsaw

Donita Manurung

UNIMED

astro physics

Universidad de Holguín

Rosita Renovita

The University Of Lampung

Jonathan Aprilio

Trisakti University

Fuji KIA

University of jambi

Interests

Uploads

Papers by Awali Priyono

Pre-stack 3D Land Seismic Data Merging in Cemara, Akasia Bagus and Tugu Areas of Northwest Java Basin, Indonesia

6th Saint Petersburg International Conference and Exhibition, 2014

Extensive exploratory effort in Cemara, Akasia Bagus and Tugu areas has produced three different ... more Extensive exploratory effort in Cemara, Akasia Bagus and Tugu areas has produced three different 3D land Survey till now. Processing and interpretation of these individual 3D Data sets have led to more exploration in the area. However, due to different acquisition and processing parameters, it was difficult to interpret subsurface image the overlap areas. The differences in Line orientation, bin size, charge size of sources, offset length, fold, amplitude, frequency contents, etc. have affected the quality of each seismic data sets. In order to solve these problems, pre-stack 3D Seismic data merging was conducted. The main goal of Pre-Stack 3D Merging is to combine these original data sets into one data sets, processed as one project at the same time using same processing software. In the end, the comparison of fold coverage distribution, composite seismic sections and time slices, after and before merging, was made and shown significantly improvement in subsurface images. These merge data sets can help interpreter to get better understanding of geological process in these area

Additional file 1 of Delineation of sedimentary basin structure beneath the Banyumas Basin, Central Java, Indonesia, using ambient seismic noise tomography

Additional file 1: Table S1. Facies analysis data of Jati-1 well (Tampubolon et al. 2014). Fig. S... more Additional file 1: Table S1. Facies analysis data of Jati-1 well (Tampubolon et al. 2014). Fig. S1. An overview of 1291 group velocity dispersion curves of Rayleigh wave in the 1 s-10 s period with a range of group velocity values of 0.4 km/s - 3.2 km/s. Fig. S2. The trade-off curve of the regularization damping and smoothing parameters in the 1 s period. A The search for damping values is varied from 5000 to 10. A damping number of 300 is obtained. B The search for smoothing values is varied from 3000 to 5. A smoothing number of 300 is obtained. Fig. S3. Inversion of the synthetic model of group velocity with Checkerboard Resolution Test (CRT) in each period with a cell size of ~ 6.5 km2, damping, and smoothing factors of 300. (A) 1 s, (B) 2 s, (C) 3 s, (D) 4 s, (E) 5 s, (F) 6 s, (G) 7 s, (H) 8 s, (I) 9 s, (J) 10 s, (K) initial synthetic model. Checkerboard patterns of synthetic models could be relatively recovered in the period 1 s - 6 s. Fig. S4. The amount of ray path in each Ra...

Seismotectonics of the 2019 Mw 6.5 Ambon Earthquake and its Aftershocks Based on Moment Tensor Inversion: Identification of Activated Fault Planes

SSRN Electronic Journal, 2022

Local earthquake tomography of the source region of the 2018 Lombok earthquake sequence, Indonesia

Geophysical Journal International, 2021

SUMMARY We develop and present a 3-D seismic velocity model of the source region of the 2018 Lomb... more SUMMARY We develop and present a 3-D seismic velocity model of the source region of the 2018 Lombok, Indonesia earthquakes by using local earthquake tomography. The data consist of 28 728 P- and 20 713 S-wave arrival times from 3259 events which were recorded by 20 local seismic stations. The results show that most of the significant earthquakes occur to the edge of high-velocity regions. We interpret these to represent coherent blocks of the Flores Oceanic Crust underthrusting Lombok. At depths shallower than the nucleation area of the largest earthquake, many triggered aftershocks are located within a low-velocity, high-Vp/Vs region which is probably a highly fractured fault zone with a large amount of fluid. This fault zone is parallel to the dip of the Flores Back Arc Thrust and probably ruptured during this earthquake sequence. A prominent low-velocity, high-Vp/Vs region is colocated with the northwest and southern flank of the Rinjani volcanic complex. This large aseismic regi...

Early Results of Eastern Indonesia P-wave Tomography Study Using Regional Events

IOP Conference Series: Earth and Environmental Science, 2021

The tectonic system of Eastern Indonesia is controlled by several major and minor plates, such as... more The tectonic system of Eastern Indonesia is controlled by several major and minor plates, such as Indo-Australian, Australian plate, and Pacific plates. This area is known for its complexity, and high seismic activity. This study tries to image the complex structures beneath this region by employing regional events data and seismic tomography methods. We used five years of regional events catalog provided by the Indonesian Agency of Meteorology, Climatology, and Geophysics. We have sorted 7336 events recorded between 120° – 136° longitude and 0° – 13°(-) latitude consisting of 46446 P and 15467 S wave arrival data. Relocated hypocenter map shows a better constrain location on seismicity along outer Bandar Arc. A dipping pattern of seismicity is seen that is going deeper to the Banda Sea. The seismicity map also images a steep angle pattern of seismicity that could be related to the subduction slab roll-back model at North of Wetar island. Interestingly, we spotted a seismicity gap i...

P-wave Attenuation Structure around the Source Region of 2018 Lombok, Indonesia, Earthquakes

AGU Fall Meeting Abstracts, Dec 1, 2020

High Resolution 3D Shallow Crustal Structure of the Lembang Fault, West Java, Indonesia, produced with Ambient Noise Tomography

The Lembang Fault is located about eight km north of Bandung City, the capital city of the West J... more The Lembang Fault is located about eight km north of Bandung City, the capital city of the West Java province of Indonesia, with about 8.79 million people (dated 2020). Several studies state that the Lembang fault is an active fault with a maximum magnitude of 6.5 to 7. In terms of earthquake disaster mitigation, distinguishing characteristics of seismic velocity structures around active faults is critical. This manuscript presents our new research on ambient seismic noise tomography around the Lembang Fault by deploying 73 local temporary seismometer networks to develop a 3D high-resolution shear wave velocity model. We conducted a design acquisition analysis to obtain optimal results using a ray density analysis approach to image the geometry of the Lembang fault with good illumination and resolution, which is essential for identifying shallow geological features. Our result showed that variations in velocity values, both laterally and vertically, were related to variations in the...

Download

Ambient Seismic Noise Cross - correlation of Ambon Island and Surrounding Area, Eastern Indonesia: Preliminary Result

IOP Conference Series: Earth and Environmental Science, 2021

The island of Ambon lies on complex tectonics, part of Banda Arc which is driven by the Australia... more The island of Ambon lies on complex tectonics, part of Banda Arc which is driven by the Australia – Eurasia collision. Historical earthquake data show that an earthquake resulting the greatest tsunami in Indonesia had occurred at Ambon Island. On 26 September 2019, Ambon was shaken by an M 6.5 earthquake at a depth of 10 km (BMKG). In this study, we use ambient noise data from 11 temporary stations deployed by ITB and 4 permanent stations owned BMKG which are recorded from October until December 2019. Here, we purely use the vertical component of seismogram to retrieve the Empirical Green’s Function of Rayleigh waves. Cross-correlations were obtained from the daily data series and stacked the day-by-day cross-correlation data into one inter-station cross-correlation. The Empirical Green’s Function is seen at the band period 1-15 s. As a part of our study, we analyze the Green’s Function with frequency-time analysis (FTAN) to get Rayleigh wave group velocity. The group velocity of Ra...

Shear Wave Splitting of the 2018 Lombok Earthquake Aftershock Area, Indonesia

Lombok is one of the islands in the transitional zone from the Sunda Arc to the Banda Arc, Indone... more Lombok is one of the islands in the transitional zone from the Sunda Arc to the Banda Arc, Indonesia. In the mid-2018, the island of Lombok was shaken by a series of strong earthquakes, started with a magnitude 6.4 earthquake on July 29, 2018 followed by earthquakes on August 5 (M 7.0), August 9 (M 5.9), and August 19 (M 6.3 and 6.9). Some researchers suggested that this phenomenon occurred due to a segmentation rupture in the northern part of Lombok Island. This study aims to obtain information on the distribution of the Lombok earthquake rupture zone 2018, through Shear Wave Splitting (SWS) study. Splitting, or S-wave separation, occurs when the S-wave passes through an anisotropic medium. The S wave is split into fast and slow S waves with almost orthogonal polarizations and has parameters such as delay time and polarization direction of the fast S wave. To determine the SWS parameters, we used a Lombok earthquake aftershock data set recorded from 4 August to 9 September 2018, us...

Download

Joint Inversion of Rayleigh-wave Dispersion Curve for Near-Surface S-Wave Velocity Estimation

Summary The most important step in surface wave method is inverting Rayleigh wave dispersion curv... more Summary The most important step in surface wave method is inverting Rayleigh wave dispersion curve for near-surface shear-wave velocity estimation. Joint inversion of Rayleigh wave dispersion curve has been developed. We combine the advantages of two different optimization methods, the genetic algorithm (GA) and Occam’s Algorithm. First, we use GA with small number generation to get the possible output solution, then we improve the solution (shear-wave velocity and thickness of layer) using Occam’s algorithm. We use several synthetic models to show performance of the joint inversion. Numerical errors are shown to appraise the value accuracy inversion model. Inversion of the synthetic dispersion curves for the four layers model provides most accurate solution models with is 0.0005 % and is 0.001%. Whereas, that for the six layers model also presents accurate result with is 1.48% and is 12.51 %. These errors indicate that the accuracy apply successfully join inversion scheme and promi...

Preliminary Results of Automatic P-Wave Regional Earthquake Arrival Time Picking using Machine Learning with Kurtosis and Skewness as The Input Parameters

IOP Conference Series: Earth and Environmental Science, 2021

The traditional method in determining first arrival time of earthquake dataset is time consuming ... more The traditional method in determining first arrival time of earthquake dataset is time consuming process due to waveform manual inspection. Additional waveform attributes can help determine events detection. One of the widely used attribute is The Short Term Averaging/Long Term Averaging (STA/LTA) which is simply division moving average of waveform amplitude between short time and longer time. Alternatively, waveform attribute can also be built using kurtosis and skewness. The kurtosis attribute is defined as sharpness of data distribution. By definition, the maximum signal should be at or close to the P wave arrival. The skewness is typically used to show normal distribution of the data. The uniqueness of this method is that it has an ability to determine whether the first P wave arrival has positive of negative number. The skewness calculation is similar to kurtosis but it uses the power of 3 instead of 4. With the objective of generating efficient scheme to pick first time arriva...

Baseline Seismic Survey for Gundih CCS Pilot Project

SSRN Electronic Journal, 2019

Gundih CCS Pilot Project will enter the new phase of implementation in 2019. Monitoring of CO 2in... more Gundih CCS Pilot Project will enter the new phase of implementation in 2019. Monitoring of CO 2injection must be planned carefully before CO 2-injection is started. One of proposed monitoring method is time-lapse reflection seismic survey in the vicinity of Jepon-1 well, i.e. the candidate of well injection, which is situated in Blora District, Central Java-Indonesia. Before injection of CO 2 is started, a pseudo 3D seismic survey was then carried out in 2016. The objective of this survey is

Download

Study on linearized inversion of surface wave dispersion for estimating near-surface shear wave velocity

IOP Conference Series: Earth and Environmental Science, 2019

Near-surface shear wave velocity is very important for geophysical and geotechnical engineering. ... more Near-surface shear wave velocity is very important for geophysical and geotechnical engineering. The use of surface wave for characterization of the near-surface subsurface involves three steps: acquisition, extracting dispersion curve, and inversion method to generate near-surface shear wave velocity. In this study, we explored a linear inversion technique to see the effect of some parameters such as the initial model and the damping factor on the inversion results. We used the Levenberg-Marquardt algorithm and written in Matlab. Synthetic data used three-layer earth model with assumed homogeneous, isotropic, flat, and half-space. Synthetic examples demonstrated that inverted result depends on the initial model. Based on our results, we have ranged from 30% to 120% of the true model in providing the initial model. It is found that there is an asymmetry between the range of the initial model which is larger and smaller than the true model. This is due to the non-linearity of the dif...

Download

A Resolution Enhancement of Rayleigh Wave Dispersive Imaging using Modified Phase-Shift Method

IOP Conference Series: Earth and Environmental Science, 2019

Multi-channel analysis of surface wave (MASW) method is an efficient tool for obtaining vertical ... more Multi-channel analysis of surface wave (MASW) method is an efficient tool for obtaining vertical shear wave profile. One of the key in the steps in the MASW method is generate a dispersion curves. To generate dispersion energy, the acquired shot gathers data have been transformed from the time-space (t-x) domain to the frequency-phase velocity (f-v) domain using 2D wavefield transformation so that the dispersion curve can be determined by picking the peak of dispersion energy. In this paper, we proposed a modified phase-shift method to increase the resolution of Rayleigh wave dispersive imaging so that the accuracy of the dispersion curve determination becomes higher. This method added processing step on the conventional phase-shift method. The algorithm is written in Matlab environment. We compared this method with the conventional phase-shift and frequency-wavenumber (F-K) method and applied to synthetic and real data. We used DISECA program in Matlab to generate synthetic wavefor...

Download

Offset aperture optimization of partial CRS stack based on AVO gradient for reservoir characterization

Proceedings of the 11th SEGJ International Symposium, Yokohama, Japan, 18-21 November 2013, 2013

Conventional seismic data processing based on the CMP stacking method can not utilize all the wav... more Conventional seismic data processing based on the CMP stacking method can not utilize all the wave energy that propagate outside the CMP line due to the complexity of the subsurface structure. CRS stack data processing method provides a solution by taking into account the wave energy outside the CMP area within Fresnel zone boundaries, especially in areas with complex structures. Partial CRS stack as a further developed method uses the local slope of the curve and reflector elements effectively to add up the amplitudes of all the seismic traces in the reflector element to produce a new super-gather. The reflector segment area is represented by two parameters, offset aperture in zero offset domain and offset aperture in offset domain. In order to maintain relative amplitude preservation, offset aperture optimization method was introduced based on AVO gradient patterns of every offset segment or angle in a CRS super-gather. The AVO gradient value on CRS super-gather, which are calculated based on offset aperture parameters, can be controlled by using the theoretical value of the AVO gradient estimated from well data. Based on the optimization results, we can collect the gathers in every segment and make a new super-gather with a gradient closest to the theoretical value. Optimization method for the partial CRS stack was tested on seismic data from sand reservoir. From the result of the offset aperture optimization showed that the AVO gradient of CRS data was deviated from the theoretical value by 3.5%, better than conventional CMP stack.

Estimation of Seismic Anisotropic Parameter using Velocity Dependent Offset Analysis

3rd Annual International Conference on Geological & Earth Sciences (GEOS 2014), 2014

Partial CRS stack seismic data characterization on AVO anomaly: case study on 2D synthetic seismic data

Proceedings of the 10th SEGJ International Symposium, Kyoto, Japan, 20-22 November 2011, 2011

Partial CRS stack is one of alternative seismic data processing methods, which does not depend on... more Partial CRS stack is one of alternative seismic data processing methods, which does not depend on macro velocity model. One of the advantages of this method compared with conventional CMP stack is the summation process in each CMPs in order to obtain better signal quality. CRS travel time formula can be applied to obtain partial stacked CRS supergathers. The CRS gathers obtained by this method are regularized on offset domain and have better S/N ratio than conventional gathers. This method is very helpful to improve the quality of land seismic data subject to the random noise and lack of fold coverage. The purpose of this research is to examine the efficiency of the partial CRS stack methods to highlight AVO/A (Amplitude Versus Offset/Angle) anomalies by using synthetic seismic data. The tested AVO anomaly is class 3 where the negative amplitude increases with increasing distance/offset. Seismic modeling is performed on a simple anticline model composed of shale and sand layers. The synthetic seismic data are provided to the processing using partial CRS stack to test the offset aperture that is considered to influence the seismic amplitude. AVO analysis was conducted on the pre-stack time migrated (PSTM) synthetic data with and without partial CRS stack for three offset apertures. The seismic amplitude responses change with the apertures. Dipping reflectors also give different amplitude response from flat reflectors. It can be inferred from the analysis that although the class of AVO anomalies on the partial CRS stack gathers are almost same, there should be an optimum offset aperture to avoid ambiguous result of the AVO attributes.

Anisotropic parameter estimation using velocity variation with offset analysis

Seismic anisotropy is defined as velocity dependent upon angle or offset. Knowledge about anisotr... more Seismic anisotropy is defined as velocity dependent upon angle or offset. Knowledge about anisotropy effect on seismic data is important in amplitude analysis, stacking process and time to depth conversion. Due to this anisotropic effect, reflector can not be flattened using single velocity based on hyperbolic moveout equation. Therefore, after normal moveout correction, there will still be residual moveout that relates to velocity information. This research aims to obtain anisotropic parameters, e and δ, using two proposed methods. The first method is called velocity variation with offset (VVO) which is based on simplification of weak anisotropy equation. In VVO method, velocity at each offset is calculated and plotted to obtain vertical velocity and parameter δ. The second method is inversion method using linear approach where vertical velocity, δ, and e is estimated simultaneously. Both methods are tested on synthetic models using ray-tracing forward modelling. Results show that δ value can be estimated appropriately using both methods. Meanwhile, inversion based method give better estimation for obtaining e value. This study shows that estimation on anisotropic parameters rely on the accuracy of normal moveout velocity, residual moveout and offset to angle transformation.

Erratum to “Improving the accuracy of the expanded anisotropic eikonal equation at larger offsets using Levin’s T-transformation”

Studia Geophysica et Geodaetica

Download

Determination of Shear Wave Splitting Parameters in 2018 Lombok Earthquake Using Rotation Correlation Method: Preliminary Result

IOP Conference Series: Earth and Environmental Science, 2021

Shear-wave splitting (SWS), or the propagation of two independent shear waves, can be used as an ... more Shear-wave splitting (SWS), or the propagation of two independent shear waves, can be used as an indicator of seismic anisotropy. In this study, we utilize this concept using aftershock data of the 2018 Lombok earthquake which had been acquired in period of August 4 – September 9, 2018. The goal of this research is to better understand the crack distribution related to the rupture zone of the 2018 Lombok earthquake. After applying instrument correction to the data, the waveform data were then windowed in each P and S arrival time. To determine the SWS parameters, we performed rotation in each horizontal seismogram components. The horizontal components were rotated from azimuth 0° to 180° with an increment of 1°. Cross-correlation coefficient (CCC) was determined for each rotation angle. The polarization direction and the SWS delay time were chosen from the parameters shown in the highest value of CCC.

Pre-stack 3D Land Seismic Data Merging in Cemara, Akasia Bagus and Tugu Areas of Northwest Java Basin, Indonesia

6th Saint Petersburg International Conference and Exhibition, 2014