Hasil Pencarian  ::  Simpan CSV :: Kembali

"Lebong, Bengkulu merupakan salah satu zona prospek panas bumi yang dilalui oleh sesar aktif, yaitu Sesar Sumatera, yang berperan penting dalam mengontrol sistem panas bumi di wilayah tersebut. Informasi terkait distribusi densitas batuan bawah permukaan di zona ini masih terbatas, terlepas dari pentingnya karakteristik densitas dalam mengidentifikasi jalur migrasi fluida panas bumi dan zona reservoir. Penelitian ini bertujuan untuk mengidentifikasi distribusi densitas bawah permukaan di daerah prospek panas bumi Lebong melalui pemodelan inversi 3D data gravitasi. Pengukuran gravitasi dilakukan dengan jarak antar titik 250 meter pada area seluas 12 × 10 km. Data yang diperoleh dikoreksi menggunakan koreksi udara bebas, Bouguer (dengan densitas rata-rata 2,5 gr/cc berdasarkan metode Parasnis dan Nettleton), serta koreksi medan, hingga dihasilkan peta anomali Bouguer lengkap. Pemodelan dilakukan dengan pendekatan unconstrained dan constrained inversion yang memanfaatkan pendekatan Cartesian Cut Cell (CCC) dan algoritma Iterative Reweighting Inversion (IRI) Interpretasi hasil pemodelan dikorelasikan dengan data magnetotellurik (MT) dan peta geologi. Hasil analisis menunjukkan adanya zona densitas rendah yang memanjang searah struktur sesar dan terpusat di area graben, yang diinterpretasikan sebagai jalur utama migrasi fluida panas bumi pada zona prospek Lebong.

---

Lebong, Bengkulu is one of the geothermal prospect zones traversed by an active fault, the Sumatra Fault, which plays a significant role in controlling the geothermal system in the area. However, information related to the subsurface rock density distribution in this zone remains limited, despite the importance of density characteristics in identifying fluid migration pathways and reservoir zones. This study aims to identify the subsurface density distribution in the Lebong geothermal prospect area through 3D gravity inversion modeling. Gravity measurements were conducted with a 250-meter station spacing across an area of 12 × 10 km. The acquired data were corrected using free-air correction, Bouguer correction (with an average density of 2.5 g/cc based on the Parasnis and Nettleton methods), and terrain correction, resulting in a complete Bouguer anomaly (CBA) map. Inversion modeling was performed using both unconstrained and constrained, which utilizes the Cartesian Cut Cell (CCC) method and the Iterative Reweighting Inversion (IRI) algorithm). Interpretation of the modeling results was supported by magnetotelluric (MT) and geological data. The analysis revealed a low-density zone trending parallel to the fault structure and centered in the graben area, which is interpreted as the main fluid migration pathway in the Lebong geothermal prospect zone."

"Penelitian mengenai tekanan pori dilakukan pada Formasi Talang Akar dan Baturaja, Sub-Cekungan Jambi, Cekungan Sumatera Selatan. Cakupan penelitian berfokuskan pada wilayah yang mengalami overpressure atau tekanan pori melebihi normal. Analisis mengenai tekanan pori dilakukan sebagai upaya mengoptimalkan proses pengeboran di wilayah yang akan dilakukan pengembangan sumur. Pengolahan data terbagi menjadi dua, yaitu pengolahan tekanan pori pada data sumur dan penyebarannya menggunakan data seismik. Dalam melakukan estimasi nilai tekanan pori pada sumur, digunakan metode dasar tekanan pori, yaitu metode Eaton. Selanjutnya, dilakukan penyebaran tekanan pori pada data seismik menggunakan neural network dengan data masukan berupa kecepatan gelombang P dan inversi impedansi elastik. Penggunaan data tersebut dipilih karena cukup efektif dalam melakukan estimasi persebaran tekanan pori yang dikontrol oleh keberadaan litologi dan fluida. Selain itu, dilakukan pemodelan substitusi fluida guna menguatkan analisis persebaran gas pasiran.Dari pengolahan data didapatkan bahwa terdapat zona overpressure di Formasi Talang Akar dengan kedalaman 6030.184 ft - 6368.5 ft. Sedangkan pada Formasi Batu Raja tidak terlihat anomali tekanan pori yang signifikan. Keberadaan fluida hidrokarbon berupa gas pasiran di antara litologi shale dan keberadaan struktur patahan diduga menjadi penyebab terjadinya anomali berupa overpressure.

