Este trabalho é motivado pelo interesse em modelar padrões espaciais em dados composicionais. A categoria de problemas de interesse envolve, por exemplo, as frações granulométricas de um solo ou composição química de uma rocha, isto é, estruturas de dados que as observações são parte de algum todo e referenciadas espacialmente. O interesse central está em predizer, conjuntamente, o padrão espacial dos componentes na área de estudo respeitando a restrição da soma das componentes corresponderem ao total. Neste sentido, combina-se a teoria de dados composicionais originalmente desenvolvida para observações independentes com métodos geoestatísitcos multivariados. A abordagem proposta declara explicitamente um modelo paramétrico que descreve a espacilização das variáveis. Em particular neste trabalho, o objetivo é propor e implementar um modelo geoestatístico bivariado para dados composicionais obtendo inferências via verossimilhaça e expressçies para predições espaciais, apresentado, discutindo e disponibilizando resultados e implementação computacional.A metodologia propostas é utilizada na análise de um conjunto de dados simulados, cujas composições são formadas por três componentes, possibilitando a obtenção de mapas de predição para cada um dos componentes.