J. of Math. and Its Appl. ISSN: 1829-605X Vol. 1, No. 1 (2004), 15–19
Metode Asimilasi Data sebagai Estimasi Penyelesaian Masalah-masalah Lingkungan

## Erna Apriliani
Jurusan Matematika Institut Teknologi Sepuluh Nopember, Surabaya

## Abstrak
Asimilasi Data merupakan metode estimasi sistem dinamik stokastik yang merupakan penggabungan model dinamik dengan data-data pengukuran. Metode Asimilasi data telah banyak digunakan dalam mengestimasi masalah-masalah meteorologi, hidrodinamika. Di sini akan dibahas tentang apa metode asimilasi data itu, bilamana dipergunakan, apa kelebih- annya dibandingkan dengan metode lainnya serta bagaimana penerapannya pada masalah lingkungan
Kata kunci : sistem dinamik stokastik, interpolasi optimum, filter Kalman

## 1 Pendahuluan
Dalam kehidupan sehari-hari seringkali kita dihadapkan pada masalah - masalah esti- masi atau prediksi, khususnya untuk ma- salah yang bergantung pada waktu/masalah dinamik. Untuk sistem dinamik determin- istik, yaitu tidak memuat faktor stokastik maka penyelesaian secara numerik dari sis- tem dinamik tersebut dapat membantu menye- lesaikannya. Sedangkan untuk sistem di- namik stokastik biasanya data-data peng- ukuran banyak dipergunakan untuk memod- elkan serta memprediksi sistem tersebut untuk waktu yang akan datang.
Pemodelan dan estimasi berdasarkan data- data pengukuran belum tentu mencerminkan sistem yang sebenarnya, karena keterbatasan data yang dapat diperoleh. Pada pemodelan dan estimasi berdasarkan sistem dinamik yang ada terjadi penumpukan error estimasi yang terjadi karena kesalahan estimasi pada tiap langkah waktu. Untuk mengantisipasi keku- rangan dan kelebihan masing-masing metode digunakan metode asimilasi data. Metode asi-
milasi data merupakan penggabungan penyele- saian sistem dinamik dan data-data peng- ukuran untuk memperbaiki hasil estimasi se- hingga penumpukan error dapat dikurangi.
Masalah-masalah lingkungan seperti penye- baran limbah cair dilaut, penyebaran limbah cair pada air tanah, pencemaran udara, distri- busi suhu udara pada suatu daerah meru- pakan sistem dinamik yang tidak terlepas dari pengaruh faktor stokastik. Oleh karena itu perlu dikaji penerapan metode asimilasi data untuk menyelesaikan masalah-masalah tersebut, khususnya masalah limbah cair di laut. Metode asimilasi data yang dipergunakan disini adalah filter Kalman, sedangkan pema- kaian kedua metode lainnya masih dalam ka- jian.

## 2 Metode Asimilasi Data
Metode asimilasi data pada awalnya dikem- bangkan untuk memprediksi cuaca dan telah banyak diterapkan dalam masalah meteorologi [12], [11]. Kemudian penerapan metode asi-

## Erna Apriliani
milasi data pada masalah-masalah yang lain- nya telah banyak dilakukan antara lain oleh Zhang, 1996 [12] untuk pemodelan polusi udara, Voorrips, 1998 [11] menerapkan pada model gelombang laut, dan pada tahun 1999, Miller [8] menuliskan metode asimilasi data untuk model-model stokastik tak linear.
Langkah sederhana dari asimilasi data adalah mengganti atau memperbaiki model yang sudah ada berdasarkan data pengukuran yang diperoleh (substitusi data). Tetapi substi- tusi data tidak selalu dapat dilakukan karena substitusi data tidak dibenarkan dalam hukum- hukum fisika akibatnya model menjadi tidak setimbang dan hasil simulasi yang diperoleh tidak bermanfaat lagi [12]. Oleh karena itu yang paling penting dalam masalah asimilasi data adalah bagaimana mengkombinasikan an- tara model dan data pengukuran sehingga diperoleh hasil estimasi yang terbaik.
Sampai saat ini ada tiga teknik yang dapat dipergunakan dalam metode asimilasi data yaitu:
– Interpolasi Optimum
– Filter Kalman
– Metode Variasional 4-D

