Table of contents

Preface

Hendra Grandis¹, Andri D. Nugraha¹, Zulfakriza^1*, Rexha V. Ry¹, David P. Sahara¹, Endra Gunawan¹, Riky Hendrawan²

Department of Geophysical Engineering, Faculty of Mining and Petroleum Engineering, Institut Teknologi Bandung, Indonesia

Pertamina, Indonesia

*email: zulfakriza@itb.ac.id

This issue of Journal IOP Conference Series: Earth and Environmental Science contains the outstanding papers presented at The 41^th Scientific Annual Meeting of Indonesian Association of Geophysicists. This is a scientific meeting which was organized by Indonesian Association of Geophysicists. It was held on August 26 - 29 September, 2016 in Lampung – Indonesia.

All papers published in this volume of IOP Conference Series: Earth and Environmental Science have been peer reviewed through processes administered by the proceedings Editors. Reviews were conducted by expert referees to the professional and scientific standards expected of a proceedings journal published by IOP Publishing.

Seismic data with wide range of frequency is needed due to its close relation to resolution and the depth of the target. Low frequency provides deeper penetration for the imaging of deep target. In addition, the wider the frequency bandwidth, the sharper the wavelet. Sharp wavelet is responsible for high-resolution imaging and is very helpful to resolve thin bed. As a result, the demand for broadband seismic data is rising and it spurs the technology development of broadband seismic in oil and gas industry. An obstacle that is frequently found on marine seismic data is the existence of ghost that affects the frequency bandwidth contained on the seismic data. Ghost alters bandwidth to bandlimited. To reduce ghost effect and to acquire broadband seismic data, lots of attempts are used, both on the acquisition and on the processing of seismic data. One of the acquisition technique applied is the multi-level streamer, where some streamers are towed on some levels of depth. Multi-level streamer will yield data with varied ghost notch shown on frequency domain. If the ghost notches are not overlapping, the summation of multi-level streamer data will reduce the ghost effect. The result of the multi-level streamer data processing shows that reduction of ghost notch on frequency domain indeed takes place.

Reservoir rock containing oil and gas generally has high porosity and permeability. High porosity is expected to accommodate hydrocarbon fluid in large quantities and high permeability is associated with the rock's ability to let hydrocarbon fluid flow optimally. Porosity and permeability measurement of a rock sample is usually performed in the laboratory. We estimate the porosity and permeability of sandstones digitally by using digital images from μCT-Scan. Advantages of the method are non-destructive and can be applied for small rock pieces also easily to construct the model. The porosity values are calculated by comparing the digital image of the pore volume to the total volume of the sandstones; while the permeability values are calculated using the Lattice Boltzmann calculations utilizing the nature of the law of conservation of mass and conservation of momentum of a particle. To determine variations of the porosity and permeability, the main sandstone samples with a dimension of 300 × 300 × 300 pixels are made into eight sub-cubes with a size of 150 × 150 × 150 pixels. Results of digital image modeling fluid flow velocity are visualized as normal velocity (streamline). Variations in value sandstone porosity vary between 0.30 to 0.38 and permeability variations in the range of 4000 mD to 6200 mD. The results of calculations show that the sandstone sample in this research is highly porous and permeable. The method combined with rock physics can be powerful tools for determining rock properties from small rock fragments.

In compliance with the resolutions of IAGA (International Association of Geomagnetism and Aeronomy), Since 1960's, every five years BMKG or Meteorology, Climatology and Geophysics Agency of Indonesia made geomagnetic field maps based on actual measurements in 53 repeat stations. It's the map for more accurate result of Geomagnetic maps Epoch 2015.0, the number of repeat stations has been increased to 68 locations. Analysis data was conducted by spatial analyses using collocated co-kriging and kriging with external drift to map the observation data in five components, such as Declination (D), Inclination (I), Vertical (Z), Horizontal (H), and Total Geomagnetic Field (F). The data reduction used one permanent observatory i.e., Kupang Geophysical Observatory, as a reference standard. The results of this Geomagnetic Maps, that the contour lines of Indonesian geomagnetic declination in range -1 to 4.5 degree, Inclination component are -5 to -37 degree, Vertical component are -4000 to -28000 nT, Horizontal component are 36000 to 42000 nT, and Total Geomagnetic Field are 39000 to 46000 nT. In conclusion, Indonesian Geomagnetic Maps for Epoch 2015.0 can be used to compute geomagnetic data around Indonesian regions until next 5 years.

