ANALISA JARINGAN PIPA LOOP-NODE DUA FASA MENGGUNAKAN METODE BEGGS AND BRILL

ANALISA JARINGAN PIPA LOOP-NODE DUA FASA MENGGUNAKAN METODE BEGGS AND BRILL Rudi Rubiandini R.S. - Insitut Teknologi Bandung Harisza Koswara – Stavanger University, Norway [email protected] PENDAHULUAN Fluida reservoir dari kepala sumur dialirkan ke proses pemisahan dan dari proses pemisahan ini fluida reservoir dari berbagai sumur dialirkan ke stasion pengumpul. Fluida reservoir tersebut dialirakan melalui pipa salur yang diletakkan di atas tanah atau ditanam di dalam tanah dan mengikuti topografi permukaan tanah. Selain itu pipa juga ditekuk untuk mengikuti belokan, disambung dengan pipa-pipa yang berbeda diameternya.

RINGKASAN Perancangan suatu jaringan pipa dibatasi oleh kebutuhan konsumen, kemampuan supply, dan bentuk jaringan. Masalah akan menjadi bertambah jika kita mencoba untuk memodifikasi suatu jaringan yang telah ada sebelumnya, karena ada batasan lain yaitu dimensi pipa lama yang sudah terpasang sebelumnya. Selain itu, kesetimbangan aliran fluida pada tiap cabang pipa dapat berubah akibat dari modifikasi tersebut. Penggunaan suatu perangkat lunak sangat diperlukan dalam menganalisa jaringan pipa untuk memberikan hasil yang lebih baik dan akurat berdasarkan pada suatu algoritma sederhana. Algoritma yang digunakan dalam analisa adalah yang bersifat step-forward. Harga tekanan pada setiap node langsung dikoreksi pada setiap iterasinya berdasarkan ketidakseimbangan laju alir massa pada node tersebut karena harga tekanan yang digunakan untuk menentukan laju alir massa dalam fluida dalam pipa pada tahap iterasi sebelumnya belum tepat. Kecepatan konvergensi sangat bergantung pada faktor akselerasi yang digunakan untuk mengkoreksi harga tekanan pada node. Algoritma ini bisa diaplikasikan pada berbagai macam bentuk jaringan pipa, baik itu berbentuk loop, non-loop, maupun parallel dan dapat diaplikasikan untuk berbagai metode perhitungan laju alir. Pada paper ini, algoritma di aplikasikan untuk jaringan pipa dua fasa, dimana perhitungan laju alirnya menggunakan metode Beggs and Brill.

Aliran Fluida Dua Fasa dalam Pipa Aliran dua fasa dalam pipa, selanjutnya disingkat aliran fluida saja, merupakan fungsi dari perbedaan tekanan inlet (upstream) dan tekanan outlet (downstream), faktor gesekan, temperatur, serta parameter-parameter sifat fisik fluida dua fasa (gas dan minyak) dan pipa itu sendiri. Persamaan aliran fluida dalam pipa dengan sendirinya merupakan persamaan yang menghubungkan antara laju alir gas dengan parameter-parameter yang disebutkan di atas. Prosedur yang dikembangkan oleh Beggs and Brill adalah untuk menghitung tekanan outlet (Poutlet) dari tekanan inlet (Pinlet) dan laju alir yang telah diketahui sebelumnya. Untuk perhitungan laju alir dengan tekanan outlet (Poutlet) dan tekanan inlet (Pinlet) diketahui.. Gambar 1 adalah diagram alir perhitungan laju alir dengan metode Beggs and Brill 2). Persamaan Beggs and Brill digunakan untuk menghitung beda tekanan (pressure drop) dalam pipa sebagai berikut:

Kata kunci : jaringan pipa, parallel, loop, non-loop, Beggs dan Brill, algoritma sederhana, simulator, node.

1