---

Research on pore pressure was carried out in the Talang Akar and the Batu raja Formation, Jambi Sub-Basin, South Sumatra Basin. The scope of research focuses on areas that experience overpressure or the pore pressure exceeds normal. Analysis of pore pressure is carried out as an effort to optimize the drilling process in areas where well development will be carried out. The data processing is divided into two, namely the processing of pore pressure in the well data and its distribution using seismic data. In estimating the value of pore pressure in drilling wells, the basic method for calculating pore pressure is used, namely the Eaton method. Next, the distribution of pore pressure estimates on seismic data is carried out using a neural network with input data in the form of P wave velocity and elastic impedance inversion. The use of input data was chosen because it is quite effective in estimating the pore pressure distribution which is controlled by the presence of lithology and fluids. In addition, fluid substitution modeling was carried out to strengthen the analysis of the distribution of sandy gas.

From the data processing carried out, it was found that there are zones overpressure in the Talang Akar Formation with a depth of 6030.184 ft to 6368.5 ft. Whereas in the Batu Raja Formation there is no significant pore pressure anomaly. The presence of hydrocarbon fluids in the form of sandy gas between lithology shale and the existence of a fault structure is thought to be the cause of the anomaly in the form of overpressure."

"Cekungan Banggai memiliki potensi hidrokarbon yang terakumulasi pada reservoir karbonat dengan porositas ganda – porositas primer dan sekunder akibat diagenesis. Hal tersebut menyebabkan variasi nilai porositas secara vertikal maupun lateral sehingga menyulitkan identifikasi zona hidrokarbon yang potensial. Untuk mengidentifikasi zona reservoir karbonat dengan porositas ganda, penelitian ini mengintegrasikan metode inversi seismik, multiatribut seismik, dan Probabilistic Neural Network (PNN), yang bertujuan untuk meningkatkan akurasi dalam pemodelan karakteristik reservoir. Hasil penelitian menunjukkan bahwa karbonat build-up Formasi Mantawa memiliki impedansi akustik rendah (6200–10000 m/s*g/cc), porositas tinggi (21%–28%), dan saturasi air rendah (30%–47%), yang mengindikasikan adanya batuan berpori yang mengandung hidrokarbon. Diharapkan dengan penerapan integrasi metode tersebut, distribusi reservoir, porositas, dan saturasi air dapat dipetakan dengan lebih akurat berdasarkan hasil dari prediksi Probabilistic Neural Network (PNN) dan membuka peluang eksplorasi lebih lanjut di zona prospektif Lapangan ZENA.

---

The Banggai Basin has hydrocarbon potential accumulated in carbonate reservoirs with dual porosity – primary and secondary porosity due to diagenesis. This causes variations in porosity values both vertically and laterally, making it difficult to identify potential hydrocarbon zones. To identify carbonate reservoir zones with dual porosity, this study integrates seismic inversion methods, multi-attribute seismic analysis, and Probabilistic Neural Network (PNN), aimed at improving the accuracy of reservoir characteristic modeling. The results show that the carbonate build-up of the Mantawa Formation has low acoustic impedance (6200–10000 m/s*g/cc), high porosity (21%–28%), and low water saturation (30%–47%), indicating the presence of porous rocks containing hydrocarbons. It is expected that with the application of the integration of these methods, the distribution of the reservoir, porosity, and water saturation can be mapped more accurately based on the results from the predictions of the Probabilistic Neural Network (PNN), thus opening up further exploration opportunities in the prospective zones of the ZENA Field."

"Cekungan Banggai memiliki potensi hidrokarbon yang terakumulasi pada reservoir karbonat dengan porositas ganda – porositas primer dan sekunder akibat diagenesis. Hal tersebut menyebabkan variasi nilai porositas secara vertikal maupun lateral sehingga menyulitkan identifikasi zona hidrokarbon yang potensial. Untuk mengidentifikasi zona reservoir karbonat dengan porositas ganda, penelitian ini mengintegrasikan metode inversi seismik, multiatribut seismik, dan Probabilistic Neural Network (PNN), yang bertujuan untuk meningkatkan akurasi dalam pemodelan karakteristik reservoir. Hasil penelitian menunjukkan bahwa karbonat build-up Formasi Mantawa memiliki impedansi akustik rendah (6200–10000 m/s*g/cc), porositas tinggi (21%–28%), dan saturasi air rendah (30%–47%), yang mengindikasikan adanya batuan berpori yang mengandung hidrokarbon. Diharapkan dengan penerapan integrasi metode tersebut, distribusi reservoir, porositas, dan saturasi air dapat dipetakan dengan lebih akurat berdasarkan hasil dari prediksi Probabilistic Neural Network (PNN) dan membuka peluang eksplorasi lebih lanjut di zona prospektif Lapangan ZENA.