Rock physical parameters value (Vp and Vs) is one of fundamental aspects in reservoir characterization as a tool to detect rock heterogenity. Its response is depend on several reservoir conditions such as lithology, pressure and reservoir fluids. The value of Vp and Vs is controlled by grain contact and contact stiffness, a function of clay mineral content and porosity also affected by mineral composition. The study about Vp and Vs response within sandstone and its relationship with petrographic parameters has become important to define anisotrophy of reservoir characteristics distribution and could give a better understanding about local diagenesis that influence clastic reservoir properties. Petrographic analysis and Vp-Vs calculation was carried out to 12 core sample which is obtained by hand-drilling of the outcrop in Sukabumi area, West Java as a part of Bayah Formation. Data processing and interpretation of sedimentary vertical succession showing that this outcrop comprises of 3 major sandstone layers indicating fluvial depositional environment. As stated before, there are 4 petrographic parameters (sorting, roundness, clay mineral content, and grain contact) which are responsible to the differences of shear wave and compressional wave value in this outcrop. Lithology with poor-sorted and well- roundness has Vp value lower than well-sorted and poor-roundness (sub-angular) grain. For the sample with high clay content, Vp value is ranging from 1681 to 2000 m/s and could be getting high until 2190 to 2714 m/s in low clay content sample even though the presence of clay minerals cannot be defined neither as matrix nor cement. The whole sample have suture grain contact indicating telogenesis regime whereas facies has no relationship with Vp and Vs value because of the different type of facies show similar petrographic parameters after diagenesis.

Banggai-Sula foreland basin in Matindok block is the product of Late Miocene to Early Pliocene collision between Banggai-Sula microcontinent and East Sulawesi Ophiolite (ESO)-magmatic arc of Sundaland. Proven petroleum system plays for Donggi, Senoro, and Matindok gas fields are Miocene build-up carbonate and platform carbonates related to wrench structures, sealed by Plio-Pleistocene Celebes Molasse consists of conglomerate, conglomeratic sandstone, shale and charged during early Pliocene related to collisional tectonic from Miocene marine clastic source rocks. Latest interpretation of 3D & 2D seismic datasets on onshore area of Matindok block shown that there is possibility of NE to SW trending channel systems existed in Plio-Pleistocene Celebes Molasse as post collision sediments. Based on biostratigraphy interpretation, Celebes Molasse interval was deposited in marine environment, generally outer sublittoral to upper bathyal. Channel filled pattern on seismic data expression shown onlap, divergent and mounded-onlap fills on 6 (six) interpreted channel bodies. Fining upward sequences is well recognized based on gamma ray log from several wells that penetrated those interpreted channel. The thicknesses of channel bodies vary from 50 m to 150 m with sand to shale ratio up to 40 % and gas composition reading from C1 to C3. Despite of there is no well testing data on these channels interval; we believed that this Pliocene deep water channel can become the new exploration play although several studies about seal effectiveness of channel bodies as stratigraphic traps, lateral and vertical migration pathways need to be carried out.

Paleogeography is one of critical points that always less considered by explorationist in the world. Almost all of the consideration is focused on trapping mechanism. Paleogeography is guidance in understanding both of physical and chemical of rock characteristic which will correlate with its depositional environment. Integration of various geological and geophysical data such as; tectonic, structural geology, stratigraphy, lithology, and biostratigraphy will lead us to a better understanding of rock characteristics. Six paleogeographic interpretations was made consist of; Early Tertiary (P5-56-55 ma), Middle Eocene (P14-41 ma), Late Oiligocene (P22-25.5 ma), Early Miocene (N7-16.5 ma), Middle Miocene (N9-14.5 ma), and Pleistocene (NN19-1.5 ma). That six paleogeographic interpretations are assumed represent the paleogeographic evolution of East Java Basin time after time. In Middle Eocene time, it would be more than hundred possibilities regarding the location where the formation deposited. This would be controlled by the existence of some local structural paleohighs and horsts which oriented NW-SE followed by their own sedimentary transportation path. With assumption that hydrocarbon generation was occurred in 15 Ma and the depth of maturation window lies on about 2,500 m depth. Therefore, the possibility of source rock maturation is high, due to almost of the clastics sediment of Ngimbang deposited into the series of grabens. The Kujung reef types simplified defines and categorize into; 1) Patch Reef 2) Berrier Reef 3) Pinnacle Reef Over Isolated Reef. Kujung Carbonates were deposited in Early Miocene when regional transgression occurred. The depositional environments were dominated by shallow marine littoral-sublittoral. Generally, the reservoir quality of this Kujung Carbonate shows fair to good quality, in range7-32% porosity, and 1-1400 mD permeability (internal SKK Migas data).

