Simulação sequencial directa
Simulação sequencial directa (SSD, ou seu acrónimo em inglês - DSS) é um tipo de simulação sequencial habitualmente utilizada no ramo da geoestatística. É utilizada para estimação[1] e especialmente simulação estocástica dos nós de uma malha (ou Grid) na qual cada um deles está condicionado aos restantes simulados anteriormente. É chamada de directa porque ao contrátio de outras metodologias, como é o caso da simulação sequencial gaussiana ou simulação sequencial da indicatriz, não necessita de qualquer transformação da variável original (Soares, 2001).[2] Foi originalmente conceptualizada por Andre Journel em 1994[3] mas só em 2001 com Amílcar Soares foi avançada um implementação que corresponda aos critérios de uma simulação sequencial. Embora menos utilizada que a sua par, simulação sequencial gaussiana, é hoje utilizada em processos nas mais diferentes áreas: caracterização de reservatórios petróliferos, ciências da agricultura e agronomia, estudos ambientais, etc..
Nota histórica
[editar | editar código-fonte]Um método de simulação sequencial directa foi originalmente conceptualizado por Andre Journel em 1994. Embora este tivesse demonstrado que, desde que as distribuições locais utilizadas para simulação fossem centradas no estimadores média e variância de krigagem simples, a reprodução do modelo de covariância (modelo de variografia) estaria assegurada nesta implementação ainda que o mesmo não poderia ser dito da distribuição de probabilidades da variável original (histograma). Esta constatação foi mais tarde denominada de teorema de Journel (i.e. Journels Theorem, Caers, 2000)[4] e acabou por ser o obstáculo mais sério à sua utilização.
Em 1999 Jef Caers[5] propôs a utilização de uma técnica de pós-processamento das simulações para transformar os valores simulados para uma nova distribuição de probabilidades. Ainda anteriormente fora proposto uma metodologia semelhante por Journel e Xu[6] no entanto verificou-se que em alguns casos esta transformação poderia destruir a reprodução do modelo de variograma (Nowak & Srivastava, 1997).[7]
Em 2001 Amílcar Soares[2] avançou com uma nova implementação que consiste em usar a média e variância de krigagem simples não para definir uma distribuição local (à semelhança da simulação sequencial gaussiana) mas para re-amostrar uma distribuição local a partir da total. Os valores são depois simulados dentro desta re-amostragem. Com esta metodologia conseguiu tanto a reprodução do modelo de variograma como a reprodução do histograma.
Em 2003 foi feita uma nova implementação deste algoritmo muito semelhante à descrita por Amílcar Soares. A sua principal diferença está na maneira com se calculam os estimadores de média e variância da krigagem simples para posterior re-amostragem da distribuição da variável original (Oz, et al., 2003).[8]
Definição
[editar | editar código-fonte]A simulação sequencial directa deverá corresponder aos passos determinados num processo de simulação sequencial no qual existem dois passos com métodos estocásticos fundamentais:
- (Passo 5) Geração de um caminho aleatório sem repetição
- (Passo 7) Escolha aleatória sobre a fdp local no nó a simular (Ripley, 1987)[9]
O seu procedimento é o seguinte:
- Definição do caminho aleatório sem repetição.
- Calcular a função de distribuição cumulativa de probabilidades local para cada nó a partir dos estimadores da krigagem simples.
- Simular um valor a partir da distribuição cumulativa de probabilidades local.
No passo 3 da geração aleatória sobre uma função de distribuição de probabilidades o método de simulação sequencial directa implementada por A.Soares (2001) segue o seguinte modelo:
![](http://178.128.105.246/cars-http-upload.wikimedia.org/wikipedia/commons/thumb/b/b9/Explica%C3%A7%C3%A3o_acerca_do_m%C3%A9todo_da_simula%C3%A7%C3%A3o_sequencial_directa.png/600px-Explica%C3%A7%C3%A3o_acerca_do_m%C3%A9todo_da_simula%C3%A7%C3%A3o_sequencial_directa.png)
O esquema acima demonstra a re-amostragem da distribuição pelos intervalos definidos pelos estimadores locais (médias e variâncias de krigagem simples) de : e . Aqui é simulado um valor na função equivalente de .
Discussão
[editar | editar código-fonte]Em 2003 foi feito um estudo comparativo[10] entre as implementações feitas na biblioteca de geoestatística de programação, GSLIB[11] para simulação sequencial directa, e a implementação feita em 2003 (Oz, et al., 2003). As conclusões revelam que as diferentes implementações têm resultados muito semelhantes.
A SSD pode ser utilizada também em co-simulação nas suas vertentes co-simulação co-localizada sequencial directa e co-simulação co-localizada sequencial directa com distribuições de probabilidade conjuntas.[12]
Foi implementada também uma versão da SSD com anisotropias locais, simulação sequencial directa com anisotropias locais[13] na qual um modelo de variografia é descritizado para cada bloco ao invés da utilização de um modelo global para todas as localizações.
Ver também
[editar | editar código-fonte]- Simulação sequencial
- Simulação sequencial gaussiana
- Co-simulação co-localizada sequencial directa
- Co-simulação co-localizada sequencial directa com distribuições de probabilidade conjuntas
- Simulação sequencial directa com anisotropias locais
Referências
- ↑ Embora tenha como base processos estocásticos é por vezes (quando fortemente condicionada) usada como processo de estimação.
- ↑ a b Soares, A. (2001), "Direct Sequential Simulation and Co-Simulation" (2006), Mathematical Geology, Vol. 33, No. 8, November 2001
- ↑ Journel, A. G., 1994, "Modeling uncertainty: some conceptual thoughts", in Dimitrakopoulos, R., ed.,Geostatistics for the Next Century: Kluwer Academic Pub., Dordrecht, The Netherlands, p. 30–43.
- ↑ Caers, J., 2000, Direct sequential indicator simulation: Proceedings of 6th International Geostatistics Congress, Cape Town, South Africa.
- ↑ Caers, J., 1999, Adding local accuracy to direct sequential simulation: Stanford Center for Reservoir Forecasting, Annual Meeting 12, v. 2.
- ↑ A.G Journel, W Xu, "Posterior identification of histograms conditional to local data" ,Mathematical Geology, 26 (6) (1994), pp. 323–359
- ↑ Nowak, M.S., Srivastava, R.M., 1997. A geological conditional simulation algorithm that exactly honours a predefined grade-tonnage curve. In: Baafi, E.Y., Schofield, N.A. (Eds.), Proceedings of the Geostatistics Wollongong 96, Vol. 2. Kluwer Academic Publishers, Dordrecht, Netherlands, pp. 669–682.
- ↑ B. Oz, C.V. Deutsch, T.T. Tran, Y. Xie,DSSIM-HR: a FORTRAN 90 program for direct sequential simulation with histogram reproduction,Computers & Geosciences, 29 (1) (2003), pp. 39–51
- ↑ Ripley, B. (1987), "Stochastic Simulation", NY: John Wiley & Sons
- ↑ Robertson, R. K., Mueller, U. A., Bloom, K. M., 2006. Direct sequential simulation with histogram reproduction: A comparison of algorithms. Mathematical Geology 32 (3), 382–395.
- ↑ C.V. Deutsch, A.G. Journel, GSLIB: Geostatistical Software Library and User's Guide,(Second Ed.)Oxford University Press, New York (1998) 369pp
- ↑ Horta A., Soares A.,2010, "Direct sequential co-simulation with joint probability distributions",Mathematical Geosciences - 42(3):269-292.
- ↑ Ana Horta, Maria Helena Caeiro, Ruben Nunes and Amílcar Soares, 2010, "Simulation of Continuous Variables at Meander Structures: Application to Contaminated Sediments of a Lagoon",Quantitative Geology and Geostatistics, 2010, Volume 16, 161-172, DOI: 10.1007/978-90-481-2322-3_15