---

The Banggai Basin has hydrocarbon potential accumulated in carbonate reservoirs with dual porosity – primary and secondary porosity due to diagenesis. This causes variations in porosity values both vertically and laterally, making it difficult to identify potential hydrocarbon zones. To identify carbonate reservoir zones with dual porosity, this study integrates seismic inversion methods, multi-attribute seismic analysis, and Probabilistic Neural Network (PNN), aimed at improving the accuracy of reservoir characteristic modeling. The results show that the carbonate build-up of the Mantawa Formation has low acoustic impedance (6200–10000 m/s*g/cc), high porosity (21%–28%), and low water saturation (30%–47%), indicating the presence of porous rocks containing hydrocarbons. It is expected that with the application of the integration of these methods, the distribution of the reservoir, porosity, and water saturation can be mapped more accurately based on the results from the predictions of the Probabilistic Neural Network (PNN), thus opening up further exploration opportunities in the prospective zones of the ZENA Field."

"Telah dilakukan penelitian guna delineasi zona prospek sistem panasbumi daerah ldquo;Z rdquo; menggunakan permodelan tiga Dimensi magnrtotellurik didukung data terpadu berupa geologi dan geokimia serta terintegrasi data gravitasi. Daerah panasbumi ldquo;Z rdquo; dalam tatanan tektoniknya termasuk pada jalur backarc Sumatera, tepat pada salah satu segmen sesar Sumatera bagian selatan, disusun oleh batuan vulkanik dan sedimen klastik yang berumur Tersier hingga Kuarter Andesit-Basalt . Gejala adanya sistem panasbumi pada daerah penelitian ditandai dengan kemunculan manifestasi permukaan berupa alterasi dan lima mata air panas bersuhu 44,4 - 92,5 oC, pH 8,19 - 9,43 dan bertipe bikarbonat, sulfat-bikarbonat, serta sulfat-klorida. Pembentukan sistem panasbumi dipengaruhi oleh aktivitas tektonik menyerong oblique antara lempeng Samudera India dan Lempeng Kontinen Eurasia searah dengan pola sesar Sumatera. Berdasarkan analisis air panasbumi temperatur reservoir diambil melalui perhitungan geothermomether SiO2 Fournier 1977 , Na-K Giggenbach 1988 , Na-K-Ca, diagram Na-K-Mg serta diagram Enthalphy - Cloride Mixing Model berkisar 145 - 155oC, termasuk dalam sistem panas bumi bertemperatur sedang. Berdasarkan inversi tiga dimensi data MT didapatkan kedalaman Top of Reservoar TOR sistem panasbumi daerah ldquo;Z rdquo; sekitar 400 m elevasi 50 mdpl sedangkan berdasarkan forward modeling data gravitasi lintasan 2 dimensi diperkirakan sumber panas berupa cooling instrusion diperkirakan batuan gabro ; resistivitas ge; 450 ?m ; densitas 2,95 - 3,15 gr/cc dan reservoar berupa batupasir resistivitas 50 - 250 ?m ; densitas 2,60 gr/cc . Sistem panasbumi daerah penelitian termasuk jenis tektonik fracture zone dengan temperatur sedang dengan luas daerah prospek sekitar 7,5 km2.

---

A study for delineating geothermal system of prospect area ldquo Z rdquo has been done by using tree dimension modeling of magnetotelluric supported unified data just like geological and goechemical and integrated gravity data. Geothermal area ldquo Z rdquo in tectonic setting included in Sumatra volcanic backarc, right on one of the southern part of Sumatra fault segment. Compodes by volcanic and clastic sendimentary rock are Tertiary to Quarternary Andesite Basalt. The existance of goethermal system in this area is indicated by the presence of thermal manifestation in form of alteration and five hot springs temperature in the ranges 44.4 ndash 92.5 oC, and pH 8.19 ndash 9.43 and type of fluida are bicarbonate, sulphate bicarbonate, and sulfate chloride. The development of geothermal system is affected by tectonic oblique between the Indian Ocean plate and the Eurasian Contenent Plate direction of the Sumatra fault patterns.