Geometry is an important parameter for the field of hydrocarbon exploration and exploitation, it has significant effect to the amount of resources or reserves, rock spreading, and risk analysis. The existence of geological structure or fault becomes one factor affecting geometry. This study is conducted as an effort to enhance seismic image quality in faults dominated area namely offshore Madura Strait. For the past 10 years, Oligo-Miocene carbonate rock has been slightly explored on Madura Strait area, the main reason because migration and trap geometry still became risks to be concern. This study tries to determine the boundary of each fault zone as subsurface image generated by converting seismic data into variance attribute. Variance attribute is a multitrace seismic attribute as the derivative result from amplitude seismic data. The result of this study shows variance section of Madura Strait area having zero (0) value for seismic continuity and one (1) value for discontinuity of seismic data. Variance section shows the boundary of RMKS fault zone with Kendeng zone distinctly. Geological structure and subsurface geometry for Oligo-Miocene carbonate rock could be identified perfectly using this method. Generally structure interpretation to identify the boundary of fault zones could be good determined by variance attribute.

Time to depth conversion is an important processof seismic interpretationtoidentify hydrocarbonprospectivity. Main objectives of this research are to minimize the risk of error in geometry and time to depth conversion. Since it's using a large amount of data and had been doing in the large scale of research areas, this research can be classified as a regional scale research. The research was focused on three horizons time interpretation: Top Kujung I, Top Ngimbang and Basement which located in the offshore and onshore areas of east Java basin. These three horizons was selected because they were assumed to be equivalent to the rock formation, which is it has always been the main objective of oil and gas exploration in the East Java Basin. As additional value, there was no previous works on velocity modeling for regional scale using geological parameters in East Java basin. Lithology and interval thickness were identified as geological factors that effected the velocity distribution in East Java Basin. Therefore, a three layer geological model was generated, which was defined by the type of lithology; carbonate (layer 1: Top Kujung I), shale (layer 2: Top Ngimbang) and Basement. A statistical method using three horizons is able to predict the velocity distribution on sparse well data in a regional scale. The average velocity range for Top Kujung I is 400 m/s - 6000 m/s, Top Ngimbang is 500 m/s - 8200 m/s and Basement is 600 m/s - 8000 m/s. Some velocity anomalies found in Madura sub-basin area, caused by geological factor which identified as thick shale deposit and high density values on shale. Result of velocity and depth modeling analysis can be used to define the volume range deterministically and to make geological models to prospect generation in details by geological concept.

The potential sources of meteorological phenomena in Nuclear Power Plant (NPP) area of interest are identified and the extreme values of the possible resulting hazards associated which such phenomena are evaluated to derive the appropriate design bases for the NPP. The appropriate design bases shall be determined according to the Nuclear Energy Regulatory Agency (Bapeten) applicable regulations, which presently do not indicate quantitative criteria for purposes of determining the design bases for meteorological hazards. These meteorological investigations are also carried out to evaluate the regional and site specific meteorological parameters which affect the transport and dispersion of radioactive effluents on the environment of the region around the NPP site. The meteorological hazards are to be monitored and assessed periodically over the lifetime of the plant to ensure that consistency with the design assumptions is maintained throughout the full lifetime of the facility.

Model based inversion was applied to inversion process of 2D seismic data in Kangean Offshore Area. Integration acoustic impedance from wells and seismic data was expected showing physical property, facies separation and reservoir quality of carbonate rock, particularly in Kangean Offshore Area. Quantitative and qualitative analysis has been conducted on the inversion results to characterize the carbonate reservoir part of Kujung and correlate it to depositional facies type. Main target exploration in Kangean Offshore Area is Kujung Formation (Oligo-Miocene Carbonate). The type of reservoir in this area generate from reef growing on the platform. Carbonate rock is a reservoir which has various type and scale of porosity. Facies determination is required to to predict reservoir quality, because each facies has its own porosity value. Acoustic impedance is used to identify and characterize Kujung carbonate facies, also could be used to predict the distribution of porosity. Low acoustic impedance correlated with packstone facies that has acoustic impedance value below 7400 gr/cc*m/s. In other situation, high acoustic impedance characterized by wackestone facies above 7400 gr/cc*m/s. The interpretation result indicated that Kujung carbonate rock dominated by packstone facies in the upper part of build-up and it has ideal porosity for hydrocarbon reservoir.

This paper describes a method to build a near-surface velocity model. The purpose of this work is to have a better understanding of the near-surface condition below area with high variation topography and high surface energy attenuation. The velocity model built then was used as input for statics calculation. The method applied to this work is first arrival tomography. This method was tested using two datasets with different acquisition geometries. And to verify the result, we compare the calculated statics correction with a conventional method. The results of this work is two similar near-surface models with detail model and the statics applied on the seismic data showed better correction compare to the conventional one. We conclude that the method is applicable accurately on the area with high variation of topography and high surface energy attenuation.