f tp Gm v m  g ∆z  ρ tp sin θ + 2g c D gc ∆p =  ρ tp v m v sg 1− gc p

Algoritma Sederhana untuk Analisa Nodal Jaringan Pipa Dua Fasa Algoritma Sederhana yang dikembangkan oleh Tian-Adewumi (1993) adalah salah satu metode pemecahan masalah jaringan yang tergolong analisa nodal, yaitu mendasarkan pemecahan atas hukum Kirchoff I. Persamaan Dasar Algoritma Seperti analisa nodal lainnya, algoritma ini memerlukan suatu perkiraan awal atau inisialisasi harga tekanan pada node-node yang tekanannya tak diketahui. Dengan harga tekanan perkiraan yang diberikan, laju alir massa gas dalam semua cabang yang berhubungan dengan suatu node (katakanlah node ke- i) dihitung dengan persamaan aliran gas dalam pipa, persamaan (2) atau (4). Kemudian error yang terjadi pada node ke- i tersebut (errori) dihitung dengan persamaan berikut:

   .........(1)

Model dan Sistem Jaringan a. Jaringan Pipa Sederhana Jaringan pipa yang yang sederhana hanya mempunyai beberapa node yang dihubungkan oleh sebuah pipa ke node yang lainnya. Analisis jaringan pipa yang sederhana ini juga dapat dikembangkan untuk jaringan yang lebih luas atau kompleks. Jaringan pipa sederhana ini terdiri dari sistem loop dan sistem biasa. Gambar 2 merupakan jaringan pipa sederhana sistem loop sedangkan Gambar 3 memperlihatkan jaringan pipa biasa. Dengan demikian diperlukan dua model penyelesaian dari problem jaringan tersebut. Dalam pembentukkan model untuk sistem jaringan loop dan biasa menggunakan asumsi sebagai berikut : 1. Temperatur fluida selalu konstan disepanjang sistem. 2. Kondisi aliran adalah steady state. 3. Komposisi fluida seragam disepanjang sistem. 4. Tidak terdapat valve, peralatan lain dan belokkan disepanjang sistem. 5. Tidak ada penurunan tekanan pada sambungan pipa 6. Tidak ada penurunan tekanan pada belokan pipa

Li

errori = ∑ qij + Qi .......................(2) j =1

dimana : li = jumlah total cabang yang terhubung dengan node ke – i N = jumah total node pada jaringan qij = jumlah total rate fluida yang mengalir dari cabang - j ke node – i. Qi = jumlah supply atau demand pada node ke -i Error yang dihitung pada node pada dasarnya merupakan ketidakseimbangan laju alir massa yang terjadi pada node tersebut karena harga tekanan yang digunakan untuk menghitung laju alir massa gas dalam pipa-pipa yang berhubungan dengan node tersebut belum tepat. Hukum Kirchoff I akan terpenuhi jika error sama dengan nol. Harga tekanan node yang tidak tepat itu kemudian dikoreksi berdasarkan harga error yang terjadi, yaitu dengan menambahkan tekanan node tersebut dengan perkalian antara error yang terjadi pada node dengan suatu faktor akselerasi (acceleration factor) β, seperti yang diberikan oleh persamaan berikut:

b. Jaringan Pipa Kompleks Gambar 4 adalah jaringan pipa kompleks yang biasanya dijumpai dalam proses pengaliran fluida produksi dari masing-masing sumur menuju stasiun pengumpul dan sistem pengolahan dan juga distribusi menuju konsumen. Jaringan ini dapat disederhanakan seperti yang diperlihatkan pada gambar 5. Hal ini untuk memudahkan dalam menganalisis sistem jaringan yang kompleks.

2

Pi

( k +1)

= Pi

(k )

dimana: P1ij = Tekanan inlet pada cabang ke - ij P2ij = Tekanan outlet pada cabang ke -ij

+ β i × (errori ) ..............(3)

