Home

DESCRIPTION ET MODE D`EMPLOI DU - Infoterre

Contents

1. a I A vm EE mue o mu en 111 2 5 Temp rature d injection 30 00 Viscosit 83 cp Estimations apr s 7 00 mois d exploitation d bit constant ae ae e Ce A me Q 50 100 125 150 175 200 225 250 m3 h dPh 5 0 10 3 13 0 15 7 18 4 21 1 23 9 26 6 kg cm2 dPe 6 2 0 3 0 4 2 5 5 7 0 8 7 10 5 kg cm2 dPap 7 0 7 0 7 0 7 0 7 0 7 0 7 0 7 0 kg cm2 lt lt e re i ana que mm mm ve mm 51 111 2 6 Temp rature d injection 35 00 C Viscosit 74 Estimations apr s 7 00 mois d exploitat
2. H SITE FICTIF FOR 3 HE 46 IE X X X X X X X X FORAGE DESIGNATION FOR5 PRODUCTION TUBAGES caidas Om 335 D 13 3 8 235 380 9 5 8 380 1599 D 7 1599 1775 D 6 R servoir TYPE COMPLETION TROU OUVERT RESERVOIR NATURE CALCAIRE NIVEAUX PRODUCTEURS 1613 0 m 1620 DEBIT 13 1627 0 m 1628 5 m DEBIT 8 1634 5 m 1637 5 m DEBIT 12 1648 5 m 1649 O m DEBIT 14 1649 0 m 1651 0 m DEBIT 10 1644 5 m 1647 8 m DEBIT 5 1665 6 m 1668 3 m DEBIT 10 1673 0 m 1675 1 m DEBIT 11 1675 1 m 1678 0 m DEBIT 5 1719 7 m 1721 0 m DEBIT 12 COTE MOYENNE 1670 m EPAISSEUR PROD CUMUL 26 3 m POROSITE MOYENNE 15 SALINITE 30 g 1 PRESSION ART OU NIV 8 5 bars ESSAI ue et PRODUCTION ET BUILD UP DATES sans ras d s rene de le 6 9 82 2h DUREE DE PRODUCTION 12 h DEBIT MOYEN 115 m3 h DATE FERMETURE 7 9 82 13h 59 35 DUREE DE FERMETURE 12 h COTE DE LA SONDE TP 1598 m TEMPERATURE RESERVOIR 58 6 C INTERPRETATION CHANTIER kh 17 0
3. E E dPh 4 3 8 9 11 2 13 5 15 9 18 2 20 6 22 9 kg cm2 dPap 7 0 7 0 7 0 7 0 750 7 0 7 0 7 0 kg cm2 an an EM Er mp E rm e 2 3 Temp rature d injection 40 00 Viscosit 67 Estimations apr s 1 00 mois d exploitation d bit constant I IILI iii dPh 4 1 8 4 10 6 12 8 15 0 17 2 19 4 21 6 2 dPc 6 1 9 2 9 4 0 5 3 6 7 8 3 10 0 kg cm2 dPap 7 0 7 0 7 0 7 0 7 0 7 0 7 0 7 0 kg cm2 OO mu OUT 50 111 2 4 Temp rature d injection 45 00 Viscosit 61 cp Estimations apr s 1 00 mois d exploitation d bit constant o E an un O
4. 11 CONTINUE ITX 0 JX 40 JY 390 DO 14 1 1 2 ITY 80 IF 14 G 2 80 IY IYX IF 2 IY IYN DO 42 1T 1 6 IX IXXC IT 4 xIDX ENCODE IBUF 4014 IX IY ITX ITY CALL ECRV24 IBUF 40 0 IF 1 60 TO 12 IP IPN IC 1T 4 IDP 10 ENCODE CIBUF 1016 JX JY IP ETX CALL ECRV24 CIBUF 22 0 12 CONTINUE 70 44 CONTINUE IF 8 1 60 TO 100 C IX IXXC ID4 ENCODE IBUF 1006 IK IYL4 CALL ECRV24 IBUF 48 0 IF 1C0D EG 1 ENCODE IBUF 1020 215 216 TEMPS ETX 1020 FORMAT LB 3A4 984 F6 1 h1 IF 1C0D E6 2 ENCODE IBUF 1020 ZL4 ZL6 TEMPS ETX CALL ECRV24 IBUF 57 0 IX 1X 240 ENCODE CIBUF 1006 IX IYLS CALL IBUF 18 0 ENCODE IBUF 1022 ZL7 LIBF ETX 4022 FORMAT 718 204 681 04 CALL ECRV24 IBUF 17 0 100 CONTINUE RETURN END RO RREO RO SUBROUTINE COURBE PH PE TE 1000 ICOUR C C TRACE DES COURBES DIMENSION PH 31 PE 31 INTEGER 1 IBUF 80 ETX IUN COHNON ETI4 IXN IXX IDC IYN IYX PN ECHPD PEN ECHPE DATA ETX 3 DATA XLIB Ti IUN C C IFL 0 IF TE GE 0 60 TO 8 1FL 1 TE TE 8 CONTINUE 00 100 1C 1 2 IXNC IXN IDC IC IXXC IXX IDCHIC ENCODE IBUF 1008 IXNC IYN IXXC IYX 1008 FORMAT 1N 16 16 16 16 CALL ECRU24 IBUF 29 0 DO 80 1P 1 31 JY 1P 4 x200 IF 27 IYD IY IF IP E0 31 IVF IY GO TO 10 12 1 10 X PHCIP PN ECHPD 60 TO 14 42 X PECIP PEN ECHPE 14 IX IXXC INTOO IF 1 60 TO 1
5. 2044 FORMAT 2X AUTRE DATE DE SORTIE 0 N 2046 FORMAT 2X SELECTION DU MODE DE CALCUL DE L INTERFERENCE 42X Hors interference le comportement propre de l ouvrage 22X est toujours calcule avec les caracteristiques specifi 32X ques fournies 42X x Pour le calcul de l interference de l autre puits S2x 3 options amp 2X option 1 k et h sont identiques 72X option 2 k et h sont distinctes et donnees option 3 k et h sont une moyenne des valeurs donnees 92 7 OPTION CHOISIE STOP END RA RR SO O ASA FUNCTION E4 U DIMENSION A 6 8 3 C 3 C FONCTION EXPONENTIELLE INTEGRALE Cossee DATA 57721566 99999193 24991055 05519968 976004E 2 1 107857 2 DATA 250421 2 334733 1 DATA 1 681534 3 330657 1 IF U 67 1 60 TO 12 Lee ell INFERIEUR 1 Et a 1 AL06 U Zei DO 10 Led 10 EA Et 4CL X RETURH C U SUPERIEUR 1 42 XN B 1 XD C 1 Zei DO 14 122 3 AN XN B 14 XD XD C I wA E1 XN XD E4 E1 UXEXP 0 RETURN END SUBROUTINE CADRE GN 0X PN PX ICOD TEMPS LIBF Causes C TRACE DES CADRES DIMENSION ZL1 4 ZL2 9 2L3 7 DP 2 IDP 2 IPN 2 2L4 3 4ZL5 3 ZL6 9 2L7 2 2L8 3 2ZL9 3 INTEGER 1 IBUF 80 ETX LIBF 6 COHNON ETI4 IXN IXX IDC IYN DO ECHPD PEN ECHPE DATA ETX 3 DATA ZL4 UDEBI T
6. au syst me MICRAL Un exemple de document obtenu par cette proc dure est donn en annexe 2 Activation WSGTIMO Le tableau des r sultats issus du pr traitement ESSOUT S peut tre inclus dans le rapport par fusion avec le fichier texte option CTL KR de WORDSTAR 8 CALCUL DES PREVISIONS D EXPLOITATION PROGRAMME PREVISS Ce traitement est ind pendant du pr c dent les donn es n cessaires tant introduites au clavier Le programme fournit les r sultats sous la forme de tableaux l cran et cr e gale ment un fichier PROUT S Ce dernier fichier peut ensuite tre incorpor directement dans le document dactylographi final paragraphe 9 Le trac simultan des courbes pression d bit et puissance lectrique d bit est optionnel Le calcul est effectu pour une exploitation fonctionnant en doublet quilibr c est dire pour un d bit d exhaure gal au d bit d injection tout instant 8 1 HYPOTHESES DU CALCUL L estimation des pressions aux deux forages est obtenue l aide de la fonction de Theis ou exponentielle int grale El et par application du principe de superposition La variation de pression chaque forage est donc la somme de la contribution propre de l ouvrage consid r et de la variation de pression due l interf rence de l autre puits Dans le cas d un forage d injection il est tenu compte du comportement la zone d extension croissante autour du puits
7. ICC 0 15 0 WRITE 7 1014 WRITE 7 1028 IF CIOPT NE LB 60 TO 62 WRITE 7 1024 ITAB J J 1 30 Bi WRITE 7 1026 DPR PWFO 60 TO 64 62 IF IOPT NE LR 60 TO 80 WRITE 7 1066 ITAB J J 1 30 CIDAT K 8 WRITE 7 1066 60 TO 64 80 IF IOPT NE LI 60 TO 82 WRITE 7 1082 CITAB J J 1 30 CIDAT K K 1 B WRITE 7 1048 PWFO 60 TO 64 82 IF IOPT NE LF 60 TO 84 WRITE 7 1086 ITAB J J 1 30 IDAT K K 1 81 WRITE 7 1088 DPR PWFO 60 TO 64 84 STOP ERCODE C C LECTURE DATA 64 CONTINUE IC 12 IDEC 0 TDEC 0 IBUF 0 TO 8 CONTINUE IC 0 WRITE 7 1014 9 IF IOPT ER LR TO 46 IF IOPT EQ LI 60 TO 66 WRITE 7 1022 WRITE 7 1016 WRITE 7 1022 60 TO 60 66 CONTINUE WRITE 7 1070 WRITE 7 1072 WRITE 7 1070 68 CONTINUE 1 7 40 CONTINUE IF FLAG NE IREP READ 9 1003 1H IPU IN 15 P IF FLAG EG IREP READ 9 1012 IPU IM IS P IF CIPU NE 60 TO 12 1C 10 1 ICC ICC 1 IF CICC GT NPHAX 60 10 22 T FLOAT IH FLOAT IH 60 FLOAT I5 3600 IF IDEC GT 0 60 TO 11 GO TO 11 IDEC 1 TDEC 24 11 T 1 TDEC END DP P PHFO 62 READ 1 1017 DXH PROGRAM PLOTEST 1017 FORMAT 13 WRITE 4 1018 Cas TRACE DES COURBES POUR INTERPRETATION DES ESSAIS 4018 FORMAT MODULES SUR Y zer Caaan UTILISE FICHIER ESSTRA S CREE PAR ESSAIS S READ 1 1012 NBHY C VERSIOH POUR TRACEUR 747
8. XLP XLP 1000 WRITE 1 1028 65 1 1022 PAP WRITE 1 2000 WRITE 1 1001 1 1003 LIBFICI 1 1 6 WRITE 1 2006 1 1022 RMI WRITE 1 2008 READ 1 1022 SI RUISRWINEXP 51 WRITE 1 2010 READ 1 1022 DITI DITI DITI 10 WRITE 1 2012 READ 1 1022 XLTI XLT1 1000 WRITE 1 0014 READ 1 1022 XLI XLI XL1 1000 WRITE 1 2036 READ 4 2022 IDE IDES 1 IF IDE NE ICOD GO TO 1 IDES 2 ENCODE 1060 ESC CALL ECRV24 IBUF 3 0 08 18 8 8 1 30 6 9 0 1 DO 5 1 1 30 J 9 Q J Q 9 IxDQ 1 CONTINUE CALL OPEN PROUT S D iesen WRITE 6 1044 WRITE 4 2024 1 100 01 0 100 DITP4 DITP 10 XLTP4 XLTP11000 XLP 1 XLP 1000 WRITE 6 2026 XKF H1 00 D1 RWP S PHI CT DITP ALTP4 PAP 1 IVERS Desses WRITE 6 2028 01 0171 10 XLTI1 XL11 1000 XL14 XL1 1000 51 WRITE 6 2030 XKI HI RWIA S1 DITIA XLTIA ALIA IPAG 0 IRC 1 1 26 3600 Los Co PUITS DE PRODUCTION 3 CONTINUE IRC IREH WRITE 1 1000 WRITE 1 2042 READ 1 1022 15 Lesen IF IDES EQ 1 GO TO 58 WRITE 1 1062 PAUSE ETAT CIBLE CALL ECRV24 CIBUF 9 0 98 CONTINUE TisT6x2 6098E 0 00 4 60 TO 30 32 34 IVERS 30 XK XKP H HP 60 TO 36 32 Xk XKI H HIx100 60 TO 36 XK
9. de une vitesse d impression nettement plus r duite et est cons quent r serv e l dition d finitive des documents apr s v rifications e OPERATOIRE Avec ce logiciel l interpr tation classique d un essai s effectue selon le mode s quentiel suivant saisie des donn es au clavier version actuelle l aide de l diteur de texte et cr ation d un fichier des donn es brutes pr traitement des donn es l aide du programme ESSAIS fournissant le tableau des donn es retenues annexe du rapport final et cr ant un fichier pour le trac trac des diff rentes courbes choisies annexes du rapport final l aide du programme PLOTEST interpr tation et d termination des param tres du r servoir d apr s les courbes et les m thodes classiques modifications et mise jour du rapport d essai type pr enregistr l aide du logiciel de traitement de texte WORDSTAR puis dition du document sur l imprimante Le document des pr visions d exploitation pressions et puissances lectriques est ensuite obtenu en deux tapes calcul des pr visions d exploitation selon diverses hypoth ses l aide du programme PREVIS avec visualisation sur l imprimante et cr ation d un fichier des r sultats mise jour du document dactylographi type pr enregistr et fusion avec le fichier des r sultats pr c dent avant dition Tous les documents obtenus peuvent
10. objectif tout aussi important en pratique 1 sortie quasiment imm diate du document de synth se dactylographi l aide du m me quipement informatique Ce document est obtenu par simples corrections d un rapport type standard pr enregistr traitement de textes Ce logiciel d aide l interpr tation des donn es d essais est exploit depuis plusieurs ann es d velopp tout d abord pour un micro ordinateur Zilog puis Micral R2E il a t ensuite am lior pour suivre l volution des besoins caract ristiques des nouveaux essais et profiter au mieux des performances des nou veaux p riph riques pr sents sur le march traitement de textes notamment L utilisation de ce logiciel relativement simple d emploi a permis un gain de temps appr ciable et une meilleure fiabilit pour l interpr tation de nombreux essais sur forages g othermiques 1 SYSTEME UTILISE Le syst me utilis est le MICRAL 21D quip des p riph riques suivants imprimante matricielle standard imprimante type traitement de texte avec accents et logiciel WORDSTAR traceur de courbes 4 HP 7470 L imprimante matricielle n est pas indispensable pour cette application elle permet toutefois de controler rap dement les r sultats d un calcul et de conserver une trace crite de l exploitation Le controle des r sultats peut galement tre effectu l cran L imprimante de type traitement de texte pos s
11. tation d essai obtenu l aide du logiciel IAO Annexe 3 em Se Sas suites gares 243 Fac simil d un rapport pr visions d ex ploitation obtenu l aide du logiciel Annexe 4 RA 199 Listings du logiciel utilis Remarques compl mentaires pour la modifi cation des programmes INTRODUCTION L interpr tation des essais hydrauliques effectu s sur les deux forages d un doublet g othermique constitue une tape im portante pour l estimation du potentiel local de la ressource productivit temp rature et la pr vision du comportement ult rieur du syst me d injection pompage durant l exploitation Les essais sont g n ralement effectu s imm diatement apr s les travaux de compl tion et de d veloppement de l ouvrage en utilisant le mat riel disponible sur le site cet instant Les mesures sont obtenues par enregistrement des variables pression et temp rature l aide d une sonde descendue sensiblement au toit de la formation r servoir sabot du casing Ces enregistrements sont r alis s par une soci t de service sp cia lis e avec un quipement identique celui utilis pour les autres diagraphies Ces donn es brutes temps pression temp ra ture sont disponibles pour interpr tation sur divers supports physiques listings cassettes disques Les test pratiqu s usuellement bien que sp cifiques cer tains gards sont assez voisins de ceux r alis s dans le domaine
12. 3 14 10 33 14 11 3 14 11 33 14 12 3 14 12 33 14 13 3 14 13 33 14 14 3 14 14 33 14 15 33 14 16 33 14 17 33 14 18 33 14 19 33 14 20 33 14 22 3 14 23 3 14 24 3 14 26 3 14 28 3 14 30 3 14 32 3 14 34 3 14 36 3 14 38 3 14 40 3 14 42 3 14 44 3 14 46 3 14 48 3 14 50 3 14 52 3 14 54 3 14 56 3 14 58 3 15 O 3 15 2 3 15 4 3 15 8 3 15 12 3 14944 15778 16611 17445 18278 19111 19944 20778 21611 22444 23278 24111 24944 26611 28278 29945 31611 33278 34944 37444 39111 40778 44111 47444 50778 54111 57444 60778 64111 67445 70778 74111 77444 80778 84111 87444 90778 94111 97444 1 00778 1 04111 1 07444 1 14111 1 20778 5 295 5 327 5 359 5 386 5 414 5 440 5 465 5 489 5 511 5 535 5 556 5 577 5 595 5 634 5 667 5 699 5 732 5 762 5 788 5 828 5 852 5 876 5 920 5 962 5 998 6 034 6 068 6 099 6 129 6 156 6 183 6 207 6 231 6 255 6 276 6 297 6 318 6 337 6 356 6 374 6 393 6 409 6 442 6 473 46 094 43 436 41 074 38 961 37 060 35 340 33 047 31 682 36 dt heures dp kg cm2 tp dt dt ou fonction S 2 27444 1 34111 1 40778 1 47444 1 54111 1 60778 1 67445 1 74111 1 80778 1 87444 1 94111 2 00778 2 10778 2 20778 2 30778 2 40778 2 50778 2 60778 2 70778 2 80778 2 90778 3 00778 3 14111 3 27444 3 40778 3 54111 3 67444 3 80778 3
13. Exp 1 lt 3 ATkh 4 k t 41 En pratique peut crire Q u B u Cr El Gene 4rrkh 4 k t La pression P2 n est donc fonction que de la viscosit du fluide initial Fonctionnement en doublet a a a ae ae La variation de pression totale chaque puits s obtient en sommant la contribution propre de l ouvrage consid r fonction du rayon efficace du puits et la contribution de l interf rence due l autre puits situ la distance D Pour le puits de production Q u D Pe CS 2nkh r Ainsi lorsque les param tres du r servoir sont homog nes et pour les temps longs le rabattement net la production est stationnaire 13 Pour le puits d injection L existence d un contraste de viscosit se traduit donc dans ce cas par une variation de pression nette fonction du temps Pertes de charges La perte de charge par kilom tre de tubage de section constante est calcul e par la relation 2 7 1 057 10 35 avec u viscosit en cp Q d bit en m3 h Di diam tre interne en mm Puissance lectrique La puissance lectrique est calcul e par la relation P kW O 04 Q Pd avec Q d bit en m h Pd pression diff rentielle de la pompe kg cm 14 8 3 SEQUENCE D EXPLOITATION EXEMPLE Activation PREVISS Exploitation SELECTION DU MODE DE CALCUL DE L INTERFERENCE OPTION CHOISIE 2 1 DONNEES POUR
14. PRODUCTION SIGLE DU FORAGE 6A1 FORS XKP D F10 3 0 805 HP M F10 3 23 00 10 3 0 5 D M F10 3 1200 RW CM F10 3 7 62 S F10 3 3 71 PHI DEC F10 3 0 15 10 3 0 100E 03 GAMME DE DEBIT 50 100 125 150 175 200 225 250 CONFIRMATION O N O TUBAGE DITP CM F10 3 15 94 TUBAGE XLTP M F10 3 1599 LONG PRODUCTION XLP M F10 3 122 PAP KG CM2 F10 3 8 5 DONNEES INJECTION SIGLE DU FORAGE 6A1 FOR6 RWI CM F10 3 7 62 S F10 3 3 7 TUBAGE DITI CM F10 3 15 94 TUBAGE XLTI M F10 3 1813 LONG INJECTION XLI M F10 3 150 DONNEES POUR PRODUCTION DATE EN MOIS F10 3 120 Sortie r sultats 15 Pause changer la feuille DONNEES INJECTION NB TEMP INJECTION MAX 9 11 4 pen LEURS CROISSANTE D O 35 40 45 ESTE VISCOSITES 9F5 2 SS 0 83 0 74 0 67 0 61 LISTE PRES ART POT 9F5 1 7 0 7 0 7 0 7 0 g DATE EN 15 10 3 2 1 Sortie r sultats Pause changer la feuille AUTRE TEMPS O N DATE EN MOIS F10 3 7 Sortie r sultats Pause changer la feuille AUTRE TEMPS O N DATE EN MOIS 120 Sortie r sultats AUTRE TEMPS O N H STOP 9 MISE A JOUR ET EDITION DU DOCUMENT DES PREVISIONS La mise jour du document s effectue de la m me mani re que pour le rapport d essai on utilise un document type pr enregis tr Le fichier des r sultats PROUT S obtenu pr cedemment est inclu directement par fusion des deux fic
15. READ 6 END 200 DT I THCI DPCI PCI 1 1 7 DO 40 K 1 7 IF DT K E4 0 60 TO 200 50 TO 50 52 54 210 58 50 1000 Cases C LOG LOG ET SPE 30 CONTINUE X DT K UMIN IF X LT 1 60 TO 54 0610 X xFLOAT NINC IX XAD INTOO IF IX BT XAF 60 TO S UMINY IF Y LT 1 TO 56 Y DYHXALOG10 V xFLOAT NINC IY YAD INT Y IF 1Y GT YAF 60 TO 56 ENCODE IBUF 444 IX I 114 FORMAT PU PA I6 16 PD PU CALL ECRV24 IBUF 24 0 60 TO 56 C C aas HORNER DEBIT CONSTANT C 52 X TH K IF X GT XHLIM 60 TO D DXM AL0610 X 40 IX XAD INT X Y P K UNINY IF CY LT 0 60 TO 56 Y DYM Y 40 DY IY YAD INT Y IF IY GT YAF 60 10 S ENCODE 114 CALL ECRV24 04 0 60 TO 56 Coros SEMI LOG C 94 X DT K UMINX 64 IF X LT 1 60 TO S X DXMXALOG10 X 4D IX X D INT X IF IX GT XAF 60 TO 56 Y P K UHENY IF Y LT 0 60 TO 56 Y DYHxY 40 DY IY YAD INT Y IF IY BT YAF 60 TO 56 ENCODE IBUF 444 IX IY CALL ECRV24 24 0 60 TO 56 Comos C HORNER DEBIT VARIABLE 58 X TH K UMINX IF X LT 0 60 TO 5 X kxDXN440 DX IX XAD INT X IF CIX GT XAF 60 TO 56 Y P K UMINY IF Y LT 0 60 TO 56 Y DYHXY 40 DY IY YAD INT Y IF IY GT YAF 60 TO 56 ENCODE IBUF 144 IX IY CALL ECRV24 CIBUF 24 0 56 CONTINUE 40 CONTINUE 60 TO 38 200 CONTINUE ENDFILE 6 60 TO 40 Css C FIN DE PROCESSUS 240 C
16. XAD 250 DXG 8 IX XAF XAD LX 9 1 00 CHOISI 77 READ 1 1012 1000 DYH 5 40 4012 FORMAT 11 DYB 7200 LY DYH IF 1 00 6 4 60 TO 210 YAD 2794DYE IF ICOD NE 6 60 TO 70 YAF YAD LY INB 10 Co GRILLE FIXE ABAQUE SPE Cicle NINC 40 NINC 1 NBHX 3 C INITIALISATION TRACEUR DXH 3060 LX 9180 74 CONTINUE XAD 810 WRITE 1 1030 XAF 9990 4030 FORMAT 1X CONTROLE ETAT TRACEUR RET OU T ET RETURN NBHY 2 PAUSE ETAT DYH 3060 ENCODE 100 LY 6120 400 FORMAT INSLT SP1 YAD 580 CALL ECRV24 IBUF 9 0 YAF 6700 Disisi INB 10 Co TRACE GRILLE 60 TO 72 Cessa AXE X LOG Ca Css GRILLE NEL HBHX 4 gei IF IC0D E8 5 GO TO 60 70 CONTINUE XX ALOG1D UMINX WRITE 1 1014 IF XX 6E 0 XX XX D 5 4014 FORMAT DONNEES GRILLE IF XX LT D XX XX D 5 41X NB MODULES SUR X 11 2 7 INE INT XX 1 READ 1 1012 NENT DO 20 I 1 NBL WRITE 1 1016 IC XAD 1 1 4DXMENINC 4016 FORMAT 1X 7LONG MODULE X EN 13 7 7 INE INE 63 ICX4 IC 420 CALL ECRV 4 IBUF 40 0 ICY1 YAD 240 IF I EG NBL 60 TO 23 2 1 1 240 00 22 4 2 10 ICY2 10Y4 420 XX FLOAT J ENCODE IBUF 102 ICX1 ICY1 INE ETX A DYHXALOG10 XX 4FLOAT NINC 102 FORKAT DI1 0 510 2 0 3 PU PA 16 16 LB 12 41 IC4 IC INT X ENCODE IBUF 104 XAD IC1 X4F IC1 CALL ECRV24 CIBUF A9 DI CALL ECRV24 CIBUF 40 0 ENCODE CIBUF 103 ICX2 1CY2 I
17. XKP XKI 2 H HP H14100 2 36 YI UD FACA 4 XK HXP TI XI PHI UDACT 4 XK PC4 1 057E 7 U0 D 21 XTI XIxDxD T4 DO 7 L 1 IDES NBQC NBQ L DO 4 1 1 NBGL J 1 L 1 8 XHP YP Q J XHI YIxG J 1 1 1 1 1 0 9678 1 1 1 79 01 4 79 PC I PU I wXLTP DD 20 RWP 9 0 2 1 1 41 1 79 79 PD I PH I PC I PAP 1 0 1 50 04 4 CONTINUE IF L E 2 60 TO IF 10 5 6 4 PN D DO 52 LL 1 3 9 50 TO 5 IF PD B LE PHAX LL 60 TO 54 52 CONTINUE PA PHAX LL CALL CADRE 011 0 8 PN PZ 1 TS LIBFP 1 56 CONTINUE IPAG IPAG 1 4 DO K 1 2 IF IRC EG 0 WRITE J 1048 WRITE J 1046 WRITE J 2040 WRITE J 1030 WRITE J 1032 WRITE 7 1030 WRITE J 1034 WRITE 3 1036 WRITE 3 1038 WRITE 3 1030 WRITE 3 1040 WRITE 3 1030 WRITE 3 1042 WRITE J 1030 6 CONTINUE 60 TO 7 IPAG TS 1 1 1 8 124 6 1 1 1 8 FAP PAP PAP PAP PAP PAP PAP PAP PD 1 1 1 8 PE 1 1 1 8 9 CONTINUE CALL COURBE PD 4 PE 1 0 4 0 7 CONTINUE WRITE 1 2044 READ 1 2022 IDE IF IDE EQ ICOD 60 TO 3 WRITE 6 1044 C PUITS D IHJECTION C IPAG 0 WRITE 1 2000 WRITE 1 1066 READ 1 1052 NBT WRITE 1 1065 READ 1 1070 TTICI I 1 HET WRITE 1 1072 READ 1 1074
18. calcul est effectu pour t 1 mois 7 mois 120 mois avec prise en compte de l effet pari tal thermique influence de la viscosit plus elev e dans la zone envahie par les eaux de r injection 46 11 2 PERTES DE CHARGE DANS LES TUBAGES dPc 1 057 10 kg cm2 di diam tre int rieur du casing Par convention le calcul des pertes de charge dans la compl tion est effectu pour un d bit moyen gal la moiti du d bit nominal 11 3 BILAN DES PRESSIONS lt Les puits tant art siens les bilans des pressions sont les suivants Production Pd dPh dPc Pap Injection Pd dPh dPc Pas 11 4 PUISSACE ELECTRIQUE DE POMPAGE La puissance l ctrique de pompage est calcul e par relation 2 4 10 kW 47 DONNEES PUITS DE PRODUCTION PERMEABILITE INTRINSEQUE 805 HAUTEUR PRODUCTIVE CUMULEE 23 000 VISCOSITE DU FLUIDE RESERVOIR 500 DISTANCE ENTRE PUITS 1200 000 RAYON PUITS DE PRODUCTION 7 620 SKIN DE PRODUCTION 3 710 POROSITE casco se meses sis 0e 58 0 0 6 150 COMPRESSIBILITE 100 DIAMETRE INTERIEUR 15 940 LONGUEUR TUBAGE DE PRODUCTION 1599 000 LONGUEUR PRODUCTION RESERVO
19. et cr e deux fichiers ESSOUT S fichier des r sultats qui peut tre visualis l cran ou sur l imprimante et inclus dans le rapport final par fusion avec le fichier texte voir paragraphe 7 fichier est formatt et son examen fournit les valeurs limites n cessaires la d finition des trac s des courbes ESSTRA S fichier n cessaire au programme de trac formatt Ces deux fichiers sont d truits chaque nouvelle exploita tion du programme ESSAIS pour les conserver il convient donc de les renommer Le traitement de la variable temp rature n cessite une exploitation suppl mentaire dans ce cas consid rer l analogie pression temp rature pour la r ponse aux questions 4 1 EXEMPLE 1 On consid re l enregistrement d une remont e de pression fichier XDATBUC S existant cons cutive 10 heures de production d bit constant 100 m3 h partir t 00 0000 Les questions sont soulign es dans la s quence d exploitation suivante Activation du programme ESSAIS NOM DU FICHIER DES DONNEES XDATBUC S T NI OPTION TYPE D INTERPRETATION RABATTEMENT R INJECTION BUILD UP B FALL OFF F CODE B DUREE DE 13 212 2 10 0000 DATE FIN DE PRODUCTION 13 212 10 0000 PRESSION PWFO F8 3 175 292 4 2 EXEMPLE 2 On consid re galement la remont e de pression suivant un
20. rieure la vitesse de Darcy Vth Vd Ca o Cf et Ca sont les capacit s calorifiques volumiques res pectives du fluide et du milieu poreux satur Le voisinage du puits est alors constitu de deux zones con centriques O lt lt Rth z6ne envahie par les eaux d injection tem p rature constante Ti avec une viscosit ui d o la solution Pl r gt Rth z ne situ e l ext rieur du volume envahi o T u uo d o la solution P2 Pression l int rieur de la z ne envahie P1 en zm pe pn e FE e E e E e EE E E E E E E E rr E e E E mn E E e EE bn E E E E E a Avec les hypoth ses ci dessus on montre que la pression Pl s exprime l aide de la fonction exponentielle int grale avec la viscosit d injection ui et d une fonction S traduisant l effet pari tal thermique Q u 8 u C r dP P1 El S 4nkh 4 k t avec A u u A A tes A Ca E u ATI et g u C Cf Q ss k Ca h 12 En pratique et pour les forages g othermiques S peut tre remplac e par la relation approch e plus simple u ga F C Dn A SS u On remarquera que S est ind pendante du temps et ne d pend que du d bit et du contraste des viscosit s Pression l ext rieur de la z ne envahie P2 Dans ce cas la pression P2 s crit B u Eeer A i dB P Esos EL eee
21. une premi re version susceptible de modifications en fonction de l volution des m thodes d analyse et des souhaits des utilisateurs L utilisateur est cens conna tre le fonctionnement g n ral du syst me utilis micro ordinateur MICRAL seuls certains points n cessaires de d tail sont rappel s 10 SOMMAIRE Introduction Syst me informatique utilis Mode op ratoire ses eee se Saisie 24 enregistrement des donn es Pr traitement programme ESSAIS 4 1 Exemple 1 4 2 Exemple 2 S E Sahaya kapaka am Trac des courbes programme PLOTEST se D termination des param tres du r servoir Mise jour et dition du rapport d essai Calcul des pr visions d exploitation programme PREVISS st ra 8 1 Hypoth ses du calcul 8 2 Rappels des formules utilis es 8 3 SE ones d exploitation exemple Mise jour du document des pr visions d exploitation Annexes qu drames esse ss EE 10 14 15 16 10 1 10 2 10 3 10 4 Diagrammes type disponibles par utilisation du programme PLOTEST Exemple de fiche caract ristiques pour interpr tation DU ere aer Fac simil d un rapport interpr
22. 0 4 23 4 26 5 29 5 kg cm2 6 2 0 3 0 4 2 5 5 7 0 8 7 10 5 kg cm2 dPap 7 0 7 0 7 0 7 0 7 0 7 0 7 0 7 0 kg cm2 ka O 3 s s E mme mm ge cm lt gt a em y s 53 111 2 10 Temp rature d injection 35 00 Viscosit 74 Estimations apr s 120 00 mois d exploitation d bit constant Q 50 100 125 150 175 200 225 250 m3 h 111 2 11 Temp rature d injection 40 00 Viscosit 67 Estimations apr s 120 00 mois d exploitation d bit constant 54 111 2 12 Temp rature d injection 45 00 Viscosit 61 cp Estimations apr s 120 00 mois d exploitation d bit constant 55 Ar ne ti n in D par 3T1311N3933310 6 q AZ A A E E 100 125 150 175 200 225 230 DEBIT EN M3 H 175 200 225 250 75 o 8 8 8 1410413 3FONVSSINd DEBIT M3 H PREVISION D EXPLOITATION t mois 120 0 PRODUCTION FORAGE FORS 56 250 Q r CO EDS S sss Q NN E Ki gt 9 3 5 225 y 225 200 173 173 130 DEBIT EN M
23. 0 ET HICRAL 210 WRITE 1 1020 C INTERFACAGE KODE HARDWARE HANDSHAKING 4020 FORMAT 1X LONG MODULE Y EN amp 12 23 Caoa NECESSITE ROUTINE ECRV24 AVEC TEST PIN 20 READ 1 1017 DYH VERSION 1 0 DU 30 01 83 rs AUTEUR A MENJOZ TRACE DES POINTS Cire Cosa INTEGER DXG DXD DYE DYM DXN DYH YAD 222272 CONTINUE INTEGER 1 LIB 80 ITERM IBLAH ET ES IBUF 80 WRITE 1 1022 DIMENSION DT 8 TH B DP 8 P 8 4022 FORMAT 1X TRACE DES POINTS DATA ESC ETX ITERM IBLAN 27 3 4H 11X VALEUR MINI SUR X E10 3 2 7 Cost READ 1 1024 UMINX SET HARDWARE HANDSHAKING 4024 FORMAT 10 3 WRITE 1 1026 ENCODE IBUF 900 ESC 4026 1X VALEUK MINI SUR Y E10 3 2 900 FORMAT 1 READ 1 1024 CALL ECRV24 IBUF 3 D IF 1 00 6 6 GO TO 74 JF ICOD NE 1 WRITE 1 1028 Cases SELECTION DU TRACE 1028 FORNA 1X VALEUR D UN MODULE Y E10 3 2 Cua IF ICOD ME 1 READ 1 1024 DY 10 CONTINUE IF ICOD EG 5 WRITE 1 1032 WRITE 1 1010 4032 FORMAT 1 VALEUR D UN MODULE X E10 3 2 7 4010 FORMAT 1X OPTION TRACE CODE IF 1 00 6 5 READ 1 1024 DX 4 Xe Elise 2 1X 06 106 14 DE LA GRILLE 3 4X HORNER Q CONST 2 Cie 4 4X SENI LOG NDH 3 LX NBNX DXH 40 5 4X FIN 4 DXD 15 40 1X HORNER Q VAR 5 DX6 10000 LX DXD 7 1X ABAQUE 4
24. 000 8 6 GO TO 78 30 J J 4 IF CICOD HE 1 60 TO 42 IF LIB J EG IBLAN GO TO 30 IF J LE 50 60 TO 24 106 105 ET SPE WRITE 1 1004 Losso 1004 FORMAT 4X LIBELLE TROP 10801 78 NBL NBMY 1 60 TO 28 XX ALOG1D UMINY gt 24 X FLOAT LY A2U FLOAT J 2 IF 0 XX XX 0 5 ILY YAD INT X IF XX LT 0 XX XX 0 5 ENCODE CIBUF 140 IX IB Fei BD INE INT XX i 410 FORMAT PUsPA 16 16 LB 5141 PU 00 22 1 1 NBL CALL ECRV24 CIBUF 75 0 IC Y4D 1 4 INE INEH Ces LIBELLES AXE X ICX1 XAD 440 Diosas ICY1 IC 60 INC 0 ICX2 10X1 240 DO 34 1 3 ICY2 10Y1 120 32 WRITE 1 1008 ENCODE CIBUF 102 ICX1 ICY1 INE ETX 4008 FORMAT 1X LIBELLE 12 AXE X 70 CAR MAX 7 CALL ECRV24 IBUF 43 0 11 70 1 ENCODE IBUF 103 ICX2 ICY2 INE ETX READ 4 1002 LIB I 1 1 71 CALL ECRV24 IBUF 44 0 ENCODE IBUF 109 ENCODE IBUF 104 XAD IC XAF IC 109 FORMAT 7014 0 510 2 0 27 CALL 15 0 ILY 279 520 INC J 71 35 Ei IF LIB J EG IBLAN 60 TO 35 IF J LE 70 60 TO 33 WRITE 1 1004 60 TO 32 33 X FLOATILX 120 FLOAT J 2 ILX XAD INT X ENCODE IBUF 112 ILX ILY 112 FORMAT 16 16 CALL ECRU24 IBUF 18 0 ENCODE IBUF 113 LIB K K 1 70j ETX 113 FORMAT LB 7461 CALL ECRV94 IBUF 73 0 INC 1NC 160 34 CONTINUE Co TRACE DES POINTS Coros 82 CONTINUE CALL OPEN 6 ESSTRA S D REWIND 6 38 CONTINUE
25. 2 PD 0 77 dt 1h tD CD 1000 CD exp 2S 10 Transmissivit 2 rr kh dP kh PD Q PD sels d o Q u u 2 kh 4045 D cm cp u avec u 0 50 cp kh 20 22 D m et h 23 m k 0 879 D Facteur de skin kh dt 5 2 0 915 10 cm3 atm C CD en B Ct h rw et avec S 29 2 7269 0 5 Ln CD exp 2S CD il vient 5 3 29 b M thode de HORNER Le report de la figure 2 repr sente l volution de la pression de fermeture Pws en fonction de la variable tp dt at o tp est la dur e de production et dt le temps coul depuis l arr t du pompage Transmissivit La pente moyenne de la m kh d o u kh et avec h 23 m k Facteur de skin droite vaut 1 632 kg cm2 d cade 0 183 Q 3702 D cm cp m 18 51 D m 0 805 D Pression extrapol e pour le calcul 1 157 533 kg cm2 30 1 15 ees ONT m u Ct rw soit 5 3 71 Cette derni re valeur de S traduit l effet de skin global de l ouvrage incluant l incidence de la d viation 29 degr s Pression statique extrapol e c M thode de MILLER DYES HUTCHINSON Pour les faibles valeurs de la dur e de fermeture dt 1 pression Pws est report e en fonction de Log dt la figure 3 Transmissivit Pente moyenne m 1 54 kg cm2 d cade kh d o 0 183 3674 D cm cp u m kh 18 37 D m k 0 798 D avec h 23 m Fa
26. 3 H PREVISION D EXPLOITATION t mois 150 DEBIT EN M3 H d d n 8 NZ a 5 10000 HITET ETET 2 Wz nuna Ho 6 441 si si Ea Lat r r ERR 29 yasa Je q R R 8 g B Ro 8 37T1311N3434110 NOISS3Ud 13 3INVSSINd 1 0 FORAGE FORG INJECTION 57 125 150 175 200 225 DEBIT EN M3 H 100 D INJECTION 30 35 40 45 TEMPERATURES 3TIZLINIHIASIO NOISS3d 500 250 225 175 125 o mz a ARES Zona A A 8 8 q va 19313 39NVSSInd DEBIT EN M3 H PREVISION D EXPLOITATION A 7 0 INJECTION FORAGE FORS PRESSION DIFFERENTIELLE Pd Kg om2 PUISSANCE ELECTRIQUE 58 TEMPERATURES D INJECTION 30 C 35 40 45 C 4 H 8 75 100 125 150 175 200 225 250 DEBIT EN M3 H TEMPERATURES D INJECTION am Ti 35 Ti 40 C 45 d PA B 100 1 CAR ARA a AA l PR RAI 055 75 100 125 150 175 200 225 250 DEBIT EN M3 H INJECTION PREVISION D EXPLOI
27. 44 7 558 1 994 158 862 AAA O CR A do soga PERS SSD AS E nn A RE A AE A A AN O AAN El E ER AS DE AECA de o As RS Wee eE ME LOG LOG FOR GRAM F A EURES TI DI H E TU FIC RE FIGURE 1 DT EN TEST DE FORMATION SIT REMONTEE DE PRESSION APRES FERME CNI IA 1553 30 NOILYIXYA 159 158 157 156 155 154 CNI IA 55394 TIF F DIAGRAMME DE HORNER DT CTIF IA REMONTEE DE PR PRESSION EN KG CM2 159 158 10 FIGURE TP DT TEST DE FORMATION SITE FIC REMONTEE DE PRESSION APRES FERMETURE D T I I A D F G FORS RAMME DE HORNER PRESSION EN KG CM2 160 158 156 154 152 150 10 TES REMONTEE DE PRE SSI Se EN HEURES APRES FERMETURE DIAGR RE MATION SITE FICTIF ME SEMI LOG 43 10 3 ANNEXE 3 Fac simil d un rapport pr visions d exploitation r alis int gralement avec le syst me informatique MICRAL 44 B R G M Orl ans le 31 Janvier 1983 DEPARTEMENT GEOTHERMIE r f 72 83 AM am SITE FICTIF FOR FORAGES FOR5 FOR6 PREVISION DES PRESSIONS DE POMPAGE ET DES PUISSANCES CONSOMMEES Les pressions de pompa
28. 45 23 16 3 146 23 32 3 147 23 40 3 148 23 48 3 149 23 56 3 150 0 4 3 151 0 12 3 152 0 20 3 153 0 28 3 154 0 36 3 155 0 44 3 156 0 52 3 157 1 0 3 158 1 8 3 159 1 16 3 160 1 24 3 161 1 32 3 162 1 40 3 163 1 48 3 164 1 56 3 6 07445 6 20778 6 34111 6 47445 6 60778 6 74111 6 87444 7 00778 7 14111 7 27444 7 40778 7 54111 7 67444 7 80778 7 94111 8 07445 8 22445 8 34111 8 47445 8 60778 8 74111 8 87444 9 00778 9 14111 9 27444 9 54111 9 67444 9 80778 9 94111 10 07444 10 20778 10 34111 10 47445 10 60778 10 74111 10 87444 11 00778 11 14111 11 27444 11 40778 11 54111 11 67444 11 80778 11 94111 7 275 7 285 7 295 7 304 7 315 7 323 7 331 7 340 7 348 7 355 7 364 7 371 7 380 7 387 7 393 7 398 7 406 7 414 7 420 7 427 7 435 7 440 7 445 7 451 7 456 7 467 7 474 7 479 7 485 7 489 7 496 7 501 7 504 7 511 7 515 7 520 7 524 7 529 7 534 7 538 7 541 7 547 7 549 7 556 158 579 158 589 158 599 158 608 158 619 158 627 158 635 158 644 158 652 158 659 158 668 158 675 158 684 158 691 158 697 158 702 158 710 158 718 158 724 158 731 158 739 158 744 158 749 158 755 158 760 158 771 158 778 158 783 158 789 158 793 158 800 158 805 158 808 158 815 158 819 158 824 158 828 158 833 158 838 158 842 158 845 158 851 158 853 158 860 38 no du t dt dp tp dt dt p ou point heures heures kg cm2 fonction 5 kg cm2 165 2 4 3 12 074
29. 6 ENCODE CIBUF 1000 IX IY 1000 FORMAT PU PA I6 16 CALL ECRV24 18UF 21 0 60 TO 18 46 ENCODE IBUF 1002 D I 1002 FORMAT 16 16 CALL ECRV24 IBUF 45 0 18 IF 1P EQ 27 IXD IX IF IP E0 31 IXF IX 80 CONTINUE IF 0 8 1 60 TO 98 RUN IXF 1XD RIS IYF IYD 22 RUN RUN RISARIS 22 6 RT 22 RUN RUN ZZ RIS RI5 Z2 ENCODE IBUF 1004 RUN RIS 4004 FORMAT DP F5 2 FS 2 810 45 0 2 CALL ECRV24 24 0 IX 1XD 40 IF CIFL EG 1 IX IXD 200 ENCODE IBUF 1006 IX IYD XLIB TE ET 1006 FORMAT 14 16 LB 44 F3 0 41 CALL ECRV24 IBUF 29 0 98 CONTINUE TEE IF 0 60 TO 99 ENCODE 1012 4012 FORMAT 7010 47 CALL ECRV94 IBUF 5 0 1X IXNC 640 ICOUR 160 IY IYN 490 ENCODE IBUF 1014 IX IY XLIB TE IUN ETX 1014 FORMAT PU PA 16 16 LB A4 F4 0 241 CALL ECRV24 IBUF 31 0 99 CONTINUE Conos ENCODE 1010 CALL ECRV24 2 0 400 CONTINUE 4010 FORMAT IW RETURN END 71 Remarques compl mentaires pour la modification des programmes Les programmes sont disponibles sous la forme d objets ex cutables Dans le cas de modifications des logiciels ceux i doivent tre compil s puis link s avec les routines syst me n cessaires biblioth que FORTRAN et module d interfa age du traceur a Cas du programme ESSAIS b F80 ESSAIS S 180 ESSAIS FORLIB S ESSAIS N E b Cas du pr
30. 6 2X POROSITE F10 3 1002 FORMAT E10 3 7 COMPRESSIBILITE 210 3 atg 1 4003 FORMAT 641 8 2X DIAMETRE INTERIEUR 10 3 1004 FORMAT 2X XKP D F1D 3 9 2X LONGUEUR TUBAGE DE PRODUCTION F10 3 4005 FORMAT 2X XKI D F10 3 1 LONGUEUR PRODUCTION RESERVOIR F10 3 1006 FORMAT 2X HP 4 F10 3 2X PRESSION ARTESIENNE POTENTIELLE F4D 3 ka cu2 1007 FORMAT 2X HI A F10 3 3 2 OPTION CALCUL INTERFERENCE 110 1008 FORMAT 2X U0 CP F40 3 2028 FORHAT 2X DONNEES PUITS D INJECTION 2X 25 4010 FORMAT 2X D K F10 3 2030 FORMAT 2X PERMEABILITE INTRINSEQUE F10 3 0 7 1012 FORMAT 2X RW CH F40 3 7 2X HAUTEUR PRODUCTIVE CUMULEE 2 F10 3 RI 4014 FORMAT 2X S F40 3 1 RAYON PUITS 0 2 10 3 68 SKIN IHJECTIOR 10 37 2X DIAMETRE INTERIEUR TUBAGE F40 3 LONGUEUR 2 5 n 7 2 LONGUEUR INJECTION RESERVOIR F10 3 2032 FORMAT 4X III 2 PUITS D INJECTION 9X 47 2036 FORMAT 2 TRACE DES COURBES 0 N 2038 FORHAT IP 2040 FORMAT 6X Estimations apres F 2 mois d exploitation 1 debit Constant 2042 FORMAT EN MOIS F10 3
31. 9 ETA 4006 FORMAT LB 984 81 CALL ECRV24 39 0 50 TO A IX IXXC 1D3 ENCODE IBUF 1006 IA 1112 CALL ECRV24 IB 198 0 ENCODE 1010 ZL3 1 1 1 7 ETX 1010 FORMAT LB 7A4 41 CALL ECRV24 31 0 6 CONTINUE ENCODE CIBUF 1012 4012 FORMAT 010 1 510 15 0 2 CALL ECRV24 IBUF 16 0 _ _ Je 1200 135 DO 10 18 1 2 ITX 80 IF 14 E0 2 ITX 60 IX IXRC IF 18 68 2 DO 8 17 1 NTIK Y IYHY IT 4 xKIDY IF CIY GT IYX IY IY ENCODE IBUF 4014 IX IY ITX ITY 1014 FORMAT PUSsPA 16 16 PD PR 16 16 CALL ECRV24 40 0 IF 18 6 1 60 TO 8 16 10N 1T 1 1D6 IF 16 6T 10X 18 16 ENCODE IBUF 1016 JX JY IQ ETX 4016 FORMAT PR 16 16 LB 13 41 CALL ECRV24 IBUF 22 0 8 CONTINUE 40 CONTINUE C IX IKXC 1D1 ENCODE IBUF 1006 IX IYL1 CALL ECRV24 CIBUF 18 0 ENCODE IBUF 4048 1 1 4 1018 FORMAT LB 4A4 A1 CALL ECRV24 19 0 IF ICOD NE 2 60 TO 11 IX IXNC 400 IY IYN 400 ENCODE IBUF 1006 IY CALL ECRV24 IBUF 48 0 ENCODE IBUF 1024 ZLB I I 1 3 ETX 1024 FORMAT LB 304 01 CALL ECRV24 15 0 IX 1X 160 ENCODE 1006 IX IY CALL ECRV94 CIBUF 18 0 ENCODE IBUF 1024 ZL9 1 1 1 3 ETX CALL ECRV24 IBUF 15 0 IX 1X 80 IYY 1Y 1080 ENCODE IBUF 1026 IX IY IX IYY 1026 FORMAT PU PA 16 16 PD PA 16 16 PU CALL ECRV24 IBUF 40 0 C
32. 9 51 14 0 6 14 0 21 14 0 36 14 0 51 14 1 6 14 1 21 14 1 36 14 1 51 14 2 18 14 2 33 14 2 48 14 3 3 14 3 18 14 3 33 14 3 48 14 4 3 14 4 18 14 4 33 14 4 48 14 5 3 14 5 18 14 5 33 14 5 48 14 6 3 14 6 18 14 6 33 14 6 48 14 7 3 14 7 33 00028 00445 00861 01278 01694 02111 02528 02944 03361 03778 04528 04945 05361 05778 06194 06611 07028 07444 07861 08278 08694 09111 09528 09944 10361 10778 11195 11611 12028 12444 13278 14111 012 081 667 2 961 3 352 3 739 3 969 4 106 4 240 4 339 4 486 4 549 4 613 4 669 4 718 4 762 4 808 4 846 4 883 4 920 4 951 4 979 5 010 5 039 5 066 5 091 5 114 5 140 5 160 5 182 5 223 5 260 43093 164 2700 616 1394 456 940 093 709 178 569 411 475 719 408 544 358 023 318 646 266 026 243 693 224 832 208 690 194 720 182 511 171 750 162 193 153 650 145 966 139 019 132 707 126 948 121 670 116 818 112 339 108 195 104 349 100 769 97 428 91 376 86 039 151 316 151 385 151 971 154 265 154 656 155 043 155 273 155 410 155 544 155 643 155 790 155 853 155 917 155 973 156 022 156 066 156 112 156 150 156 187 156 224 156 255 156 283 156 314 156 343 156 370 156 395 156 418 156 444 156 464 156 486 156 527 156 564 35 dt heures dp kg cm2 tp dt dt ou fonction S t heures 14 8 33 14 9 3 14 9 33 14 10
33. 94111 4 07445 4 20778 4 34111 4 47445 4 60778 4 74111 4 87444 5 00778 5 14111 5 27444 5 40778 5 54111 5 67444 5 80778 5 94111 6 501 6 528 6 554 6 580 6 602 6 624 6 645 6 666 6 685 6 704 6 724 6 739 6 765 6 788 6 811 6 831 6 853 6 873 6 892 6 910 6 927 6 944 6 966 6 987 7 006 7 024 7 042 7 060 7 076 7 092 7 107 7 122 7 135 7 150 7 163 7 175 7 189 7 201 7 213 7 223 7 235 7 246 7 256 7 267 6 276 9 948 9 524 9 139 8 787 8 464 8 167 7 892 7 638 7 402 7 182 6 977 6 693 6 435 6 200 5 984 157 805 157 832 157 858 157 884 157 906 157 928 157 949 157 970 157 989 158 008 158 028 158 043 158 069 158 092 158 115 158 135 158 157 158 177 158 196 158 214 158 231 158 248 158 270 158 291 158 310 158 328 158 346 158 364 158 380 158 396 158 411 158 426 158 439 158 454 158 467 158 479 158 493 158 505 158 517 158 527 158 539 158 550 158 560 158 571 mm mm o E EE mm mm mm mme 37 dt heures dp kg cm2 tp dt dt ou fonction S no du t point heures 121 20 4 3 122 20 12 3 123 20 20 3 124 20 28 3 125 20 36 3 126 20 44 3 127 20 52 3 128 21 O 3 129 21 8 3 130 21 16 3 131 21 24 3 132 21 32 3 133 21 40 3 134 21 48 3 135 21 56 3 136 22 4 3 137 22 13 3 138 22 20 3 139 22 28 3 140 22 36 3 141 22 44 3 142 22 52 3 143 23 O 3 144 23 8 3 1
34. AT 1X 3041 5 D INJECTION DU BA1 1086 FORMAT 1 3041 14X FALL OFF TEST DU D i 1088 FORMAT 1X DUREE D INJECTION FB 3 Heures 1 1 F8 3 KG CM 1090 FORMAT 2X OPTION TRAITEMENT TEMPERATURE O N 2 CALL OPEN 6 ESSTRA S Di CALL OPEN 7 55007 5 0 5 DE L ESSAI C NPHAX 300 WRITE 1 1036 READ 1 1038 B 1 1 1 11 ENCODE 4 1038 CALL OPEN 9 4 0 WRITE 1 1090 READ 1 1042 FLAG WRITE 1 1060 READ 1 1042 IF IOPT EG LR 60 TO 60 IF 10 8 11 60 TO 60 C BUILD UP ET FALL OFF IF 10 8 18 WRITE 4 4008 IF IOPT EG LF WRITE 1 1009 READ 1 1002 1H1 1M1 151 DPR FLOAT 1H1 FLOAT 1M1 60 FLOAT 1S1 3600 IF IOPT EQ LB WRITE 1 1000 IF CIOPT EG LF WRITE 1 1001 READ 1 1002 1 2 1 2 152 TPSFLOAT IH2 FLOAT IH2 60 FLOAT 152 3600 1 WRITE 1 1004 READ 1 1005 PNFO WRITE 1 1048 READ 1 1042 IP IF IR EQ IRR 60 TO 40 WRITE 1 1050 READ 1 1002 IF NBPAL GT 15 60 TO 22 NBPAL4 NBPAL 1 TPAL 13 0 GPAL 1 0 00 24 N 1 NBPAL WRITE 1 1052 N READ 1 1054 PAL M WRITE 4 1056 READ 1 1054 TPAL H 24 CONTINUE 60 TO 40 C RABATTEMENT ET INJECTION 60 CONTINUE IF IOPT EQ LR WRITE 1 1062 IF IOPT EG LI WRITE 1 1063 READ 1 1002 2 2 152 TP FLOAT IH2 FLOAT IH2 60 FLOAT I52 3600 WRITE 1 1064 IF 1O
35. BUREAU DE RECHERCHES G OLOGIQUES ET MINI RES SERVICE G OLOGIQUE NATIONAL B P 6009 45060 Orl ans Cedex T l 38 63 80 01 DESCRIPTION ET MODE D EMPLOI DU LOGICIEL IA0 POUR LE TRAITEMENT DES DONN ES D ESSAIS SUR FORAGES G OTHERMIQUES VERSION 1 0 POUR SYST ME R2E MICRAL 21D par A MENJOZ y3 R G M T 125 1983 BIBLIOTHEQUE INSTITUT MIXTE DE RECHERCHES G OTHERMIQUES B P 6009 45060 Orl ans Cedex T l 38 63 80 01 Rapport du B R G M 83 SGN 578 IRG Ao t 1983 Le logiciel IAO Interpr tation Assist e par Ordinateur dont le mode d emploi est d crit dans ce rapport est destin faciliter l interpr tation des donn es d essais sur forages r alis s de mani re syst matique et r p titive 11 se compose d un ensemble de programmes facilitant la tache de l op rateur depuis la saisie des donn es jusqu la sortie du document final dactylographi 11 s agit essentiellement d un logiciel d assistance effectuant le pr traitement des donn es la visualisation sous forme de courbes puis la mise en forme des r sultats obtenus au cours de l interpr tation Bien qu il soit possible d effectuer une interpr tation automatique par l ordinateur la d termination des param tres du r servoir est laiss e aux soins de l op rateur en raison des nombreux cas d esp ce pr voir pour la mise au point d une proc dure fiable et g n rale Le pr sent document est
36. D m 5 3 14 REMARQUES D b art max 115 D viation dans 6 22 CHEF DE PROJET ING FORAGE ING ESSAIS 23 10 2 ANNEXE 2 Fac simil d un rapport interpr tation d essais r alis int gralement avec le syst me informatique MICRAL 24 B R G M DEPARTEMENT Orl ans le 31 Janvier 1983 R f YY 83 XXX AM am SITE FICTIF FORAGE FOR INTERPRETATION DU TEST DE FORMATION DU 30 01 83 SOMMAIRE I RESUME ET RAPPEL DES CARACTERISTIQUES II INTERPRETATION DU TEST DE FORMATION 11 1 Caract ristiques de l essai 11 2 Interpr tation de l essai III CONCLUSIONS IV ANNEXES 25 1 RESUME ET RAPPEL DES CARACTERISTIQUES L interpr tation de l volution de la remont e de pression cons cutive l essai de production art sienne montre un r servoir de transmissivit moyenne de l ordre de 18 D m et un ouvrage d velopp facteur de skin n gatif Rappel des caract ristiques Niveaux producteurs entre 1613 m et 1721 m d v Hauteur productrice cumul e h 26 3 m d v 23 m vert Porosit moyenne 15 Temp rature stabilis e 1598 m d au d bit moyen de 115 m3 h Tf 58 6 C Viscosit moyenne de l eau u 0 5 Salinit moyenne 30 g l 4 1 Compressibilit totale Ct 10 atm Pression art sienne mesur e Par 8 5 bars potentielle fin de build up 2 Pression statique extrapol e Pe 159 35 kg cm 1598
37. ENT M3 H DATA ZL2 PRES SION NNN ON ta 69 DATA ZL3 PUIS GANC E EL ECTR IQUE DATA IDC ID1 102 103 1D4 4600 360 425 740 600 DATA IYL1 IYL2 IYL4 1YLS 3494 599 1649 3449 DATA IXN IXX IYN IYX 4150 150 1075 7079 DATA 214 2157 PROD 7 DATA ZL PRE 151 0 0177 4 15 2 DATA ZL7 FORA GE DATA ZL8 TEMP ERAT URES DATA ZL9 D INJEC Cissi PON PN DE 6X AN TQN INT QN 106 25 IQX INT 0X Xs FLOAT IDQ Y D8 X 0 5 NTIK INT Y 1 Y 450 XXx40 D0 IDY INT Y DP 4 PX PN S IPN T INT PN IDP 4 INT DP 1 IDX 800 2 QXXPX QNEPN XO 04 5 10P 2 INTCDP 2 PEH QNXPNx0O 04 IPN 2 INT PEN ECHPD 4000 PX PN ECHPE 4000 0 04 C DO 100 IC 1 2 IXNC IXN IDCX IXXC IXX IDCXIC ENCODE IBUF 1000 IXHC IYX IXXC IYX IXXC IYN 1000 FORMAT PUSPA 16 16 sPD PA 16 16 16 16 CALL ECRV24 IBUF 54 0 ENCODE IBUF 1002 DON DON IXNC IYX 1002 FORMAT PA 16 16 16 16 PU CALL ECRV24 1BUF 32 0 ENCODE 4004 1004 FORMAT 01 1 0 910 15 0 27 CALL ECRV24 CIBUF 17 0 60 TO 2 4 10 2 IX IXAC IDO ENCODE IBUF 1006 IX IYL2 1006 FORMAT Beb 16 16 CALL ECRV24 IBUF 18 0 ENCODE IBUF 1008 ZL2 1 1 1
38. I est la longueur r elle d vi e la perte de charge dans la compl tion est calcul e pour un d bit moyen gal la moiti du d bit nominal La longueur prendre en compte XLP ou XLI est la longueur d vi e entre le sabot et la bas du dernier niveau producteur Remarques On notera que les estimations des pressions et des puissances lectriques n cessaires ne tiennent pas compte 1 pression du fluide souhait e dans le r seau de surface Cette derni re donn e fonction notamment des caract ristiques des pompes de la valeur du point de bulle et de l installation de surface en g n ral peut tre facilement introduite dans le pro gramme Pour le puits de production la pression de refoulement en t te est suppos e nulle la puissance lectrique d termin e est donc un minimum majorer de la valeur correspondante la sur pression d sir e en surface Pour le puits d injection la pression d aspiration la pompe est suppos e nulle la puissance lectrique d termin e est donc dans ce cas un maximum minorer de la valeur correspondante la pression du fluide souhait e dans le r seau Ainsi aux pertes de charge pr s dans le r seau de surface seule la puissance lectrique totale de pompage a une signification pratique imm diate 10 8 2 RAPPEL DES FORMULES UTILISEES Les hypoth ses classiques utilis es pour la d termination du comportement hydrodynamique d un forage en exploi
39. IR 122 000 PRESSION ARTESIENNE POTENTIELLE 8 500 OPTION CALCUL INTERFERENCE 1 DONNEES PUITS PERMEABILITE INTRINSEQUE 805 HAUTEUR PRODUCTIVE CUMULEE 23 000 RAYON PUITS D INJECTION 7 620 SKIN 3 700 DIAMETRE INTERIEUR TUBAGE 15 940 LONGUEUR 1813 000 LONGUEUR INJECTION RESERVOIR 150 000 atm 1 cm m kg cm2 48 111 1 1 PUITS DE PRODUCTION Q 50 100 125 150 175 200 225 250 m3 h 3 7 7 3 9 2 11 0 12 9 14 7 16 5 18 4 kg cm2 dPc 5 1 6 2 4 3 3 4 4 5 6 6 9 8 3 kg cm2 111 2 PUITS D INJECTION 111 2 1 Temp rature d injection 30 00 Viscosit 83 Estimations apr s 1 00 mois d exploitation d bit constant Q 50 100 125 150 175 200 225 250 m3 h dPh 4 6 9 5 12 0 14 5 17 0 19 5 22 1 24 6 kg cm2 dPc 6 2 0 3 0 4 2 5 5 7 0 8 7 10 5 kg cm2 dPap 7 0 7 0 7 0 7 0 7 0 7 0 7 0 7 0 kg cm2 em ne am ae ue ae ne O ee ee 49 111 2 2 Temp rature d injection 35 00 Viscosit 74 Estimations apr s 1 00 mois d exploitation d bit constant mm LIL
40. LOG LOG HORNER Q CONST SEMI LOG MDH EIN HORNER Q VAR ABAQUE SPE CODE CHOISI 7 2 DONNEES GRILLE NB MODULES SUR X 11 2 3 LONG MODULE X EN MM 13 2 70 CONTROLE ETAT TRACEUR taper return LIBELLE AXE Y 50 CAR 2 LIBELLE 1 AXE X 70 CAR MAX 7 FIGURE 2 TP DT 7 DT LIBELLE 2 AXE X 70 CAR 2 TEST DE PRODUCTION SITE FICTIF CAS1 LIBELLE 3 AXE X 70 CAR 2 REMONTEE DE PRESSION APRES FERMETURE DIAGRAMME DE HORNER Trac du dessin puis retour au menu de d part Remarque Pour chaque nouveau dessin le programme re demande les param tres de la grille et du trac des points ce qui permet d effectuer un agrandissement de certaines z nes d un graphique ou un cadrage La taille maximale de la grille est de 210 mm selon X et de 160 mm selon Y Les figures 1 5 repr sentent un exemple des divers types de courbes disponibles selon le code de trac choisi annexe 1 6 DETERMINATION DES PARAMETRES DU RESERVOIR Les relations utilis es pour la d termination des param tres du r servoir sont rappel es dans l exemple de traitement figurant en annexe Le syst me d unit s employ est le syst me CGS DARCY 7 MISE A JOUR ET EDITION DU RAPPORT D ESSAI La mise jour du rapport s effectue par correction d un rapport type pr enregistr l aide du logiciel de traitement de textes WORDSTAR classique adapt
41. NE 22 CONTINUE 403 FORMAT D11 0 510 15 0 2 PU PA 16 16 sLB 12 41 23 CONTINUE IF 1 00 6 6 60 TO 82 CALL ECRY24 44 0 50 TO 28 ENCODE 104 IC YAD IC YAF Co 104 FORMAT 6 7 I6 16 L6 PU Casas AXE Y KORNER OU HDH CALL ECRU24 IBUF 40 0 IF 1 E HBL 60 0 21 42 ININY INT UMINY 00 20 J 2 10 IDY INT DY XX FLOAT INE IMINY IDY X ALOG10 XX NBL NBHY 1 X Y xDXM FLOAT DO 44 154 881 IC41 IC INT X IC YAD I 1 XDYMx40 ENCODE 104 1014 YAD IC4 YAF INE INE IDY CALL ECRV24 CIBUF 40 0 ICX1 XAD 440 20 CONTINUE ICY4 IC 60 60 TO 21 ENCODE IBUF 116 ICX4 ICY4 INE ETX 116 FORMAT DI1 0 S10 2 0 3 PUsPA 16 16 LB 13 61 Co AXE X LINEAIRE MII CALL ECRV 4 IBUF 44 0 60 CONTINUE ENCODE IBUF 104 XAD IC XAF IC XNE UMINX DX CALL ECRV 4 IBUF 40 0 DO 42 I 1 NBL 44 CONTINUE IC XAD 1 1 DXHx40 Doses XNE XNE DX C LIBELLE AXE Y ICX4 IC 240 Decais ICY1 YAD 160 28 WRITE 1 1000 ENCODE IBUF 118 ICX1 ICY1 XNE 4000 FORMAT 1X LIBELLE AXE Y 50 CAR MAX 2577 118 FORMAT D11 0 510 2 0 3 PU PA 16 16 LB F4 1 41 41X 50 1Hx READ 1 1002 LIB K K 1 70 CALL ECRU24 CIBUF A8 Di 1002 FORMAT 7141 ENCODE IBUF 104 IC YAD IC YAF ENCODE IBUF 108 CALL ECRUZ4 IBUF 40 0 108 FORMAT 010 1 510 2 0 3 62 CONTINUE CALL ECRV24 CIBUF 45 0 ILX 10X1 280 21 CONTINUE ii 1
42. ONTINUE STOP END PROGRAM PREVIS C C CALCUL DES PREVISIONS D EXPLOITATION D UN DOUBLET C ESTIMATION DES PRESSIONS D EXPLOITATION ET DES C PUISGANCES ELECTRIQUES CONSOMMEES Cao VERSION 1 0 DU 30 01 83 Ceana AUTEUR HENJOZ C DIHENSION 8 39 PH 31 PC 31 01 39 PD 31 PE 31 1 N80 2 PHAX 3 TTIC9 THUC9 9 INTEGER 1 ICOD IR IDE IBUF 80 ESC 111 LIBFP 6 6 DATA 8 50 100 125 150 175 200 225 250 DATA P1 3 1415927 DATA 1C0D 1H0 DATA 25 50 7 DATA 8 31 5 97 Coon Co LECTURE DONNEES C WRITE 1 2046 READ 1 1052 IVERS WRITE 1 1000 WRITE 1 1001 READ 1 1003 LIBFP I 1 1 6 WRITE 1 1004 READ 1 1022 XKP XKI XKP IF CIVERS NE 1 WRITE 1 1005 IF IVERS NE 1 READ 1 1022 WRITE 1 1006 READ 1 1022 HF HP HP 100 IF IVERS NE 4 WRITE 1 1007 IF IVERS NE 1 READ 1 1022 WRITE 1 1008 READ 1 1022 UD WRITE 1 1010 READ 1 1022 D D Dx100 WRITE 1 1042 READ 1 1022 WRITE 1 1014 READ 1 1022 5 RHA RWPXEXP 5 WRITE 1 1010 READ 1 1022 PHI WRITE 1 1018 READ 1 1002 C1 WRITE 1 1050 READ 1 2022 IF III E 1C0D 60 TO 20 WRITE 1 1054 READ 1 1056 1 1 1 8 20 CONTINUE WRITE 1 1020 READ 1 1022 DITF DITP DITP 10 WRITE 41 1004 READ 1 1022 XLTP XLTP XLTP 1000 WRITE 1 1026 READ 1 1022 XLF
43. PT EG LR OR JOPT EG LF DP DP DT T TP VTH 0 IF CIOPT E LI 60 TO 52 IF IOPT EG LF 60 TO 50 IF IOPT EQ LR 60 TO 52 0 lt IF IR EG IRR 60 5 DO 26 M 2 NBPAL1 H N 1 UV DPR DT UU VU TPAL M VU TPAL N 0640 UTH UTH VUXQPAL H 26 CONTINUE 60 TO 52 50 VTH DPR DT DT C TRANSFERT BUFFER DISQUE Cod XUTH ICC UTH XPCICC P DT ILL jet IBUF GE 7 GO TO 34 32 IBUF IBUF 1 DTB IBUF DT THB IBUF UTH DPB IBUF DP PB IBUF P 60 TO 36 34 IBUF 0 WRITE 6 DTB K THB K DPB K PB K K 1 7 60 0 32 36 CONTINUE Coca OUTPUT IMPRIMANTE IF IOPT EG LB OR IOPT E LF ANRITE 7 1010 ICC IH 15 DT DP NTH P IF CIOPT EG LR OR IOPT E LI ANRITE 7 1074 ICC IH 1M 15 DT OP P IF IC LE 50 GO TO 10 IF TOPT EG LB OR IOPT EQG LF WRITE 7 1022 IF IOPT EG LR OR IOPT EG LI WRITE 7 1070 60 TO 8 22 CONTINUE WRITE 1 1046 STOP 42 CONTINUE IF I0FT EG LB OR IOPT EG LF WRITE 7 1022 IF IOPT EG LR OR LOPT EG LI WRITE 7 1070 WRITE 7 1014 C TRANSFERT DERNIER BUFFER DISQUE IF CIBUF E6 7 60 TO 30 18UF1 1BUF 1 DO I IBUF1 7 078 1 50 THB 1950 DPD I D PB I zO 38 CONTINUE WRITE 6 DTB K THB K DPB K PB K K 1 7 30 CONTINUE ENDFILE ENDFILE 9 REAL 1 1007 40 CONTINUE C C IDENTIFICATION ESSAI READ 9 1042 ITAB J J 1 30 READ 9 4042 CIDAT K K 1 8
44. RMAT 3 6 8 F7 0 2X m3 h 1034 FORMAT 1X dPh 1X 8 F5 1 2X ko cm2 gt IF L EG 2 50 TO 2 4036 FORMAT 4X dPc 1X B F5 1 2X ka cawo IF IDES EQ 1 0R ICOUR 61 1 60 10 68 1038 FORMAT 4X dPap 8 F5 1 2X ka ca2 PN 0 1040 FORMAT 2X Pd 1X 8 F35 1 2X k0 cm2 DO 64 LL 1 3 1042 FORMAT 2X Pe 1X 8 F5 0 2X x84 IF PD B LE PHAX LL GO TO 66 1044 FORMAT 1H1 44 CONTINUE 4046 FORMAT 1 111 1 11 PUITS DE PRODUCTION 14X 19 66 PX PHAXILL 1048 FORMAT 1 111 BILAN DES PRESSIONS ET PUISSANCE CALL CADRE 9 1 0 8 PN PX 2 TS LIBFI 1 4 ELECTRIQUES 7X 44 68 CONTINUE 1050 FORMAT 2 DE 50 100 125 150 175 200 2205 050 Costas 12X CONFIRMATION 0 N 27 J 4 4052 FORMAT 11 00 12 1 2 1054 FORMAT 1 CHOISIE 8 VALEURS 8F4 0 237 4056 FORMAT 8F4 0 IF IPAG EQ 1 WRITE 2032 4060 FORMAT 41 WRITE J 20148 IPAG TIN UI 4062 FORMAT 4X CONTROLE ETAT TRACEUR OU T ET WRITE J 2040 TS 1064 FORMAT INSLT SP1 WRITE J 1030 1066 FORHAT 1 NB TEMP INJECTION MAX 9 11 7 7 WRITE 1032 0 1 1 1 8 1068 FORMAT 4 7115 TEMPERATURES 9F4 0 2 WRITE 4 1030 1411X VALEURS CROISSANTES WRITE J 1034 PH I 1 1 8 4070 FORMAT 9F4 0 1036 PC I 1 1 8 4072 FORMAT 4X LISTE VISCOSITES 9F5 2 WRITE J 1036 PA PA PA PA PA PA PA PA 1074 FORMA
45. T 9F5 2 WRITE 1030 1076 FORMAT 1X LISTE PRES ART POT 9 5 1 2 77 3 1040 PD 1 1 1 8 1078 FORMAT 9 5 1 WRITE 3 1030 2000 FORHAT 2X DONNEES WRITE J 1042 1 1 8 2002 FORMAT 2 INJECTION F10 3 WRITE 1030 2004 FORMAT 2X UI CP F10 3 42 CONTINUE 2006 FORMAT 2X RWI CM F10 3 60 TO 43 2008 FORMAT 2X S F40 3 62 CONTINUE 2010 FORMAT 2X TUBAGE DITI CH F10 3 17 1 2012 FORMAT 2X TUBAGE XLTI M F10 3 IF ICOUR EQ 4 1722 2014 FORMAT LONG INJECTION XLI H F10 3 IF ICOUR NE NBT GO 63 2016 FORMAT 2X PA KG CH2 F10 3 1722 2018 FORMAT 1 111 2 12 TIN TIN 12X Temperature d iniection F3 2 C 63 CONTINUE 212X Viscosite pA CALL COURBE PD 1 PE 1 TIN IZ ICOUR 2020 FORMAT 2X AUTRES TEMPS 0 N 25 43 CONTINUE 2022 FORMAT A1 2024 FORMAT 2X DONNEES PUITS DE PRODUCTION 2X 27 WRITE 1 2020 2026 FORMAT 2X PERMEABILITE INTRINSEQUE F40 3 0 READ 1 2022 IR 1 2X HAUTEUR PRODUCTIVE CUMULEE F40 3 m IF IR EG 1C0D 60 TO 8 2 2X VISCOSITE DU FLUIDE RESERVOIR F10 3 ENDFILE 6 3 2x DISTANCE ENTRE 5 10 3 4 RAYON PUITS DE PRODUCTION F10 3 cu 1000 FORMAT 2X DONNEES POUR PRODUCTION 5 2X SKIN DE PRODUCTION 1 F10 3 1001 FORMAT 2X SISLE DU FORAGE 6A1 227
46. TATION t mois 120 0 FORAGE FORS 59 10 4 ANNEXE 4 Listings des logiciels utilis s fichier ESSAIS S fichier PLOTEST S fichier PREVISS S Remarques compl mentaires pour la modification des programmes 60 PROGRAM ESSAIS PRETRAITEMENT POUR INTERPRETATION DES ESSAIS C SORTIE DU TABLEAU RESULTATS IMPRIMANTE C CREATION DU FICHIER DESSIN ESSTRA S PLOTEST CREATION DU FICHIER SORTIE ESSOUT S IMPRIMANTE VERSION 1 0 DO 30 01 83 AUTEUR HERJOZ DIMENSION 300 300 4 4PAL 15 16 8 1 8 8 PBIB INTEGER oi 1748 30 1081 IR IRR IPU IPUN 1 IDRAW IREP IOPT LR LI LB LF ESC H1 N2 N3 N4 FLAG DATA IRR IPUN IREP 1HX 1H 1H0 DATA LR LI LB LF 1HR 1H1 1HB 1HF DATA ESC 27 N1 54 N2 48 N3 59 52 1000 FORMAT 2X DATE FIM DE PRODUCTION 13 212 1001 FORMAT 2X DATE FIN D INJECTION 13 212 1002 FORMAT 13 1X 212 F8 9 1003 FORMAT I3 41 212 F8 3 1004 FORMAT 2X PRESSION PWFO F8 3 25 1006 FORMAT F8 3 4008 FORMAT 2X DUREE DE PONPAGE 13 212 1009 FORMAT 2X DUREE 0 INJECTION 13 212 1010 FORMAT 1 15 2X 13 1H 12 1H 12 2X F0 5 0X FB 3 12 210 3 2 28 3 1012 FORMAT 13 41 212 8X F5 1 1044 FORMAT 181 1046 FORMAT 1X 6X t 9X dt 3X tp dt dt 244 ou 3 point h
47. THUCI 1 1 HBT WRITE 1 1076 READ 1 1078 1 1 8 CONTINUE WRITE 1 2042 READ 1 1022 TS IF IDES EG 1 60 TO 60 WRITE 1 1062 PAUSE ETAT ENCODE 1064 CALL ECRV24 IBUF 9 0 60 CONTINUE Ca as DO 13 ICOUR 1 NBT PA TPACICOUR TIN TTICICOUR UI THUX ICOUR IPAG IPAG C 14 5 2 298 GO TO 40 42 40 IVERS 40 X2 X1 GO TO 44 Ad X41 YIwXKIxHIx 400 XKP HP X2 XI XKI XKP 44 CONTINUE X3 UIxFACO 4 4PIXXKIXHIX100 XA PHIXULECT 4 XXR1 PC2 1 057E 7 UI x0 21 GAH PHI 1 PHI 0 5 G1 PHIAUOXCT XKIKGANXHI X100 RU UI UD RUI 1 RU XTI X4 RWIXRUI T1 1 DO 13 L 1 IDES NBQC HBQ L DO 10 1 1 NBQC JEI L 1 XB XMI X3 0 J XM X1 8 1 1 1 1 XTE G 61 0 J FACQ FACS ALOG 6 1 9338 1 RUI FACS FACS ALOG RUI FACS FACSEXNI PHCI PHOD FACO 67 1270 9678 1016 FORMAT 2X PHI DEC F40 3 BEI ef n 6 3 179 791 4018 FORMAT 2 1 10 3 PC I PC I wgXLTI 1020 FORMAT 2 10 3 DD 20 4RWI1 1022 FORMAT F10 3 6 Q J 1024 FORMAT 2X TUBAGE XLTP M F10 3 PC I PC I PC2XXLIx 00xx1 79 DDKx4 79 1026 FORMAT 2X LONG PRODUCTION XLP M F1D 3 PD ID PH I J PC I PA 4028 FORMAT 10 3 PE I PD I w 8 J z0 04 4030 FORMAT 1 66 1 40 CONTINUE 4032 FO
48. cteur de skin Pression extrapol e pour le calcul P 1h 157 9 kg cm2 il vient S 3 32 31 III CONCLUSIONS L interpr tation de l volution de la remont e de pression cons cutive l essai de production art sienne montre un r servoir de transmissivit moyenne de l ordre de 18 D m et un ouvrage d velopp facteur de skin n gatif voisin 3 7 32 TABLEAU RECAPITULATIF Niveaux producteurs entre 1613 m et 1721 m d v Hauteur productrice cumul e h 26 3 m d v soit 23 m v Porosit moyenne 15 Temp rature stabilis e 1598 m d au d bit de 115 m3 h Tf 58 6 C Viscosit moyenne de l eau u 0 50 cp Salinit moyenne 30 g 1 4 1 Compressibilit totale Ct 10 atm Pression art sienne mesur e Par 8 5 bars potentielle Pression statique extrapol e 159 35 kg cm2 abs 1598 m d v moyen d apr s retenu Log Log Horner MDH Transmissivit relative kh u 37 02 40 45 37 02 36 74 D m cp Transmissivit intrins que kh 18 51 20 22 18 51 18 37 D m Perm abilit intrins que k 0 805 0 880 0 805 0 800 D Facteur de skin 5 3 71 3 29 3 71 3 32 33 IV ANNEXES SITE FICTIF FORS 34 REMONTEE DE PRESSION DU 30 01 83 DUREE DE PRODUCTION 12 000 Heures 151 304 E E E E FE E lt lt lt lt lt lt lt lt lt lt lt lt i a dt heures dp kg cm2 tp dt dt ou fonction S t heures 13 59 36 13 5
49. deux forages du doublet et sont destin s contr ler la continuit hydraulique du r servoir ou ses anomalies entre les puits Selon la localisation g om triques des deux t tes de puits il peut tre pratiqu soit un test d interf rence grande distance entre les t tes de puits avec pompage ou injection dans l un des forages et enregistrement de l interf rence hydraulique dans l autre puits soit un test en boucle cas des t tes de puits voisi nes avec pompage dans l un des puits et r injection dans l autre puits Les travaux d interpr tation consistent tout d abord s lectionner l int rieur de l enregistrement global certains groupes de donn es correspondant aux s quences analyser puis effectuer diverses op rations sur chaque donn e soustraction de valeurs initiales calcul de nouvelles variables Ayant r alis es ces transformations les points doivent tre report s sur divers graphiques log log semi log etc pour permettre l interpr tation La constitution d un fichier num rique et l utilisation d une cha ne de programmes de traitement permettent donc une grande souplesse d utilisation en limitant les risques d erreur surtout lorsque les essais sont nombreux et syst matiques comme dans le cas de l ing ni rie g othermique Enfin hormis la recherche d une m thode de traitement des donn es fiable et rapi de le logiciel pr sent permet galement de satisfaire un autre
50. e p riode de production de 10 heures d bit variable 4 paliers S quence d exploitation NOM DU FICHIER DES DONNEES XDATBUV S F v OPTION TYPE D INTERPRETATION RABATTEMENT R INJECTION 1 BUILD UP B FALL OFF F CODE B E DUREE DE POMPAGE I3 2I2 10 0000 DATE FIN DE PRODUCTION 13 212 10 0000 PRESSION PWFO F8 3 175 290 met DONNEES PALIER 1 F10 3 2 130 DEBUT 0 000 FIN F10 3 2 2 DONNEES PALIER 2 F10 3 2 70 DEBUT 2 000 FIN F10 3 5 DONNEES PALIER 3 DEBIT 10 3 2 115 DEBUT 5 000 FIN F10 3 7 DONNEES PALIER 4 DEBIT F10 3 2 100 DEBUT 7 000 FIN F10 3 2 10 Remarque On remarquera que dans la r ponse aux questions relatives la variation de d bit les temps sont fournir sous forme d cimale heures avec l origine au d but de l essai 5 TRACE DES COURBES PROGRAMME PLOTEST Ce programme exploite le fichier des r sultats cr par le programme ESSAIS ESSTRA S Les donn es introduites au cours de l ex cution concernent le choix des param tres pour le trac de la grille et des points exp rimentaux Les diverses options possibles sont fonction des choix effectu s lors de l ex cution du programme de pr traitement ESSAIS activation 4 PLOTEST S quence d exploitation exemple 1 pr c dent OPTION TRACE CODE
51. ensuite tre archiv s sur disquettes pour reprise ult rieure 3 SAISIE ET ENREGISTREMENT DES DONNEES Les documents n cessaires sont la fiche technique d essai voir annexes et les mesures fournies par la soci t de service logging Cr ation du fichier et formats EDIT XXX5BUC S par exemple SITE FICTIF 5 commentaire d identification 30 carac 30 01 83 commentaire date 8 caract res 12 0107 175 200 64 2 1 ligne par point temps pression 12 0108 175 215 63 9 temp rature 13 212 F8 3 F5 1 etre es SS die S Sis ee l heure est introduite sous la forme de esa heure minute seconde ligne blanche E sortie de l diteur L enregistrement de la temp rature n est pas indispensable mais l on peut par la suite visualiser galement l volution de cette variable cas des essais d injection par exemple La fin du fichier d entr e sera d tect e par la ligne blanche c est dire l absence de point en colonne 4 d une ligne Remarques Avec quelques modifications la phase de saisie des donn es peut tre remplac e par deux proc dures plus rapides la transmission directe des donn es par le r seau t l pho nique TRANSPAC le syst me tant quip du logiciel correspondant Cette m thode suppose que le v hicule de logging soit muni du syst me de communication n cessaire la lecture d une disquette ou d une cassette enregistr e sur le site il convie
52. envahie par les eaux de r injection effet pari tal thermique Le calcul est effectu d bit d exploitation constant et pour un doublet isol En pratique les param tres principaux du r servoir perm a bilit et paisseur productive ne sont pas homog nes et d autre part cette variation n est caract ris e que par deux d terminations soit une valeur chaque forage Le calcul n cessite donc des hypoth ses suppl mentaires le comportement propre d un ouvrage sera d termin par ses caract ristiques sp cifiques d duites des essais pour le calcul de l interf rence de l autre puits 3 options sont offertes option 1 k et h sont identiques option 2 k et h sont distinctes et donn es option 3 k et h sont une moyenne des valeurs donn es La pression art sienne introduire doit tre calcul e au pr alable version actuelle elle est d duite de la valeur de la pression statique fond extrapol e r f essais et de temp rature moyenne la colonne de fluide 11 s agit donc de la pression art sienne potentielle correspondant l tat d un forage en quilibre thermique stationnaire La valeur introduire est positive pour un forage art sien La chute de pression due aux pertes de charge comprend deux composantes la perte de charge dans le casing un seul diam tre interne est consid r depuis la surface du sol jusqu au sabot La longueur prendre en compte XLTP ou XLT
53. eures heures ka cm2 2X fonction S 2X 4 3 1022 FORMAT 1 5B 1H 3 1024 FORMAT 1 30 1 11X REMONTEE DE PRESSION DU 841 1026 FORMAT 1 DUREE DE PRODUCTION F8 3 Heures 1 Pwfo F8 3 77 1027 FORMAT 1X DUREE D INJECTION 28 3 Heures 1 4 F0 3 2 1028 FORMAT 77 1036 FORHAT 2X NOM DU FICHIER DES DONNEES 1038 FORMAT 1141 1042 FORMAT 3041 4046 FORMAT 2X DEPASSENENT 1048 FORMAT 2X OPTION DEBIT VARIABLE 0 N 7 1050 FORMAT PALIERS DE DEBIT I3 1052 FORMAT ON DONNEES PALIER 13 2X DEBIT F10 3 i 1054 FORMAT F40 3 1056 FORNAT 2X DEBUT F40 3 FIN F10 5 2 7 1060 FORMAT 2X OPTION TYPE D INTERPRETATION 12X RABATTEMENT R INJECTION 32X BUILD UP B 40X FALL OFF F CODE 2 1062 FORNAT 2X DATE DEBUT PRODUCTION I3 212 1063 FORMAT DATE DEBUT INJECTION 13 212 9 7 1064 FORMAT 27 PRESSION PWS0 F8 3 227 1066 FORMAT 1X 3041 145 RABATTEMENT EN DEBIT DU 7 651 1068 FORMAT 1 FB 3 kg cu2 4070 FORMAT 1 46 1 1072 FORMAT 1X 1 UND TM BA 217 3 point 3x Heures 2X Heures 2x kua cg2 2X 4 1074 FORMAT 47 15 2 13 12 46 12 27 28 5 2 28 2 12X FB 3 1080 FORMAT 1X 6841 1082 FORM
54. ge sont calcul es pour une installation fonctionnant en boucle avec un puits de pompage et un puits de r injection op rant m me d bit Les donn es utilis es sont issues des r sultats des essais sur le forage FORS seul I CARACTERISTIQUES DE LA BOUCLE GEOTHERMALE Les caract ristiques moyennes retenues sont les suivantes Hauteur productive cumul e Se dre se ere 23 m Porosit moyenne CUA RAS 15 Transmissivit intrins que dass se eee 18 5 D m Perm abilit intrins que Soie este 805 mD Temp rature moyenne du r servoir en production 58 6 C Viscosit du fluide acesa cie rr 0 50 cp Pression art sienne potentielle exhaure 8 5 kg cm2 45 1 2 DISTANCE ENTRE OUVRAGES Conform ment au programme d op ration la distance entre les puits au niveau du r servoir est prise gale 1200 m toit CARACTERISTIQUES PARTICULIERES DES FORAGES Tubage 7 jusqu 1599 m d vi s simplification R servoir 6 facteur de skin 3 71 Puits de r injection Tubage 7 jusqu 1813 m d vi s R servoir 6 facteur de skin 3 7 hypoth se II RELATIONS UTILISEES POUR LES CALCULS Gus ou mue lt lt lt lt lt lt lt lt lt lt rm II 1 VARIATION DE PRESSION HYDRODYNAMIQUE dPh iS dPh Ln 8 atm Le
55. h Porosit moyenne g 2 2 Salinit moyenne La salinit pr visionnelle est de 30 g 1 II 2 3 Compressibilit totale moyenne lt lt lt lt lt lt lt lt lt lt lt lt Les zones productrices tant r parties entre 1613 et 1721 m d v la compressibilit est valu e la profondeur verticale moyenne de 1550 m D apr s les abaques Dodson et Standing la compressibilit de l eau est estim e 5 1 Cw 4 18 10 atm pour une pression de confinement de 155 atm et une temp rature de 58 6 C D apr e l abaque de Hall et pour une porosit 15 1 compressibilit effective de la roche vaut 5 1 Cr 5 7 10 atm d o la compressibilit totale 5 4 1 Ct 9 88 10 10 atm 11 2 4 Viscosit de l eau de la formation La viscosit estim e d apr s l abaque de Chesnut est de l ordre de u 0 50 cp 28 11 2 5 Perm abilit et facteur skin La perm abilit et le facteur de skin sont d termin s partir de l enregistrement de la remont e de pression finale avec la sonde TP place 1598 m d v Les points exp rimentaux sont donn s au tableau I e Courbe type diagramme log log La premi re estimation des param tres est effectu e l aide de l abaque Flopetrol SPE 8205 1979 et du report des points exp rimentaux en graphique bi log figure 1 Points de coincidence des graphiques dP 1 kg cm
56. hiers La proc dure est identique celle du paragraphe 7 cr ation d un nouveau fichier temporaire TAMPON S par ex insertion du fichier texte type XXXRAP S par exemple insertion des r sultats du calcul PROUT S corrections et mise en page du document global Un exemple de document complet est fourni en annexe 3 16 10 ANNEXES 10 1 ANNEXE 1 Diagrammes disponibles par utilisation PLOTEST Exemple de fiche caract ristiques pour interpr tation 17 0 EWNI IA NI 1553 30 NOILVISVA CAS 1 AGRAMME LOG LOG N HEURES E FICTIF TURE DI REMONTYE DE PRE 185 184 N 183 CNI IA NI 155388 181 10 3 180 CAS 1 IAGRAMME DE HORNER REMONTEE 19 185 184 183 182 181 180 1553 1 3 DT EN HEURES TION SITE FICTIF CAS PRES FERMETURE DIAGRAMME SEMI LOG E REMONTEE PRESSION EN KG CM2 183 p 182 181 18070 1 0 5 0 FIGURE 5 FONCTION 5 TEST DE FORMATION SITE FICTIF CAS d REMONTEE DE PRESSION APRES FERMETURE DIAGRAMME DE HORNER Q VARIABLE 22 FICHE CARACTERISTIQUES ESSAI DONNEES POUR INTERPRETATIONS
57. ion d bit constant dPh 4 6 9 5 11 9 14 4 16 9 19 4 21 9 24 4 kg cm2 dPc 6 2 0 3 0 4 1 5 4 6 9 8 5 10 2 kg cm2 dPap 7 0 7 0 7 0 7 0 7 0 7 0 7 0 7 0 kg cm2 111 2 7 Temp rature d injection 40 00 Viscosit 67 Estimations apr s 7 00 mois d exploitation d bit constant iii me 52 111 2 8 Temp rature d injection 45 00 C Viscosit 61 Estimations apr s 7 00 mois d exploitation d bit constant O e e e dPh 4 1 8 3 10 4 12 6 14 7 16 8 19 0 21 1 kg cm2 aD O DD bm e mm 111 2 9 Temp rature d injection 30 00 i Viscosit 83 cp Estimations apr s 120 00 mois d exploitation d bit constant em en srfnWmY e A Ce eue gt E o ou 5 6 11 4 14 4 17 4 2
58. m d Param tres moyens Transmissivit relative kh u 37 02 D m cp Transmissivit intrins que kh 18 51 D m Perm abilit intrins que k 0 805 D Facteur e skin 6 3 71 26 II INTERPRETATION DE L ESSAI DU 30 01 1983 Le test de formation a t r alis le 30 01 83 apr s les travaux de d veloppement de l ouvrage Il se compose d une p riode de production d bit constant art sien suivie d une remont e de pression apr s fermeture L enregistrement utilis pour l interpr tation concerne cette derni re p riode il est obtenu l aide d une sonde de pr cision positionn e sensiblement au toit du r servoir 11 1 CARACTERISTIQUES DE L ESSAI d but du pompage 30 01 83 t 02 00 00 arr t du pompage 30 01 83 t 13 59 36 dur e du pompage tp 12 h fin de remont e de pression 31 01 83 t 02 06 00 d bit nominal D 115 m3 h temp rature 1598 m d v au d bit nominal 58 6 C temp rature de m lange durant la production localisation de la sonde 1598 m d rayon du forage dans le r servoir 7 62 cm 27 11 2 INTERPRETATION DE L ESSAI Ces deux param tres sont d termin s par l examen des diagraphies et du flowm tre Les valeurs retenues sont les suivantes 23 m vertical 26 3 m d vi 15 Hauteur productive
59. nt alors de disposer du lecteur appropri connect au syst me de traitement En version standard le nombre de points est limit 300 points message d passement capacit l cran dans cas augmenter la dimension des tableaux 4 PRETRAITEMENT PROGRAMME ESSAIS a A e rn mm me eue mue Les diff rents programmes sont crits en mode conversationnel l exploitation consiste par cons quent r pondre une s rie de questions Activation du programme ESSAIS Quatre options sont possibles quant au type d essai traiter rabattement injection build up fall off Le programme peut tenir compte d un test d bit variable pour les deux derni res options A partir des donn es lues le programme d termine pour chaque point dt temps depuis le d but de la s quence s lectionn e dP variation de pression tp dt dt variable pour la repr sentation de Horner o tp est la dur e de production ou d injection Dans le cas d un test d bit variable de type build up re mont e de pression apr s production ou fall off chute de pression apr s injection la variable de Horner pr c dente est remplac e par la fonction S suivante 9 tp t dt 5 somme 2 Log HR tp t dt Le programme fournit le tableau des variables interpr t es
60. ogramme PLOTEST F80 PLOTEST S 180 PLOTEST 5 242 FORLIB S PLOTEST N E c Cas du programme PREVISS F80 PREVISS S 180 PREVISS ASV242 FORLIB S PREVISS N E
61. p trolier ou celui des forages d eau Trois groupes d essais peu vent tre identifi s les essais sur forages de production 115 se caract risent g n ralement par une p riode de produc tion soit par art sianisme soit par air lift suivie d une re mont e de pression apr s fermeture Build up Les essais de ce type fournissent l identification des param tres conditionnant la productivit et la nature des traitements ult rieurs ventuels pratiquer d veloppement etc les essais sur forages d injection On distingue tout d abord un test de production identique au pr cedent destin contr ler la productivit et estimer la variabilit des param tres Le forage est test pour sa fonction propre par essai d injection isotherme ou non isotherme Ce test peut tre suivi ventuellement de l enregistrement de chute de pression accompagnant le retour l quilibre apr s fer meture Fall Off type d essai fournit les l ments n ces saires pour appr cier le caract re r versible ou non r versible des param tres d injectivit vis vis de ceux de productivit les essais de type interf rence Les deux groupes d essais pr c dents ne concernent qu un seul forage et fournissent une estimation moyenne des param tres et des variations lat rales de faci s l int rieur d un certain rayon d investigation autour de chaque puits Les essais de type interf rence concernent les
62. tation sont les suivantes puits de faible rayon th oriquement infinit simal soit la solution la ligne source p n tration totale dans le r servoir horizontal ou consi d r comme tel aquif re homog ne isotrope et d extension infinie Dans ces conditions la variation de pression cons cutive au d bit constant produit Q s exprime l aide de la fonction de Theis ou exponentielle int grale El dP El avec les variables suivantes dans le syst me CGS DARCY pression atm 3 d bit de production ou d injection s perm abilit intrins que Darcy paisseur productive cm viscosit du fluide porosit 1 compressibilit totale atm gale la somme des compressibilit s de la roche et du fluide dans r servoir distance au puits o est valu e la pression cm temps 5 COOC gt x O Pour les valeurs suffisamment faibles de l argument de El soit guCr kt lt 0 1 environ peut utiliser la solution asymptotique dP Ln t 0 80907 11 Injection temp rature diff rente dans le r servoir Pour l estimation du comportement hydrodynamique admet l existence d un front thermique abrupt caract ris par un rayon thermique croissant avec la dur e d injection La vitesse thermique ou vitesse de d placement des isothermes est inf rieure la vitesse r elle du fluide et sup

Download Pdf Manuals

Related Search

Related Contents

Sunways SM 60PP 02-13 V2 FR  Sony VAIO VGN-AR71ZU notebook  installation and operation manual urban multi air conditioner    manual de usuario  USER MANUAL  IT - Rasor 1946  File  Seraclone® Anti  Manual London 8  

Copyright © All rights reserved. Failed to retrieve file