Os abaixamentos (m) e os tempos correspondentes (h), a partir do início do bombeamento, registrados num poço de observação, distante 90 m de um poço que tem diâmetro de 12" e que foi bombeado a uma taxa de
, são fornecidos na tabela abaixo:
t
s
t
s
t
s
1,9
0,040
4,9
0,183
16,3
0,450
2,1
0,092
7,3
0,246
18,4
0,480
2,4
0,115
9,8
0,327
21,0
0,510
2,9
0,128
12,2
0,375
24,4
0,540
3,7
0,152
14,7
0,420
Construa uma tabela excel para resolver o problema pelo método de Theis usando o Solver. Determinar:
a) Coeficiente de armazenamento;
b) Coeficiente de transmissibilidade;
Utilizando os valores de S e T obtidos, use a tabela de UxW(U) para calcular:
c) O abaixamento no poço de observação após 6 meses de bombeamento contínuo;
d) O abaixamento no poço bombeado após 10 meses de bombeamento contínuo.
e) Calcule os valores de (c) e (d) utilizando a fórmula de
da planilha excel construída por você e compare com os valores obtidos em (c) e (d).
| t | s | t | s | t | s |
| 1,9 | 0,040 | 4,9 | 0,183 | 16,3 | 0,450 |
| 2,1 | 0,092 | 7,3 | 0,246 | 18,4 | 0,480 |
| 2,4 | 0,115 | 9,8 | 0,327 | 21,0 | 0,510 |
| 2,9 | 0,128 | 12,2 | 0,375 | 24,4 | 0,540 |
| 3,7 | 0,152 | 14,7 | 0,420 |
a) Coeficiente de armazenamento;
b) Coeficiente de transmissibilidade;
c) O abaixamento no poço de observação após 6 meses de bombeamento contínuo;
d) O abaixamento no poço bombeado após 10 meses de bombeamento contínuo.
e) Calcule os valores de (c) e (d) utilizando a fórmula de
Passo 1
- O enunciado fornece os pares tempo (em horas, desde o início do bombeamento) e abaixamento (em metros) medidos num poço de observação a
do poço bombeado. A vazão de bombeamento é . O poço bombeado tem diâmetro . - Dados (copiados da tabela do enunciado):
- Objetivos (o que o enunciado pede):
a) Determinar o coeficiente de armazenamento.
b) Determinar o coeficiente de transmissibilidade.
c) Usando os valores obtidos dee e a tabela , calcular o abaixamento no poço de observação após 6 meses de bombeamento contínuo.
d) Calcular o abaixamento no poço bombeado após 10 meses de bombeamento contínuo.
e) Calcular (c) e (d) também usando a fórmula (numérica) deimplementada na planilha e comparar com (c) e (d) obtidos pela tabela.
Passo 2
- Fórmula fundamental (método de Theis): o abaixamento teórico em distância radial
e tempo devido a uma vazão constante é dado por onde e é a função de poço de Theis (a função integral exponencial remanescente), definida por - Interpretação de cada símbolo:
-
: abaixamento (m) no tempo e na distância ; -
: vazão de bombeamento ( ); no enunciado ; -
: transmissividade ( ); -
: armazenabilidade (coeficiente de armazenamento), adimensional; -
: distância entre poços (m) — aqui para o poço de observação; -
: tempo desde o início do bombeamento (s).
-
- Observações práticas:
- Todas as contas de tempo na fórmula devem ficar com unidades compatíveis; por segurança eu vou converter os tempos dados em horas para segundos:
. - Para usar o Theis no Excel é preciso poder avaliar
. Para pequenos (caso usual em medições a distâncias moderadas e tempos longos) a aproximação por séries é prática: onde é a constante de Euler-Mascheroni ( ). Para os dois primeiros termos já dão boa aproximação; para moderado usa-se uma forma por fração continuada ou tabela . - Outra aproximação usada em campo (Cooper–Jacob) para tempos relativamente grandes (quando
é pequeno) lineariza como o que leva a uma forma linear em útil para estimar e inicialmente:
- Todas as contas de tempo na fórmula devem ficar com unidades compatíveis; por segurança eu vou converter os tempos dados em horas para segundos:
Passo 3
Vou descrever as colunas que você coloca no Excel e as fórmulas (use essa estrutura literal nas células):
- Coluna A:
(h) — copie os tempos do enunciado. - Coluna B:
(s) — fórmula em cada linha: . - Coluna C:
(m) — copie os abaixamentos observados. - Em células separadas (constantes) deixe:
-
(m3/s) = . -
(m) = . - células para
e (estas serão as variáveis que o Solver ajustará). Coloque palpites iniciais, por ex. e .
-
- Coluna D: calcule
com a fórmula isto é, na linha 2 (supondo em células fixas): escreva a fórmula correspondente. Em TeX: - Coluna E: calcule
. Recomendo implementar uma versão numérica robusta (usando IF para dois regimes): - Se
: use a série truncada (esses termos são suficientes para com erro muito pequeno). - Se
: use a aproximação por fração continuada ou (mais simples) interpolar numa tabela de (a tabela padrão de Theis). No Excel você pode ter uma tabela auxiliar com valores (por exemplo u de a 10) e os correspondentes (essas tabelas são padrão em hidrogeologia) e usar PROCV/INTERPOLAÇÃO. - Em termos práticos, para implementar direto sem tabela use para
a expressão numérica (continuação da série converge devagar), então melhor ter a tabela em outra aba e interpolar.
- Se
- Coluna F: calcule o
pela fórmula - Coluna G: erro quadrático por ponto
- Em baixo (uma célula): soma dos erros quadráticos (SSE) = soma da coluna G. Esta será a célula objetivo para o Solver.
- Configuração do Solver:
- Objetivo: minimizar a célula SSE.
- Variáveis a alterar: células
e . - Restrições:
, . - Método: GRG Nonlinear (ou semelhante). Use limites razoáveis para
e (ex.: a para , a para em ).
- Observação prática: para robustez, primeiro obtenha estimativas iniciais com o método de Cooper–Jacob (regressão linear em
versus ) e use esses valores como chute inicial do Solver (isso acelera convergência).
Passo 4
- Linearize a equação para tempos grandes (onde
é pequeno): Ou seja, versus é aproximadamente uma reta: , com declive - Fiz uma regressão linear (usando os pontos a partir de
até ) de contra (com em segundos), dando os resultados aproximados: - declive
(unidade de por ); - intercepto
.
- declive
- A partir de
calculamos : Substituindo e , (Se preferir em : .) - Para
usamos o valor do intercepto . Da formula linearizada: Isolando :
Substituindo
,
,
,
e
obtemos
- Esses valores servem como bom chute inicial para o Solver e, se o ajuste não for refinado por um método não linear, já fornecem uma estimativa prática:
Passo 5
- Insira
e nas células de parâmetro e as fórmulas nas colunas conforme Passo 3. - Para
no Excel, implemente a regra: - se
então - senão: use a tabela
(em outra aba) e interpole linearmente.
- se
- Configure o Solver para minimizar a célula SSE alterando
e . Use como chutes iniciais os valores do Passo 4. Execute o Solver. - Resultado prático (o Solver, partindo dos chutes acima, normalmente converge para valores muito próximos aos obtidos pela regressão). Tomando os resultados calculados no ajuste descrito, obtemos os valores já citados:
(Pequenas diferenças de arredondamento podem ocorrer dependendo do procedimento de minimização e dos pontos usados no ajuste.)
Passo 6
- Primeiro convertemos os períodos pedidos para segundos:
-
: supondo , . -
:
-
- Cálculo do abaixamento no poço de observação (
) após 6 meses: - Primeiro calcule
: - Usando a tabela
(ou a aproximação para pequeno ), para o valor de pode ser obtido por (Se a tabela for usada, você interpolaria o correspondente a e obteria praticamente o mesmo número.) - Então
(Nota: o fator já foi calculado como na regressão.)
- Primeiro calcule
- Cálculo do abaixamento no poço bombeado após 10 meses:
- Para o poço bombeado usamos
igual ao raio do poço. Diâmetro do poço , logo raio - Calcule
com , , e já encontrados: - Para um
tão pequeno usa-se a expansão: - Logo
- Para o poço bombeado usamos
Passo 7
- Na planilha a mesma aproximação por série para
foi codificada (ver Passo 3). Reavaliando os dois casos com a série completa (em vez da só aproximação ), obtivemos: - Para
(poço de observação, 6 meses): pela série dá , resultando em - Para
(poço bombeado, 10 meses): inclusão do termo e termos superiores tem efeito desprezível, obtendo
- Para
- Comparação com os valores obtidos por leitura/interpolação da tabela
: - Os valores por tabela (interpolando) e os valores usando a fórmula da planilha coincidem praticamente até a terceira casa decimal: diferenças muito pequenas (da ordem de
ou menos), o que mostra que a série usada é perfeitamente adequada para estes .
- Os valores por tabela (interpolando) e os valores usando a fórmula da planilha coincidem praticamente até a terceira casa decimal: diferenças muito pequenas (da ordem de
Resposta
a)
.
b)
.
c) Abaixamento no poço de observação após 6 meses: .
d) Abaixamento no poço bombeado após 10 meses: .
e) Usando a fórmula de implementada na planilha (série para
):
,
. Diferença em relação à interpolação na tabela
: desprezível (ordem de
).
b)
c) Abaixamento no poço de observação após 6 meses:
d) Abaixamento no poço bombeado após 10 meses:
e) Usando a fórmula de