dimana: P(k+1) = Tekanan node ke-i setelah dikoreksi, = Tekanan node ke-i sebelum P(k) dikoreksi, β = Faktor akselerasi node ke - i.

Garis Pedoman untuk Analisa Nodal Jaringan Pipa Dua Fasa Untuk sistem jaringan pipa yang sederhana, ada dua variabel yaitu tekanan dan supply atau demand. Laju alir pada pipa tergantung tekanan dan diameter pipa. Jika parameter-parameter pipa diketahui, maka terdapat 2N variabel bebas pada sistem yang mempunyai N node. Dengan demikian dapat dibuktikan bahwa spesifikasi beberapa N variabel untuk menyelesaikan N variabel tersisa adalah sebuah kondisi yang diperlukan. Gambar 2 adalah contoh sistem yang tidak mempunyai solusi khusus. Contoh pertama pada Gambar 2 mempunyai 5 node dengan spesifikasi 5 variabel P2, P3, P4, Q2 dan Q3. Solusi untuk 5 variabel yang tidak diketahui P1, P2, Q1, Q4 dan Q5 adalah tidak tentu. Terdapat kombinasi yang tidak terbatas dari nilai-nilai variabel yang tidak diketahui yang memenuhi kesetimbangan massa pada node dan persamaan aliran pada tiap-tiap pipa. Contoh kedua pada Gambar 3 merupakan sistem sederhana tanpa loop dengan 5 node. Jika P1, P2, Q1, Q2, dan Q3 diketahui solusi untuk variabelvariabel P3, P4, P5, Q4 dan Q5 adalah tidak tentu, karena terdapat kombinasi yang tidak terbatas dari nilai-nilai variabel yang tidak diketahui yang memenuhi persamaan kesetimbangan massa pada node dan persamaan aliran pada tiap-tiap pipa. Dua contoh tersebut merupakan contoh problem yang tidak terdefinisi dengan baik. Meskipun demikian dapat dibuktikan bahwa N variabel yang terdistribusi sepanjang N node mempunyai solusi yang tidak pasti dalam sistem tersebut. Dengan kata lain, jika sebuah sistem memenuhi 2 hal yaitu 1. Hanya terdapat satu variabel yang tak diketahui pada setiap node, 2. Sekurang-kurang satu harga tekanan diketahui pada suatu node dalam jaringan.

Jika node ke-i merupakan node referensi dengan supply atau demand (Qi) tidak diketahui, maka supply atau demandnya yang dihitung, yaitu dengan persamaan berikut li

Qi = −∑ qij .....................................(4) j =1

Kemudian perhitungan dilanjutkan pada node berikutnya (i = i + 1) sampai semua node dihitung atau i = N. Setelah semua node dihitung (akhir iterasi ke- k), dihitung jumlah error total absolut pada semua nodenya. Jika error total absolut tersebut lebih besar dan error yang ditoleransi (ε) atau persamaan (5) belum terpernuhi, maka proses perhitungan diulangi lagi (iterasi ke- k + 1). Jika tidak berarti telah konvergensi telah tercapai. errori( k ) < ε i = 1,2,3,..., N ........(5) Faktor Akselerasi Prinsip dasar yang dipakai pada Algoritma Sederhana adalah mengoreksi harga tekanan pada suatu node sedemikian rupa sehingga harga tekanan yang telah dikoreksi dapat memberikan keseimbangan massa pada node tersebut. Jadi faktor akselerasi pada dasarnya merupakan suatu faktor konversi yang digunakan untuk mengubah ketidakseimbangan laju alir massa ke dalam harga penambahan (increment) tekanan yang diperlukan untuk tercapainya keseimbangan massa. Faktor akselerasi dapat dihitung dengan persamaan berikut: P1ij − P2ij β i = min ; j = 1,2,3,...., li ....(6) qij

3

