******************************************** * * * FAILLE * * * * Ecrit par * * Arnaud Etchecopar et Guy Vasseur * * Porte sur PC par * * Chistophe Tourneret et Alfredo Taboada * * puis sur Macintosh par * * Bernard Celerier * * * ******************************************** Revisions de cette documentation: 02-04-96;06-03-95;29-02-95;24-02-97;23-05-97;08-04-99 OBJET: ****** Trouver un tenseur minimisant pour un certain pourcentage de failles la somme des carres des ecarts entre les stries observees et les stries calculees. A : MESURES *********** 1. L'azimut sera un entier compris entre 0 et 360 ,le pendage un entier compris entre 0 et 90. 2. Les donnees seront lues sur l'Unite Logique 4. Le titre comprend jusqu'a 80 caracteres sur une ligne. Les mesures sont rentrees a raison d'une par ligne (un plan, une strie, nature et numero). Si un meme plan comporte plusieurs stries, il faut le repeter. Col. 3-4-5: Direction (0-180) Plan Col. 10-11: Pendage (0- 90) Col. 17: Secteur geographique du pendage(N,S,E ou W) i) Col. 21-22: Pitch (0-90) Strie Col. 26: Secteur geographique de la demi-droite de reference ii) Col. 29-30-31: Azimut(0-180) (NB: une strie d'azimuth 0 est assimilee avec une strie de pitch 0: ceci est un bug du a la convention du fichier de donnees. Eviter donc d'entrer un azimuth de strie zero) Nature Col. 35: Nature de la faille(I,D,S ou N) de la En cas d'indetermination on peut repeter le plan faille et la strie. Col. 70-71-72: Le numero en trois chiffres de la strie. Fin du Sur la ligne qui suit la derniere mesure, fichier il faut mettre en colonne 3-4-5 le chiffre 450 qui indique la fin des donnees. Exemple de fichier de donnees: NICE STATION 3 JUR. SUP. DEFILE DU CHAUDAN 170 73 W 35 N S 001 020 90 22 N S 008 170 74 E 45 N D 026 162 90 30 S D 027 450 Les terminer par 450. Voir aussi le fichier joint: fstest_data.txt B: PARAMETRES ************* Les parametres d'execution sont en format libre et ils seront lus sur l'UL 3.Il faudra TOUS les indiquer en les separant par un ou plusieurs blancs et en mettant un 0 (Zero) pour ceux qui ne seront pas utilises. Ceux qui commenceront par une lettre comprise entre I et N seront des entiers et les autres des reels. 1. ITRO IWRI *ITRO = 1: Toutes les failles sont prises en compte. 2: On ecarte les failles dont les numeros suivent. 3: On retient les failles dont les numeros suivent. *IWRI = 1: Pour avoir l' impression des mesures. 0: Pour ne pas les avoir. Si ITRO = 2 ou 3, on met une ou plusieurs lignes qui contiennent les numeros des failles separes par blancs en terminant par 999. 2. INDX IPCNT NTIR IX IPRINT * INDX = 1: Tout sera tire au hasard. 2: Sigma1 sera impose en azimut et pendage. 3: Sigma3 sera impose en azimut et pendage. 4: Sigma1 et Sigma3 seront imposes. 5: Sigma1 , Sigma3 et R seront imposes. (N.B. certaines des options 2,3,4 ou 5 sont boguees: les restrictions ne sont pas appliquees) * INPCT= Pourcentage de la mesure que l'on desire prendre en compte. * NTIR = Nombre de tirages au hasard pour fixer le tenseur de depart (100 a 200 est une valeur habituelle). * IX = Nombre servant a varier le tirage au hasard que l'on incrementera de 1000 en 1000 par exemple. * IPRINT = 0. Si INDX = 2,3,4 ou 5,il faut mettre une ligne supplementaire: 3. IDXA IPXA IDXC IPXC R * IDXA et IPXA : Azimut et pendage de Sigma1 * IDXC et IPXC : Azimut et pendage de Sigma3 * R : Valeur du rapport R Exemples de fichier de parametres: 1) 2 1 12 45 35 7 26 999 1 60 200 0 0 2) 1 1 2 55 100 1000 0 127 42 0 0 0 Voir aussi le fichier joint: fstest_params.txt