Based on the analysis of geothermal water reservoir temperature are taken through the calculation geothermometer SiO2 Fournier 1977, Na K Giggenbach 1988 , Na K Ca, Na K Mg diagram and Enthalpi Mixing Cloride Model range 145 ndash 155 oC, classified as intermediate temperature. Base on a three dimensional inversion of the magnetotelluric data obtained depth Top of Reservoir TOR geothermal system area ldquo Z rdquo about 400 m elevation 50 meters above sea leavel , while based on the two dimensional of the gravity data predicted heat sources such as cooling instrusion estimated gabbro density 2,95 ndash 3,15 gr cc and reservoar such as sandstone resistivity 50 ndash 250 m density 2,60 gr cc . The Geothermal systems of research area classified as the type of intermediate temperature tectonic fracture zone with prospect area about 7,5 km2."

"

Daerah X merupakan salah satu daerah yang memliki potensi panas bumi di Jawa Barat ditandai dengan keberadaan manifestasi mata air panas dan hangat di permukaan, serta berada di kawasan vulkanik. Penelitian ini bertujuan untuk mengidentifikasi basemen di daerah panas bumi X dengan menggunakan metode inversi tiga dimensi data gravitasi dan hasilnya diintegrasikan dengan informasi geologi, manifestasi, hasil analisis First Horizontal Derivative (FHD). Complete Bouguer Anomaly (CBA) digunakan sebagai data input untuk inversi tiga dimensi. Hasil analisis CBA menunjukkan kontras anomali tinggi dan rendah berarah tenggara – barat laut, yang mana sesuai dengan strike struktur yang berada di kawasan tersebut, dengan anomali tinggi berada di selatan (diduga berasal dari batuan berumur Miosen - Pliosen) dan anomali rendah dominan di utara (diduga berasal dari batuan vulkanik muda yang belum terlitifikasi sempurna). Peta kedalaman basemen diperoleh dengan mengekstraksi nilai densitas  dari hasil inversi tiga dimensi CBA. Pemilihan nilai densitas ini didasarkan pada estimasi rata-rata densitas batuan dari Formasi Bentang yang diasosiasikan sebagai basemen di daerah panas bumi X. Berdasarkan analisis peta kedalaman basemen, pengangkatan basemen dominan berada di kawasan Gunung X (zona prospek) diduga akibat keberadaan heat source yang mendorong basemen dengan rentang kedalaman 1000 hingga 4500 meter. Basemen dengan kedalaman yang dalam dominan berada di timur dari pengangkatan basemen Gunung X, diduga akibat batuan vulkanik muda yang belum terlitifikasi sempurna cukup tebal di permukaan. Manifestasi panas bumi berada di tepi pengangkatan basemen yang juga berkorelasi dengan struktur geologi, serta hasil analisis kelurusan FHD. Batas tepi pengangkatan basemen dan struktur saling berasosiasi dengan kontras densitas dan diduga membentuk zona permeabel sebagai pathway fluida panas bumi.

 

---

Area X is one of the geothermal prospect regions in West Java, as indicated by the presence of hot and warm spring manifestations at the surface and its location within a volcanic zone. This study aims to identify the basement structure in the X geothermal area using three-dimensional gravity data inversion, integrated with geological information, surface manifestations, and the results of First Horizontal Derivative (FHD) analysis. Complete Bouguer Anomaly (CBA) data was used as the input for the 3D inversion. The CBA analysis revealed a southeast-northwest trending contrast between high and low gravity anomalies, which aligns with the structural strike direction in the region. High anomalies are located in the southern part, presumably associated with Miocene–Pliocene-aged rocks, while low anomalies are dominant in the north, likely caused by the presence of relatively thick, unconsolidated young volcanic deposits. The basement depth map was obtained by extracting the 2.6 g/cm³ isosurface from the 3D CBA inversion results. This density value was selected based on the estimated average density of the Bentang Formation, which is interpreted as the basement in the geothermal area of X. The basement uplift is primarily observed in the Mount X area (the prospect zone), presumably due to the presence of a heat source that drives the basement upwards, with depths ranging from approximately 1000 to 4500 meters. Deeper basement is observed to the east of the uplift, likely due to the thickness of unconsolidated young volcanic rocks in that area. Geothermal surface manifestations are found along the margin of the basement uplift, which also correlates with geological structures and lineaments identified from the FHD analysis. The boundary of the basement uplift and associated structures correspond with density contrasts and are interpreted to form a permeable zone that may act as a pathway for geothermal fluid migration.

 