Jika suatu sistem telah memenuhi syarat dia atas maka sistem ini mempunyai problem yang terdefinisikan dengan baik. Karena sederhananya petunjuk-petunjuk diatas akan membantu pengguna untuk menentukkan data minimum yang diperlukan untuk menganalisis sistem tertentu. Pada Gambar 4 memperlihatkan sistem jaringan pipa yang kompleks. Jika tekanan pada node 1 dan laju alir masuk atau keluar diketahui, maka ini merupakan petunjuk untuk dapat menyelesaikan sistem jaringan tersebut. Sistem jaringan ini dapat disederhanakan seperti pada Gambar 5. Penyederhanaan menjadi sistem yang terpisah-pisah menjadi sistem jaringan sederhana. Penyelesaian masingmasing sistem jaringan yang sederhana telah dibahas diatas. Penyelesaiaan problem tersebut adalah sama dengan atau tanpa penyederhanaan.

Jaringan Pipa pada kasus yang pertama ini merupakan jaringan pipa parallel dengan 2 titik node dan 2 cabang pipa, seperti terlihat pada gambar di bawah ini. Fluida masuk ke dalam jaringan melalui titik node 1 dan keluar melalui titik node 2 : Hasil perhitungan software didapatkan hasil seperti pada Tabel 1a dan 1b. Kasus 2 : Jaringan Non Loop dengan 8 Titk Node dan 7 Cabang Bentuk jaringan pipa pada kasus 2 berupa jaringan tanpa loop dengan 8 titik Node dan 7 cabang pipa seperti pada gambar dibawah ini. Node 1 sebagai node demand dengan harga rate sebesar 150000 BBL dan harga tekanan diketahui sebesar 1000 psia. Untuk node 2 dan 3 harga tekanan tidak diketahui dan tidak ada massa fluida yang masuk atau keluar dari node ini. Node 4,5,6,7, dan 8 merupakan node supply dengan harga rate masingmasing sebesar 30000 BBL.

Data Masukkan untuk Analisa Data masukan (Input data) yang diperlukan untuk analisa jaringan gas terbagi menjadi 1. Data sifat fisik fluida (SGgas, SGoil, Trata, Pb, GOR, µ, Rs, Bo) 2. Data node (data jumlah node, nomor node, P node (referensi), laju alir (supply atau demand) pada node-node tertentu) 3. Data cabang, (data jumlah cabang pada jaringan, nomor cabang, nomor node pada ujung-ujung setiap cabang, parameter-parameter pipa seperti panjang pipa, diameter dalam pipa dan sudut elevasi pipa).

Hasil perhitungan software didapatkan hasil seperti pada Tabel 2a dan 2b. Kasus 3 : Jaringan Loop dengan 5 cabang dan 5 titik Node Jaringan pipa pada kasus yang kedua ini, berbentuk loop dengan 5 titik node dan 5 cabang pipa. Fluida masuk ke dalam jaringan pipa melalui node 1 dan 2 dengan rate pada node 1 sebesar 60000 BBL dan pada node 2 sebesar 40000 BBL. sedangkan pada titik 3, 4, dan 5 fluida meninggalkan jaringan dengan rate pada node 3 sebesar 30000 BBL, node 4 40000 BBL dan node 5 sebesar 30000 BBL. Pada jaringan berbentuk loop ini, arah aliran fluida sebelumnya tidak diketahui. Harga tekanan hanya diketahui pada titik node 1, yaitu sebesar 1000 psia. Hasil perhitungan software didapatkan hasil seperti pada Tabel 3a dan 3b.

Diagram Alir Proses Analisa Jaringan Pipa Diagram alir (flow chart) proses analisa jaringan pipa dua fasa dengan Algoritma Sederhana yang telah dijelaskan sebelumnya dapat dilihat pada Gambar 6. STUDI KASUS DAN PEMBAHASAN Kasus 1 : Jaringan Pipa Parallel dengan 2 titik Node dan 2 Cabang

4

