Story Transcript
�-
♦
� Diputación
Instituto Tecnológico
Provincial
GeoMinero de España
de Al icante
L
L NORMAS U N= D A D
DE DE
EXPI OTAC 1QN
DE
LA
J 2 J O N A- C A R R A S O U E T A Junio de 1992.
i L
L
L
p
i
y L
MINISTERIO DE INDUSTRIA, COMERCIO Y TURISMO
........ -Dio. . . .
... 3unio de 1992.
C>
542/E
Nº
AaM SUBTERRANEAS
Na
PROYEM AGRUMM TITUM PiWYECM
ESTUDIOS DE EVAMA=0N Y PROTE=N DE 11ECURSOS REDRICOS SUBITERRAMOS DE IA PROVINCIA DE ALICANTE
N- DIVISICK AMAS, G.
Nº PIANIFICP=CN Fma mm=cet
incmo
27-XI--91
A- T 26-VI-92
INOME (Título)
NGMAS DE EXPW=0N DE LA WIDAD JI~CAMUMM92h
C=M(S)
wica(s)
JUCAR
a»C«U»D(ES) ADTC~(S)
VAIENCIANA
PROVINCIAS
ALICANTE
2
Este
estudio
Tecnológico Provincial cooperación
y de
ha
sido
Geominero Alicante
existente
colaboración de TEYGE,
realizado de
en entre
S.A.,
España el
por y
marco
ambos
del
el la
Diputación convenio
organismos,
Director del estudio.
DPA.Luis Rodríguez Hernández.
Director del estudio.
TEYG8,S.A. Joaquin Barba-Romero Muñoz. Emilio Orejudo Ramírez. Ernesto García Sánchez. Enrique Vázquez Suñé. Enrique Martí Requena. José Haría Catalán Alonso.
con
como empresa contratista.
ITGU.Ramón Aragón Rueda.
Instituto
de la
3
INDICE MEMORIA. Pág. 1 SITUACION GEOGRAFICA,
7
2.
9
CLIMATOLOGIA,
2.1.
TEMPERATURA.
10
2.2.
PLUVIOMETRIA.
11
2.3.
EVAPOTRANSPIRACION.
11
3.
HIDROLOGIA,
12
4.
GEOLOGIAe
13
4.1.
ESTRATIGRAFIA.
14
4.2.
TECTONICA.
18
HIDROGEOLOGIA,
20
5.
5.1.
CARACTERISTICAS LITOLOGICAS.
21
5.2.
CARACTERISTICAS GEOMETRICAS.
22
5.2.1.
Aculfero de Jijona.
23
5.2.2.
Aculfero de Carrasqueta.
24
5.3.
CARACTERISTICAS PIEZOMETRICAS.
25
5.3.1.
Aculfero de Jijona.
25
5.3.2.
Aculfero de Carrasqueta.
27
5.4.
PARAMETROS HIDRODINAMICOS
28
5.4.1.
Aculfero de Jijona
28
5.4.2.
Aculfero de Carrasqueta
28
CAPTACIONES EXISTENTES.
29
5.5.1.
Aculfero de Jijona.
29
5.4.2.
Aculfero de Carrasqueta.
30
5.5.
5.6.
FUNCIONAMIENTO HIDRAULICO
5.6.1.
Aculfero de Jijona.
30 30
4
5.6.2.
Acultero de Carrasqueta.
32
5.7. BALANCE HIDRICO.
34
5.7.1.
Aculfero de Jijona.
34
5.7.2.
Aculfero de Carrasqueta.
34
5.8. GRADO DE EXPLOTACION DE LA UNIDAD.
35
5.8.1.
Aculfero de Jíjona.
35
5.8.2.
Aculfero de Carrasqueta.
35
5.9. RESERVAS EXPLOTABLES.
36
5.9.1.
Acultero de Jijona.
36
5.9.2.
Aculfero de Carrasqueta.
36
6. HIDROOUIMICA,
37
6.1. ACUIFERO DE JIJONA.
38
5.2.
ACUIFERO DE CARRASºUETA.
38
7.-EVALUACION DE VOLUMENES DE AGUAS RESIDUALES,
39
8. ANALISIS DEL USO DEL AGUA. VOLUMENES DISPONIBLES Y DEMANDAS.
41
8.1.
USO DEL AGUA.
42
8.2.
ESTINACION DE LA DEMANDA.
44
8.3.
RELACION APORTACION-DEMANDA.
44
9. ANALISIS DE INFRAESTRUCTURAS.
45
10. NORMAS DE EXPLOTACION,
47
10.1.
NORMATIVA GENERAL.
48
10.2.
ZONIPICACION.
so
10.3.
PERIMETRO DE PROTECCION.
51
11. MODELO ESTOCASTICO.
54
METODOLOGIA
55
11.1.
11.2.
MODELIZACION LINEAL.
62
11.3.
IDENTIFICACION DEL MODELO LINEAL.
85
11.4.
PASE DE CALIBRACION DEL MODELO LINEAL.
98
11.5.
ANALISIS DE RESIDUOS Y VERIFICACION.
107
11.6.
EXPLOTACION DEL MODELO.
112
11.7.
CONCLUSIONES AL MODELO.
120
12.
CONCLUSIONES GENERALES.
126
6
INDICE DE PLANOS.
Plano nº
1:
Plano nº 2:
Cartografía hidrogeológica.
Perímetro de protección.
IL
8
La
unidad
de
Norte-Noroeste de
de
la Sierra de
Jijona-Carrasqueta la población
está
de Jijona,
La Carrasqueta.
situada
ocupando
Su práctica
al
parte
totalidad
se
encuentra incluida en el término municipal de Jijona.
La
cota
topográfica
materiales
permeables
parte alta
de
de
oscilan
la Sierra de
los entre
afloramientos 1.000
m.s.n.m.
Carrasqueta y 450
las inmediaciones de la población de Jijona.
de
los
en
la
m.s.n.m.
en
or C> T-a C> OG 3: p,
10
Los
datos
climáticos
parámetros
climatológicos
realizadas
periódicamente
Meteorológico
Zonal
termopluviométrica utilizados
han
por
el
para
basado
calcular
en
las
colaboradores
Valencia,
Jijona. período
pluviometría el período
2.1.
se
de
de
abarcan
utilizados
datos
Los
medidas
del
la
en
termométricos
(1944-1976)
y
los
de
(1980-1991).
TEMPERATURA.
Los
meses
en
que
alcanzan
se
temperatura
oscilación días
incluso
anual que
40ºC
16,75ºC. a
alcanzarse
(7-VIII-1976)
negativos
El
y
riesgo sólo
En
el
se
los
de
de
alcanzan
de
de
es en
extremos
bajo, los
Diciembre-Enero y ni siquiera todos los años.
La
Agosto
extremos
mínimos
heladas
Enero.
meses
valores
valores
altos
mientras que el de
diferencia,
con
de
es
llegan
(2-1-1991).
valores
es,
más
valores
temperatura son los de Julio y Agosto, menor
14,8 ºC;
se muestra en la figura 1.
su distribución anual
-7ºC
Centro
Estación
La temperatura media anual de la zona es de
hay
los
ya
meses
de de que de
Temperuturu-j
(valores en
ESTACION METEOROLCG4CA DE J.JONA En 7,02
Tcny-vrntum
�m 23.7
20,13
1
12,44
10.2
P3
Jn
my
At.
mz
Ft>
se
Ag
20,7
2.3.76
TEMPERAITURAS, so 28 26 24 22 20 18 id 14 12 10
4
£m
Fb
mz
Ab
11
my
Ag
Oc
Se
MIRSES
FIG.
:L
No
Di
CIC, 6
Nc1.2H
Di 7.9
2.2.
PLUVIONETRIA.
La pluviometría período
anual es de 382 mm para
(1980-91).
Su
distribución
la media del
mensual
puede
observarse en la figura 2.
El
mes
mientras f inal
que
del
mínima
los
de
verano
Noviembre) los
de
y
en
últimos
de
otro
del
primavera
1991
en que
otoño
(Mayo) . 1990
y
el
se
de
concentran
(Septiembre
Hay
que
1991
han
Julio, al
Octubre
destacar sido
que
secos,
se han registrado cantidades de
netamente
mm,
es
pluviosidad
inicios
1989,
200
interanual.
mayor
años
especialmente agua
e
pluviosidad
a
inferiores
la
media
Sin embargo este período ha sido precedido de
1986-88
interanual
pluviometría
con
(461
mm)
y
el
año
superior
muy 1982
en
que
a
se
la
media
alcanzaron
548 mm/año.
2.3.
EVAPOTRANSPIRACION.
Los
valores
la fórmula de 7hornthwaite
calculado por medio de media del mm.,
período.
figura
3.
Se
ha obtenido
distribución
cuya
Se
corresponde
observa al
mes
potencial
evapotranspiración
de
de
mensual
claramente Julio
produce en el mes de Enero.
y
un
un
valor
puede un
anual
observarse
máximo
mínimo
de
de 15
235
se
han
para la de
1150
en mm
mm que
la que se
PluAomctrfc (Vdoms en mm)
EGW10N METEOR=GICA DE JJONA En
Fb
27.t
Plv%iomotr(c
Ab
mz
24
40,1
27,2
Jn
my 53,7
-A 5,74
17.75
Ag 16.2
PLUVIOMETRIA 70 so -
44 -
30 -
Y
0 Fin
Fb
M&
Ab
My
ín
j1
kg
Se
MESES
FIG.
2
oc
No
se 42,5
C542
Ng 57.1
Di 16.6
12
La
evapotranspiración
evapotranspiración entre
Noviembre
de
campo
a
la
de
25
real
potencial
y Febrero. mm
la
para
Si
se
es
equivalente
los
meses
evapotranspiración
potencial
una
real
entre
la
comprendidos
considera
evapotranspiración
a
capacidad
sería
los
igual
meses
de
Noviembre a Marzo y parte de Abril.
La
figura
3
muestra
el
balance
la
para
media
del
período en la estación de Jijona.
3.- HIDROLOGIA.
La en
las
red
hidrográfica
sierras
confluye
de
una
épocas de la
red
fuertes
fuente
variaciones
estación que
se
y
inicia
Carrasqueta
y
atravesando
Castalla,
únicamente
se
acuífero.
conduce
agua
Se en
Solamente existe un manantial,
cuyo y
caudal responde
sufre
importantes
directamente
a
la
la época.
hace
un
de
un
(1991)
año
hidrométricaen carece
que
lluvias.
estacionales
pluviosidad de
Desde
hídrica
Nuches,
de
zona
la
la superficie aflorante del
transversalmente trata
de
Barranco
en
Madroñal
Peñarroya,
de el
en
existente
el
se
ha
Manantial
de
registro
histórico
anual de agua drenada por la fuente.
instalado Nuches, de
la
por
una lo
cuantía
T~
Axo.
P (mm)
igay-91
1600
11,29
cet.
HM.
Dic.
5,k 20
57.10
iraso
7.90
No0
1420
IZ44
Ene.
F;tt�
m3f.
Abr.
27,10
Zk no
27,20
7, co
lan
7Z 76
M70
M30
23170
Mi/.
jun-
JUL
4alo
5370
17.75
6,74
1623
4260
3nig
%81
Se>
Totol
FV~N~C" Pr~ »40R~RE fx 82
14.3
200
1,66
2134
294
.197
f� sú
1134
1(156
IQW
aw
1,90
E
8635
"46
2ZO2
17,99
2117
3473
5aw
87.41
13ffh W
181,60
lez67
14(164
d
31.co
mcn
31.00
31.0D
21LM
31.00
mm
31.co
3Q03
31,ci3
31.00
Mm
N ET.P.(~
0~
11,10
ík7!S
!k lo
9.45
insa
11,95
1445
DUS
IN 3D
15km
1-590
iz5o
C247
12
17.73
14.64
lag2
Pleo
enoo
iiciz7
17280
2xez
21elo5
14«W
occ
Nm
-2ft27
meo
(k. -1.13
Em. t240
Feb, isce
-Meo
Abr.
mov-
an
id
-fiago
-%57
-MOS
_=ce
BALANCE MEDIO
mo leo 160 140 120 100 en m 40 co 0 11
3
FIG.
3
A9_=45
sal -10400
IM72
Totd -7%W
13
14
La del
zona
dominio
Prebético
ocupada
por
Prebético,
Interno,
Castalla-Ibi
y
que
que
el
acuífero
se
concretamente
en
se
sitúa
podría
al
Sur
enmarca el
de
corresponder
dentro
denominado alineación
la a
una
unidad
intermedia entre el Prebético y el Subbético.
4.1.
se
ESTRATIGRAFIA.
La
serie
sitúa
la
base a techo ,
estratigráfica unidad
que
aflora
Jijona-Carrasqueta
el
en está
área donde
definida,
de
amarillas
en
por los siguientes materiales:
4. 1. 1.CRETACICO SUPERIOR
4.1.1.1 Cenomaniense-Turoniense.
Está constituido bancos
gruesos
calcáreas. estar
muy
con
Estas
por
juntas
calizas,
diaclasadas
250
(C211
m.
de
margosas llenas y
a
calizas
de
margas
radiolarios,
suelen
pasadas
y
de
veces
apareciendo con frecuencia estilolitos.
karstificadas,
15
4.1.1.2.
Senoniense- IC3-J-.
Formado
por
m
200
de
y
blancas
margocalizas
rosadas
con abundantes Globotruncanas.
4.1.2.
Terciario.
4.1.2.1.
Eoceno Inf~Medio
Está constituido tonalidad
verde,
biolititos.
4.1.2.2.
por
(C-E).
Boceno Medio-Sup
duros
de
lentejones
dolomitizadas.
Se
bioclásticas fauna
encuentra
masivas
calizas
de
M.
calcarenitas
4.1.2.3.
de
(E).
pararrecifales,
Alveolinas,
margosas
fauna de Globorotalias.
100
de
arcillas
de
veces
a
presentando
Existe abundante
Compuesto
m.
100
de
y
biomicritas Nummulites,
etc.
Olicoceno
Constituido calizas margosas
(QZL.-
por a
550
veces
m.
de
margas,
detríticas,
margocalizas,
areniscas
de
matriz
16
arcillosa
y
cemento
calizo
algún nivel de conglomerados
A
de
techo
la
serie
incluso
e
se
localizado
ha
intraformacionales.
afloran
180
de
m.
calcarenitas,
calizas margosas y margocalizas blancas.
4.1.2.4.
constituyen
Lo
fauna
máximo
un
abundante
muy
Moluscos,
Corales,
cambia muy
de
cuenca
rápidamente
esa
en
de
100
de
e
m.
y
potencia
Algas.
dirección.
La
de
la
posición
calizas
arcillosas
Briozoarios,
Esta
hacia el
aumento
y por
de
incluso
Amphisteginas,
Equinodermos
su carácter arrecifal la
(M
a veces levemente detríticas
claras, con
Mioceno Inferior
formación debido
SE,
profundidad
a de
estratigráfica
puede atribuirse al Aquitaniense.
4.1.2.5.
Burdicaliense
Entre
0
y
40
"TAP
m.
de
lateral,
En
el
techo
un
máximo
del de
30
(M21-1-
margas
blancas
intercalados,
niveles de calcarenitas y continuidad
1"
escasos
con
de pequeña potencia
con una microfauna característica.
"tap
V'
m.
de
puede
distinguirse
calcarenitas
que
en
ocasiones
intercalan
nivel margoso de potencia no superior a los 20 m.
un
17
4.1.2.6.
Mioceno Medio
Constituido
por
bioclásticas,
se
contiene
ha
4.1.2.7.
De alcance
de
m.
distinguido Algas
un nivel
a
con
siltosas, de
veces
delgadas
calizas calizas
y Anphisteginas.
La
margosas grises edad
y
que
puede
al Serravaliense.
Tortoniense-W5-J-
potencia varios
cientos
de
continuidad lateral.
carácter
cabalgamiento,
probablemente
aunque
indeterminada,
constituido
metros,
con escasos niveles detríticos
De
calcarenitas,
de
finas
margas
abundantes
atribuirse
200
areniscas
intercalaciones margas,
(M31-
o carbonatados
por margas limitada
de
Esta formación se denomina "Tap 2".
claramente
conservada
bien
aparece
manto
post
0
post-
en
fosas
tectónicas.
4.1.3.
El escaso
Cuaternario.
Cuaternario desarrollo,
dentro
de
la
tratándose
zona de
estudiada
tiene
terrenos
con
granulometría heterogénea y origen aluvial o coluvial.
18
4.2.- TECTONICA.
4.2.1.- Estructuras tectónicas.
La
región
plegado
y
fracturado,
longitudinales que
están
a
comprende
siendo
por los
un
las
gran
fallas
conjunto
principales
las direcciones de plegamiento,
genéticamente
accidentes obstante,
estudiada
una
misma
límites
están determinados
relacionados
NE
dirección y
SO
de
debido
ambos de
la
a
tipos
de
esfuerzos.
No
unidad
por importantes accidentes
geológica
de
carácter
transversal y dirección NO-SE.
Las estructuras importantes son:
Estructuras las
sierras
anticlinal,
de
de la Carrasqueta:
Cuartel
cuyo
flanco
importantes fallas un
sinclinal
sector
y
de
Ocupan el núcleo
Carrasqueta.
septentrional
Se
está
trata
de
afectado
de un por
directas que hacia el Sur evoluciona a
ocupado
meridional
por
está
las
calcarenitas
cabalgado
por
del los
M3.
En
el
materiales
eocenos.
Cabalgamiento NE, los
pone
en contacto
materiales
cabalgamiento
está
Jijona-Penáguila. los materiales impermeables interrumpido
Con
dirección
calcáreos del al
eocenos
Mioceno. NE
por
SO-
la
con Este
falla
Cocentaina-Penáguila y al oeste por la falla del Tibi.
19
4.2.2.
Pases de deformación.
Dentro dos
de
Carrasqueta
la Unidad
pueden
distinguirse
fases tectónicas mayores:
La inició
primera
durante
plegamiento sucesivos.
de
fase
el
es
de
Terciario
fondo
que
comprensivo
forma
en
origina
Esta fase alcanza
Inferior-Medio,
carácter
de
surcos
y
su clímax durante
desarrollándose
un
y
se
suave
umbrales el Mioceno pliegues
grandes
apretados cuyos flancos aparecen fuertemente cabalgados.
La condiciona Mioceno
segunda cuencas
Medio.
fase
es
miocenas
Los bordes
de
de muy
carácter subsidentes
estas
cuencas
sinformales limitados por importantes los
que
frecuentemente
halocinéticas de Triásico.
se
fallas
producen
distensivo a
y
partir del
coinciden directas,
con en
inyecciones
20
21
5.1.
CARACTERISTICAS LITOLOGICAS.
5.1.1.
Formaciones permeables.
susceptibles
formaciones
Las transmitir
cantidades
de
agua
de
almacenar
y
hidrogeológicamente
interesantes son:
Dolomías Estas
(C2).
calizas
y
calizas
son
las
del que
Cenomaniense-Turoniense constituyen
el
acuífero
Negre.
Calizas Sup.
masivas
pararrecifales
del
Eoceno
Medio-
(E).
Calizas pararrecifales
del Mioceno
Inferior
(M1).
Constituyen el acuífero de Jijona.
en ocasiones bioclásticas
Calcarenitas Medio
( M3)
5.1.2.
forman el acuífero de Carrasqueta.
Formaciones
Las
del Mioceno
impermeables.
formaciones
impermeables descritas son:
Margocalizas blancas del Senoniense
(C3).
22
margosas
Arcillas Inferior y Medio
margocalizas
y
el
perarrecifales
tonalidad
del
Eoceno
Oligoceno
(02)-
verde
(C-E).
Margas Constituyen
de
impermeable del
Mioceno
de
del base
Inferior
las
de
(M,)
calizas
que
constituyen
(M2)
constituyen
el aculfero de Jijona.
Margas el
del
impermeable
Míoceno
"Tap
de techo
Inferior
calcarenítas
V',
(M,)
del
Burdigaliense
de y
las
impermeable
el
Míoceno
calizas pararrecifales del
Medio
de
base
que
(H3)
de
albergan
las el
aculfero de Carrasqueta.
Margas de
techo
de
del
las
"Tap
2",
(M5).
Tortoníense
calcarenitas
que
Impermeable
constituyen
el
acuífero
de Carrasqueta.
5.2.
CARACTERISTICAS GEOMETRICAS.
La dos
unidad
de
acuíferos
Carrasqueta.
de
Se
del anticlinal
Jijona-Carrasqueta
sitúa
diferente
está
litología:
estructuralmente
de La Carrasqueta,
cabalgamiento de Jijona-Penáguila.
compuesta
estando
en
Jíjona el
flanco
por y SE
limitada por el
23
5.2.1.
Aculfero de Jijona.
La
roca
acuífera
Mioceno
Inferior,
de base
y la
En
sector
el
formada
actuando
formación
calcarenitas estas
está
el
"Tap
vértice
1"
por
Oligoceno como
las
calizas
como
impermeable
impermeable
Peña-rroja-Manantial
del
Serravaliense
calizas,
existiendo
están
(M3)
del
de
techo.
Rosset
las
conectadas
probablemente
a
cierta
comunicación hidráulica entre ellas.
Los
límites
hidrogeológicos
de
la
formación
M,
pueden establecerse del modo siguiente:
Al 0,
entre
del Machet,
Desde
el vértice Peñarroja y la zona N del Cabezo
por el afloramiento del
ahí
hacia
caracteristicas, encontrarse
el
pero
NE no
impermeable de base.
continúa es
observable
por
recubierto
límite
el
en
con
estas
superficie
por del
calcarenitas
las
Serravalliense(M3).
Todo
límite
el
condicionado
el
por
meridional cabalgamiento
Torremanzanas,
que
acuífera
formaciones
con
pone
y
suroriental
de
la
lateralmente
unidad
en contacto
arcillosas
del
está
Jijonala
roca
Mesozoico
y
Terciario.
El Km2.
aculfero
así
definido
tiene
una
superficie
de
35
24
5.2.2.
Acuifero de La Carrasqueta.
La
roca
acuífera
está
calizas
del
actúa
fundamentalmente
formada
por
Como
Serravaliense(M3) . la
calcarenitas
impermeable
serie
de
y
base
arcillosa
del
Burdigaliense.
Los modo
límites
hidrogeológicos
establecerse
del
siguiente:
Al NO por el afloramiento arcillosa
del
septentrional sector
de
tramo
lo
por
de
las
base,
en
por el
del
que
aunque
suficientemente
y E
reforzado
margas
límite
para poder actuar como
S
o subafloramiento
Burdigaliense,
este
impermeable
Al
pueden
en
Oligoceno.
no
se
de la serie
es
en
arcilloso
zona
Existe
un
cartografiable
el
podido
ha
la
el
constatar
un
Serravalliense
límite.
afloramiento
de
la
serie
arcillosa del
Burdigaliense.
Al
NE
el
hidrogeológico. anticlinal la
del
Carrasqueta,
flujo
Este eje que
subterráneo.
subterráneas 5004,
límite
drenan
mientras
que
corresponde
se
umbral
ocasionado
está
transversal origina
al
una
sinclinal
clara
Hacia
el
SO
hacia
el
manantial
hacia
el
NE
con
del
la
umbral
un por
un
suave
principal
de
divergencia
en
umbral
las
aguas
Nuches
2933~
de
infiltración
que
el
se
25
produce el
en
la calcarenita serravallienses
manantial
definido
5.3.
del
Molinar,
2932-5002.
se El
dirige
hacia
aculfero
así
tiene una superficie de 7,7 Km2.
CARACTERISTICAS PIEZOMETRICAS.
5.3.1.
Aculfero de Jijona.
Los Rosset
puntos 28338008
28338033,
representativos
son
y
28338011,
sondeos
Ayuntamiento
manantial
de
los
sondeos
que
fueron
construidos
por
el
para
utilizarlos
en
el
Jijona
de
el
28338028
y
abastecimiento público.
El
ITGE
controlando Diciembre
las
A
del mismo
de
evolución totalmente
este
consecuencia el
suavemente
para
cota
cota
orden de Por
sondeo
del
la
lo
Alicante
de
del
2833-8028
agua
se
vienen desde
situaba
la cota del manantial tanto,
piezómetro
el
análisis
puede
de
a de la
considerarse
representativa del acuífero Jijona.
acuífero,
la
cuando
m. s. n. m.
(600
Provincial
piezometrías
1976,
de
607 m.s.n.m., Rosset
Diputación
y
del
nivel
de
la
del
situarse
manantial.
explotación agua
en
que
que
se
sondeo
el
a mediados de Puesto
a
somete
al
desciende
1978 por debajo
de
circunstancia
se
esta
26
mantiene
durante
manantial
comienzos
de
Sin
inicia
un
1980
seco, Las
m),
de
se estabilicen
aprecia
muy
precipitaciones origina
un
descendente. otoño
de
Sin
ese
1988
el
recuperación,
del
fuertes lluvias recuperación
del
orden
de
80
período
del
orden
para 1991.
de
30
recuperar A
esta
año
1990
m
entre
nuevamente
el
nivel
Julio el
recuperación
en
como
ha
escasez
de
la
que
evolución
registradas
en
recuperación
una
fuerte
respuesta
a
las
1987.
conducido
al
La
nivel
1980.
piezométrico
de
nivel la
caídas
100 m).
m.
1986-88
el
Septiembre),
experimenta
agua
hasta
(60 m).
de
fortísima
(casi
se
se llega a
acaecidas a partir de Octubre de
del
el
e incluso
a
1980
lluvias
efecto
piezométrico al máximo histórico de
Durante
las
lluvias
una
de
a
del
vuelve
descenso
(hasta
las
originan
nivel
ascenso
primavera
pendiente
piezométrica en ese mismo año
En
el
1986 la
producen
manantial
cuando
embargo,
año
la
claramente
de
se
del nivel del agua
durante
aumento
el
los niveles
producir un fuerte aumento
Se
y
el
aportaciones
que
ininterrumpido 1982,
período
espectacular
partir de
casi
y
50
este sus
lluvias
(unos a
en
con
un
embargo,
nuevo
año,
originan
invierno-primavera hacen que
un
reguladas.
piezométrico
surgir.
de
permaneció
completamente
nivel
más
1990
y
Enero
anterior sucede
desciende
en un
de
1991
Mayo
de
descenso
27
progresivo m.s.n.m.
del en
alcanza los de 70
m.
nivel
Mayo
piezométrico
de
1991
521 m.s.n.m.,
Durante
hasta
lo que
desde
Mayo
de
la
cota
1992
en
supone un descenso
el mes de Mayo de
1992,
593 que anual
se ha registrado
un ascenso de 36 m en el nivel piezométrico.
La sondeo el
observación 2833-8028
manantial
en que
de
todavía
de
la
permite Rosset
vierte
conocer
niveles
los
permanece agua.
cota del nivel piezométrico
Los
evolución
Esto
piezométrica
períodos
seco
y
último
en
en
los
sucede
los
del que
períodos cuando
la
es superior a 600 m.s.n.m.
piezométricos
en
fecha
27-V-92
se
indican a continuación:.
SONDEO
PNP(E)
COTA(
28338011
81,77
518,23
28338028
78
557
28338033
101,5
m.s.n.w.)
553,5
Sondeo mal nivelado.
5.3.2.- Acuifero de la Carrasqueta.
En
este
representativo (29335004)
aculfero
el
corresponde
de cota 610 m.s.n.m.
único al
dato
manantial
piezométrico de
Nuches
28
5.4.- PARAMETROS HIDRODINAMICOS.
5.4.1.- Aculfero de Jijona.
El
cálculo
de
basado en un bombeo sondeo 1978
Sereña
nivel
de ensayo
(28338028)
un caudal
con
un
de
58
hidrodinámicos
se
realizado por el ITGE
entre
29-11-1978
el
l/seg.
y nivel
en
el
sondeo
valor
de
transmisividad
estático
deducido
parámetros
los
era
de
76
final m.
en
De
torno
y
el
ha
en el 2-12-
a 95
m.
El
aquí
se
ha
a
los
300
el valor de porosidad eficaz aunque
por
M2/día.
Se desconoce similitud
con
otros
acuíferos
próximos
se
estima
próxima
al 2%.
5.4.2.- Aculfero de Carrasqueta.
Se debido el
de
desconocen a
la
sus
similitud
Jijona,
se
parámetros de
podrían
hidrodinámicos.
características considerar
para
Aunque
litológicas él
válidos
parámetros estimados para el mencionado aculfero.
con los
29
S.S.- CAPTACIONES EXISTENTES.-
5.5.1.
Aculfero de Jijona.
Las el
captaciones
volumen
de
que
explotan
extracciones
(año
el
aculfero,
1991)
y
el
así
uso
a
como
que
se
destina el agua se expresan en el cuadro adjunto.
Nº
INVENTARIO
(»3/año)
VOLUMEN
uso
28338011
43.356
Abastecimiento urbano
28338028
815.670
Abastecimiento urbano
28338033
201.200
Abastecimiento urbano
El Sondeo Rosset 210
m.
Norte
y
está
situado
la población
de
reserva,
por
(28338011)
tanto
en
el
es
Castalla,
al
utiliza
como
de
de
una
Barranco
de Jijona. no
una profundidad de
tiene
Se
de
explotado
pozo
manera
continua.
El 570
m.
en el
Sondeo capta
Sereña las
Paraje de
(28338028)
calizas
del
con
Mioceno
Sereña y sobre él
una
profundidad
Inferior,
recae
se
de
sitúa
la mayor parte de
la explotación.
Por
último
la Casa de Los
el
sondeo
Pinetes,
Pinetes
tiene
(28338033)
una profundidad
situado de
300
en
mts.
30
Atraviesa
la
formación
tras perforar
acuífera
las margas
del
de
"Tap
Carrasqueta
1"
alcanza
en
la
seco
y
formación
aculfera de Jijona.
5.5.2.
Aculfero de Carrasqueta.
No punto
existen de
sondeos
que
exploten
representativo
agua
manantial de Nuches
este
acuífero.
únicamente
Como
existe
el
de
la
(29335004).
5.6.- FUNCIONAMIENTO HIDRAULICO.
5.6.1.
Aculfero de Jijona.
La
alimentación
infiltración permeables
Dado régimen Roset, se
hará
relación bombeos.
lluvia
acuífero
sobre
los
procede
afloramientos
de
rocas
(9,5 Km2).
que
no la
de
del
no
existen
influenciado estimación
mediante funcional
un
de
datos por
los
modelo
entre
la
la
descarga,
en
del
manantial
de
sobre
bombeos, recursos
medios
estocástico piezometría,
del
acuífero
basado
en
pluviometría
la y
31
Este epígrafe
modelo 11
valor de
de
desarrollado este
con
informe,
los recursos
todo
ha
detalle
permitido
medios para el período
en
el
estimar
el
1980-1991
en
0,983 hm3/año.
Este
alimentación la
correspondería
valor
directa
transferencia
por
fundamentalmente
infiltración
vertical
de
lluvia,
por
directa
el
a
la
ya
que
acuífero
de
del
de
Carrasqueta no parece probada.
En
este
lluvia,
caso
caída
tasa
la
sobre
sus
permeables,
materiales
infiltración km2
9,5
ascendería
período(1989-1991).
del
de
El
él
de al
se
agua
afloramiento para
25%
la
de
media
recogieron
en
la
estación meteorológica de Jijona 382 mm.
La
descarga
del
aculfero
se
de
realiza
dos
modos
diferentes:
a)
A
través
2833-8008;
manantial
del
del
que
se
dispone
solamente de dos medidas de aforo.
- 30-11-76.
Caudal
- 31-10-86.
32
Antes
del
lls.
inicio
caudal
fluctuaba
llegar
a
alcanzar
considerablemente
inapreciable.
de
mucho 50
la
explotación
entre lls)
y
del
la
primavera
el
verano
sin llegar a agotarse.
que
acuífero, (donde se
el
podía
reducía
32
Con la
tan
escasos
descarga
explotación.
del
datos
manantial
Igualmente
descarga a partir de caudal
son muy
resulta antes
resulta
1986,
acusadas
imposible
inicio
imposible
puesto
y el
del
que
las
período
cuantificar de
la
cuantificar
la
oscilaciones
de
existente
entre
los
dos aforos disponibles es demasiado dilatado.
b)
A
través
de
la
explotación
2833-8028 y 2833-8033, hm3
destinados
que
de
en el
íntegramente
los
sondeos
1991
año para
2833-8011,
extrajeron
1,06
abastecimiento
de
Jijona.
5.6.2.
Aculfero de Carrasqueta.
La
alimentación
infiltración
de
del
lluvia
la
Suponiendo una tasa de aculfero de Jijona,
25%
la
este
valor'medio.
sobre
procede
sus
de
la
afloramientos
de
(7,7 Km2).
rocas permeables
alimentación
acuífero
de
Este
de
con
la
un
se
similar a
lluvia total,
acuífero
volumen
Nuches(2933-5004),
infiltración
es
drenará
caudal
de por
medio
se estima
que
hm3/año
de
manantial
de
0,73 el
la del
continuo
de
unos
23 l/seg.
La Nuches
descarga
del
(2933-5004),
aculfero
se
controla
a
través
del
manantial
recientemente.
Se
de
dispone
34
5.7.
BALANCE HIDRICO.
5.7.1.
El
Acuífero Jijona.
balance
(1981- 91)
del
acuífero
para
la
media
(h®3 /año) .
*Infiltración de
Salidas
0,983.
lluvia
( hm3 /año )
( Año 19911 .
1,06.
* Bombeo
- 0,077
-Variación de almacenamiento
Admitiendo lluvia,
las
existiría
cifras un
5.7.2.
por
entradas
de
desequilibrio
considerando las extracciones
�-
periodo
puede expresarse del modo siguiente:
Entradas
de
del
de
de
infiltración
-0,077
hm3/año
1991.
Acuífero de Carrasqueta.
Para equilibrado
este para
acuífero el
el
periodo
balance estudiado,
está
claramente
puesto
que
las
35
entradas
por
infiltración
de
lluvia
deben
las descargas a través del manantial
La del
recarga
período
medias valor
del en
es
del
orden
analizado. manantial
principio
En
de
de
función
coherente
iguales
a
de Nuches.
HM3/año,
0,73
Nuches
ser
de
se
esto
estiman
con
los
como
media
las
salidas
23
l/seg.,
en
aforos
de
que
se
de
Jijona
dispone.
5.8.
GRADO DE EXPLOTACION DE LA UNIDAD.
5.8.1.
Aculfero de Jijona.
Las superan
extracciones a
las
actuales
entradas
por
del
infiltración
de lluvia en un valor próximo a 0,077 cabría
considerar
que
éste
aculfero
está
directa
Hm3/año.
actualmente
del
agua
Por tanto, ligeramente
desequilibrado.
5.8.2.
Aculfero de Carrasqueta.
Este
presenta
aculfero
excedentarios,
que
se
drenan
actualmente de
forma
cuantificados con actualmente con detalle.
recursos
natural,
no
36
5.9.
RESERVAS EXPLOTABLES.
5.9.1.
Aculfero de Jijona.
Las
reservas
actuales del
se
muro
han
del
movilizables
deducido
nivel
a
desde
partir
aculfero
y
del de
las mapa
la
captaciones de
isohipsas
posición
de
la
superficie piezométrica.
Suponiendo y una porosidad de
potencia
eficaz
del
del
1%,
nivel
el
acuífero
de
70
m.
volumen movilizable
es
10 Hm3.
5.9.2.
Aculfero de Carrasqueta.
Del se
una
han
mismo
modo
calculado
1,15 Hm3.
descrito
las
para
reservas
en
el
aculfero
éste,
de
Jijona,
estimandolas
en
700
900
800
l�p -CP qFP
sol¡ solo 0009
C?
8008
-500
0028 8033
.000
Jijona
ESBOZO ISOHIPSAS DE LA BASE DEL MIOCENO INFERIOR
40 0
1
Escala: 1150.000
2
3Km.
ACUIFERO DE JIJONA
-
�o
0
2 Km.
Escala: 1/25.000
e0o
100
400
500
5004 6OO
ESBOZO ISOHIPSAS DE LA BASE DEL M3 SERRAVALIENSE ACUIFERO CARRASQUETA
37
38
6.1.
ACUIFERO DE JIJONA.
La
Diputación
periódico
en
obtenidos
de
el
entre
225
y
Alicante
sondeo
este
bicarbonatada
de
300
entre 325 y 375 regadío C2-S1.
2833-8028.
control
cálcica, mg/1
la
ohmios/cm.
un
Según
hidroquimico,
el y
realiza
residuo
conductividad
los
la
seco
muestreo datos
facies
es
tiene
valores
alcanza
valores
La clasificación de aptitud al
Desde un punto de
vista sanitario,
el
agua
es potable con respecto a los parámetros analizados.
6.2.
ACUIFERO DE CARRASQUETA.
El sólidos
agua
disueltos
bicarbonatada C2-S1.
manantial
del
Desde
el
cálcica un
punto
de
Nuches
5-11-86. y de
una
Es
tenia un
agua
clasificación
vista
sanitario
412
se
para
mg/1 de
de
tipo
regadío
trata
agua potable atendiendo a los parámetros analizados.
de
un
39
E:"SJPLI-Oupkc Mom W10
^C�U^
OE:
IVOICaumIsmno
40
cauce
vertidas
son
Jijona
residuales
aguas
Las
del
sin y
Cosc6
río
al
ningún
tipo
Barranco
parcialmente
son aprovechadas
en
generadas
aguas
del
población
de
tratamiento
al
la
de
Purgatorio.
Estas
abajo para regadío de
pequeñas huertas.
M3/año
850.000 pequeño
cuya
incremento
Diciembre
como
anualmente
generado
volumen
El
en
meses
consecuencia
del
de
octubre,
mayor
orden
de
muestra
un
del
mensual
distribución los
es
Noviembre
consumo
y
producido
por las empresas en la elaboración de turrones.
Existe
planta
una
funcionamiento;
se
físico únicamente.
trata
depuradora de
una
no
que
planta
de
está
en
tratamiento
.-1
H
z
z
1
1
1
1
i
J
1
1
1
42
8.1.
USO DEL AGUA.
El
agua
los sondeos
extraida
del
28338011,
acuífero
28338028
Jijona
de
y 28338033
abastecimiento urbano
de
la población de
en
de
parte
el
abastecimiento
través
a
de
se emplea en el así como
Jijona,
del
consumo
es
de
industrial
de
la ciudad.
La durante en
población todo
que
es
que
de
del
350
por
distribución industrial
orden
dos lado
un
las
de
los
meses
muy
y
Diciembre,
en
que
en
perdidas
a
aunque
el
de
este
red
la
otro,
elaboración
llega
dotación
alta
concentra principalmente en los meses
Noviembre
verano
provocan
que
y por
de
La
10%.
dotación
30%)
habitantes,
8 .000
un
factores
en un
empresas
las
de
de
1/hab/día,
(cifrados
de
excepción
a
considerar
incremento;
éste se
año,
disminuye
aplicada hay
el
abastecida
de
consumo
turrones,
de Octubre,
incrementarse
un
25% o 30% el consumo mensual de la población.
El
agua
riego del valle
La -
Nuches, es:
manantial
del
del
superficie o
en
su
de
Nuches
es
empleada
para
río Coscón.
regada,
defecto
con
con los
agua sondeos
del del
manantial
de
Ayuntamiento
43
Cultivos
has.
Frutales
10
Alimendro
so
Hortalizas
15
Total
El
volumen
Privadas de 350.000
M3�
regadio
de
anual
75
inscrito
en
Registro
el
la Confederación Hidrográfica del que 168
atienden has,
de
una
las
superficie
Aguas
Júcar es de
escriturada
actualmente
que
de
se
de
ríegan
75
has.
Existe
un
acuerdo
el Ayuntamiento por el se
seca,
entre
Comunidad
la
de
Regantes
cual cuando el Manantial
el agua necesaria
para
el
regadio
es
y
de Nuches donada por
el Ayuntamiento a los regantes.
Cuando conducidos Diputación Jijona.
existen al de
sondeo
excedentes Sereña,
Alicante
para
en
donde
el
manantial
son
recargar
estos
utilizados el
por
acuífero
son la de
44 8.2.~ ESTIMACION DE-LA DEMANDA.
Teniendo m3/año
para
Jijona,
y
en el
cuenta
la
la
población
consumo
red,
consumo
abastecimiento
dotación real es de
El
el
180
de
demanda de
525.600
urbana
de
hab.,
la
8.000
1/hab/día.
de
cifrado
éste
concentrándose
la
existente
industrial,
está
de
estimado
las
en
básicamente
empresas
unos en
conectadas
M3/año,
204.400 último
el
a
trimestre
del año.
Las
pérdidas
estimado
en
un
en
30
la
red
de
lo
que
supone
%,
distribución un
han
se
volumen
anual
cercano a los 300.000 m3.
8.3.
RELACION APORTACION-DEMANDA.
En base aculfero embargo
a los
de
Jijona aculfero
el
regulados
recursos
no
Alicante
está
valor
de
podrían
servir
ligeramente de
en
excedentes para
se debe considerar
actualmente la
la
recarga
utilizados
compensar
del acuífero
parece
disponer
de
Diputación
de
primero.
El
del en
sobradamente
de Jijona,
el Sin
desequilibrado.
Carrasqueta
que
empleando
estos
déficit actual Hm3/año.
datos existentes
la el
estimados
recarga ligero en 0,077
45
46
La
distribución
del
agua
del
manantial
cuenta con una primera balsa situada en con capacidad de
Desde
ahí
capacidades ,
de
Nuches
sus proximidades,
1.100 m3.
vierte
se
una
de
1770
1.200 m3 y tres de 450 m3.
a
otras,
m3,
dos
con de
siguientes
las
m3,
2.000
una
de
es de 9.420
La capacidad total
m3.
Actualmente Alicante, con
los
la
se
recarga
sobrantes
del
para ello el sondeo de
realiza, artificial manantial la Sereña.
por del de
la
Diputación
acuífero Nuches ,
de
de
Jijona
utilizandose
47
IL 0
14001:tmplo
]DE:
1EXE>T_a«Dor^c JECYm
48
10.1.
NORMATIVA GENERAL.
10.1.1.
Sectorización.
La los
normativa
dos
que
acuíferos,
se
adopta
Jijona
mismo
uso
población
es
el
de
agua de
El
directamente
por
captaciones
mientras que
en el
de
de
fuente
las
informe
primero
de
ambos se
ellos
de
son
el
propiedad
aprovechados
como
destina a
urbano
es
aculfero de Carrasqueta
Nuches
incluye
tratándolos
abastecimiento
Jijona.
la
de
que
este
y Carrasqueta
un sólo conjunto puesto que el un
en
a
la
explotado municipal,
los excedentes
para
recargar
el
aculfero de Jijona.
10.1.2.
El los
Volumen máximo de explotación.
volumen
recursos
hm3/año,
más
hm3/año.
Dado
máximo
estimados los que
de
para
del el
explotación
manantial
el de
aculfero
para
de
de
de
Jijona,
Carrasqueta
de
recargar
los el
comprender
Carrasqueta,
recursos,
aprovechamiento Nuches
acuífero
aculfero
de capacidad para regular sus continuar
el
deberá
sería
no
0,73
dispone
conveniente
sobrantes
aculfero
0,98
de
del
Jijona.
49
ligero desequilibrio
Esta recarga mitigaría su
recursos regulados en la unidad.
tiempo aumentaría los
unidad,
las
aumentar
Cabría
el
en
si
previamente
abastecida
actualmente todos
se
depuradas, por
recursos
los
hídricas
disponibilidades
futuro
residuales,
Entonces,
y al mismo
utilizasen
las
en el riego
de
el
manantial
de
Carrasqueta
de
la
aguas la
zona
Nuches.
de
podrían
ser
regulados en Jijona.
10.1.3.
Captaciones de escasa importancia.
Se considera definido
por
los
deberán autorizarse hasta el
que
recursos
renovables
y
límite
explotación
la
actual.
10.1.4.
Prescripciones
generales de orden técnico.
establece
de
algunos
elementos
constructivos
permitan
Estas
adecuado
el
ejercer
concesiones,
forma
obligada
Se
especialmente deberían
de
inclusión
rango
menor,
que
sobre
las
control
en sondeos y pozos.
equiparse
con
los
siguientes
elementos:
- Tubo piezométrico. - Contador. - Espita de toma de
de
la
agua a la salida de la tubería.
so
Los
facilitar,
parte
perseguidos
objetivos
posteriores
desarrollo
el
de
dos
por
tipos,
estudios
de
otra,
por
y,
son
y
una
trabajos los
verificar
poder
condicionantes impuestos a la concesión. orden de cosas y para el
En otro
sondeos,
y
pozos
general
obligatoriedad
la
actuaciones
Si
la
presentase
especiales
oportunas de
la
fuese
negativa
y
no se
hidrogeológicos,
terreno
con
las
hermético
problemas
relleno
su
a
a
cierre
perforación
captación.
procederá
el
requerirán
que
proceder
de
de
carácter
con
también
imponer
deberá
se
abandono
caso de
cementando
natural,
la misma.
solamente el tramo superior de
ZONIFICACION.
10.2.
Dado acuífero
el
carácter
de
Jijona
la
concesión
desequilibrio
debiéndose
concesionales,
salvo
para
propone
como ZONA NO AUTORIZADA,
ella
en
permitir
no
autorizaciones
nuevas
las correspondientes
el
las normas se
explotaciones
nuevas
clasificar toda su superficie
exclusivo
establecer
a la hora de
de
presenta
que
casi
uso
su
y
abastecimiento urbano, para
de
a explotaciones
de escasa importancia.
Existen Estas
son:
industrial
otras que
en
razones la
insatisfecha
que
motivan
actualidad en
la
esta
existe
población
de
propuesta.
una
demanda
Jijona
que
51
obliga a
importar agua del acuífero
prevista población la
de
demanda
los
ejecución
la
Jijona,
de
lo
un
polígono
cual
población.
la
escasos
de
recursos
de Maigmá,
industrial
producirá Todo
ello
disponibles
y que
un
está
en
la
incremento
en
obliga
a
reservar
abastecer
para
la
población de Jijona.
10.3.
PERIMETRO DE PROTECCION.
Existe
un
abastecimiento
Jijona,
a
partirá
(torre),
en
del
alineación
punto
intersección
minutos
Este
treinta
y
geodésico
Murcia,
por
(torre)".
8'
(38º
recta punto
E),
el
sondeos
B.O.E.
nº
de 263
en
Iglesia
de
el
vértice
geodésico
hasta
33' el
SO"
y
tres
N) ,
vértice
alineación desde
este
punto, de
114
Alicante,
a Valencia,
y desde
al
vértice
geodésico
ocho
paralelo cincuenta punto
en
Zurdo;
el
vértice
alineación
en la
el
Pozo
hasta
kilométrico
recta,
este
geodésico
recta,
hasta
el
minutos
desde
Jijona
grados
con
hito
el
recta,
tres
Madrid,
de
treinta
hasta
Montagur;
hasta
alineación
los
geodésico
meridiano
del
grados
Norte
este
recta,
en
de
en alineación
este punto
(3º
ocho
alineación desde
publicados
recta,
desde
segundos
protección
vértice
Peñarroja; de
de
definido de la forma siguiente:
de 3/11/77,
"Se
perímetro
carretera
este
Iglesia
punto, de
de en
Jijona
52
Las
coordenadas
del
perímetro
oficial
quedan
definidas del siguiente del modo:
PUNTO
x
y
Z(W)
A'
720,450
4276.950
1182
B*
714.450
4268.200
1226
C,
713.050
4271.500
1020
D'
722.700
4275.850
1081
El
717.350
4269.800
424
for
717.400
4273.000
860
El
propuesto
perímetro
anterior
e
incluye,
afloramiento de
las
calizas
parte Sur
la zona bajo
acuíferas
y
margas
sobre
blancas
en
las
del
en su del
la cual que
se
este
informe
parte
Norte
Mioceno se
Burdigaliense.
definido es el siguiente:
El
parte
Inferior y en
encuentran
disponen
engloba
las
al del la
calizas
discordantes perímetro
las así
53
PUNM
x
y
z(2)
A
717.350
4269.800
424
B
714.500
4267.350
740
c
713.200
4269.700
1000
D
713.050
4271.500
1020
E
715.250
4274.700
1040
p
718.950
4276.550
1010
G
720.450
4276.950
1182
H
723.900
4278.000
980
1
722.700
4275.850
1081
54
2-2--
MC>E>1011.C>
loorrc)clksrrjccc>-
55
11.1.- METODOLOGIA.
El de
objetivo
métodos
del
estocáticos,
comportamiento del conocidos sus
trabajo
de
de
la
un
y
renovables
obtención,
modelo
aculfero ante
pluviometría
recursos
es
que
por
medio
reproduzca
el
la influencia de valores
bombeos,
ante
y
permita
distintas
estimar
hipótesis
de
recarga por lluvia.
Matemáticamente, una
caja
salidas
negra de
con
agua,
funcional
si
unos
el
que
se
considera
valores
problema
explique
piezométrico como
acuífero
conocidos
será la
respuesta a
el
de
encontrar
como
entradas la
variación
relación nivel
del
las variables
y
de entradas
y
salidas.
En el
caso presente
estocástico
entre
rechazándose debido
a
las
por
va
piezometrías,
cualquier
que
obtención,
se
otra
que
buscar
un modelo
pluviometrías variable
utilizadas
lo
a
son
de
permitirán
de
y
lineal
bombeos,
influencia,
fácil
y
rápida
la
realización
se
ajusta
de
rápidas predicciones.
Hay
que
únicamente calidad
señalar
que
acuerdo
con
de
de
conclusiones
la lo
misma, será
y en
el la que
modelo
información la
función
si
sirve
disponible
universalidad de
y
los
y
la
de
las
datos
del
56
período
elegido
mismos
todas
acuífero.
para
la
las
Es
contienen
características
decir,
representativa de
modelización
particulares
muestra
la
en
es
sí del
suficientemente
la población.
11.1.1.- Pases de desarrollo.
Para
la
del
realización
presente
trabajo
se
en la modelización
seguido la metodología siguiente
ha
de
la
piezometría en función de la pluviometría y los bombeos:
Evaluación y análisis de la información disponible.
1
La estadística la (no
información
se
del
todo,
deterministas)
ello,
primer
en
básica,
ha
disponible,
de
basa
que
lugar,
en
llegar
y
de
evaluarse
y
porque
principios
a
y
expliquen
lo
a partir
una
lo
ella
se
Por
importante
la
y
información
a construir
va
de
generales
predigan.
muy
forma
analizarse de
una parte
cogiendo
que
el
modelo.
En
esta
detectan las
los
analizan
las
valores
extraños
(outlayers),
características
establece
particulares de
necesidad
la
transformación
series
se
fase
que
ajusta
los
de
cada
de se
a
la
se
observan
serie,
aplicarles
datos
datos,
y
se
aquella
distribución
57
normal, el
antes
ajuste
de
proceder
de
un
características ser aplicable
que
a
las
fases
modelo
los
siguientes,
lineal
datos
han
de
dicho modelo y ser útil
ya
exige
cumplir
que
ciertas para
poder
en su explotación.
2 - Identificación del modelo lineal.
Una
vez
realizados
evaluados sobre
transformaciones
precisas
lineal,
identificación
de
de
las
variables
y
en
tiempo
el
el
analizados
ellos
modelización
análisis
y
las
adecuación
ha
de
proceder
lineal,
que
establecimiento
existir
de
posibles
los tener
y a
la
a
consistente
puedan
puedan
y
depuraciones
su
dicho modelo
que
datos
para
se
relaciones
los
la en
entre
el las
retrasos variables
las
independientes.
El por
concepto
ejemplo,
manifestación se produce el
mismo
retrasos
que
puede meses,
retrasos
ante
como
un mes,
que
de
un
mes, otro
se
en
dos ejemplo
el
en
el
mes
tiempo
lluvias los
más
etc ...
quiere
un
mes
de
que
etc.. ante
dado
0
bombeo
dado
piezométricos
posterior,
o
cabo
de
meses
los niveles
al
su
o bién en
entonces 0
un
decir,
piezometría
tarde,
hablando
meses, es
en
niveles
produce,
afección sobre
manifestar etc ...
de
dos meses,
en
su
unas
aumento
mes
respectivamente; en un mes,
de
de
se 2
58
Para lineal
ello
se
en
la
fase
selecciona
características
y
de
el
identificación
que
más
peculiaridades
se
del
modelo
ajusta
a
las
observadas
con
la
el
paso
es
información disponible.
3 - Calibración del modelo.
Una vez mediante
identificado el
un
procedimiento
mínimos cuadrados, lineal con
que
las
A
del
de
modelo,
en
fase
de
los parámetros
de
dependiente
piezometría
consistente,
anterior
por
la ecuación
cálculo
de
ajuste
de
independientes,
parámetros la
siguiente
algorítmico
variable
proceso
los
seleccionada
la
variables este
calibración cálculo
calcular
relaciona
dos
bombeos.
modelo,
se en
la
pluviometría le
denomina
resumen, ecuación
de
y
en
el
lineal
identificación
del
modelo.
4 - Verificación del modelo y análisis de residuos.
Por último, proceder
a
metodología grado de
una vez
calibrado
la verificación del
ajuste
análisis de
la
de
mismo
en
modelo,
es menester
mediante
resíduos,
información
con los datos observados construcción del mismo.
del
el
para
la
llamada
conocer
predicha por
el
el
modelo
la realidad y empleados en
la
59
Como se
toma
la
consecuencia
la
vuelta
decisión a
la
predicciones
que
y por sea
lo
más
tanto,
de
11.1.2.-
Parámetros
en
evaluación
en
mejor
se
a
ajuste
en
de
calibración
ajusten
y útil
proceso,
depuración
de
fase se
este
la
el modelo
válido
la
de
algunos
ciertos
para
los
a
datos
obtener
datos
mejor
casos y
unas
observados,
la
realidad y
la explotación.
estadísticos
univariantes
análisis
y
la
de
utilizados información
disponible..
Para
la
disponible,
evaluación así
transformaciones ajustar
los
parámetros
necesarias a
las
lineal,
se
estadísticos
análisis la
como
datos
modelización
y
de
la
información
determinación
que
se
le
han
condiciones utilizan
de
de
las
aplicar
exigidas los
en
para la
siguientes
univariantes:
Media muestral 1 n
X
donde
n
es
el
valores
de
la
n E xi i-1
número
de
variable
en
los
datos,
cuestión
y
xi,
cada
observada.
uno
de
los
60
Varianza muestral
insesgada:
1 S2
n
-
-
n-1
Donde
n
es
muestral,
el
y xi,
E
(Xi
-
X)2
i-1
número
de
cada uno de
los
datos,
x
los valores de
es
la
media
la variable en
cuestión observada.
* Desviación típica
S
=
�«S2
* Asimetria n 3 E (xi - x) i-l
1 a n
s^3
a
es
el
valor
muestral
de
Pisher,
y
Donde
función
distribución
de
cuando
la moda se
menor
de
0
generales
se
(totalmente
coeficiente
del
valor
0
presentara
el
totalmente
simétrica,
encuentre a la izquierda de
cuando puede
lo
ocurra
la
aceptar
simétrica)
si
a,
sobrepasa el valor de 2-,/(6/n),
contrario. hipótesis en
valor
para el 95%
asimetría
de
de
para mayor
una de
0
la media,
y
En
líneas
normalidad
absoluto, de confianza.
no
61 Kurtosis
1
n E (xi i=l
k n
Donde
k
es
-
el
s^4
valor
del
4 - x)
coeficiente
kurtosis muestral de Fisher, distribución normal,
3
de
y presentará
será mayor que 0
apuntamiento
o
el valor 0 para
para distribuciones
apuntadas y menor de cero para distribuciones más planas. Se
acepta
si
k
en
en
valor
líneas
generales
absoluto
es
la
menor
hipótesis de
de
2N«24/n)
normalidad
para
el
95%
de confianza.
11.1.3.- Modelo lineal multivariante.
En
el
caso
concreto
presente,
indentificaci6n proporciona un modelo multiple, término
en
el
que
la
dependiente,
y
variable las
de
la
de
regresion
piezometría
variables
fase
de
de
lineal
sería
pluviometría
el y
bombeos los términos independientes.
La
forma
general
de
una
ecuación
de
regresion
multiple es la siguiente
y � cO + c1-x1 + c2'x2 +
---
+ Cn*xn
lineal
62
Donde presente
y
es
es
la
la
variable
dependiente
piezometría)
y
x1,
x2,
diferentes variables independientes bombeos
y
lluvia),
coeficientes
o
co,
c1,
parámetros
(en
...
(en el
son
%
las
son
Cn
la
caso
caso presente,
c2,,
de
el
los
ecuación
del
modelo
se
tres
series
lineal.
11.2.- MODELIZACION LINEAL.
11.2.1.- Información disponible.
Para
el
temporales, bombeos, las
ajuste
correspondientes
pluviometrías
fechas
piezómetro
este
de
1-1-86
y
y
los
a
de
de
observados cada
series
Estas
y
toman
valores
piezometrías
1-6-90.
28-33/8/0028,
modelo
los
mes
de
entre
provienen
del
mensuales
datos
procedentes de la estación pluviométrica de Jijona. Las cota
medidas
sobre
el
pluviometrías estación bombeos
piezométricas
nivel en
del
litros/m2
pluviométrica en m3/mes.
mar
Se
se
en
Alicante
de
durante
agua un
u observaciones.
La t de
54
callos
de
1
Y
aproximadamente 2.
un
mes
incluyen en el
54 casos -
expresan
95%
en
(m.s.n.m. ), recogida (mm/mes),
estudio
student
intérvalo
metros
de
de las
en
la
y
los
un total
calculada
de
Para
confianza
es
63
La hip6tesis de normalidad grupo
de
54
casos,
si
la
se puede
Asimetría
encuentran dentro de los siguientes
-0.66 -1.33
Es si
"a"
decir, en
valor
series
la
serie
de
absoluto
es
menor
o
de
la
Kurtosis
se
s a s 0.66 s k s 1.33
en valor absoluto es menor o igual
Las
y
para un
límites
admite
se
aceptar,
datos
que
se
se incluyen en la siguiente tabla:
como
datos igual
a
normal,
0.66,
y
"k"
a 1.33.
emplearán
en
el
estudio
64
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54
FECHA 01-86 02-86 03-86 04-86 05-86 06-86 07-86 08-86 09-86 10-86 11-86 12-86 01-87 02-87 03-87 04-87 05-87 06-87 07-87 08-87 09-87 10-87 11-87 12-87 01-88 02-88 03-88 04-88 05-88 06-88 07-88 08-88 09-88 10-88 11-88 12-88 01-89 02-89 03-89 04-89 05-89 06-89 07-89 08-89 09-89 10-89 11-89 12-89 01-90 02-90 03-90 04-90 05-90 06-90
PZ 532.70 523.30 517.40 524.80 519.70 504.40 508.00 492.50 502.10 531.30 569.80 559.90 580.80 583.85 577.90 573.10 583.10 570.60 575.80 559.30 564.80 563.10 587.90 615.80 637.20 643.55 639.50 645.20 648.80 653.20 640.40 621.90 615.70 606.10 593.30 577.70 589.30 603.40 621.30 616.60 611.40 609.22 601.80 608.30 634.40 630.50 639.30 648.80 659.20 658.93 656.40 660.30 657.70 646.37
PV 9.50 6.70 8.20 43.80 7.00 7.90 40.50 6.50 153.60 135.50 18.10 1.30 42.40 54.10 .30 1.20 45.30 .00 20.60 .00 36.20 90.40 191.40 62.10 40.40 30.30 10.00 60.70 78.50 86.70 3.40 5.00 54.00 61.00 45.00 .00 24.00 20.00 81.00 34.00 52.00 36.00 .00 33.00 162.00 12.00 59.00 58.00 43.00 .00 30.00 73.00 41.00 1.00
BM 63000 63000 63000 63000 63000 63000 64000 64000 64000 127000 127000 127000 65000 65000 65000 65000 65000 65000 65000 65000 65000 130000 130000 130000 67800 67800 67800 67800 67800 67800 67800 67800 67800 135600 135600 135600 71200 71200 71200 71200 71200 71200 71200 71200 71200 142400 142400 142400 74760 74760 74760 74760 74760 74760
Tabla 11.1.1 - Datos básicos empleados
Donde
PZ
es
Nº28-33/8/0028,
en
bombeos
la
en
m.s.n.m.
de
todos
los
en
litros/m2
pluviométrica
piezómetro
Alicante,
recogidas
Jijona,
de
pozos
del
que
y
extraen
PV
en el
BM
los
agua
del
son
los
expresados en m3/mes.
aculfero,
Los
piezométrica
expresadas en
estación
totales
cota
expresados
las pluviometrías mes
la
estadísticos
de
las
series
anteriores
siguientes:
Número de observaciones válidas Variable PZ PV BM
Media
Desv.
594.40 41.05 82584.44
tip.
54 Kurtosis
48.24 42.97 27965.59
Asimetría -.47 1.68** 1.36**
-.80 3.14** -.01
de la hipótesis de normalidad)
Superan el límite
Tabla 11.1.2 - Estadísticos de los datos brutos empleados
Los parámetros que sobrepasan
los límites
por la hipótesis de normalidad expuestos
requeridos
anteriormente
se
han marcado con dos asteriscos.
Se puede observar en la
serie
de
la
asimetría.
bombeos Por
pluviometría
variable
gran apuntamiento variable
la tabla anterior (PV)
y una cierta asimetría, (BM)
ello,
presenta
mientras
presenta
asimismo
series
incumplen
ambas
(11.1.2),
una
que
que un la
ligera
la hipótesis
65
66
de
normalidad
necesaria
para
su
inclusión
directa
en
los
cálculos del modelo lineal multivariante.
11.2.2.-
Evaluación
análisis
y
de
la
información
disponible
PIEZONETRIAS
Las observan
piezometrías en
el
piezómetro
esquemáticamente 4.1.
se
como
PZ,
expresan
m.s.n.m.,
en
Nº28-33/8/0028, representandose
y
se en
se
denominan la
figura
67
ACUIFERO DE JIJONA 670 660 650 64-0 630 620
Piez!lmetro Ni¡ 28-33/8/002a -
610 600 E
590
¿
560 570 -
N
IL
560 550 54-0 530 520 sio 500 490
~ J ..........
ol-ae
01-87
Figura 4.1
Como se a utilizar será
las
como DPZ,
la
construcción
diferencias
de
es
la variable
del
presente
piezometría,
que
dependiente
modelo
lineal
representaremos
siendo:
DPZ(t)
siendo
01-90
- Gráfica de piezometrías
verá a continuación,
en
ol-ag
ol-a8 t [Mescal
PZ(t)
la
piezometría
= PZ(t)
del
mes
la piezometría del mes anterior.
- PZ(t-1)
presente,
y
PZ(t-1)
68
En
la
figura
4.2
se
representa
la
variable
diferencia de piezometría:
ACUIFERO DE JIJONA Diferencias piez!lmetro NI¡ 28-33/8/0028
50 4-0 30 ~ 20 E
10 -
C>
0 -10
-20 -30 01-86
......... 01-89 01-90
...........01-88 01-87 t (Messel
Figura 4.2 - Diferencia de piezometría
De trata picos
la figura anterior
de de
Octubre,
una
máxima y
observación los
serie
el
menor
de
sufre
amplitud
anterior
períodos
que
la
(4.2)
en
valor no
se
se puede
oscilacione los en
tal
periódicas
se con
meses
de
Septiembre-
mes
de
Agosto.
el
ajusta
muestra,
observar que
estrictamente y
como
puede
a
La
todos verse,
69
si
aunque
sirve para explicar,
En la tabla 11.2.1 piezometría
piezometría de
la
(DPZ) ,
DPZ
en
sean
mayores
95%
de
como
la
que
(la
de
1986
en
la
de
variable
se
valor
para
54
de
indicó absoluto
-
1
casos
considerados
dicha
la
SDPZ,
diferencia
en
son
de
neperiano lado
como
que
en
Al
de
student
Obsérvese
al
para
ajusta la
lado
de
más
tabla,
las
en
el
ha
y
como
que
lineal,
de
la
un
valor
de
LNDPZ, sido
modelo,
condiciones
estudentización
existe
variable
adelante,
utilizar a
la
modelización
neperiano que
"logaritmo
y
confianza),
disponibles
en
se produce un valor extraño positivo.
SLNDPZ,
(logaritmo
tal
de
justificar
exigidas
Obsérvese
t
va
se
forma
diferencia
(LNDPZ).
valores
lado,
a
del
variable
otro
mejor
del 86,
2
variable
variable
la donde,
extraños.
Por se
de
datos
de
representado
aquellos
margen
Octubre de
como
ha
los datos
valores
piezometría"
de
anteriormente,
los los
de
estudentización
valores
y
se
piezometría,
se muestran
(PZ),
diferencia
variable
la
la serie.
y variación de
de
esquemáticamente,
la
diferencia extraño
de se
por
elegida ser
la
normalidad observa
variable de
que,
la
LNDPZ,
piezometría).
positivo
y otro también positivo en Agosto del 88.
en
Nov.
70
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 (**
FECHA DPZ 01-86 02-86 -9.40 03-86 -5.90 04-86 7.40 05-86 -5.10 06-86 -15.30 07-86 3.60 08-86 -15.50 09-86 9.60 10-86 29.20 11-86 38.50 12-86 -9.90 01-87 20.90 02-87 3.05 03-87 -5.95 04-87 -4.80 05-87 10.00 06-87 -12.50 07-87 5.20 08-87 -16.50 09-87 5.50 10-87 -1.70 11-87 24.80 12-87 27.90 01-88 21.40 02-88 6.35 03-88 -4.05 04-88 5.70 05-88 3.60 06-88 4.40 07-88 -12.80 08-88 -18.50 09-88 -6.20 10-88 -9.60 11-88 -12.80 12-88 -15.60 01-89 11.60 02-89 14.10 03-89 17.90 04-89 -4.70 05-89 -5.20 06-89 -2.18 07-89 -7.42 08-89 6.50 09-89 26.10 10-89 -3.90 11-89 8.80 12-89 9.50 01-90 10.40 02-90 -.27 03~90 -2.53 04-90 3.90 05-90 -2.60 06-90 -11.33 Valores extraños)
SDPZ
LNDPZ
-.88 -.61 .40 -.55 -1.32 .11 -1.34 .57 2.05** 2.76** -.91 1.42 .07 -.61 -.53 .60 -1.11 .23 -1.42 .26 -.29 1.72 1.96 1.46 .32 -.47 .27 .11 .17 -1.13 -1.57 -.63 -.89 -1.13 -1.35 .72 .91 1.20 -.52 -.56 -.33 -.73 .33 1.82 -.46 .51 .56 .63 -.18 -.35 .13 -.36 -1.02
3.42 3.53 3.86 3.55 3.21 3.78 3.20 3.90 4.24 4.36 3.40 4.11 3.76 3.53 3.56 3.91 3.31 3.81 3.16 3.82 3.65 4.17 4.22 4.12 3.84 3.58 3.82 3.78 3.79 3.30 3.07 3.52 3.41 3.30 3.19 3.94 3.99 4.06 3.56 3.55 3.63 3.48 3.84 4.19 3.59 3.89 3.90 3.92 3.68 3.62 3.78 3.62 3.36
SIBDPZ .87 -.52 .54 -.44 -1.56 .27 -1.58 .69 1.76 2.17** -.92 1.35 .23 -.52 -.42 .72 -1.21 .39 -1.72 .41 -.14 1.55 1.70 1.38 .47 -.35 .43 .27 .33 -1.25 -2.01** -.55 -.89 -1.25 -1.60 .82 .97 1.19 -.41 -.45 -.18 -.67 .48 1.62 -.33 .64 .68 .74 -.03 -.21 .30 -.22 -1.08
Tabla 11.2.1 - Transformaciones de piezometrías
71
En
la
siguiente
estadísticos calculados
tabla
univariantes partir
a
de
de
54
se
muestran
variables
las
los
DPZ
y
LNDPZ,
disponibles
datos
los
desde
1990:
1986 a Junio de
Enero de
11.2.2.
Número de observaciones válidas = 54.00 Variable
esta
observar
puede
se
tabla
se
logarítmica
transformación
(Desv.
variable la
ha
se
ajuste
se
asimetría
0.66.
y
Tal
(11.2.2), ser de
no
la
0.70,
con
indicado a
ser,
debe
como
DPZ
por ser
aunque
no
el
en
cumple
este
la diferencia
la
valor
de
requisito pequeña
tan
a
anterior
tabla
la
que
superior
absoluto,
valor
observa
se
variable
en
que
Tip.).
normal,
distribución
la
más
para
anterioridad,
la
con
que
disminuido
ha
considerablemente la desviación típica
Como
.08 -.57
- Estadísticos de las series de piezometría
Tabla 11.2.2
En
.70 .04
13.17 .31
2.14 3.69
DPZ LNDPZ
Asimetría Kurtosis
Tip.
Desv.
Media
por no
se rechazará.
Respecto observa 0.08,
que que
distribución
a
la
presenta, como
se
normal,
Kurtosis, para
la
comentó no
en
la
variable
11.2.2.
tabla DPZ,
el
anteriormente,
sobrepasa
en
valor
se
valor
de
para
la
absoluto
el
72
valor
de
1.33,
lo
por
que
se
acepta
hipótesis
la
de
normalidad.
Aunque la
la
hipótesis
valores
de
de
del
la
positivos
o
serie
presentar podrían
decide el
de
ésta
transformar
efecto
diversas
posibles,
y
elimina
de
la
estos
máximo
la
series de pluviometrías y bombeos.
algunos los
serie
a
valores.
transformaciones,
el
maximiza
cumple
distorsionar
se elige una transformación
normal,
como
extraños
se
preparación
al
que
minimizar
tras varias pruebas, la
incrementos de piezometría
modelo,
eliminar
Mediante
de
normalidad,
extraños
resultados fin
serie
que
de
y
deja
valores
correlación
con
las
Esta es:
LNDPZ - ln(DPZ + 40).
donde
1n
es
el
logaritmo
neperiano,
y
DPZ
la
variable
diferencia de piezometría.
En nueva
la
tabla
serie,
estadísticos. de
asimetría
inferiores la
en
asimetría
acepta
El
y
11.2.1,
se
en
la
tabla
Como
se
puede
y
Kurtosis
valor y
1.33
muestran 11.2.2. ver,
de
absoluto para
a
la
los se
los
esta
nueva
los
límites
Kurtosis,
transformada
de
indica
los que
valores aunque
se
de
esta
representan
valores
la nueva variable como normalmente
análisis
datos
sus
calculados
variable de
por
0.66 lo
son para
que
se
distribuida.
extraños observan
en dos
la
serie
valores
73
que
son
ligeramente
extraños,
no
es
necesario
depurarlos
para poder trabajar con los datos.
PLUVIOMETRIAS.
Las recogida
pluviometrías por
desde Enero
M2 de
en
1986
la
se
expresan
estación
hasta Julio
en
litros
pluviométrica de 1992,
y se
de
de
agua
Jijona,
representan
en la figura S.
ACUIFERO DE JIJONA PluviometrOca, Est. de Jijona 260 24-0 220 200 180 160 140 120 100 80 #o 00
50 -
20 0 01-86
01-87
01~88
01-89
t [Meses]
Figura 5
- Gráfica de pluviometrías
01-90
74
En una
esta
serie
los
y
el
se
para
88,
el
igualados,
En
la tabla
y
razones
junto
(SPV,
de a
de
mismo
que
régimen
del
los
del
sobre
en
todo
y
año.
medidas
Este
niveles
en de
de
En
los
la
años
avenidas,
muy
régimen
mensuales
es
diente
mínimo.
típicamente
grandes,
resto
Su
piezometrías,
continuación
la
(LNPV),
se
entrada
ellas a
las fin
presentan
serie que
estadísticas,
SLNPV),
extraños.
apariencia
11.3.1
(PV),
pluviometría"
de datos
Agosto.
los mínimos
bajas se
o
rompe
están
más
y los picos no son tan fuertes.
pluviometría
por
los
-
Julio
es
en
-
un
donde
situados y
así
el
máximos
Noviembre,
extremadamente
nulas
pluviometría
-
a
sigue
la
Octubre
máximo
se
que
con
diferencia
observa
1987
picos
casi en
la de
con
gráfica
con
-
meses de Junio
sierra
1986
anual,
Septiembre
similar a
observa
se
periódica
meses de en
gráfica
de
como ha
"logaritmo se
sido
verá
con
apreciar
los la
en
brutos
de
neperiano
de
más
elegida
pluviometría
series de
del
datos
los
el
datos
adelante,
como modelo
la
y
serie
lineal,
estudentizados
existencia
de
valores
75
FECHA 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54
PV
01-86 9.50 02-86 6.70 03-86 8.20 04-86 43.80 05-86 7.00 06-86 7.90 07-86 40.50 08-86 6.50 09-86 153.60 10-86 135.50 11-86 18.10 12-86 1.30 01-87 42.40 02-87 54.10 03-87 .30 1.20 04-87 05-87 45.30 06-87 .00 07-87 20.60 08-87 .00 09-87 36.20 10-87 90.40 191.40 11-87 12-87 62.10 40.40 01-88 02-88 30.30 03-88 10.00 60.70 04-88 05-88 78.50 06-88 86.70 07-88 3.40 08-88 5.00 54.00 09-88 10-88 61.00 11-88 45.00 12-88 .00 01-89 24.00 02-89 20.00 03-89 81.00 34.00 04-89 05-89 52.00 36.00 06-89 07-89 .00 08-89 33.00 09-89 162.00 10-89 12.00 59.00 11-89 58.00 12-89 01-90 43.00 02-90 .00 03-90 30.00 73.00 04-90 05-90 41.00 06-90 1.00 Valores extraños)
SpV -.73 -.80 -.76 .06 -.79 -.77 -.01 -.80 2.62** 2.20** -.53 -.93 .03 .30 -.95 -.93 .10 -.96 -.48 -.96 -.11 1.15 3.50** .49 -.02 -.25 -.72 .46 .87 1.06 -.88 -.84 .30 .46 .09 -.96 -.40 -.49 .93 -.16 .25 -.12 -.96 -.19 2.81** -.68 .42 .39 .05 -.96 -.26 .74 .00 -.93
LNPV 2.97 2.82 2.90 3.99 2.83 2.88 3.92 2.80 5.10 4.98 3.34 2.42 3.96 4.16 2.33 2.42 4.01 2.30 3.42 2.30 3.83 4.61 5.31 4.28 3.92 3.70 3.00 4.26 4.48 4.57 2.60 2.71 4.16 4.26 4.01 2.30 3.53 3.40 4.51 3.78 4.13 3.83 2.30 3.76 5.15 3.09 4.23 4.22 3.97 2.30 3.69 4.42 3.93 2.40
SLNPV -.74 -.92 -.82 .45 -.90 -.84 .38 -.94 1.76 1.62 -.31 -1.38 .42 .66 -1.49 -1.39 .49 -1.53 -.21 -1.53 .27 1.19 2.01** .80 .38 .11 -.71 .77 1.04 1.14 -1.18 -1.05 .66 .78 .48 -1.53 -.09 -.23 1.07 .22 .62 .27 -1.53 .19 1.82 -.60 .75 .73 .44 -1.53 .10 .96 .39 -1.41
76
Tabla 11.3.1 - Datos brutos y transformados de pluviometrías
Como extraños el
se
indicó
aquellos
valor
de
intérvalo
2
de
de
(PV)
Septiembre del de
ligeramente
en
la
89.
Se
Octubre Se
extraño
observa
transformación
en
4
estudentizada
superan
54
tan
ha
datos
que
valores 86,
que
Noviembre
logarítmica
valores
observa
del
pluviometrías
consideran
para
presenta
86,
se
serie
estudent
confianza).
del
neperiano
que
(t
pluviometrías Septiembre
anteriormente,
la
la
y
95%
de
serie
de
extraños,
en
Noviembre
del
serie
logaritmo
sólo del
del
presenta 87,
mejorado
es
un
87
y
valor
decir,
la
considerablemente
la disposición de los valores extraños de pluviometría.
En
la
estadísticos
tabla de
transformación
11.3.2
las
se
series
presentan de
los
parámetros
pluviometrías,
y
Kurtosis
PV LNPV
41.05 3.60
Tabla 11.3.2
- Estadísticos básicos de pluviometrías
indica
y la
su
logarítmica:
Número de observaciones válidas - 54.00 Variable Media Desv. Tip. Asimetría
Tal
de
como
42.97 .85
se
tendencia
indicó de
1.68 -.06
3.14 -.95
anteriormente, la
mayoría
de
la los
asimetría datos
a
77
separarse
de
la
separación
es
hacia
los
cuando
es
hacia
los
lo
normalmente en
la
para
media,
54
menor
datos
de
y
un
0.66.
anterior,
valores más
distribuida,
modelización,
la
siendo más
positiva
cuando
pequeños,
y negativa
grandes.
y por
Se
tanto
útil
aquella
cuyo
95%
confianza,
de
Como
serie
se
valor
puede
de
considera para
como
empleo
de
asimetría
en
valor
apreciar
datos
su
brutos
en
de
la
sea,
absoluto la
tabla
pluviometría
incumplen este requisito de normalidad.
La
kurtosis,
tendencia los de
de
los
la moda,
mayor de
la
y negativa
que
54
está
normalmente
casos
valor
tabla
y
siendo se
pluviometría
a
esperaría
95%
de
valores
similares
en
cuando una
intérvalo
1. 33.
un
si
la
esta
Se
la
tendencia
distribución
de
su
Como
(11.3.2),
presenta
indica
lo contrario.
distribuida a
muestral,
presentar
positiva
ocurre
absoluto
anterior
apuntamiento
datos
cuando
para
en
o
a es
normal,
considera
que,
confianza,
la
serie
apuntamiento
es
menor
se
puede
serie
de
apuntamiento
observa
datos de
3. 14,
en
la
brutos
de
superior
al de una serie normal.
Para y
eliminar
mediante un
normalizar
la
método
siguiente:
los
esta
valores
preparación de
serie
tanteo,
de
extraños de la
series
datos se de
de
ha
pluviometría seleccionado,
prueba
transformación
y
siguiendo
logarítmica
78
LNPV = ln(PV + 10)
por presentar parámetros normales o apuntamiento, desviación extranos,
presentar
maximizar
la
diferencia de piezometría
Por empleada
todo en
asimetría y kurtosis
disminuir la relación entre
típica, y
de
lo
el
menor
correlación
la media y
número con
de
la
valores
la
serie
de
(LNPV)
será
la
(DPZ).
expuesto
la modelización
esta
serie
lineal
como
datos
de
entrada
de pluviometrías.
BOMBEOS
Los mes,
bombeos
se
expresan
en
m3
y se representan en la figura 6.
de
agua
extraídos
al
79
ACUIFERO DE JIJONA Bombsos totales en el acuOfero 150 1 40 130 120 110 100 90 ao 70 (630 50 40 30 20 1 C) 0 01-86
01-90
01-89
01-88
01-87
t [meses]
gráfica de bombeos
Figura 6
gráfica
observa
se
En
la
función
de
Dirac
determinista,
-
0/11)
se
("golpe"
Noviembre-Diciembre, agua
en
produce debido
meses
estos
en
en en
la
que meses
los
al empleo
la
a
corresponde
que
el de
una
"impulso" Octubre-
de mayor volumen
industria
turronera
de
de la
población.
El
cero
74.760 m3� de donde
de
Dirac
está
en
un
y el uno en un bombeo de
bombeo 127.000
se deduce que existe una deriva de
de
63.000
a
a 142.400 m3�
80
en el
0
74.760 m3 -63.000 m3)14,5 años
de Dirac=(
2.613,3 m3/año. *
en
el
1
de
Dirac-
m3
(142.400
-
M3)/4,5
127.000
años=3.850,0 m3/año.
Se
trata,
antrópica o Es
de
las
tanto
de
tecnológica e
constante
como
por
de
la
diferente para el 0 y para el
El porcentaje
de
desarrollo
del
necesidades
variable
independiente
depende
y
una
determinista, la naturaleza.
industrial,
población,
con
una
así
deriva
1 de Dirac.
anual de agua empleada en
la industria
anualmente,
entre
el
con respecto
al consumo total varía,
18 y el 20%.
La diferencia en la deriva entre el 0 y el
de
Dirac
indica
que
el
crecimiento
del
el abastecimento urbano es mayor que industrial,
decir,
es
las
agua
en el
necesidades
empleada
1 en
abastecimiento
de
la
población
crecen más deprisa que las de la industria.
En la tabla 11.4.1. bombeos de
los
aplicada la
(BM),
su
valores (LNBM),
seleccionada
modelo
lineal,
transformación:
estudentización extraños, que
como
para
su
así
los datos brutos
se presentan
como
(SBM)
para
la
transformacíon
se
indicará
empleo y
la
en
la
más
el
de
análisis
logarítmica
adelante
construcción
estudentización
de
será del la
81
NUM FECHA 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 (**
BM
01-86 63000 63000 02-86 03-86 63000 04-86 63000 63000 05-86 06-86 63000 07-86 64000 64000 08-86 64000 09-86 10-86 127000 11-86 127000 12-86 127000 01-87 65000 02-87 65000 03-87 65000 04-87 65000 05-87 65000 06-87 65000 65000 07-87 65000 08-87 65000 09-87 10-87 130000 11-87 130000 12-87 130000 01-88 67800 02-88 67800 03-88 67800 04-88 67800 67800 05-88 06-88 67800 67800 07-88 08-88 67800 67800 09-88 10-88 135600 11-88 135600 12-88 135600 01-89 71200 02-89 71200 71200 03-89 04-89 71200 05-89 71200 06-89 71200 07-89 71200 08-89 71200 71200 09-89 10-89 142400 11-89 142400 12-89 142400 74760 01-90 02-90 74760 74760 03-90 04-90 74760 05-90 74760 74760 06-90 Valores extraños)
SBM -.70 -.70 -.70 -.70 -.70 -.70 -.66 -.66 -.66 1.59 1.59 1.59 -.63 -.63 -.63 -.63 -.63 -.63 -.63 -.63 -.63 1.70 1.70 1.70 -.53 -.53 -.53 -.53 -.53 -.53 -.53 -.53 -.53 1.90 1.90 1.90 -.41 -.41 -.41 -.41 -.41 -.41 -.41 -.41 -.41 2.14** 2.14** 2.14** -.28 -.28 -.28 -.28 -.28 -.28
LNBM
SLNBM
-4.61 -4.61 -4.61 -4.61 -4.61 -4.61 -4.07 -4.07 -4.07 -.78 -.78 -.78 -3.73 -3.73 -3.73 -3.73 -3.73 -3.73 -3.73 -3.73 -3.73 -.73 -.73 -.73 -3.13 -3.13 -3.13 -3.13 -3.13 -3.13 -3.13 -3.13 -3.13 -.65 -.65 -.65 -2.69 -2.69 -2.69 -2.69 -2.69 -2.69 -2.69 -2.69 -2.69 -.57 -.57 -.57 -2.38 -2.38 -2.38 -2.38 -2.38 -2.38
-1.44 -1.44 -1.44 -1.44 -1.44 -1.44 -1.03 -1.03 -1.03 1.53 1.53 1.53 -.76 -.76 -.76 -.76 -.76 -.76 -.76 -.76 -.76 1.56 1.56 1.56 -.29 -.29 -.29 -.29 -.29 -.29 -.29 -.29 -.29 1.62 1.62 1.62 .04 .04 .04 .04 .04 .04 .04 .04 .04 1.69 1.69 1.69 .29 .29 .29 .29 .29 .29
Tabla 11.4.1 - Datos de bombeos
82
Se observa extraños,
la existencia
de
Octubre
correspondientes Asimismo
se
aplicada,
al
ha
que
eliminado
tres valores
a
período
observa
de
ligeramente
Diciembre
de
la
bombeo
de
para
la
transformación
por
completo
1989,
industria. logarítmica
esos
valores
(la
serie SLNBM no presenta valores extraños).
En
la
tabla
estadísticos de
la
11.4.2. serie
y de su transformación
se
de
presentan
datos
logarítmica
los
brutos
de
parámetros
bombeos
(BM)
(LNBM):
Número de observaciones válidas - 54.00 Variable BM LNBM
Media Desv. 82584.44 -2.75
Tip.
Asimetría Kurtosis
27965.59 1.29
1.36 .47
Tabla 11.4.2 - Estadísticos de
Como
se
puede
bombeos
(BM),
superior
(1.36).
presentan
As¡. mismo, serie
de
kurtosis
se
datos casi
observar,
observa
brutos
igual
a
de 0,
las series de bombeos
la
un
en
serie valor
la
lo
de
datos
de
tabla
bombeos por
-.01 -.81
asimetría
anterior
(BM), que,
brutos
con
que
presenta respecto
de muy
la una al
83
apuntamiento,
la
serie
se
considera
como
normalmente
distribuida.
A
fin
de
los valores
extraños,
mediante una de
series
mejorar
función
de
el
valor
de
asimetría
se ha decidido logarítmica.
pruebas
y
y
transformar
Mediante
aplicando
eliminar la
serie
la preparación
procedimiento
un
de
se ha elegido la siguiente transformación:
tanteo,
LNBM = ln(0.000007-BM - 0.43)
por ser
maximiza la
correlación
de piezometría
Como serie
se
debajo
que
la ha
puede
de
serie visto,
de
la
requisitos con
serie
de normalidad,
la
de
diferencia
(DPZ).
calculada
por
se
manteniendo los
la que,
presenta los
pueda al
serie
calculada
ello,
la
en
ver
serie LNBM
de
comentados
la
valores
se
como
que
extraños.
variable
la modelización lineal del acuífero.
y
normal.
11.4.l.,
tabla
la
de
nueva
kurtosis
anteriormente
como
presenta elige
11.4.2.,
asimetria
considerada
comentar no
tabla
valores
límites ser
la
para
También la Por
nueva todo
entrada
en
84
11.2.3.-
Conclusiones
a la
fase
de evaluación
y análisis
de la información disiponible.
A fase
de
Junio las de
la
vista
lo
evaluación,
de
1990
en
variables las
de
a
en
el
el
utilizar
y
en
como
bombeos
anteriormente
período
aculfero
piezometrías
pluviometrías
expresado
de el
Enero
de
Jijona, ajuste
variable
se
del
de
esta
1986
concluye modelo
a
que
lineal
dependiente,
variables
como
en
y
las
independientes
son:
Variable dependiente:
LNDPZ -
ln(DPZ + 40)
Variable independiente PV:
LNPV - ln(PV + 10)
independiente BM:
Variable
LNBM = ln(0.000007-BM - 0.43)
Donde PV
y
BM,
1n
es
las
la
función
variables
diferencia de piezometrías,
Se ha apreciado en
la
forma
diferencias
de
de
logaritmo explicadas
y
anteriormente
DPZ, de
pluviometrías y bombeos.
la existencia de las
neperiano,
semejanza
correspondientes
curvas
piezometrías
una clara
y
pluviometrías,
a
las
con
una
85
forma cíclica periódica anual Octubre - Noviembre, siendo
en Junio -
y mínimos
presenta un
impulso en
Diciembre,
que
corresponde
a
- Noviembre
utilización
su
que
Dirac
de
de Octubre
los meses
-
Julio - Agosto,
función
una
bombeos
variable
la
en Septiembre
con máximos
en
la
industria turronera típica de la zona.
11.3.- IDENTIFICACION DEL MODELO LINEAL
En
esta
variables,
sobre
entre
correlación,
identifican
se
las
entre
relaciones
establecer las que pudieran tener más
a fin de
influencia desfase
fase
la las
mismas
que
cual
y
proporciona
es
el
mayor
la grado
el
establecer
como
así
dependiente
variable
de
independencia entre las variables independientes.
Todo ello se realiza de la siguiente manera:
lº)
Análisis
variable
de
las
dependiente
correlaciones piezometría
variables independientes,
2º)
y
entre
cruzadas cada
una
de
la las
pluviometría y bombeos.
Análisis de la matriz de correlacion bivariante.
Para
que
significativa
la se
correlación requiere,
entre
para
54
dos datos
variables
sea
el
de
y
95%
86
intérvalo de
confianza,
que
el
coeficiente
de
correlación
sea superio a:
Rho = V(t2/(n-2+t2))
11.3.1.-
Análisis
de
= V(4/(54-2+4)
- 0,26
correlaciones
cruzadas
entre
diferencia de Diezoigetrías y r>luvionetría.
La de
las
gráfica series
de
de
la
las
función
correlaciones
de
diferencias
de
piezometría
cruzadas con
las
pluviometrías se observa en la figura 7.1.
Correlaciones cruzadas : LIDPI LUY Corr. Error Dif. Cruz. Std. -1 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7
-.129 -AG7 -391 .181 .032 .012 .120 .672 -366 -.109 -341 .030 -.168 -.028 -.106
-.75
-A -.25
0
.25
.5
.75
1
.147 .146 .144 .143 .141 .140 .139 .137 .139 .140 � 141 .143 .144 .146 .147
Siabolos : Autocorrelaciones
Litites de error standard
Figura 7.1 - Funcion de correlaciones cruzadas, dif. pluviometrías
piezometrías
87
Tal aprecia
y
como
una
0
con
observa
correlación
desplazamiento pluviometría
se
en
el
influye
retrasos
piezométrico,
en
esta
positiva
tiempo,
importante
lo
que
instantáneamente
inferiores
figura
a
(7.1), (0,67)
indicaría (en
éste,
se sin
que
la
el
mismo
mes),
en
el
nivel
es decir,
"LNDPZ(t)
es
linealmente
dependiente
de una forma
significativa con LNPV(t)"
o
dicho
de
manifiesta de
un
En
modo,
en
incremento
piezometría),
descenso
de
otro
y
en
el
del el
mes
nivel mes
en
en
que
llueve,
(diferencia que
no
se
positiva
llueve,
un
(diferencia negativa).
figura
7.2.
las dos variables:
se
presenta
la
gráfica
de
dispersión
88
4.05L N D p Z
1 11
2
2 11
11 3.6-
1
11
2
1
21
11 1
1 11
2 1
1
1 11
3.15-
2 1 2.975 2.55
3.825 3.4
4.675 4.25
5.1
LNPV Figura 7.2 - Gráfica de LNDPZ con LNPV
En agrupa
esta en
inclinado la
figura un
(alrededor
cada
y
vez
grande
al más
de
y
de
los
Noviembre de
que
45º),
se
grandes.
Sin
constituyendo
(outlayer),
con
decrecer
baja,
se
valor
corresponde
1986,
de
la
y que
nube
un
eje
tal
forma de
producen
produce
puntos
el
que
la
al
crecer
una
la
observación
son
también
incremento
más
pluviometría
fuera de la nube en
se
dominantemente
piezometría
con
extraño a
de
decrementos
embargo,
saliéndose
un
que
incrementos
piezometría
relativamente
aprecia
elipsoide
pluviometría
mayores,
se
de puntos
correlación del
sería necesario depurar.
mes
de
89
Este
valor extraño
interpretar que
un
corresponde
propia
valor
a
un
naturaleza), anteriores
meses subidas 4. 2. ) , en
como
de
la
que
estos
en Noviembre defectuoso
valor
ya
que
en
piezometría
corresponden lugares
al a
año,
los
fuertes 4.1.
que
producidas
riadas
la
en
figuras
natural
sino
por
observan
las
puede
medida,
mismo
se
f enómeno
debidos
la
se
(producido
este
(véanse
no
1986
en
anómalo
posteriores,
y
de
y
acontece por
gotas
frías.
Se
observa
asimismo
grupo de puntos que corresponden,
muy
bajas
piezometría
no
observa
existen
87,
que 7.2.),
Abril
del
la
figura
estas con
extraños 5
lluvia
seguramente
negativas
de
esperarse.
Se
8
puntos
hay
Marzo
4
del
(ver
(Marzo 90)
el
87,
seco.
de
año
sin
año
bastante
pluviometrías), a
meses
de
el
altas. uno
de
formas, modelo
se
sin
importancia,
todas
embarrarán
un
a
Es
puede
decir,
sequía, aunque se
que
se
estos En
ver que
pluviometría, los
lineal,
del
que de
relativamente
De
y
a
la mayor acumulación
precedido
precipitación.
89
los cuales
siendo
corresponden
pluviometrías
pluviometrías de
en
unos
claramente
del
pequeño
anterior,
debiera
grupo
un
la nube,
diferencias
este
Julio
de
caso
como
cuales
(gráfica
medidas
con
en
los
87,
al
bajas
la nube,
separan de valores
tan
de
existencia
se sitúa fuera de
igualmente
pluviometrías
figura
la
o
siguen
un
periodo
sin
apenas
estos
datos
puede
pensar
90
en
su
utilización,
explotación, garantizando
para
valores que
ante
forzarle
desde un
el
lado
período
produzca la sobreexplotación,
a
de
predecir, de
la
sequía
no deseada,
11.3.2.- Análisis de las correlaciones
en
su
seguridad,
similar
no
se
del acuífero.
cruzadas entre
las
diferencias de i>iezonetría y los bombeos.
La de
las
gráfica
de
variables
la de
función
de
diferencia
correlaciones de
piezometrías
variable de bombeos se observa en la figura 8.1.
cruzadas con
la
91
Correlaciones cruzadas
Dif. -9 ~8 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6
:LNDPZ LNBM
Corr. Error Cruz. Std. -1 .151 -.263 -.421 .149 .147 -.218 -.174 .146 .144 -.066 .143 .010 .141 .199 .140 .217 .139 .323 .137 .208 .139 .294 .140 -.054 .141 -.124 .143 -.325 .144 -.146 .146 -.110
Simbolos
-.5 -.25
Autocorrelaciones
Figura 8.1
En
75
Correlaciones
esta
figura
negativa
correlación
se
0
y
correlación +l
meses.
observa
casual,
posiblemente
con
series
y
la
a
coincidencia
un
de
Asimismo,
esta
que
a
la
en
la
una
desplazamiento se observa
con desplazamientos
señalar
debida
existencia
la
significativa
que
.75
cruzadas LNDPZ con LNBM
positiva menor Hay
.5
de
última
periodicidad evolución
-1 es
de
las
dentro
del
período.
Como
se
ha
principio
bombeos,
en
a
misma,
ella
antrópica
y
comentado
no
y es
tecnológica.
por una El
1
Limites de error standard
de 8 meses atrás y 4 meses adelante. una
.25
anteriormente, propia
la
naturaleza
variable hombre,
variable intrínseca
aleatoria,
para
su
sino
desarrollo
92
económico este
y
de
su
propia
aculfero
para
dos
y prioritario
para
el
cual
un
implica
Enero
de
1986
segundo
Diciembre
hasta
mes
de
de
1989.
la
cada
tendencia
principales:
de
m31mes
en
74.760 la
año,
en
Se
observa
de
crecimiento
los
varía
que
en
o
población,
lo
M3/mes
en
1990;
en
de
de
muy
Octubre m3
127.000 el
en
lugar
turronera,
mes
consumos
lineal,
agua
primer
Junio
desde
ambos
el
63.000
meses
142. 200 m3
a
la
industria
que
86,
en
de
estimado
zona,
Octubre del
utiliza
abastecimiento
para
en de
usos
consumo
lugar,
importante
actividad,
en
a el
de Octubre
presentan
incluso
una
ligeramente
exponencial.
Esto de
8
implica
meses
de
piezometría, ciclo
que
la
correlación
retraso
sea
y
4
simplemente
natural
de
piezometría y el
una
aguas
fenómeno
de
incremento
produce ciclo los
de
simultáneamente
natural
datos
causal
brusco
del
del
directa
los
uso
no
inmediata
sobre entre
bajas
industrial
de
del
agua
Es de observar que,
bombeos
En
observada
coincidencia
en
la época de aguas
aculfero.
presentados
adelanto
altaslaguas
los meses de Octubre a Diciembre. el
negativa
estos altas
consecuencia
se
puede
entre
los
meses,
observar bombeos
y
la el la en con se
debida al
de
esto,
en
correlación los
niveles
plezométricos.
Por
las
estocástica
necesidades de
incluir
impuestas a
esta
en
variable
la
modelización
bombeos
en
la
93
construcción grado
de
del
modelo
influencia
acuífero,
y
a
en
fín
el
de
poder
funcionamiento
teniendo
en
anteriormente
sobre
simultaneidad
cíclica anual
de las
la
modelización con
retraso
la
cuenta
tres series,
la
variable
respecto
a
determinar
del
todo
lo
expuesto
en
la
evolución
se decide
bombeos
la
general
el
introducir
con
variable
8
en
meses
de
diferencias
de
piezometría,
En
la
relación
de
figura estas
8.2. dos
se
presenta
variables,
la
gráfica
desplazando
8
variable de bombeos.
1
2
4.05L N D p Z
1 2 1 1
3.6-
1
1 1
1 3 1 1
1 2 1 1 3
1 2
1 2 1 1 21
3.15
1
-4.2 -4.8
1
-3 -3.6
-1.8 -2.4
LNBM8
Figura 8.2
- Gráfica de LNDPZ con LNBM8
-.6 -1.2
de
meses
la la
94
Se observa en esta figura, nube
de
puntos,
estos
alrededor
decir,
a
menos
con de
una una
una gran dispersión de
cierta linea
bombeos
inclinada
incrementos
bombeos,
decrementos más grandes.
11.3.3.
Análisis
de
tendencia
la
matríz
el
análisis
más
a
agruparse
descendente, altos,
y
a
correlaciones
de
la
es más
del
modelo.
Una
vez
cruzadas,
expuesto
entre
en el
con
de
modelo,
ajusten
modo
de
las
correlaciones
y aunque se realizó a la vez el análisis
correlaciones
se
de
para
resumen
correlaciones
el el
y/o
cada una fin de objetivo
de
las
variables
encontrar
aquellas
propuesto,
conclusión,
entre
las
el
se
análisis
diferentes
de las
utilizadas que
mejor
presenta, de
la matriz
variables
que
han considerado(Tabla 11.5).
MATRIZ DE CORRELACION ENTRE LAS VARIABLES DPZ LNDPZ PV LNPV BM8 LMBMB
1.00000 .98415 1.00000 .56227 .56990 1.00000 .89584 1.00000 .60260 .64195 -.40524 -.43315 -.32285 -.40098 1.00000 .92738 1.00000 -.45831 -.46943 -.32120 -.35848 BM8 LNBMB PV LNPV LNDPZ DPZ
Tabla 11.5
- Matriz de correlacíon
a
se
95
Debido
al
desplazamiento
de
la
número de casos ha descendido de para
46-1
casos
correlación
es
mínimo
de
54 a 46.
2,01,
y
requerido
de
0,29,
con
lo
que
para el caso en estudio
LNDPZ y LNPV
las
de
bombeos,
el
La t de studend
el
por
correlación significativa con 46 es
serie
coeficiente
la
de
hipótesis
casos y 95%
correlacioens
de
de
confianza
significativas
son
es de
0.64
LNDPZ Y LNBM8 es de -0.47
Que son las más juego, la
con
lo que
hipótesis
entre
y LNBMB,
de
puede
solamente
las variables puestas son
normalidad, entre
correlación es
no
de
correlacionadas
altas entre
las
sino
ellas. dos
variables
-0.36,
que
aunque
suficientemente
que
verifican las
más
lado,
la
independientes,
LNPV
Por
las
considerar
que más
en
otro
siendo baja
son
significativa
como
para
se
poderlas
introducir en la estimación del modelo.
12.3.4.- Conclusiones a la identificación del modelo
a)
Relación
de
pluviometrías.
la
diferencia
de
piezometrías
con
las
96
Se entre
observa
las dos
una
correlación
series,
positiva
sin desplazamiento
indica la influencia inmediata de nivel
piezométrico
correlación es lo que en en
del
aculfero
superior cuando
acuífero
medio
de
son
casi
períodos
cantidad
importante
se
Fenómeno
que
con
lluvias
en
de
que
siempre
agua
temporal,
lo
de
las
series,
lluvias
fuertes
más
torrencial
perdiéndose
escorrentía
menores
Esta
las
caracter
secos,
que
en el
Jijona.
transforman
bastante
de
grande
las pluviometrías
se explica por el hecho de
este
bastante
una
superficial.
repartidas
en
el
tiempo no se produce.
Se ha decidido, entre de
por todo esto,
la siguiente
forma
+ 40)
-> ln(pv(t)
+ 10)
Relación diferencias de piezometría con los bombeos.
La retraso de
la relación
las diferencias de piezometrías y las pluviometrías
ln(dpz(t)
b)
modelízar
correlaclon
negativa
es
un
evolución
altos
del
a
debida
cíclica
período
período
de
que
necesidades
las
aguas
de
altas de
encontrada
fenómeno anual, uso del agua
que
natural
acuífero.
de que
hace
industrial
son
con
8
meses
coincidencia los
coincidan Es
decir,
periódicas
de
al
bombeos con
el
debido igual
a
que
97
la evolución se
decide
de las
diferencias de piezometrías.
modelizar
datos de la siguiente
ln(dpz(t)
+ 40)
la
relación
entre
ambas
Por ello series
forma:
ln(0,000007-ba(t)
0,43).
de
98
11.4.- FASE DE CALIBRACION DEL MODELO LINEAL
Con
todo
seleccionado
lo
un
expuesto
modelo
de
anteriormente,
regresion
lineal
se
ha
multiple
para
representar al acuifero
ln(pz(t)
-
pz(t-1)
a2-1n(0.000007-bm(t-8)
+
40)
-
ao
+
al.ln(pv(t)+10)
+
- 0.43)
donde despejando se tiene
pz(t)
-
pz(t-1)
a2-1n(0.000007*ba(t-8)
exp(a0
- 0.43»
en la que:
ln:
exp:
es el logaritmo neperiano.
exponencial en base e.
pz:
piezometría.
pv:
pluviometría.
ba:
bombeos.
- 40
+
al-ln(pv(t)+10)
+
99
ao,
al y a2:
Para la suma
los coeficientes a calcular.
ello
de
se
las
ha
empleado
distancias
al
puntos
(ln(dpz(t)+40),
0.43))
al plano de regresion.
En resultado
la de
algoritmo de
tabla los
un
algoritmo
cuadrado de
ln(pv(t)+10),
siguiente coeficientes
minimiza
cada uno
de
los
ln(0,000007-bm(t-8)-
11.7.1,
se
calculados,
ajuste por mínimos cuadrados
cálculo del modelo
que
a
muestra
el
partir
del
empleado,
en el
lineal:
DESVIACION CORR COEFICIENTE ERROR STD. VAR. MEDIA TIPICA X / Y COEFICIENTE STANDARIZ. DEL COEF. F VALOR P LNPV 3.68217 .87267 .64195 .19238 .54352 .04178 21.20007 .00004 LNBM8 -2.88848 1.30762 -.46943 -.06486 -.27459 .02788 5.41092 .02479 DEPEND. LNDPZ 3.72239 .30888 TERMINO INDEPENDIENTE 2.82666 ERROR STD. DEL TERM. IND. .14969 CORRELACION MULTIPLE (R) .69124 COEFF. DE DETERMINACION .47781 ERROR STD. DE LA ESTIMACION .22834 DETERMINANTE DE LA MATRIZ .87150 NUMERO DE OBSERVACIONES 46
Tabla
11.7.1
- Resultado del algoritmo de
la regresion
VALOR
100
En
dicha
independientes series que
tabla
se
LNPV,
originales
observa,
LNBM8,
de
series
las
variables
transformadas
pluviometría
los coeficientes calculados
(PV)
y
de
bombeos
las
(BM),
son:
(Término independiente) aO - 2.83 t 0.15
2 - 2.83 t 0.30
(Pluviometrías)
al = 0.19 ± 0.04
2 - 0.19 ± 0.08
(Bombeos)
a2 -0.06 ± 0.03
2 -0.06 ± 0.06
El
intérvalo
correspondiente propio
valor
a
incertidumbre
de la
variable
coeficiente,
del
lo cual
de esta última variable
es más
intérvalo de incertidumbre
bombeos cosa
variable pluviometrías,
al
para
que
del
coeficiente
es no
parecido ocurre
con
al la
indica que el coeficiente significativo
que el
con respecto
coeficiente
relativo
a la variable bombeos.
Por valor de
de
otro 2.83,
y
Fisher
deduce para
probabilidad
las
es
el
independiente
término
significativo
respecto
tiene
al
un
intérvalo
al presentar éste un valor de solo 0.3.
incertidumbre,
Se
la
lado,
también
en
variables
dicha en
correspondiente
significación
de
dicha
el
tabla modelo,
dejado
al
variable
11. 7. 1 y
el
azar,
la
F
de
P
de
valor que
incluida
indica en
la
101
regresión. variable indica
A
la
vista
pluviometría que
es
de
esta
tiene
tabla
una F
altamente
de
se
observa
del
5%,
trabajando, confianza. superior
nivel
o
lo
Si
la
al
de
21.2,
que
significativa,
y
que
la
más
allá
del
5.41,
de
2.48%,
que
que
lo
variable
la regresión, de
se
el
95%
observa
estandar
final
a
una
la
desviación
la
rechazada
predecir.
valor
del
en
al
esta
estimaci6n
la
análisis de
de
la
Así mismo
coeficiente
regresión es de 0.69,
En
azar
requerido.
5%
de
fuera como
Respecto
dejada
requerido
al
es
siguiente
variable se observa
de
que
de
azar
por
el
error que
al
de confianza nos valor
similar
a
dependiente
que
se
en
que
el
de
la
realizado
el
correlación
11.7.2
el
0.23,
la
se
la varianza de la regresión:
tabla
múltiple
valor suficientemente
tabla
al
significativa.
tabla
incertidumbre de 0.46,
típica
quiere
grado
sería
probabilidad
inferior
además
intérvalo de
estamos
al
multiplicarlo por la t de student al 95% da un
de
al
al estar esta probabilidad
confianza
ser
en
que
dejada
grado de confianza demuestra ser altamente
Se
inferior
observa que presenta una F de Pisher
lugar al
muy
con
mismo,
probabilidad
grado
da
es
la
valor de bombeos, de
incertidumbre
que
5%
significativa en
de
la
Pisher
probabilidad dejada al azar es de 0.004%, margen
que
ha
alto.
102 ANALISIS DE LA VARIANZA DE LA REGRESION FUENTE DE VARIACION
GRADOS SUMA DE MEDIA DE LIBERTAD CUADRADOS CUADRAD. VALOR F
VALOR P ATRIBUIBLE A LA REGRESION .00000 DESVIACION DE LA REGRESION TOTAL
2
2.051
1.026
43 45
2.242 4.293
.052
19.67297
Tabla 11.7.2 - Analisis de la varianza de la regresión
En esta tabla se observa que Fisher al azar
de
19.67,
que
da
lugar
inferior a 0.001%,
a
lo que
se ha obtenido una F de una
probabilidad
dejada
indica que la regresión
es altamente significativa.
De esta forma,
(i)
el modelo resultante es el siguiente:
LNDPZ - 2.83 + 0.19
donde E es el error,
Aplicando transformaciones
-
LKPV - 0.64
que varía entre ±
a
la
realizadas
despejando el valor de
-
LNBM + E
2xO.23
anterior sobre
0.46.
ecuación las
las
variables,
la piezometría se tiene:
(¡¡)pz(t) - pz(t-1) = erp[(2.82666+0.19238zln(pv(t)+10)-0.06476xln(O.000007zba(t-8)-0.43)�E)J-40
y
103
la
pluviometría,
en
piezometría
la
pz(t)
neperiano,
base
en
exponencial
la
es
exp
donde
el
1n
el
instante
bombeos,
los
bm(t-8)
e,
y
E
logaritmo t,
pv(t)
el
error
en la estimación.
cometido
siguiente tabla
En la
la
de
(DPZ),
transformación
cálculo
diferencias resIduos predichos
calculados
o diferencias entre (RES
= LNDPZ*
los resIduos
(SRES),
extraños
esta
en
empleada
modelo
de
(LNDPV),
los
observados y los
y la estudentización
cuyo
variación normal del resto de los datos).
el las
poder observar
(residuos
en de
valores
los valores
- LNDPZ),
a fin
serie
el
por
los
piezometrías
de
logarítmica
(LNDPZ),
modelo
del
los valores
se muestran
diferencia
variable
observados la
11.8.1
valor
los
de
valores
excede
la
104 NUM FECHA 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 (**
DPZ
01-86 02-86 -9.40 03-86 -5.90 04-86 7.40 05-86 -5.10 06-86 ~15.30 07-86 3.60 08-86 -15.50 09-86 9.60 10-86 29.20 11-86 38.50 12-86 -9.90 01-87 20.90 02-87 3.05 03-87 -5.95 04-87 -4.80 05-87 10.00 06-87 -12.50 07-87 5.20 08-87 -16.50 09-87 5.50 10-87 -1.70 11-87 24.80 12-87 27.90 01-88 21.40 02-88 6.35 03-88 -4.05 04-88 5.70 05-88 3.60 06-88 4.40 07-88 -12.80 08-88 -18.50 09-88 -6.20 10-88 -9.60 11-88 -12.80 12-88 -15.60 01-89 11.60 02-89 14.10 17.90 03-89 04-89 -4.70 05-89 -5.20 06-89 -2.18 07-89 -7.42 08-89 6.50 09-89 26.10 10-89 -3.90 11-89 8.80 12-89 9.50 01-90 10.40 02-90 -.27 03-90 -2.53 04-90 3.90 05-90 -2.60 06-90 -11.33 Valores extraños)
Tabla 11.8
LNDPZ
LNDPZ*
3.42 3.53 3.86 3.55 3.21 3.78 3.20 3.90 4.24 4.36 3.40 4.11 3.76 3.53 3.56 3.91 3.31 3.81 3.16 3.82 3.65 4.17 4.22 4.12 3.84 3.58 3.82 3.78 3.79 3.30 3.07 3.52 3.41 3.30 3.19 3.94 3.99 4.06 3.56 3.55 3.63 3.48 3.84 4.19 3.59 3.89 3.90 3.92 3.68 3.62 3.78 3.62 3.36
4.11 4.08 3.77 3.59 3.89 3.93 3.54 3.56 3.86 3.32 3.54 3.32 3.81 3.96 4.09 3.89 3.82 3.78 3.64 3.89 3.93 3.75 3.37 3.40 3.83 3.85 3.80 3.47 3.71 3.68 3.90 3.76 3.82 3.61 3.31 3.59 3.99 3.60 3.82 3.81 3.77 3.44 3.71 3.85 3.76 3.32
Resultados del modelo
RES
SRES
.20 -.15 -.60 .19 -.22 .16 .01 -.01 -.05 .01 -.28 .16 -.01 .31 -.08 -.33 -.29 -.06 .06 .07 .16 -.04 .07 .33 .31 .44 .50 .28 -.24 -.31 -.16 .19 .27 -.03 -.17 -.25 -.20 .01 -.07 -.09 -.15 -.24 .09 .07 .14 -.03
.92 -.70 -2.71** .85 -1.01 .74 .05 -.02 -.22 .03 -1.25 .74 -.05 1.41 -.37 -1.48 -1.34 -.26 .28 .30 .71 -.18 .32 1.49 1.41 1.98 2.26** 1.26 -1.07 -1.39 -.73 .88 1.25 -.12 -.78 -1.12 -.91 .04 -.33 -.40 -.70 -1.08 .40 .32 .62 -.14
105
Se observa en esta tabla de piezometría predichas son, los mínimos máximos
lo
cual
sobre
se
interpreta
de
modelo. los la
identificación, Así
mismo,
valores
al
grandes,
seguridad,
que
se
que
las
máximo
equilibrio
y
que lo
sacaría vaciaría.
valores
existencia
en
al
más
de
han encontrado
los
en
resto
sobre
encuentra
todo
del
predicciones
cálculo del bombeo de equilibrio no valor
superiores
los
siempre, se
diferencias
inferiores en los
forzado
han
modelo
asegurando
la
por
ser menor el
e
para
que
las
siempre,
todo,
valores extraños y embarrados fase
casi
que
a los valores observados,
observados,
grandes,
11.8.1,
la del
para
lado
de
para
el
sobrepasarán nunca el
al
acuífero
Lo
que
se
de
su
estado
traduciría
en
de una
ligera infravaloración de la recarga.
En
la
figura
9
se
representa
observados con los predichos.
la
gráfica
de
los
valores
106
4.05L N D p z
1 1 1 2 1 1 1 1 12
11 1 1 1
1
3.6-
1 1
1
2
1
1 2 2
1 2 3.15
1
1
1
1
1
1 1
1
1
1 3.375
3.25
1
1 3.625
3.5
1
1 3.875
1 1 4.125
1 4
3.75 LNDPZ*
Figura 9 - Grafica de valores observados predichos
En alrededor se
ella de
se una
encuentran
inferior
y
observa recta
dentro
superior
de de
que
los
ascendente las
rectas
y
puntos que que
incertidumbre,
corresponde al valor extraño observado de piezometría.
1 valores
se
todos marcan excepto
en los
agrupan los el
pares límite
uno,
que
incrementos
107
11.5.-
ANALISIS
DE
RESIDUOS
Y
VERIFICACION
DEL
MODELO-
LINEAL
En la
la
figura 9 del apartado
representación
piezometría
gráfica
frente
a
los
de
anterior,
los
valores
se ha mostrado
valores
observados
predichos,
destacar ninguna particularidad especial
no
de
siendo de
y pareciendo que
el ajuste es correcto.
En de los
los
la
figura
residuos
primeros
10.1
y
frente
a
la diferencia
los observados.
10.2 los
se
presentan
valores
entre
los
RES - LNPDZ*
- LNDPZ
gráficas
predichos,
valores
Es decir:
las
siendo
predichos
y
108
ACUIFERO DE JIJONA 4,3
ResOduos del modelo
4,21 4,143.93.83,7 L
3.6-
9 1(23 11 [311 45 22 El 21 11 241 j9 2811 17 43 34 55[] "40 Va CIC] 03 5 11 49 540 El ¿�711"T . 5100 3 El 29 El r 44 12
3,5-
36
50 31
3.4-
20
435 g
3,3-
52 11
3.23,113 0. 8 - 0. 6 - 0,4 - 0,2
0
0,2
0,4 0,6
LNDPZ* - I-NIDPZ Figura 10.1 - gráfica de residuos predichos,
frente a valores
con numeración
0,8
1
109
ACUIFERO DE JIJONA ResOduos de¡ modelo 4.34.24.143.93.83.7L
09-89 87 11 w=:W. 11 o,,12&Wkj 11 01 -88 Ib1 11-86 El u. job %3-190 09-1 o%-s-%8 11 08-8W~-86
3,63.53,43,3-
12-88 02-90 11 OV-8,09-1313 07-CWW-87[3 [Di
3.23.13 1 -0,8-0,6-0,4-0,2
0
0,2
0,4
0,6
LN13PZ* - I-NIDPZ Figura 10.2 - gráfica de residuos predichos,
frente a valores
con fechas
0,8
1
110
En
la
figura
10.1
el
número
que
aparece
en
cada
punto se
corresponde a la numeración secuencial
del
mientras
que
colocando
la
figura
10.2
es
la
misma,
pero
dato,
encima el mes y el año correspondiente.
En
estas
figuras
no
se
configuración
especial
de
los
puede afirmar
que
el
ajuste
correcto.
Unicamente
Noviembre
de
corresponde
punto
respecto
forma
con
que
se
lineal
es
punto
11,
de
la
nube,
al
que
aparece anómalo
ninguna
residuos,
destacar
1986,
al
observa
modelo
el
aparte
de
descrito
lo
en
la
primera
o
y
fase
de análisis y evaluación de la información disponible.
Se
observa que para los otros valores
en
la
fase
(Marzo han
de
-
identificación
Abril
predicho
formación
del
87,
valores
de
un
que
Julio
aparente
han
del
bastante
extraños
89
forzado y Marzo
ajustados,
grupo
de
En
la
presentado los
residuos
valores
y
que
ya
se
del
del
90),
obligando
a
se la
anómalos
a Diciembre
apartado
predichos
al
Se
observan
punto
ha hablado,
11, y el
anterior
por el
estudentización,
extraños.
correspondientes del
11.8
valores su
modelo
del
88
que sin embargo son valores correctos.
tabla los
el
valores
correspondientes a los casos de Agosto (números 32 a 36),
observados
a
fín
dos
valor punto
modelo, de
35
han
así
como
analizar
valores de
se
extraños,
Noviembre que
los
del
86
pertenece
al
grupo
descrito
de
correctos
valores
embarrados
por
el
grupo defectuoso descrito en el párrafo anterior.
Se
que
observa
(observación
valor
el
sólo
representa
11)
predicha,
siendo
este
valor
del
residuo
valor
predicho.
el
Es
decir,
el
del
residuo
la
variable
desfavorable,
y mínimo
el
caso más
(observación
máximo
21)
15%
de
solo
es
del
predice
modelo
los
observados con una exactitud que va del 2.6% al
En
la
estadísticos
tabla
11.9
que
definen
por
presenta,
se el
test
de
2.6%
del
valores
15%.
último,
los
de
normalidad
los
residuos:
Número de observaciones válidas - 46.00 Variable
Media Desv.
RES
.00
Tip.
Asimetría Kurtosis .10
-.03
.22
Tabla 11.9 - Estadísticos de los residuos
Se
en
esta
0,
es
decir,
una media no
se
de
encuentra
ni
negativa).
Presentan
hablado
en
estandar el
que
curva
los de
residuos
valores (media
por encima
de
debajo
(medía
también
una
desviación
de
apartado
la
regresión
anterior),
y
del una
poseen
predichos
de
por
siempre
ni
(error
la
siempre
positiva),
0.22,
tabla,
aprecia
residuos residuos
típica
que
ya
se
asimetría
de ha y
112
kurtosis
inferiores
distribución Por ello
la
normal serie
los
a con
de
46
límites casos
residuos
se
y
requiere
que
un
95%
acepta
confianza.
de
como
una
normalmente
distribuida.
11.6 - EXPLOTACION DEL MODELO.
En valores el
punto
la
figura
observados de
control
por el modelo:
11.1 de NQ
se
las
presenta diferencias
283380028,
la
gráfica
de
junto
de
piezometría a
los
los en
predichos
113
ACUIFERO DE JIJONA RESULTADO MODELO, DIF. PIEZOMETRIAS 50
4-0
30
20
10
E N
0 -10
A
V
-20
-30 01-86 -
Di -87
........ . . . . . .Di. -88
Dif pz (2833a0028)
+
Figura
11.1
-
Diferencias
01-89
01-90
Dif. predichas
observadas
y diferencias
predichas
En
esta
figura
del
lado
sitúa,
simultáneamente a
igualar
En
los
la
a
la
picos
se
puede de
serie más
figura
11.2
piezometria medida
junto
observar la
original
que
el
seguridad, pero
sin
modelo
se
oscilando llegar
nunca
altos.
se a
presenta
la predicha.
la
gráfica
de
la
114
ACUIFERO DE JIJONA 670 660 650 640 630 620
E
RESULTADOS FINALES DEL MODELO -
600 590 sao 570 560 550 54-0 530 520 510 500 490 01-8e
01-87
01-90
01-89
DI-a8 meses
Valores observados
+
Valores predichos
Figura 11.2 - Piezometría observada y piezometría predicha
En
ella
siempre
por
se
debajo
puntos mínimos, del
lado
gráfica, lo que
de que
observa de
que la
la piezometría medida,
seguridad.
Se
aprecia
la forma de ambas curvas es
la estructura de
la variable
ha sido calculada por el modelo.
llegando
se observa que
con lo que también
la
pero
predicha
también
a
está los
estamos en
muy similar,
la por
observada piezometría
La
ecuación
estimar
el
bombeo
equivalente considere.
El
obtenida
a
de
los
recursos
posible de
de
el
del
del
medios
modelo
permite
acuífero,
del
período
valor que
se
se define a continuación:
equilibrio
bombeo en
determinado
ajuste
equilibrio
Este concepto
bombeo
del
significa
el aculfero,
incermento
de
aquel
para que
piezonetría
valor
en un período
entre
el
origen
y final del período sea cero.
Este
bombeo
ecuación
(ii)
(apartado
11.4),
puede obtenida
de
pz(n)
la
que
bmeq,
es
el
partiendo
fase
para
pz inicial y pz
equilibrio,
que satisface
calculado
en
haciendo
quiere considerar, bombeo
ser
el
de
de
la
calibración
período
que
se
final sean iguales.
El
valor
de
bm
1n
el
constante
las condiciones siguientes:
- pz(0)
- 0
Y
(M)
0 =
donde
t-n exp 12.826666+0.19238-1n(pv(t)+10)-ln(0.000007-baeq(t-8)-0.43)]-40
1
exp
es
la
función
logaritmo neperiano,
Los
bombeos
de
pv(t)
exponencial
de
base
e,
la pluviometría en el mes t.
equilibrio
calculados
ecuación anterior son los siguientes:
mediante
la
116
periodo 86-90 (calibracida) : 87.401 J/xes = 1148 H§3,año
81.921 131tes = V983 H§3/aio
periodo 81-91
Con
el
fin
de
determinar
el
períodos más amplios de tiempo, piezométricas
medidas
(28338028)
el
Alicante, junto
el
que
con
estación estos
por
las
de
valor
se del
piezometría éste
abarcan
y
coincide
valor final.
se
ese
de
observa
año
mismo
1981
de
con
el
que
la
con
al
año
recogidas
período.
agua
A
modelo,
de
1991, en
la
partir
de
utilizando
del
período,
la
evolución
predicha
por
equilibrio
sensiblemente
para
las medidas
punto
mensuales
reconstruido
bombeo
modelo
han tomado mismo
el
pluviometrías
ha
del
la Diputación Provincial
y
desde
para
se
este
I.T.G.E.
Jijona
datos
en
ajuste
la
real,
coincidiendo
el
117
ACUIFERO DE JIJONA EVOLUCION PIEZOMETRICA SEGUN BOMBEOS 700-
680 eco64-0 620COC) 5150 E
CE
:I
560 54-0 - 0 520-
12
500 4430 -
13
+60 4-4-0 4»20 4-001981
19e2
1ga3
1984
E]
I.T.C.E.
1985
1986
1957
1988
1989
1990
1991
0.ga3 hM3
+
Figura 16 - Predicciones en la evolución piezométrica
Los valores se
de piezometría y pluviometría utilizados
incluyen en la tabla 9 del Anejo.
En
la
figura
piezométricos
del
17
adjunta
ITGE-DPA
Jijona para el mismo punto,
y
se los
representan del
(28338028).
los
datos
ayuntamiento
de
ACUIFERO DE JIJONA PIEZOMETRIAS
283340026
700 680 660 640 620 600 580 E e
560 540 520 500 480 460 440 420 400 1........ ...... 1981 1982 1983
1954
1955
1986
1987
1988
1989
1990
1991
t [Mes-es] Calibraci5n
12
I.T.G.E.
+
Figura 17 - Comparación de medidas piezométricas
En esta figura se puede apreciar
como
las medidas
el ITGE-DPA presentan una diferencia de cota del unos
metros
90
construcción
con
del
respecto
a
modelo,
presente
orden de
empleadas
las
teniendo
de
en
ambas
la
series
de datos una evolución bastante similar.
En
siguiente
la
evoluciones suponiendo 1981-1991
figura
piezométricas que mes
que
el
modelo
pluviometrías
se
repiten
las
a
mes,
manteniendo
y
presentan
se
18
un
las
predice,
del bombeo
periodo anual
119
constante
equilibrio, último año
igual
en
el
primer
caso
y en el segundo a lo bombeado
al
bombeo
de
realmente en el
(1991):
ACUIFERO DE JIJONA EVOLUCION PIEZOMETRICA SEGUN BOMBEOS 550 540 5.30 520 510 500 490 480 E d Ld
470 -
-É
450-
450 440430420 410400390 380 370
........... .................................................. 1992 l�;3 1997 19941995 1996
1998
1999
......... 2000 2001
t [meses] 0,9133 Hm3
Figura 18
0Predicción
I,W HM3
de evolución piezométrica
2002
120
11.7.- CONCLUSIONES AL MODELO.
Las conclusiones y recomendaciones trabajos
surgidas de
los
de modelizaci6n lineal de las piezometrías,
función de las pluviometrías y los bombeos,
en
en el
acuífero de Jijona han sido las siguientes:
Res¡>ecto a la calidad de la información
lº)
disí>onible
1.1 Variable Piezometría.
a)
Los
datos
calidad otro
piezómetro
razonable,
punto
bombeos
del
y
y
acuífero a
suficiente,
la
vez
dado que
28330028
número
que
esté
no
existe
menos
disponga
de
en
este
influenciado
una
serie
una
momento por
los
piezométrica
recomienda realizar un control
se
tienen
sistemático
con medidas mensuales.
b) la
Se
ha
encontrado
variable
que
diferencias
transformación
la
de
piezometría
período Enero-1986 a Junio de
ln(dpz(t)
donde
dpz(t)
piezometría
es en
la el
normaliza
es,
para
el
la siguiente:
1990,
+ 40)
diferencia momento
que
t y
entre el
la
observación
anterior
t-1,
de
por
la ser
121
la que
más
normaliza
la
pluviometrías y bombeos,
serie
y
siendo
más
se
correlaciona
la variable
con
a utilizar
en
técnicas y modelos estadísticos.
c) La variable piezometría del
acuífero
correlativa con tendencia
Jijona,
de
de punto utilizado muestra
los períodos
ascendente
de
durante
una
periodicidad
lluvia y de período
el
(283380028) anual
sequía, de
y una
ajuste
del
modelo.
1.2.
Variable pluvionetrlas.
Los
datos
estación
de
máximos de de es,
cada
de
pluviometrías
Jijona, lluvia
año,
y
son
en
la
los
de
correspondientes calidad.
buena
meses
transformación
para el período
Presentan
Septiembre
de
que
la
unos
a Noviembre
normaliza
estudiado de Enero de
a
la
serie
1986 a Junio de
1990.
ln(pv(t)
donde
1n
es
pluviometría
en
periodos
producen Agosto,
el
y
que
en
+ 10)
neperiano
logaritmo
un
mes de
dado.
sequía
los meses
de
Es en
de
los
máxima,
presenta con un caracter torrencial.
y
pv(t)
observar meses la
de
que
la se
Abril
a
pluviometría
se
122
1.3.
Variable bombeos
Los
datos
explotados
de
en
considerado,
bombeos
el
de
corresponden
acuífero
Enero
de
para
1986
todo
a
Junio
totalmente antrópicos y tecnológicos, claramente Dirac de
cada
para el
con
0
la
determinista el
"golpe"
año,
que
Dirac,
meses
corresponden
industria
de
los
turronera,
que
a
de
corresponde
valores
el
período
de
a
1990.
una
Octubre
la
y una
los
función
deriva
de
agua
exponencial
consumo
de
a Diciembre
utilización
al
son
siendo una variable
correspondiente
en
a
en
de
agua
para
que
para
este
considerado,
la
abastecimiento urbano.
Tras
muchas
acuífero,
mayor
dentro
y
transformación
pruebas
que
correlación
más con
se del
ha
encontrado
período
normaliza la
la
variable
serie de
y establece diferencias
la de
piezometrías es:
ln(0.000007
donde
1n es el
un mes dado.
-
bm(t)
- 0.43)
logaritmo neperiano y bm(t)
los bombeos
en
123
2.
lineal aiustado
Respecto al model
ha
Se
diferencias
de
del
dentro
que
la
piezometría
se
observado
en
pluviometrías, estudiado
de
1990,
de
transformada
de
período
el
en
Junio
a
variable
la
con
64%,
del
1986
de
transformada
variable
la
de
positivamente
correlaciona
grado
un
Enero
de
negativamente
con
mes
mismo
transformada
variable
y
bombeos,
con un retraso de 8 meses y un grado del 47%.
en
Igualmente,
análisis
el
de
correlaciones
las
cruzadas se ha visto que las máximas correlaciones
con
8
de
meses
retraso
pluviometrías
transformada de
el mismo mes con la serie
con
variable
la
son en
transformada
y de
bombeos.
cuadrados
de mínimos
Como resultado del ajuste
se ha
llegado a la siguiente ecuación:
pz(t) - pz(t-1) - exp[(2.82 + 0.19-la(py(t) + 10) - 0.06-la(0.000007-bl(t-8) - 0.43) + E)] - 40
donde
exp
neperiano,
es
la
las
la
bm(t)
estimación,
predicciones
en
piezometría
la
pz(t)
la pluviometría, en
exponencial
en
los bombeos,
siendo del
base
el
modelo
error de
por el análisis de residuos que
e,
1n
instante
el
y E el
y
logaritmo t,
pv(t)
error cometido
estandar
0.23,
el
estimado
habiéndose
de
visto
el modelo es correcto.
124
3.
Explotación
El
del
modelo
modelo.
estocástico
realizado
permite
reproducir
la
piezometría en el aculfero en función de las variables pluviometrías y bombeos.
En consecuencia conociendo dos de
las tres variables
se
puede estimar la tercera.
Por otra parte permite calcular el bombeo de equilibrio, equivalente
a
considerado,
los recursos si
renovables
valor
del acuífero para el período
se establece que el incremento de piezometría entre
el origen y final del período es nulo, por lo que:
Pzn-PzO=0
(Y) 0 z
donde
Z exp[(2,82 + 0.19.1n(py(t) + 10) - 0.06-la(0.000007-ba(t-8) - 0.43) + E)] - 40 t=1
exp
neperiano,
es
la
pv(t)
función es
la
exponencial pluviometría
de en
base la
e,
1n
unidad
el
de
logaritmo
tiempo
t
y
bmeq(t-8) es el bombeo de equilibrio mensual.
Para conocer el bombeo de equilibrio en el período considerado, 1981 a 1991, la ecuación
se
introducen los valores de pluviometría mes a mes en
anterior
(v),
distribuidos mensualmente
y de
se
tantea
acuerdo
un
con el
valor
del
sistema
bombeo de
anual,
explotación
125
de
aculfero,
hasta
que
pz(n)-pz(o)
sea
igual
a
cero.
El
valor
obtenido para el período (1981-1991) es de:
bme = 0'983 HM3/año.
Otra
posible
píezometría
a
pluvionetrías. series
de
simular ello
la
se
aplicación partir
de
Partiendo
bombeo
y
de
reservas
una
captaciones
pueden
la
movilizar
de
si
inicial
prefijadas, en el
el
de
bombeos
piezometría
resultante
para
estimación
serie
una
determinar
suficientes
riesgo de secarse.
la
pluviometrías
piezometría podría
es
se
acuífero.
aculfero
regulación suficientes
y/o
la y y
puede Con tiene
si
las
reservas
sin
126
127 La unidad de Jijona-Carrasqueta acuíferos:
Jijona
y
está
Carrasqueta,
formada por dos
ambos
de
diferente
litología.
está
El acuífero de Jijona pararrecifales
del
Carrasqueta
por
las
Zonalmente
ambas aunque
contacto, se
Mioceno
que
descartarse
el
y
se
de
éste entre
apreciable
hidráulica
en
ponen
través
a
de
Medio(M3).
Mioceno
del
acuíferas
formaciones
comunicación
produzca
Inferior(M1),
calcarenitas
parece
por las calizas
integrado
ambas.
acuífero
El
hm3/año, año
1991
que
se
estimó
un
con
valor
el mismo
periodo ;
regadíos
del
valle
del de
forma natural por el manantial
de
de
estimado agua
este
de
del
0,983
de
acuífero
medio
parte
recursos
explotación
la
a
río
El
hm3/año.
1,06
en
Carrasqueta se drena de Nuches
inferior
ligeramente
valor
unos
de
1980-1991
periodo
el
para
renovables
medios
dispone
Jijona
de
Coscón
aprovecha
se
parte
y
hm3/año
0,73
en
para los
recientemente
en la recarga del acuífero de Jijona.
aguas
Las
prácticamente hm3/año, donde
al
son
sin cauce
depurar, del
en
río
parcialmente
pequeñas huertas.
de
residuales
un
Coscó
vertidas,
Jijona
son
volumen
próximo
y
del
Bco.
aprovechadas
en
el
a
0,85
Purgatorio, riego
de
128
Se
considera
depurasen
y
se
actualmente has,
se
si
utilizasen
abastecida
podrían
de Carrasqueta, acuífero
que
de
las en
con
el
el
liberar unos
aguas
residuales
regadío
manantial
de
de
m3/año
480.000
aumentando
sus
la
zona
Nuches, del
que se podrían utilizar en la Jijona
se
75
aculfero
recarga del
disponibilidades
hidráulicas.
Se de
ha
realizado
Jijona,
entre
la
que
un
modelo
reproduce
pluvíometría
y
estocástico
una
función
bombeos
y
de
del
acuífero
transferencia
piezometría.
la
Este
modelo podría ajustarse más aún en el futuro mejorando el control de
Dado
los datos de bombeo y piezometría.
que
los
Ji j ona-C arras queta que
existe
plazo unidad
de
un
están
aumento
carácter
debe
recursos
de
aprovechados la
demanda
industrial,
clasificarse
nuevas concesiones,
regulados
como
por lo que
se ZONA
de en
su
la
unidad
totalidad,
previsible
corto
que
esta
AUTORIZADA
para
se propone el
protección que se adjunta en el epígrafe
y
a
considera NO
de
10.
perímetro
de
ANEJO TABLAS Y DOCUMENTOS
ACUIFERO DE JIJONA - CARRASQUETA Series empleadas
NUM MES 01 02 03 04 05 06 07 08 09 10 11 12 01 02 03 04 05 06 07 08 09 10 11 12 01 02 03 04 05 06 07 08 09 10 11 12 01 02 03 04 05 06 07 08 09 10
AÑO 1981 1981 1981 1981 1981 1981 1982 1982 1982 1982 1982 1982 1982 1982 1982 1982 1982 1982 1983 1983 1983 1983 1983 1983 1983 1983 1983 1983 1983 1983 1984 1984 1984 1984 1984 1984 1984 1984 1984 1984 1984 1984 1985 1985 1985 1985
AQUAGES 8028
I.T.G.E. 8028
541.26 541.5 533.96 527.42 521.98 516.27 483.9 483.25
485.7 542.9 533.3
525.1
513.8 506
494.3
480.1
PV
26.7 37.3 0 30.6 11 2.9 0 6.7 56 30.9 48.8 80.7 113.1 0 0 2.5 28.9 171.5 16.5 0 0 10.7 5.9 13.8 11 14 0 63 0 10.6 80.8 12.4 5.7 48.3 31.4 14.5 94.6 6.7 0 25.5 31.1 2.3
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40
11 12 01 02 03 04 05 06 07 08 09 10 11 12 01 02 03 04 05 06 07 08 09 110 11 12 01 02 03 04 05 06 07 08 09 10 11 12 01 02 03 04 05 06 07 08 09 10 11 12 01 02 03 04
1985 1985 1985 1985 1985 1985 1985 1985 1986 1986 1986 1986 1986 1986 1986 1986 1986 1986 1986 1986 1987 1987 1987 1987 1987 1987 1987 1987 1987 1987 1987 1987 1988 1988 1988 1988 1988 1988 1988 1988 1988 1988 1988 1988 1989 1989 1989 1989 1989 1989 1989 1989 1989 1989
471.7
462.95 532.70 523.30 517.40 524.80 519.70 504.40 508.00 492.50 502.10 531.30 569.80 559.90 580.80 583.85 577.90 573.10 583.10 570.60 575.80 559.30 564.80 563.10 587.90 615.80 637.20 643.55 639.50 645.20 648.80 653.20 640.40 621.90 615.70 606.10 593.30 577.70 589.30 603.40 621.30 616.60
444.8
406.1
503.85
495.8
483.1
563.55
519.7
94.2 0.3 13.9 17.9 10.9 42.7 101.4 3 2 0 13.6 60 89.1 26.9 9.5 6.7 8.2 43.8 7 7.9 40.5 6.5 153.6 135.5 18.1 1.3 42.4 54.1 0.3 1.2 45.3 0 20.6 0 36.2 90.4 191.4 62.1 40.4 30.3 10 60.7 78.5 86.7 3.4 11.6 56.3 61.4 45.1 0 23.7 19.8 80.6 34.1
41 42 43 44 45 46 47 48 49 50 51 52 53 54
05 06 07 08 09 10 11 12 01 02 03 04 05 06 07 08 09 10 11 12 01 02 03 04 05 06 07 08 09
1989 1989 1990 1990 1990 1990 1990 1990 1990 1990 1990 1990 1990 1990 1991 1991 1991 1991 1991 1991 1991 1991 1991 1991 1991 1991 1992 1992 1992
611.40 609.22 601.80 608.30 634.40 630.50 639.30 648.80 659.20 658.93 656.40 660.30 657.70 646.37
513.22
534.5
562.93
550.37 537.9
554.42
549.2
528.1
51.6 35.7 0 32.8 161.8 12.4 59.4 58.5 43.2 0 24.9 72.5 41.3 0.5 1.4 18.8 10.9 34.2 8.2 6.1 57.1 29.3 59.5 5.2 5.3 7.3 0.7 0 9.5
TABLA 9 - PIEZOMETRIAS DEL Y PLUVIOMETRIAS EMPLEADAS
C QID L IP
IP ] CHAS
TJ
] EAD
D
]:
cD]2-
QTJE PA.
ARCHIVO DE PUNTOS ACUIFEROS
Coordenod-as º*graficos
13
N* de registro .............. instituto Tecnoicq1co l0 ".' ��, ',
r
o
FM119 H,24 11 125
1
Fecha
Alíturo de¡ agua
Caudal
(esp e CtO a 1.0 referencia
M3/
Cotaobsoluec
Aktodo
de¡
.............
agua
meld.
....................... ...... .............. ..........................
mg 1 131
w
L3)
138
¡42
t5Q.
(54
W
159
...... .
................................... ... ...............................
149
148
143
............. ENSAYOS
BOMBEO
DE
................
1 horas
Dvraci6n U bombeo
1 163
1
.......................
.. ............
............. ...... ........................ .. . .. 1 .............
............... ................... ............................ . ..................... ................. ................................................. ............... : .................................................. 1 ....... 1 ..........
Fecha Caudal extroido (m3/h)
........ ...
................................................ .... .... .
1
........................................... 1 ......................... .............. . ................................................................. ....... .......... . .................... .......... ..............
minu.
. Depresión en M.
..................... .......................
Tronsmisividad (mYseJ -rROVECIIAI.II1:-14'1*0 ITMPORAI. DE AGIJIJ,� PRIVIIIii
D.ANTONIO BERNABEU GOMEZ, en representación del Ayto. Jijona ................................. . .............. ............ P- 0308300-C « J con D.N.I. 6 C.I.F . ........................... , domiciliado en ........ 03100 Avda. Constitución, 6 calle .................................................... C. 1« .......... Don
Dentro del plazo establecido por las disposi.c,�*.ortes
trari.,-,itori¿I-,i
como titular de un aprovechamiento de aguas pr . das con las* siguientes car acteristi cas:
y 3§ de la Ley de Aguas,
-
Lugar o paraje, término municipal y pvovincia en que je en ESTaW
DI A
S~A LECT. NIVEL
IÉ tos
Hoila nºw>*4: AÑO:
Pozos
PINETA NIVEL LECT.
EL ROSE T NIVEL LECT,
-2
Lkl
ic
x li
z
3 (y-)
2b J01)
(e)
11
laz
si
3s2xfi9c
SS
W25
c,'� 4 5 Z-tv(>
s(4
qQvo (ps
(11) tjb2-CAc> 65 (-k�
b(>-Z� 11D
SS
ski' u. s
131 -Z1 A i
11
(P) 3*�*Z9-23c -33
4'JJ 1-,
S' Zitm 5006-7 0 10 S 5 1-P 2
N,
12,9 -y q
3
(t-�
ll�-o
Gis- (e) siZ15,10
s�
9,15
(MURVO
ESTADO
=A
S~A NIVEL LECT.
ggj
]DÉ rOS
PINET.A LECT. NIVEL
Gi�i- (e�
CSOIOI-70
G615(P)
Hoja nº AÑO: J A t�
POZOS EL ROSET LECT. NIVEL
(MSERV*
zL
'�33�79c> -4-31 to
0[1
71'2,c,
C 03 17-0
G4 (P)
is�>"(0
GL
-11,yo
-t
--im-s 2¿
11
4,011
cl
Gis9fo
U2010
1
54 140
-10'5
7.V95 3.3 9 1,17í; o
q� 138
V3 ip 4-
-z /kf 2'�
2,
401 �
3(-,t«
f(210 1,11,7s-
eo
2 z,
IA- m- a,16 c 190* m - cq, li/
o60 a
6- A,
(P
q 350 0
6 o, i;o (P)
w Y-XA
=A 21-G
'IC
SEWA LECT. NIVEL 67.6.4(M)
6 Giq:S
2 b
6 i.J 1
6-
6 5, -¿1 40- 4-
PINETA LECT. NIVEL
j,&s.3:>a 4t�4,010
26.-568
4?S-SOC2
-52,63
0
bpc.jso
EL ROSET IECT. NIVEL lo,-To(p)
6'3,2
e),9
52
Hoja,iiº AÑO:
POZOS
IS LOS
ESTADO
6 0,-2
Z
6
¿3 o j) o
3 3 G.1 1
0
)1� 0
-2
6
6514Z V 46-6
10) D¿3.¿/0
0 yjj.j3 0
ci
Jo
61..1,0
¿3_5
t; 5 c, o
6ci
510
e,
A2# jil 6
nio p
ip
Ji-', soc
c
'4 .2, cl
oj 0
260
2 u -2 o
15
60
CIO Ejo
q0
00
Cio. 3-2o
6
o
q0
800
�7
so, ti
C)
o
p
z,� 7 l� 30-
3.5 -4 qA, 8 -7
C750 2�8 .230
82
--o
q0 i�6
09,30
�7 /
o ir
1-ps tl
cu
52 J12
OBSERV.
ESTAW
DI A
Poz2i
PINETA NIVEL LECT.
S~A NIVEL LECT.
6.11,116
320 - tih t7
14- 1 Z
iz Los
Hoja nº A AÑO:
'91733 EL ROSE T NIVEI, IE()T. q q 3
92..430
01 3 1
138,40
)-390.100
pi
J.lt-20.400
100
Á, -12 Í 8
AA
-1 t,
320
0
11.20
130
L3 2
14
4 � 9,3 3 -3 ?>0
Af 62 :520
)02
-57 los
J ¿?0
103.-
no IP
1/1-2
7,9,50
J-3
m 74-1
98.-Vilo
P 1033-30
ql S
p
650
81' p-
loss
»h -- ? 0
10
270
-9
150,
27
¿ 90-61.p 11¿(2,go. s>o-?q-p p t; 8,� q 70 j�cos -0) D
q o, .4 0
-4. 0) 0
po. 0
iz g -0. � 3 �.615 220
ti g
q 5
61,25 m 100,1(0
Ll
Q CI, 15 0
/y c
o o
82
).01.670
q q
o
jo 1o,
Ao m
so
J67 72v
cLí
OBSERV*
NOR TE
í-� -Ocio~
17
7,
S. E.V,
Sondeos
exsterites
Sondeos recome-ndadíos
Sondeo
posible
Perímetro
Cortes
protección
Hidrogeo!,,-,gic-,s
<
REAU,
REALIZ.ADO PARA FECHA
19-4-84 E
1/150.000
AY. J 1 JON: A JOSE A. HERVAS
DIBUJADO
JOSE FUSTER
COMPROBADO PLANO N.*
PLANO
1
E
SI TUACION -E
PROTE-1-
y
DECLARACION DE LOS APROVECIJAIMIrNTOS DE AGIJAS PRIVAI)t,�.; PA-qA SI] INZ EN EL REGISTRO DE AGUAS COIMO "AI:>RC)V§-Ciltlt'ITI':141'0
DE AGO/J; PnIWI
AIMNIO BE~EU GONIEZ , en representaci6n del Ayto. Jijona. .............................. P~ 0308300-C JIJONA con D.N.I. 6 C.I.F . ............................. domiciliado en ....... 03100 6 calle. ....... .............................. o C. 13 ......... Don
Dentro del plazo establecido por las disposiciones trannitoria,y 32 de la Ley de Aguas',
como titular de un ¿ií)rOveckllliii(-,nto, de
das con las* siguientes caracteristicas:
-
-
-
Lugar o paraje, térinino municipal y provincia en que Le eju.L tra la captación de agua. Par-aje de Sereña-Jijona (Alicante). Destino de las'aguas, (Usos domésticos; abastecimiento (le P(ci6n, especifícando el número de babitantes; regadío, cando la superfic.ic y el tipo de cultívo; etc) Abastecimiento de poblaci6n de 8.000 habitantes:
-
Profundidad del pozc, o longitud de la galería. 434 m. Mo del comienzo de la explotación (si no se conoc^e, aproximadamente) ., afto: 1.972
-
Caudal máximo y volumen anual
estinar
(si se con0ue).
45 lit./seg. S 0 L I C
T A:
Su inscripción en el Registro de Aguas como aprovecliairiiento'te:r:: ral de aguas privadas,
para lo que se acompaña la dc)cLii;ifi.il:,ici6ti --cret-;�:c,
va del título legítimo del dereclio sobre el aprovechamiento de aíotias das que se venían utilizando con anterioridad al. 1 de enero de 1986*, tente en (Registro de la Propiedad, nas,
Documento Público,
interesado.
Fecila y firma. EN L r 1',0,,, �k
-2t:
L
r,elii