Sulawesi, one of the biggest island in Indonesia, located on the convergence of two macro plate that is Eurasia and Pacific. NOAA and Novosibirsk Tsunami Laboratory show more than 20 tsunami data recorded in Sulawesi since 1820. Based on this data, determination of correlation between tsunami and earthquake parameter need to be done to proved all event in the past. Complete data of magnitudes, fault sizes and tsunami heights on this study sourced from NOAA and Novosibirsk Tsunami database, completed with Pacific Tsunami Warning Center (PTWC) catalog. This study aims to find correlation between moment magnitude, fault size and tsunami height by simple regression. The step of this research are data collecting, processing, and regression analysis. Result shows moment magnitude, fault size and tsunami heights strongly correlated. This analysis is enough to proved the accuracy of historical tsunami database in Sulawesi on NOAA, Novosibirsk Tsunami Laboratory and PTWC.

Advances in science and technology showed that ground-based GPS receiver was able to detect ionospheric Total Electron Content (TEC) disturbances caused by various natural phenomena such as earthquakes. One study of Tohoku (Japan) earthquake, March 11, 2011, magnitude M 9.0 showed TEC fluctuations observed from GPS observation network spread around the disaster area. This paper discussed the ionospheric earthquake effects detection using TEC GPS data. The case studies taken were Kebumen earthquake, January 25, 2014, magnitude M 6.2, Sumba earthquake, February 12, 2016, M 6.2 and Halmahera earthquake, February 17, 2016, M 6.1. TEC-GIM (Global Ionosphere Map) correlation methods for 31 days were used to monitor TEC anomaly in ionosphere. To ensure the geomagnetic disturbances due to solar activity, we also compare with Dst index in the same time window. The results showed anomalous ratio of correlation coefficient deviation to its standard deviation upon occurrences of Kebumen and Sumba earthquake, but not detected a similar anomaly for the Halmahera earthquake. It was needed a continous monitoring of TEC GPS data to detect the earthquake effects in ionosphere. This study giving hope in strengthening the earthquake effect early warning system using TEC GPS data. The method development of continuous TEC GPS observation derived from GPS observation network that already exists in Indonesia is needed to support earthquake effects early warning systems.

Overpressure is an abnormal high subsurface pressure of any fluids which exceeds the hydrostatic pressure of column of water or formation brine. In Offshore Nova Scotia Canada, the values and depth of overpressure zone are determined using the eaton and equivalent depth method, based on well data and the normal compaction trend analysis. Since equivalent depth method is using effective vertical stress principle and Eaton method considers physical property ratio (velocity). In this research, pressure evaluation only applicable on Penobscot L-30 well. An abnormal pressure is detected at depth 11804 feet as possibly overpressure zone, based on pressure gradient curve and calculation between the Eaton method (7241.3 psi) and Equivalent Depth method (6619.4 psi). Shales within Abenaki formation especially Baccaro Member is estimated as possible overpressure zone due to hydrocarbon generation mechanism.

Sulphide minerals are any member of a group of compounds of sulphur with one or more metals. Some of these sulphide minerals are economically important. This study used induced polarization method to identify distribution and mineralized zone of sulphide mineral (Pyrite), in Libureng, Bone Regency, South Sulawesi. The data processing yielded resistivity value, percent frequency effect (PPE) value, and metal factor (MF) value which were then used to produce 2-D and 3-D section model. Based on the data interpretation, an anomaly linked to pyrite deposits was seen in four trajectories with resistivity value < = 50, PFE = > 3%, and MF > = 150, deposited in hydrothermal alteration zone, sericite zone.

The resistivity of Magnetotelluric (MT) data show the resistivity mapping in the volcanic reservoir zone and the geochemistry information for confirm the reservoir and source rock formation. In this research, we used 132 data points divided with two line at exploration area. We used several steps to make the resistivity mapping. There are time series correction, crosspower correction, then inversion of Magnetotelluric (MT) data. Line-2 and line-3 show anomaly geological condition with Gabon fault. The geology structure from the resistivity mapping show the fault and the geological formation with the geological rock data mapping distribution. The geochemistry information show the maturity of source rock formation. According to core sample analysis information, we get the visual porosity for reservoir rock formation in several geological structure. Based on that, we make the geological modelling where the potential reservoir and the source rock around our interest area.