Kasus 4: Jaringan Loop Dengan Cabang Baru Antara Titik Node 1 dan 3 Bentuk jaringan pada kasus 4 mirip dengan pada kasus 3, tapi antara node 1 dan 3 ditambahkan satu buah cabang pipa dengan diameter yang sama dengan pipa lainnya, tapi memiliki panjang yang lebih, yaitu sebesar 8000 feet. Pada kasus ini, akan dilihat pengaruh penambahan cabang pipa pada jaringan berbentuk loop yang telah ada sebelumnya pada kasus 3

Kasus 8: Jaringan Loop Dengan Cabang Baru Antara Titik Node 3 dan 5 Pada kasus 8, ditambahkan cabang pipa baru antara node 3 dan node 5 dengan panjang 8000 ft, diameter 10 in, dan sudut elevasi pipa sebesar 00. Hasil perhitungan software didapatkan hasil seperti pada Tabel 8a dan 8b. KESIMPULAN 1. Algoritma sederhana ini dikembangkan untuk menganalisa dan mendesain jaringan pipa dua fasa gas dan minyak, algoritma ini tidak memerlukan matriks aljabar ataupun tebakan awal yang baik pada inisialisasinya. Lebih jauh lagi, tidak memerlukan penentuan arah aliran. 2. Algoritma yang dikembangkan ini, bisa menangani jaringan pipa dengan aliran satu fasa ataupun aliran dua fasa. Algoritma ini juga dapat diaplikasikan pada berbagai korelasi dua fasa. 3. Program komputer yang dibuat berdasarkan algoritma sederhana dapat digunakan untuk menghitung berbagai jenis jaringan pipa baik itu jaringan berbentuk parallel, non-loop dan loop. 4. Arah aliran pada jaringan ditentukan berdasarkan hasil iterasi tekanan pada titik node, dimana fluida mengalir dari tekanan tinggi ke tekanan yang lebih rendah. 5. Banyaknya node mempengaruhi jumlah iterasi. 6. Penambahan cabang pipa baru pada suatu jaringan yang telah ada dapat merubah keseimbangan laju alir pada jaringan tersebut yang mengakibatkan berubahnya harga tekanan pada tiap node dan harga laju alir pada tiap cabang. 7. Penempatan posisi dari pipa baru dalam suatu jaringan yang telah ada berpengaruh terhadap jumlah iterasi yang dilakukan.

Hasil perhitungan software didapatkan hasil seperti pada Tabel 4a dan 4b. Kasus 5: Jaringan Loop Dengan Cabang Baru Antara Titik Node 1 dan 4 Pada Kasus 5, satu buah cabang pipa ditambahkan antara node 1 dan 4, dengan dimensi pipa yang sama pada kasus 3 yaitu panjang 8000 ft, diameter 10 inci dan sudut elevasi pipa 00. Hasil perhitungan tekanan dari kasus 5 ini dibandingkan dengan hasil dari kasus 3.

Hasil perhitungan software didapatkan hasil seperti pada Tabel 5a dan 5b. Kasus 6: Jaringan Loop Dengan Cabang Baru Antara Titik Node 2 dan 4 Pada kasus 6, jumlah aliran yang masuk dan keluar pada titik-titik node sama dengan kasus 3, tetapi ada tambahan cabang pipa dimana satu buah cabang pipa ditambahkan antara node 2 dan 4, dengan dimensi pipa baru sama pada kasus 4 dan 5, yaitu diameter 10 inci, panjang 8000 ft dan sudut elevasi pipa 00

Hasil perhitungan software didapatkan hasil seperti pada Tabel 6a dan 6b. Kasus 7: Jaringan Loop Dengan Cabang Baru Antara Titik Node 2 dan 5 Pada Kasus 7, satu buah cabang pipa ditambahkan antara node 2 dan 5, dengan dimensi pipa baru yang Sama pada kasus 3, 4, dan 5.

