P
petrolman
Guest
Salut a tous,
un probleme existenciel me pose des difficultes et je pense que l'on peut le resoudre par excel (et par fortran) mais comme je suis ''debutant''.....
Le but final est de trouver du petrole dans le sous-sol...avis aux scientifiques et aux curieux, et je sais qu'il y en a ds ce forum:
Je pose le probleme assez difficile et qqe peu long ds l'explication...
Imaginons 5 couches (5 lignes d'une colonne A avec des vitesses Vi correspondant a des profondeurs Zi inscrits ds une colonne B). Je mets un puits droit (abcsisse=0, ordonne=0,la normale) a travers ces 5 couches.
Probleme direct:
Je balance un rayon avec un angle theta5 (cad depuis le fond du sos-sol, ce theta5 est le seul que je peux fixer et je le mets en C5) depuis cette ordonne x=0. Si je calcule bien (et si je connais mes formules trigo), le rayon arrivera sur la couche 4 avec un decalage (depuis la normale) de x5 (colonne D5)=profondeur entre couche 5 et 4 (Z5) *tan(theta5)=Z5*tan(theta5). (Excel: D5=B5*tan(C5)).
Mon but est de connaitre x1!!!
Pour cela, je connais une relation:sin (theta5)/V5=sin(theta4)/V4 et plus generalement sin(theta Z+1)/(V+1)=sin(theta Z)/V donc si je remplace, je trouve dans la couche 4, x4=Z4*tan[asin(V4*sin(theta5)/V5] et ainsi de suite jusqu'a obtenir x1.TB
Mon but : Probleme inverse : Me fixant X1 (et en utilisant les relations precedentes), je desire connaitre theta5!!! Ligne par ligne je sais le faire en utilisant le solver, D5=SOMME((TAN(ASIN(SIN($C$5*PI()/180)*($A$1:$A5)/$A5)))*$B$1:$B5) ce qui me donne x1 (c inverse evidemment), je le passe dans le solver en affichant valeur a atteindre 2000(par ex) et il me change mon C5 (theta5). Je suis ensuite obliger de retourner en D1 pour changer $C$5 en $C$4 et ainsi de suite...). Comment puis je le generer directement sans retourner par le case de depart (D1) et abaisser constamment mon angle de reference (cad C5 puis C4 puis C3) surtout que j'ai bcp plus de 5 couches!!!!
Je paye l'apero (Paris ou Pau) a tous ce qui auront lu mon message en entier et l'auront assez compris (cad auront essayer de comprendre) pour me donner des semblants d'explications ou des pistes.
Merci et a+ au bar je l'espere
un probleme existenciel me pose des difficultes et je pense que l'on peut le resoudre par excel (et par fortran) mais comme je suis ''debutant''.....
Le but final est de trouver du petrole dans le sous-sol...avis aux scientifiques et aux curieux, et je sais qu'il y en a ds ce forum:
Je pose le probleme assez difficile et qqe peu long ds l'explication...
Imaginons 5 couches (5 lignes d'une colonne A avec des vitesses Vi correspondant a des profondeurs Zi inscrits ds une colonne B). Je mets un puits droit (abcsisse=0, ordonne=0,la normale) a travers ces 5 couches.
Probleme direct:
Je balance un rayon avec un angle theta5 (cad depuis le fond du sos-sol, ce theta5 est le seul que je peux fixer et je le mets en C5) depuis cette ordonne x=0. Si je calcule bien (et si je connais mes formules trigo), le rayon arrivera sur la couche 4 avec un decalage (depuis la normale) de x5 (colonne D5)=profondeur entre couche 5 et 4 (Z5) *tan(theta5)=Z5*tan(theta5). (Excel: D5=B5*tan(C5)).
Mon but est de connaitre x1!!!
Pour cela, je connais une relation:sin (theta5)/V5=sin(theta4)/V4 et plus generalement sin(theta Z+1)/(V+1)=sin(theta Z)/V donc si je remplace, je trouve dans la couche 4, x4=Z4*tan[asin(V4*sin(theta5)/V5] et ainsi de suite jusqu'a obtenir x1.TB
Mon but : Probleme inverse : Me fixant X1 (et en utilisant les relations precedentes), je desire connaitre theta5!!! Ligne par ligne je sais le faire en utilisant le solver, D5=SOMME((TAN(ASIN(SIN($C$5*PI()/180)*($A$1:$A5)/$A5)))*$B$1:$B5) ce qui me donne x1 (c inverse evidemment), je le passe dans le solver en affichant valeur a atteindre 2000(par ex) et il me change mon C5 (theta5). Je suis ensuite obliger de retourner en D1 pour changer $C$5 en $C$4 et ainsi de suite...). Comment puis je le generer directement sans retourner par le case de depart (D1) et abaisser constamment mon angle de reference (cad C5 puis C4 puis C3) surtout que j'ai bcp plus de 5 couches!!!!
Je paye l'apero (Paris ou Pau) a tous ce qui auront lu mon message en entier et l'auront assez compris (cad auront essayer de comprendre) pour me donner des semblants d'explications ou des pistes.
Merci et a+ au bar je l'espere