East Cepu High is a part of horst – graben series which formed by extensional tectonic processes during Paleogene in North East Java Basin. Due to excellent paleogeography position, the carbonate build-up was growth very well and as the main reservoir in East Cepu High. Sea level change have important factor to provide variation of facies in each carbonate buildup, one of emerging facies is detrital carbonate. Detrital carbonate indicated by onlap horizon featured with carbonate build up body. Based on paleogeography, fluctuation of sea level change and sediment source, detrital carbonate formed in leeward area in lowstand or highstand phases. Distinguish between detrital carbonate facies with other facies, advanced seismic processing performed by using continuous wavelet transform (CWT) and seismic inversion. CWT is one method of spectral decomposition used to find the frequency that represent a facies. The result from seismic inversion will support the interpretation for facies distribution. As the result, seismic data which have interval frequency 10 – 45 Hz and Acoustic Impedance (AI) value above 35000 (from cross plot between acoustic impedance and gamma ray) can be interpreted as detrital carbonate. Based on seismic interpretation, detrital carbonate facies distributed along leeward area with geometrical spreading. The lateral facies change from detrital carbonate to shale was identified which causing this facies become potential as hydrocarbon reservoir with stratigraphic trap. Based on the earlier studies, North East Java Basin have a strong hydrocarbon migration to fill the reservoir, therefore the detrital carbonate have high chance to be a new hydrocarbon prospect in this area.

Exploration activities with standard methods have already encountered many obstacles in the field. Geological survey is often difficult to find outcrop because they are covered by vegetation, alluvial layer or as a result of urban development and housing. Seismic method requires a large expense and licensing in the use of dynamite is complicated. Method of gravity requires the operator to go back (looping) to the starting point. Given some of these constraints, therefore it needs a solution in the form of new method that can work more efficiently with less cost. Several studies in various countries have shown a correlation between the presence of hydrocarbons and Radon gas concentration in the earth surface. By utilizing the properties of Radon that can migrate to the surface, the value of Radon concentration in the surface is suggested to provide information about the subsurface structure condition. Radon is the only radioactive substance that gas-phased at atmospheric temperature. It is very abundant in the earth mantle. The vast differences of temperatures and pressures between the mantle and the earth crust cause the convection flow toward earth surface. Radon in gas phase will be carried by convection flow to the surface. The quantity of convection currents depend on the porosity and permeability of rocks where Radon travels within, so that Radon concentration in the earth surface delineates the porosity and permeability of subsurface rock layers. Some measurements were carried out at several locations with various subsurface geological conditions, including proven oil fields, proven geothermal field, and frontier area as a comparison. These measurements show that the average and the background concentration threshold in the proven oil field (11,200 Bq/m³) and proven geothermal field (7,820 Bq/m³) is much higher than the quantity in frontier area (329 and 1,620 Bq/m³). Radon concentration in the earth surface is correlated with the presence of geological faults. Peak concentrations of Radon takes place along the fault.

The self-potential (SP) method is one of the oldest geophysical methods that are still available for today's application. Since its early days SP data interpretation has been done qualitatively until the emerging of the fixed geometry analysis that was used to characterize the orientation and the electric-dipole properties of a mineral ore structure. Through the expansion of fundamental theories, computational methods, field-and-lab experiments in the last fifteen years, SP method has emerge from its low-class reputation to become more respectable. It became a complementary package alongside electric-resistivity tomography (ERT) for detecting groundwater flow in the subsurface, and extends to the hydrothermal flow in geothermal areas. As the analysis of SP data becomes more quantitative, its potential applications become more diverse. In this paper, we will show examples of our current SP studies such as the groundwater flow characterization inside a fault area. Lastly we will introduce the application of the "active" SP method - that is the seismoelectric method - which can be used for 4D real-time monitoring systems.

Anelastic attenuation process during seismic wave propagation is the trigger of seismic non-stationary characteristic. An absorption and a scattering of energy are causing the seismic energy loss as the depth increasing. A series of thin reservoir layers found in the study area is located within Talang Akar Fm. Level, showing an indication of interpretation pitfall due to attenuation effect commonly occurred in deeper level seismic data. Attenuation effect greatly influences the seismic images of deeper target level, creating pitfalls in several aspect. Seismic amplitude in deeper target level often could not represent its real subsurface character due to a low amplitude value or a chaotic event nearing the Basement. Frequency wise, the decaying could be seen as the frequency content diminishing in deeper target. Meanwhile, seismic amplitude is the simple tool to point out Direct Hydrocarbon Indicator (DHI) in preliminary Geophysical study before a further advanced interpretation method applied. A quick-look of Post-Stack Seismic Data shows the reservoir associated with a bright spot DHI while another bigger bright spot body detected in the North East area near the field edge. A horizon slice confirms a possibility that the other bright spot zone has smaller delineation; an interpretation pitfall commonly occurs in deeper level of seismic. We evaluates this pitfall by applying Gabor Deconvolution to address the attenuation problem. Gabor Deconvolution forms a Partition of Unity to factorize the trace into smaller convolution window that could be processed as stationary packets. Gabor Deconvolution estimates both the magnitudes of source signature alongside its attenuation function. The enhanced seismic shows a better imaging in the pitfall area that previously detected as a vast bright spot zone. When the enhanced seismic is used for further advanced reprocessing process, the Seismic Impedance and Vp/Vs Ratio slices show a better reservoir delineation, in which the pitfall area is reduced and some morphed as background lithology. Gabor Deconvolution removes the attenuation by performing Gabor Domain spectral division, which in extension also reduces interpretation pitfall in deeper target seismic.