Hasil perhitungan software didapatkan hasil seperti pada Tabel 7a dan 7b.

5

gc : Percepatan gravitasi GOR :Perbandingan Gas dan Minyak, SCF/STB HL(0) :Liquid hold up pada kondisi horizontal, fraksi HL(θ) :Liquid hold up pada kondisi inklinasi, fraksi λ : No slip liquid hold up, fraksi Nre : Reynold number NLV : Liquid velocity number NFR : Froude Number D : Diameter pipa, ft A : Area pipa, sq ft L : Panjang pipa, ft θ : Sudut pipa ψ : Faktor koreksi liquid hold up pada kondisi inklinasi z : ketinggian Z : Faktor deviasi gas

SARAN Saran-saran untuk pengembangan lebih lanjut pada model dan algoritma ini adalah sebagai berikut : 1. Penambahan data temperatur sebagai variabel pada data node. Dengan memasukkan temperatur sebagai variabel pada node, selain tekanan dan supply atau demand, range temperatur yang kenyataannya bisa sangat besar mulai dari wellhead (titik supply) sampai titik demand dapat dimasukkan ke dalam model. 2. Penambahan pesan-pesan kesalahan pada program yang dibuat, sehingga memudahkan di dalam melakukan analisa dan desain suatu jaringan yang kompleks. DAFTAR SIMBOL Pb : Tekanan bubble point, psi Rs : Solution Gas Liquid Ratio, Scf/STB γ g : Spesific gravity gas API : Derajat kualitas minyak, oAPI T : Temperatur, oF γ o : Spesific gravity minyak β : Faktor Akselerasi Bo :Faktor volume formasi minyak, STB/Bbl Bg : Faktor volume formasi gas, cuft/SCF ε : Toleransi kesalahan ρo : Densitas minyak, lb/cuft ρg : Densitas air, lb/cuft ρtp : Densitas total, lb/cuft µo : Viskositas minyak, cP µg : Viskositas gas, cP µm : Viskositas campuran, cP σo :Tegangan permukaan minyak, dyne/cm2 fns : No-slip friction factor ftp : Total friction factor qo : Laju alir minyak, bbl/d qg : Laju alir gas, MMScf/d q : Laju alir fluida, STB qij : Laju alir pada pipa/cabang, STB Qi : Laju alir pada titik node, STB g : Percepatan gravitasi

DAFTAR PUSTAKA 1. Ahmed, Tarek., “ Hydrocarbon Phase Behavior”, Gulf Pub. Co., Houston – Texas, 1989. 2. Brown, Kermit E., “The Technology of Artificial Lift Methods “, Volume 1, PennWell Pub. Co., Tulsa, Oklahoma, 1984. 3. Brown, Kermit E., “The Technology of Artificial Lift Methods “, Volume 4, PennWell Pub. Co., Tulsa, Oklahoma, 1984. 4. Budiawan, Yusman., “Aplikasi Modeel Jaringan Nodal Untuk Evaluasi Aliran Gas dalam Jaringan Pipa yang Kompleks”, Tugas Akhir, ITB, 1995. 5. Ikoku, C.U., “Natural Gas Production Engineering”, John Willey & sons Inc., New York, 1984. 6. Kadir, Abdul., “Dasar Pemrograman Delphi 5”, Jilid 1 dan 2, Penerbit Andi, Yogyakarta, 2002. 7. Stewart, Maurice., “CPI Gas Handling Facilities Design, Operation and Maintenance”, LDI Training. 8. Rubiandini R.S, Rudi., “ Pipeline Technlogy and Maintenance”., LDI Training.

6

9. Tian, S and Adewumi, M.A., “ A Simple Algorithm for Analyzing Gas Pipeline Network”, SPE Paper 25475 presented at the Production Operation Symposium held in Oklahoma, March 21 –23, 1993. 10. Tian, S and Adewumi, M.A., “ A New Algorithm for Analyzing and Designing Two Phase Flow Pipeline Networks”, SPE Paper 28177, November, 1993