"

"ABSTRAKTujuan utama dari eksplorasi adalah penentuan lokasi pemboran. Kriteria kesuksesan target pemboran adalah area yang memiliki temperatur dan permeabilitas yang tinggi. Temperatur berasosiasi dengan keberadaan sumber panas dan jumlah energi termal yang tersimpan, sedangkan permeabilitas berhubungan dengan keberadaan struktur geologi baik patahan maupun kekar yang terisi fluida yag dapat menjadi media perpindahan energi panas. Identifikasi struktur geologi harus dilakukan melalui kombinasi pengamatan struktur geologi melalui citra penginderaan jauh remote sensing maupun pemetaan geologi sebagai metode identifikasi struktur geologi di permukaan dan interpretasi hasil survei geofisika sebagai metode identifikasi kemenerusan struktur geologi di bawah permukaan. Analisis kerapatan kelurusan lineaments density analysis adalah salah satu metode yang dapat digunakan untuk mengidentifikasi kerapatan struktur geologi di permukaan melalui citra satelit. Area dengan kerapatan struktur atau patahan yang tinggi diasumsikan memiliki permebilitas yang tinggi. Untuk konfirmasi kemenerusan struktur geologi di bawah permukaan digunakan analisis data survei gravitasi. Dari deliniasi area yang memiliki permeabilitas tinggi dan didukung data hasil pengukuran MT dan geokimia dapat dilakukan rekonstruksi model konseptual daerah penelitian

---

ABSTRACTThe primary objective of exploration is well targeting determination. Success criteria on well targeting is high temperature and permeabilities. Temperture is associated with the heat source existence and the amount of thermal energy stored within, while permeability is associated with the presence of the geological structure wether fault or joint which have fluid in it as heat energy tranfer medium. Identification of geological structures must be made through a combination of surface observation and subsurface interpretation. Surface observation is carried out through remote sensing imagery interpretation and geological mapping while subsurface continuities is obtained from interpretation of geophysical survey data. Lineaments density analysis is one of the method that can be used to identify the density of geological structures on the surface through satellite imagery. Area where contained high structure density on the surface and confirmed by the result of gravity data interpretation is assumed to have a high permeabilities. Geothermal conceptual model could be defined by deliniation of high permeability area supported by geochemical and magnetotelluric data."

"Lapangan geotermal X berada di area gunung A yangmana berdasarkan data geologi ditemukan adanya manifestasi berupa hot spring dan fumarole. Pengukuran MT dilakukan untuk mengetahui persebaran resistivity batuan di bawah permukaan. Pengolahan data MT dilakukan dari analisis time series dan filtering noise kemudian dilakukan Transformasi Fourier dan Robust Processing. Setelah itu baru dilakukan crosspower untuk menyeleksi data sehingga output dari proses ini berupa kurva MT. Setelah didapatkan kurva MT dilakukan koreksi statik dikarenakan kurva TE dan TM terjadi shifting. Untuk proses akhirnya baru dilakukan inversi 2D dan inversi 3D. setelah itu dilakukan perbandingan antara 2D dan 3D. Wilayah interest lapangan X berada di lintasan AA dan lintasan AB. Berdasarkan analisis 3D diidentifikasi bahwa zona alterasi menipis di wilayah upflow dan menebal ke arah outflow yangmana sesuai dengan teori. Wilayah upflow dapat diketahui dengan melihat manifestasi berupa fumarole.

---

The geothermal field X is located in the area of Mount A which based on geological data found the presence of hot spring and fumarole manifestations. MT measurements were carried out to determine the distribution of rock resistivity in the subsurface. MT data processing is starts from time series analysis and noise filtering then Fourier Transform and Robust Processing are performed. After that, crosspower is done to select data so that the output of this process is an MT curve. After got the MT curve then a static correction is done because the TE and TM curves are shifting. For the final process are 2D inversion and 3D inversion. After that make a comparison between 2D and 3D. The area of interest in field X is on the line AA and line AB. Based on the 3D analysis, it was identified that alteration zones thinned in the upflow region and thickened towards the outflow which is make sense with the theory."