Very Low Frequency (VLF) measurement has been done at Pagerkandang Volcanic, Dieng Volcanic Complex (DVC) to examine the possible existence of conductive zones that related with geothermal manifestation. VLF – EM survey used tilt mode with T-VLF BRGM Iris Instrument operated with two frequencies, they are 22200 Hz from Japan (JJI) and 19800 Hz from Australia (NWC). There are five lines with distance between lines is 50 m, and distance between measure points is 20 m. The parameters measured from VLF method are tilt angle (%) and elliptisity (%). Data processed by tilt angle value with fraser and Karous – Hjelt filter used WinVLF program. Karous – Hjelt filter resulted current density contour to estimate lateral location from conductive and resistive zones. The conductive zone is interpreted as the area which have high current density value. This area located at eastern dan western of Pagerkandang Volcanic. The conductive zone related to geothermal manifestation as like as fumarol that appeared because presenced of normal fault. Whereas the resistive zone is interpreted as the area which have low current density value. This area spread almost in the middle of the Pagerkandang Volcanic. The resistive zone was caused by the high weathering in claystone.

The need of ERT method in nickel latentes exploration can't be deny. This method have capability to make exploration more effective and efficient. In reality this method still remain ambiguity in its application, especially for geologist and mining expert. These ambiguity related with layer zone determination (limonite, saprolite and bedrock), the same resistivity values in the different zones and determination of bedrock. This paper try to expose interesting fact to overcome this ambiguity by using ERT data and drill data tie. This tie will show characteristic of nickel lateric based on resisitivity value and the contribution of chemistry element for resistivity value. Data ERT was collected by using gradient configuration and well data consist of mayor element and minor element. Tie result showed difference resistivity value in limonite layer influence by Fe, H₂O and Ni, where resistivity value from saprolite layer influenced by Fe, H₂O, Si0₂, MgO, Al, Cr, and Ni in certain accumulation. In bedrock layer, almost all drill data did not reach bedrock but only reached the boulder alone, it is supported by the value of the resistivity of rock unserpentinized peridotite which should show a relatively large resistivity value

Geologically, Ketaling area consists of a local high considered as flexure margin of Tempino-Kenali Asam Deep in west part and graben in east part also known as East Ketaling Deep. Numerous proven plays were established in Ketaling area with reservoir in early Miocene carbonate and middle Miocene sand. This area underwent several major deformations. Faults are developed widely, yet their geometrical features and mechanisms of formation remained so far indistinct, which limited exploration activities. With new three-dimensional seismic data acquired in 2014, this area evidently interpreted as having strike-slip mechanism. The objective of this study is to examine characteristic of strike slip fault and its affect to hydrocarbon trapping in Ketaling Area. Structural pattern and characteristic of strike slip fault deformation was examined with integration of normal seismic with variance seismic attribute analysis and the mapping of Syn-rift to Post-rift horizon. Seismic flattening on 2D seismic cross section with NW-SE direction is done to see the structural pattern related to horst (paleohigh) and graben. Typical flower structure, branching strike-slip fault system and normal fault in synrift sediment clearly showed in section. An echelon pattern identified from map view as the result of strike slip mechanism. Detail structural geology analysis show the normal fault development which has main border fault in the southern of Ketaling area dipping to the Southeast-East with NE-SW lineament. These faults related to rift system in Ketaling area. NW-SE folds with reactive NE-SW fault which act as hydrocarbon trapping in the shallow zone. This polyphase tectonic formed local graben, horst and inverted structure developed a good kitchen area (graben) and traps (horst, inverted structure). Subsequently, hydrocarbon accumulation potentials such as basement fractures, inverted syn-rift deposit and shallow zone are very interesting to explore in this area.