7

INPUT DATA Hitung Harga Tekanan Rata-Rata Hitung Z, Rs, Bo, ρg, ρo, µg, µo, λ, dan µm

Asumsikan harga ftp Hitung harga perkiraan awal ρtp

Hitung NRen, NFR, NLV, L1,L2, L3, dan L4

Hitung harga q' awal

Tentukan Pola Aliran

q' = q''

Hitung HL(0), C, Y, HL(t), ρtp, fns dan ftp Hitung harga q''

TIDAK

q'' - q' < ε YA

q=q

Gambar 1. Diagram Alir Perhitungan Laju Alir dengan Metode Beggs and Brill 2)

8

3 2

3

2

4

1

1

4

5 5

Gambar 2. Sistem loop pada jaringan pipa

1

4 3

1

3 2

4

2

5

Gambar 3. Jaringan pipa sederhana

Gambar 4. Sistem jaringan kompleks 10)

9

Gambar 5. Penyederhanaan dari sistem jaringan kompleks 10) INPUT DATA dan Inisialisasi Tekanan

LAKUKAN LOOP PADA NODE

LAKUKAN LOOP PADA CABANG

HITUNG qij PADA CABANG DENGAN BEGGS AND BRILL

HITUNG FAKTOR AKSELERASI

AKHIR DARI LOOP PADA CABANG

TIDAK

Pi DIKETAHUI

TIDAK

YA

AKHIR DARI LOOP PADA NODE

APAKAH KONVERGEN

UPDATE HARGA TEKANAN (Pi)

YA

HITUNG Q PADA NODE (Qi)

HASIL

Gambar 6. Diagram Alir Algoritma Jaringan Pipa 10) 10

1

1

2

2

Gambar 7. Jaringan Pipa Paralel

4

3

5

4 6

2

1

1

5

2 6

7

3

7 8

Gambar 8. Jaringan Pipa Non-Loop Kasus 2

3 2

3

2

4

1

1

4

5 5

Gambar 9. Jaringan Pipa Loop Kasus 3

11

3 2

3

2 6

4

1

1

4

5 5


3 2

3

2

4

1

6

1

4

5 5

Gambar 11. Jaringan Pipa Loop Kasus 5 3 2

3

2 6

4

1

1

4

5 5


12

3 2

3

2

4

1

6

1

4

5 5

Gambar 13. Jaringan Pipa Loop Kasus 7 3 2

3

2

1

4

6

1

4

5 5

Gambar 14. Jaringan Pipa Loop Kasus 8 Tabel 1a. Output Data Node Kasus (Jumlah iterasi : 10) Node 1 2

Tekanan(psia) Q (BBL) 1000 40000 976.118 -40000

Betha 0 0.00119

Error 0 0.00001

Tabel 1b. Output Data Pipa Kasus 1 Cabang 1 2

Q(BBL) Dari Node Ke Node Diameter(in) Panjang(ft) 20000 1 2 10 9000 20000 1 2 10 9000

13

Theta 0 0

Tabel 2a. Output Data Node Kasus 2 (Jumlah iterasi : 348) Node 1 2 3 4 5 6 7 8

Tekanan(psia) Q (BBL) 1000 -150000 1258.447 0 1110.07 0 1281.281 30000 1281.281 30000 1281.281 30000 1135.658 30000 1135.658 30000

Betha 0 0.00076 0.00085 0.00076 0.00076 0.00076 0.00085 0.00085

Error 0 0.00002 0 0.00002 0.00002 0.00002 0 0

Tabel 2b. Output Data Pipa Kasus 2 Cabang 1 2 3 4 5 6 7

