= Procédures préliminaires:
== Convertir les poids d'oxydes en nombres d'atomes
== Apatite: Ca = Ca - 5/3 P
== FM = Fe+Mg+Mn
= Plagioclase
== Anorthite = Ca; Albite = Na; An% = Ca/ (Na+Ca)
== IF An% > An%_Maxi (30%, valeur généralement
observée sur les plagioclases détritiques)
== THEN
=== on calcule l'anorthite à partir de l'Albite: Albite=(1-An).Plagio
> Plagio=Albite/(1-An)
=== > Anorthite=Albite.An/(1-An); Calcite = Ca - Anorthite
= Si et Al pris par l'Albite > Si = Si - 3.Albite; Al = Al - Albite
= Al suffisant pour fabriquer l'Anorthite ?
== IF Al < 2.Anorthite
== THEN repasser du Ca sous forme Calcite : diminuer An% jusqu'à
ce que (Anorthite<Al/2) ou Anorthite=0
Calcul des Argiles:
MU = Muscovite = K.Al3.Si3 ____K2.Al4.(Si6.Al2)
CH= Chlorite = FM5.Al2.Si3
IL = Illite = K.Al2.FM/2.Si7/2
BIO= Biotite = K. Mg3. Si3. Al
KA = Kaolinite = Al4Si4
GI = Gibbsite = Al
Suivant le degré de saturation de K + Fe + Mg par l'alumine,
on calcule différents assemblages;
les cas sont examinés à saturation décroissante:
IF Az_>((15*K_ +2*FM_)/5) THEN Kaolinite-Muscovite-Chlorite-(Qtz
ou Gibbsite)
ELSE
IF Az_>((9*K_ +2*FM_)/5) THEN Illite-Muscovite-Chlorite-Qtz
ELSE
IF Az_>((5*K_ +2*FM_)/5) THEN
Illite-KFelspar-Chlorite-Qtz
ELSE Biotite-KFelspar-Chlorite-Qtz
On doit obtenir une distribution entre associations mine'ralogiques
comparables à celle obtenue graphiquement suir une projection triangulaire
Al-2.Ca-Na / Fe+Mg / K.
Al>(15.K +2.FM)/5 | KA-MU-CH-(QZ ou GI) | |||
Al> (9.K +2.FM)/5 | IL-MU-(CH ou KF)-QZ | |||
Al> (5.K +2.FM)/5 | IL-KF-CH-QZ | |||
Al> (5.K +2.FM)/5 | IL-KF-CH-QZ | |||
Al< (5.K +2.FM)/5 | BI-KF-CH-QZ |
= KA-MU-CH-(QZ ou GI):
MU = K
CH = FM/5
Al = Al - 2.CH - 3.MU
Si = Si - 3.CH - 3.MU
IF Al <Si THEN KA = Al/4; QZ =
Si - KA {KA-MU-CH-QZ}
ELSE QZ = 0; KA = Si/4; GI= Al - Si {KA-MU-CH-GI}
= IL-MU-(CH ou KF)-QZ
IF Al <(3.K - 2.FM) THEN
| = IL-MU-KF-QZ
| IL = 2.FM
| K = K - IL; Al = Al - 2.IL
| MU = (Al - K)/2; KF = (3.K - Al)/2;
| QZ = Si - 3.KF - 3.MU - 7.IL/2;
ELSE
| = IL-MU-CH-QZ
| CH = (Al +2.FM - 3.K)/12;
| IL = (15.K + 2.FM - 5.AL)/6;
| MU = (5.AL - 9.K -2.FM)/6;
| QZ = Si - 3.CH - 7.IL/2 - 3.MU;
PROCEDURE ILKFCHQZ;{Az_>((5*K_ +2*FM_)/5)}
BEGIN
IF Az_>(K_+2*FM_) THEN ILMUKFQZ
ELSE
BEGIN {IL-KF-CH-QZ}
w[9]{CH}:= (K_ + 2*FM_ - Az_)/8;
w[10]{IL}:= (-5*K_ - 2*FM_ + 5*Az_)/4;
w[5]{KF}:= (9*K_ +2*FM_ -5*Az_)/4;
w[1]{QZ}:= Sz_ -3*w[9]{CH} -7*w[10]{IL}/2 -3*w[5]{KF};
END
END;{ILKFCHQZ}
PROCEDURE BIKFCHQZ; {BI-KF-CH-QZ}
BEGIN
w[9]{CH}:= (Az_ - K_)/2;
w[11]{BI}:= (5*K_ +2*FM_ -5*Az_)/6;
w[5]{KF}:= (5*Az_ +K_ -2*FM_)/6;
w[1]{QZ}:= Sz_ -3*(w[9]{CH} +w[11]{BI} +w[5]{KF});
END;{BIKFCHQZ}