We present a magnitude of completeness (Mc) quantification based on BMKG improved earthquake catalog which generated from Ina-TEWS seismograph network. The Mc quantification can help us determine the lowest magnitude which can be recorded perfectly as a function of space and time. We use the BMKG improved earthquake catalog from 2008 to 2016 which has been converted to moment magnitude (Mw) and declustered. The value of Mc is computed by determining the initial point of deviation patterns in Frequency Magnitude Distribution (FMD) chart following the Gutenberg-Richter equations. In the next step, we calculate the temporal variation of Mc and b-value using maximum likelihood method annually. We found that the Mc value is decreasing and produced a varying b-value. It indicates that the development of seismograph network from 2008 to 2016 can affect the value of Mc although it is not significant. We analyze temporal variation of Mc value, and correlate it with the spatial distribution of seismograph in Indonesia. The spatial distribution of seismograph installation shows that the western part of Indonesia has more dense seismograph compared to the eastern region. However, the eastern part of Indonesia has a high level of seismicity compared to the western region. Based upon the results, additional seismograph installation in the eastern part of Indonesia should be taken into consideration.

The Cimandiri fault which is running in the direction from Pelabuhan Ratu to Padalarang is the longest fault in West Java with several previous shallow earthquakes in the last 20 years. By using continues and campaign GPS observation from 2006-2016, we obtain the deformation pattern along the fault through the variation of strain tensor. We use the velocity vector of GPS station which is fixed in stable International Terrestrial Reference Frame 2008 to calculate horizontal strain tensor. Least Square Collocation is applied to produce widely dense distributed velocity vector and optimum scale factor for the Least Square Weighting matrix. We find that the strain tensor tend to change from dominantly contraction in the west to dominantly extension to the east of fault. Both the maximum shear strain and dilatation show positive value along the fault and increasing from the west to the east. The findings of strain tensor variation along Cimandiri Fault indicate the post seismic effect of the 2006 Java Earthquake.

The common usage of a vertical seismic profile is to capture the reflection wavefield (upgoing wavefield) so that it can be used for further well tie or other interpretations. Borehole Seismic (VSP) receivers capture the reflection from below the well trajectory, traditionally no seismic image information above trajectory. The non-traditional way of processing the VSP multiple can be used to expand the imaging above the well trajectory. This paper presents the case study of using VSP downgoing multiples for further non-traditional imaging applications. In general, VSP processing, upgoing and downgoing arrivals are separated during processing. The up-going wavefield is used for subsurface illumination, whereas the downgoing wavefield and multiples are normally excluded from the processing. In a situation where the downgoing wavefield passes the reflectors several times (multiple), the downgoing wavefield carries reflection information. Its benefit is that it can be used for seismic tie up to seabed, and possibility for shallow hazards identifications. One of the concepts of downgoing imaging is widely known as mirror-imaging technique. This paper presents a case study from deep water offshore Vietnam. The case study is presented to demonstrate the robustness of the technique, and the limitations encountered during its processing.

Start your abstract Batu Dinding Kilo Tiga or famously called Batu Dinding is one of tourist destinations in Kilo Tiga village, Amurang district, South Minahasa regency, North Sulawesi. In terms of Geology, Batu Dinding is a rock cliff with columnar structure and constituted in the form of andesitic thrakit rocks which is one of intermediate igneous rocks of volcanic aged Late Tertiary until Early Quarter, precisely during the Pleistocene. The structure and appearance of Batu Dinding was like a ladder, make Batu Dinding in demand as a destination for wall climbing or cliff hanger. The area around the Batu Dinding, there are Pamurapa River that at a time when it was high tide, can be used for rafting. Unfortunately, the accommodation and accessibility of the Amurang Batu Dinding are poor, causing Batu Dinding rarely to visit. The aim of the paper is to introduce as well as media publications of Batu Dinding tourism in order to increase local income and in addition it can be as a geological site for learning purposes.

Coal Bed Methane (CBM) as a newly developed resource in Indonesia is one of the alternatives to relieve Indonesia's dependencies on conventional energies. Coal resource of Muara Enim Formation is known as one of the prolific reservoirs in South Sumatra Basin. Seismic inversion and well analysis are done to determine the coal seam characteristics of Muara Enim Formation. This research uses three inversion methods, which are: model base hard- constrain, bandlimited, and sparse-spike inversion. Each type of seismic inversion has its own advantages to display the coal seam and its characteristic. Interpretation result from the analysis data shows that the Muara Enim coal seam has 20 (API) gamma ray value, 1 (gr/cc) – 1.4 (gr/cc) from density log, and low AI cutoff value range between 5000-6400 (m/s)*(g/cc). The distribution of coal seam is laterally thinning northwest to southeast. Coal seam is seen biasedly on model base hard constraint inversion and discontinued on band-limited inversion which isn't similar to the geological model. The appropriate AI inversion is sparse spike inversion which has 0.884757 value from cross plot inversion as the best correlation value among the chosen inversion methods. Sparse Spike inversion its self-has high amplitude as a proper tool to identify coal seam continuity which commonly appears as a thin layer. Cross-sectional sparse spike inversion shows that there are possible new boreholes in CDP 3662-3722, CDP 3586-3622, and CDP 4004-4148 which is seen in seismic data as a thick coal seam.