Q(BBL) Dari Node Ke Node Diameter(in) Panjang(ft) 90000 2 1 10 5000 60000 3 1 10 5000 30000 4 2 10 5000 30000 5 2 10 5000 30000 6 2 10 5000 30000 7 3 10 5000 30000 8 3 10 5000

Theta 0 0 0 0 0 0 0

Tabel 3a. Output Data Node Kasus 3 (Jumlah Iterasi : 103) Node Tekanan(psia) 1 1000 2 996.446 3 912.545 4 898.272 5 913.582

Q (BBL) 60000 40000 -30000 -40000 -30000

Betha 0 0.00037 0.00073 0.00073 0.00075

Error 0 0.00001 0.00003 0.00003 0.00002

Tabel 3b. Output Data Pipa Kasus 3 Cabang 1 2 3 4 5

Q(BBL) 9620.724 49620.724 19620.724 20379.276 50379.276

Dari Node Ke Node Diameter(in) Panjang(ft) Theta 1 2 10 5000 0 2 3 10 5000 0 3 4 10 5000 0 5 4 10 5000 0 1 5 10 5000 0

14

Tabel 4.a. Output Data Node Kasus 4 (Jumlah iterasi : 78) Node 1 2 3 4 5

Tekanan(psia) 1000 1001.735 966.126 938.659 943.812

Q (BBL) 60000 40000 -30000 -40000 -30000

Betha 0 0.00026 0.00096 0.00045 0.00045

Error 0 0.00001 0.00001 0.00004 0.00004

Tabel 4b. Output Data Pipa Kasus 4 Cabang 1 2 3 4 5 6

Q(BBL) 6780.952 33219.048 28522.833 11477.167 41477.167 25303.785

Dari Node Ke Node Diameter(in) Panjang(ft) 2 1 10 5000 2 3 10 5000 3 4 10 5000 5 4 10 5000 1 5 10 5000 1 3 10 8000

Theta 0 0 0 0 0 0

Tabel 5a. Output Data Node Kasus 5 (Jumlah iterasi: 135) Node 1 2 3 4 5

Tekanan(psia) 1000 1000.7 958.998 957.583 958.887

Q (BBL) 60000 40000 -30000 -40000 -30000

Betha 0 0.00017 0.00024 0.00023 0.00023

Error 0 0.00001 0.00003 0.00003 0.00002


Q(BBL) Dari Node Ke Node Diameter(in) Panjang(ft) 4078.572 2 1 10 5000 35921.43 2 3 10 5000 5921.428 3 4 10 5000 5655.249 5 4 10 5000 35655.25 1 5 10 5000 28423.32 1 4 10 8000


Tekanan(psia) 1000 985.412 947.475 946.776 950.184

Q (BBL) 60000 40000 -30000 -40000 -30000

15

Betha 0 0.0007 0.00018 0.00018 0.00037

Error 0 0.00001 0.00004 0.00004 0.00001

Theta 0 0 0 0 0 0



Theta 0 0 0 0 0 0


Tekanan(psia) 1000 982.79 928.988 924.766 955.081

Q (BBL) 60000 40000 -30000 -40000 -30000

Betha 0 0.00076 0.00041 0.00041 0.00102

Error 0 0.00001 0.00004 0.00004 0.00001



Theta 0 0 0 0 0 0

Tabel 8a. Hasil Output Node untuk Kasus 8 (Jumlah iterasi: 857) Node 1 2 3 4 5

Tekanan(psia) 1000 996.523 913.087 898.315 913.115

Q (BBL) 60000 40000 -30000 -40000 -30000

Betha 0 0.00037 0.00006 0.00074 0.00006

Error 0 0.00001 0.00024 0.00002 0.00024

Tabel 8b. Hasil Output Pipa untuk Kasus 8 Cabang 1 2 3 4 5 6


16

Theta 0 0 0 0 0 0

ANALISA JARINGAN PIPA LOOP-NODE DUA FASA MENGGUNAKAN METODE BEGGS AND BRILL

Recommend Documents