## 2.1 Interpolasi Optimum
Metode Interpolasi Optimum merupakan metode yang paling sering dipergunakan untuk memprediksi cuaca atau masalah meteorologi. Metode ini merupakan metode asimilasi data dengan pendekatan statistik yaitu teknik akar kuadrat terkecil (least squares technique). Metode ini dikembangkan oleh Lev Gandin (1963). Misalkan diberikan sistem dinamik
x k = f k ( x k − 1 , u k − 1 ) (1)
dengan pengukuran
z k = H k x k (2)
maka algoritma Interpolasi Optimum [11] adalah
(i) variabel keadaan x k diprediksi berdasar- kan
x f k = E ( x k ) = ¯ x k (3)
dengan kovariansi kesalahan estimasi
p f k = E (( x k − ¯ x k )( x k − ¯ x k ) t ) (4)
dimana E ( . ) adalah nilai ekspektasi.
(ii) Bobot koreksi kesalahan estimasi
K k = P f k H t k [ H k P f k H t k + R k ] ( − 1) (5)
(iii) Estimasi optimum dari variabel keadaan setelah data pengukuran dimasukkan
x k = x f k + K k ( z k − H k x f k ) (6)
dengan R k kovariansi kesalahan peng- ukuran.

## 2.2 Filter Kalman
Filter Kalman merupakan perumuman dari in- terpolasi optimum. Perbedaan antara Filter Kalman dan Interpolasi Optimum antara lain [9]:
• memasukkan model prediksi linear yang tidak tampak pada interpolasi optimum
• kovariansi kesalahan prediksi diperbarui se- tiap satuan waktu
• diterapkan secara global yaitu semua data pengamatan/pengukuran digunakan untuk mengestimasi semua variabel secara simul- tas, sedangkan interpolasi optimum hanya menganalisa setiap titik/variabel yang dip- ilih disekitar pengamatan/pengukuran.
Misalkan diberikan sistem dinamik sto- kastik
x k +1 = A k x k + B k u k + G k w k (7)
dengan pengukuran
z k +1 = H k +1 x k +1 + v k (8)
Untuk sistem dinamik stokastik (7),(8) tidak dapat dilakukan estimasi variabel keadaan x k dengan cara sistem observer biasa karena kesalahan estimasinya ˜ x k tidak kon- vergen ke nol [3]. Estimasi variabel keadaan dilakukan dengan menggunakan filter Kalman.
Filter Kalman merupakan estimator linear terbaik untuk x o , w k dan v k berdistribusi se- barang dan merupakan estimator optimum untuk x o , w k dan v k berdistribusi normal. Pada filter Kalman estimasi variabel keadaan dilakukan dengan cara memprediksi varibel tersebut dan selanjutnya dilakukan peng- ukuran untuk mengkoreksi prediksi variabel se- cara rekursi dengan cara meminimumkan kova- rian kesalahan penaksiran ( x k − ˆ x k ) ; x k meru- pakan variabel keadaan sebenarnya dan ˆ x k merupakan penaksiran dari variabel keadaan.
Secara lengkap, algoritma dari filter
Kalman adalah sebagai berikut [1]
Limits Journal of Mathematics and Its Applications, Volume 1, Number 1, May 2004
(i) Tahap Prediksi (Time Update)
ˆ x k +1 ( − ) = A k ˆ x k + B k u k P k +1 ( − ) = A k P k A t k + G k Q k G t k (ii) Tahap Koreksi (Measurement Update)
K k = P k ( − ) H t k ( H k P k ( − ) H t k + R k ) − 1 P k = ( I − K k H k ) P k ( − ) ˆ x k = ˆ x k ( − ) + K k ( z k − H k b x k ( − ))
Diberikan taksiran awal ˆ x 0 = x 0 dengan kovarian kesalahan penaksiran awal P 0 (+) = P 0 > 0