Java is the most populous island in Indonesia with 50 million people live there. This island geologically formed at the Eurasia plate margin by the subduction of the Australian oceanic crust. At the south part of Java, beside the occurrence of 2-plate convergence earthquake (interplate), there are also the activities of the intraplate earthquake. Research for distinguish this 2 different earthquake type is necessary for estimating the behavior of the earthquake that may occur. The aim of this research is to map the b-value in the south of Java using earthquake data from 1963 until 2008. The research area are divided into clusters based on the epicenter mapping results with magnitude more than 4 and three different depth (0-30 km, 30-60 km, 60-100 km). This location clustering indicate group of earthquakes occurred by the same structure or mechanism. On some cluster in the south of Java, b-value obtained are between 0.8 and 1.25. This range of b-value indicates the region was intraplate earthquake zone, with 0.72-1.2 b-value range is the indication of intraplate earthquake zone. The final validation is to determine the mechanism of a segment done by correlating the epicenter and b-value plot with the available structural geology data. Based on this research, we discover that the earthquakes occur in Java not only the interplate earthquake, the intraplate earthquake also occurred here. By identifying the mechanism of a segment in the south of Java, earthquake characterization that may occur can be done for developing the accurate earthquake disaster mitigation system.

Well logging data provide many geological information and its trends resemble nonlinear or non-stationary signals. As long well log data recorded, there will be external factors can interfere or influence its signal resolution. A sensitive signal analysis is required to improve the accuracy of logging interpretation which it becomes an important thing to determine sequence stratigraphy. Complete Ensemble Empirical Mode Decomposition (CEEMD) is one of nonlinear and non-stationary signal analysis method which decomposes complex signal into a series of intrinsic mode function (IMF). Gamma Ray and Spontaneous Potential well log parameters decomposed into IMF-1 up to IMF-10 and each of its combination and correlation makes physical meaning identification. It identifies the stratigraphy and cycle sequence and provides an effective signal treatment method for sequence interface. This method was applied to BRK- 30 and BRK-13 well logging data. The result shows that the combination of IMF-5, IMF-6, and IMF-7 pattern represent short-term and middle-term while IMF-9 and IMF-10 represent the long-term sedimentation which describe distal front and delta front facies, and inter-distributary mouth bar facies, respectively. Thus, CEEMD clearly can determine the different sedimentary layer interface and better identification of the cycle of stratigraphic base level.

Upper Red Bed, Menggala Formation, Bangko Formation, Bekasap Formation and Duri Formationare considered as the major reservoirs in Central Sumatra Basin (CSB). However, Telisa Formation which is well-known as seal within CSB also has potential as reservoir rock. Field study discovered that lenses and layers which has low to high permeability sandstone enclosed inside low permeability shale of Telisa Formation. This matter is very distinctive and giving a new perspective and information related to the invention of hydrocarbon potential in reservoir sandstone that isolated inside low permeability shale. This study has been conducted by integrating seismic data, well logs, and petrophysical data throughly. Facies and static model are constructed to estimate hydrocarbon potential resource. Facies model shows that Telisa Formation was deposited in deltaic system while the potential reservoir was deposited in distributary mouth bar sandstone but would be discontinued bedding among shale mud-flat. Besides, well log data shows crossover between RHOB and NPHI, indicated that distributary mouth bar sandstone is potentially saturated by hydrocarbon. Target area has permeability ranging from 0.01-1000 mD, whereas porosity varies from 1-30% and water saturation varies from 30-70%. The hydrocarbon resource calculation approximates 36.723 MSTB.

In between 5 – 10 years, in our country, passive seismic method became more popular to finding hydrocarbon. Low price, nondestructive acquisition and easy to mobilization is the best reason for choose the method. But in the other part, some people are pessimistically to deal with the result. Instrument specification, data condition and processing methods is several points which influence characteristic and interpretation passive seismic result. In 2010 one prospect in East Java Basin has been measurement constist of 112 objective points and several calibration points. Data measurement results indicate a positive response. Furthermore, in 2013 exploration drliing conducted on the prospect. Drill steam test showes 22 MMCFD in objective zone, upper – late oligocene. In 2015, remeasurement taken in objective area and show consistent responses with previous measurement. Passive seismic is unique method, sometimes will have difference results on dry, gas and oil area, in field production and also temporary suspend area with hidrocarbon content.