"Daerah Prospek panasbumi "B" terletak di Pesawaran, Kabupaten Lampung Selatan, Lampung. Dari data remote sensing diketahui bahwa arah utama dari kelurusan-kelurusan pada daerah panasbumi prospek "B" adalah Baratlaut-Tenggara yang sesuai dengan pola struktur geologi utama dan berhubungan dengan kehadiran manifestasi permukaan. Dari data geokimia diketahui bahwa zona outflow prospek panasbumi "B" berada pada daerah manifestasi mata air panas dan dari plotting ternary diagram Na-K-Mg menunjukkan temperatur reservoar sebesar 220 C. Analisis geofisika dari data gravitasi sebanyak 163 titik pengukuran dan dari data Magnetotellurik sebanyak 58 titik pengukuran menunjukan bahwa lapisan clay cap dengan densitas 2.2 gr/cc memiliki nilai resistivitas sebesar.

---

Area prospect of B geothermal area is located in Pesawaran, South Lampung District, Lampung. From remote sensing data is known that the main direction of the lineaments in the area of geothermal prospect B is Northwest Southeast in accordance with the pattern of major geological structures and associated with the presence of surface manifestations. From the geochemical data known that the prospects for geothermal outflow zone B in the region of hot springs and the manifestation of plotting Ternary Diagram Na K Mg shows a reservoir temperature of 220 C. Geophysical analysis from gravity 163 data and magnetotelluric 58 data measuring point indicate that the clay cap layer with a density of 2.2 g cc and resistivity of "

"ABSTRAKPenelitian yang dilakukan di cekungan Sumatera Utara bertujuan untuk mengetahhui bagaimana kondisi batuan sedimen pra-tersier di bawah permukaan. Hal ini dilakukan karena adanya isu dimana batuan sedimen yang telah berumur pra-tersier diisukan belum semua nya termertamorfkan. Hal ini menjadi menarik karena ada kemungkinan cadangan minyak dan gas bumi yang tersimpan pada batuan sedimen pra-tersier. Metode gravitasi merupakan metode geofisika yang digunakan untuk memetakan struktur bawah permukaan. Karena nilai yang terukur oleh metode ini adalah nilai anomali gravitasi dari batuan yang ada di bawah permukaan. Metode inversi menggunakan persamaan Fast Fourier Transform pada aplikasi Matlab yang mana bertujuan untuk dapat membuat pemodelan dari bawah permukaan. Data masukan yang digunakan pada pemograman Matlab berupa nilai data anomali graitasi, jumlah kolom dan baris atau koordinat titik pengambilan data, nilai kedalaman rata-rata batuan atau lapisan targer, nilai kontras densitas batuan dan nilai kriteria konvergensi. Hasil pemograman menampilkan horizon bawah permukaan pada nilai kedalaman dan horizon hasil perhitungan nilai kedalaman dengan nilai anomali gravitasi. Proses inversi dilakukan pada saat melakukan proses iterasi dan akan berhenti jika nilai kesalahan RMS yang didapat dari proses iterasi telah mendekati nilai yang ditentukan sebagai kriteria konvergensi atau sampai iterasi maksimum tercapai. Hasil iterasi berhenti saat nilai RMS error sekitar 0.00016879 % pada iterasi kedua. Horizon yang terbentuk sesuai dengan keadaan geologi yang terbentuk di bawah permukaan dimana batuan sedimen pra-tersier berada di kedalaman rata-rata 5.8 km di bawah permukaan dengan ketebalan dari 5 km hingga 6 km.

---

ABSTRACTThe research conducted in the North Sumatra basin aims to find out how the conditions of pre-tertiary sedimentary rocks below the surface. This is done because there are issues where sedimentary rocks that have been pre-tertiary are not yet included. This is interesting because it is possible that oil and gas reserves are stored in pre-tertiary sedimentary rocks. The gravity method is a geophysical method used to map the subsurface structure. Because the value measured by this method is the value of the gravitational anomaly of the rock below the surface. The inversion method uses the Fast Fourier Transform equation in the Matlab application, which aims to be able to make modeling from the subsurface. Input data used in Matlab programming is in the form of gravity anomalous data values, column number and row or coordinates of data retrieval points, average rock depth or targer layer values, rock density contrast values and convergence criterion values. The programming results show the subsurface horizon at the depth and horizon values resulting from the calculation of the depth value with the value of gravity anomalies. The inversion process is carried out when carrying out the iteration process and will stop if the RMS error value obtained from the iteration process has approached the value specified as the convergence criterion or until the maximum iteration is reached. The iteration results stop when the RMS error value is around 0.00016879% in the second iteration. Horizon is formed in accordance with the geological conditions formed below the surface where pre-tertiary sedimentary rocks are at an average depth of 5.8 km below the surface with a thickness of 5 km to 6 km."