## 2.3 Metode Variasional 4-D
Metode lain yang juga sering digunakan dalam metode asimilasi data adalah Metode Varia- sional Empat Dimensi. Ide dasar metode ini adalah meminimumkan jarak antara trayek- tori penyelesaian model dengan pengamatan pada interval waktu tertentu [12]. Metode yang digunakan untuk mendapatkan jarak min- imal adalah metode adjoint . Pada metode ini diperlukan data untuk sepanjang interval waktu untuk memperbaiki model sehingga san- gat sulit diterapkan untuk peramalan secara on-line .

## 3 Asimilasi Data pada Masalah Lingkungan
Masalah-masalah lingkungan yang dimaksud- kan di sini adalah masalah-masalah yang meru- pakan sistem dinamik stokastik, baik time- variant maupun time-invariant, seperti keting- gian air sungai, ketinggian gelombang air laut, proses salinasi, pencemaran air tanah, pence- maran udara dan masalah lainnya. Metode asi- milasi data ini telah diterapkan pada pemod- elan polusi udara [12], model-model kimia di atmosfir [10], model gelombang laut [11].
Penulis beserta tim penelitiannya telah menggunakan filter Kalman sebagai metode asimilasi data pada masalah estimasi keting- gian air sungai [4], estimasi penyebaran limbah cair di laut [7], estimasi aliran air tanah [2]. Pada saat ini sedang dikaji metode asimilasi data lain untuk masalah estimasi ketinggian air tanah, estimasi penyebaran limbah cair pada air tanah dan estimasi penyebaran limbah cair
di laut dengan metode variasional empat di- mensi. Pengembangan algoritma pada metode asimilasi data juga dikaji dan dikembangkan misalnya modifikasi filter Kalman [5, 6].
Dalam melakukan penerapan Filter Kalman sebagai metode asimilasi data dilakukan langkah-langkah sebagai berikut:
(i) Pemodelan: Model dari masalah terse- but dibangun berdasarkan hukum-hukum fisika dan fenomena alam yang berlaku, misalnya model aliran air sungai mem- punyai model matematika yang dikenal dengan persamaan Saint Venant.
∂h ∂t + D ∂u ∂x = 0 (9) ∂u ∂t + g ∂h ∂x + C f u = 0 (10)
dengan syarat batas dan syarat awal
h (0 , t ) = ψ b ( t ) h ( x, 0) = 1 u ( x, 0) = 0
u ( L, t ) = u N ( t )
dimana h ( x, t ) ketinggian air diatas D , D kedalaman sungai, u ( x, t ) rata-rata ke- cepatan aliran sepanjang irisan penam- pang melintang, t menyatakan waktu, x posisi sepanjang sungai, g kecepatan gravitasi, C f konstanta gesekan, ψ b ket- inggian air pada posisi x o dan u ( L, t ) ke- cepatan air pada batas x N .
(ii) Pendiskritan: perlu dikaji metode pendiskritan yang sesuai dengan masalah yang dikaji, misalnya untuk persamaan Saint Venant tersebut digunakan skema Pressmann. Dari pendiskritan ini diper- oleh bentuk sistem ruang keadaan ( state space system )
X k +1 = A k X k + B k U k (11)
(iii) Pengukuran yang dapat dilakukan:
Pada tahap ini harus dikonstruksi persa- maan pengukuran yang mengkaitan data- data pengukuran yang diperoleh dengan variabel keadaan yang akan diestimasi. Perlu diperhatikan bahwa model sistem dan pengukuran haruslah terobservasi ( observable ).
(iv) Pembentukan Sistem Stokastik: Di sini perlu dikaji faktor-faktor stokastik
Limits Journal of Mathematics and Its Applications, Volume 1, Number 1, May 2004
yang mempengaruhi sistem dan meru- pakan noise pada sistem serta noise pada pengukuran.
X k +1 = A k X k + B k U k + G k W k (12)
dengan W k gangguan sistem yang di- asumsikan berdistribusi normal Gauss dengan mean nol dan kovarians Q k .
Misalkan pengukuran yang diperoleh mempunyai persamaan
z k +1 = H k +1 X k +1 + V k +1 (13)
dengan V k gangguan pengukuran yang diasumsikan berdistribusi normal Gauss dengan mean nol dan kovarians R k .
(v) Simulasi penerapan filter Kalman : dibuat program komputer untuk mensi- mulasikan hasil estimasi berdasarkan data-data pengukuran.
Pada langkah-langkah tersebut di atas, fak- tor stokastik dapat muncul pada saat pe- modelan sehingga sistemnya dinamakan sistem persamaan diferensial stokastik atau dibentuk setelah pendiskritan dan pembentukan persa- maan pengukuran untuk memudahkan ma- salah.
Pada penelitian-penelitian yang telah dilakukan penulis dan tim, data-data peng- ukuran dibangkitkan dari simulasi yaitu menambahkan variabel random berdistribusi normal pada sistem deterministiknya (model dan pengukuran). Hasil estimasi dengan meng- gunakan filter Kalman adalah sebagai berikut:
0 5 10 15 20 25 30 35 40 0 5 10 15 20 25 30 35 40 45 ketinggian air o- : ketinggian hasil +- : ketinggian real estimasi posisi ke x k
Gambar 1: Estimasi Ketinggian Air
0 5 10 15 20 25 30 35 40 0 5 10 15 20 25 30 35 40 45 posisi ke x k kecepatan aliran +- : kecepatan real o- : kecepatan hasil estimasi

## Gambar 2: Estimasi Kecepatan Aliran Air

## 4 Penutup
Berdasarkan kajian pustaka dan beberapa penelitian yang telah dilakukan dapat disim- pulkan bahwa metode asimilasi data dapat dit- erapkan untuk mengestimasi masalah-masalah lingkungan; sedangkan efisiensi waktu kom- putasi dan tingkat akurasi masing-masing dari ketiga metode masih perlu dikaji. Modifikasi algoritma masih perlu dikembangkan untuk mendapatkan hasil yang lebih efisien dan aku- rat.

## Daftar Pustaka
[1] Anderson, R., Optimal Filtering , Prentice Hall, 1979.
[2] Handayani,A., Estimasi Filter Kalman terhadap Distribusi Konsentrasi Polusi Air Tanah , Tugas Akhir S1, ITS, 2003.
[3] Apriliani, E., “Perbandingan antara Ob- server Deterministik dan Filter Kalman”,
Majalah Ilmiah Himpunan Matematika Indonesia , Vol. 6, no. 1, 2000.
[4] Apriliani, E., “The Application of the Modified RRSQRT Filter on One Dimen- sional Shallow Water Problem”, Proceed- ings ITB , Vol 34 no 1, ITB, 2001.
[5] Apriliani, E., “The Reduced Rank SRIF”, Majalah Ilmiah Himpunan Matematika Indonesia , vol. 7 no. 2, 2001.
[6] Apriliani, E., “The Estimation Of The Water Level By The Reduced Rank
Limits Journal of Mathematics and Its Applications, Volume 1, Number 1, May 2004
Square Root Information Filter”, Proceed- ings of the Asia - Pasific Vibration Con- ference , vol. II, Jilin Science and Technol- ogy Press, China, 2001.
[7] Endang S., Estimasi Penyebaran Limbah Cair di Laut dengan Filter Kalman dan Modifikasi Filter RRSQR , Tugas Akhir S1, ITS, 2002.
[8] Miller, R.N., Carter, E.F. and Blue, S.T, Data Assimilation into Nonlinear Stochas- tic Models , Tellus 51A, 1999.
[9] Schlatter, T.W., -, Four-Dimensional Data Assimilation , internet.
[10] Seger, A., Data Assimilation in Atmo- spheric Chemestry Models Using Kalman Filtering , Desertation TU Delft, 2002.
[11] Voorrips,A., Sequential Data Assimilation Methods for Ocean Wave Models , Diserta- tion, The Netherlands, 1998.
[12] Zhang, Xue-Fen, Data Assimilation in Air Pollutin Modelling , Disertation, TU Delft,
1996.