Mémoire
Mémoire
Mémoire
Elaboré pour l’obtention du diplôme Master2 en :
MammeriWalid,AllaouiAmhamed,SayahDjeddi,Mohamed Hocine
Thème
OPTIMISATIONDESPARAMETRES
D’UNFORAGEHORIZONTAL
Soutenupubliquementle,25Juin2012,devantlejurycomposéde:
juin2012
Remerciements
Remerciements
Nous voudrions ici témoigner toute notre gratitude à Monsieur Belghit Chafik, notre
directeur de mémoire, pour sa confiance, sa disponibilité pendant l’élaboration de ce
travail. Son soutien nous a permis de mener à bien ce travail dans le Département de Forage et
Exploitation.
Nos profonds remerciements à messieurs les membres de jury, d’avoir accepter de participer
au jury, et d’évaluer notre travail.
XI
Dédicace
Dédicace
Je dédie ce travail à :
Mes parents.
Mes frères et sœurs.
La petite rose Mohamed, Said et Hichem,Mariem et Khadija,Sarah et Lyna
Mes tentes et oncle.
La famille Mammeri
La famille Barir
La famille Karima.
La famille Boukchache
Mes amis et collègues : Allaoui Amhamed Sayah, Djeddi Mohamed Hocine, Belkhira
Hichem,Boumada Athman Ben afan ,Ouis Omar Motiaa.Ben maza Hocem.
Walid Salhi, Baboukha Fathi, Ben Azza Abd Araouf, Boungab Abd Allah, Hamrouni
Mohamed, Derouiche Yakoub,Hichem Aribi,Bedouda Djalal,Mohamed Abd Eldjwad.
Mes amis et frères :Khedir Mohamed:, Ben Seddik Abd El MONIM,Homi Homia.
Walid
xii
Résumé
Abstract
The objective of this work is the optimization of horizontal drilling parameters. For this, an
optimization model has been developed by theoretical and empirical formulas, adopted for the
horizontal drilling. We used also mathematical methods to solve the problem.
Résumé
L’objectif de ce travail est l’optimisation des paramètres de forage horizontal. Pour cela, un
modèle d’optimisation a été développé par des formules théoriques et empiriques, adoptées au
forage horizontal. On a également utilisé des méthodes mathématiques pour résoudre le
problème.
II
Sommaire
Résumé…………………………………………………………………………………………... II
Sommaire……………………………………………………………………………………….. III
Liste des figures………………………………………………………………………………….VI
Liste des tableaux………………………………………………………………………………VIII
Nomenclatures…………………………………………………………………………………...IX
Remerciements…………………………………………………………………………………..XI
Dédicace ………………………………………………………………………………………..XII
Introduction générale…………………………………………………………………………………….....1
III
Sommaire
IV
Sommaire
Conclusion générale……………………………….……………………………………………..37
VI
Sommaire
VII
Liste des tableaux
8
Résumé
9
Liste des figures
VI
Liste des figures
VII
Liste des tableaux
VIII
Nomenclatur
es
Nomenclatures
BPX BP Exploration
DB B Diamètre de l’outil
N Vitesse de rotation
Pression ECD pression exercée par un fluide de forage sous condition de circulation (densité
équivalente de circulation)
T Torque (ft-lbf)
IX
Nomenclatur
es
φ Porosité effective.
X
Introduction générale
Introduction générale
Par conséquent, l’efficacité du forage a une influence sur le succès économique de l’industrie du
pétrole et du gaz.
Les performances de l’outil de forage, qui est habituellement mesurées en termes de coût par
pied pendant un intervalle d'un trou, sont un facteur important dans les coûts de forage globaux. Les
performances de l’outil de forage dépendent de beaucoup de paramètres, telles que la vitesse
d’avancement de l’outil, à quelle profondeur on fore, et le coût des outils de forage, et les autres
ressources nécessaires pour l’opération de forage. D’autres paramètres entrent aussi en jeu, tels
que : type de l’outil, condition d’opération de forage, type de la formation, profondeur du puits,
etc. L’un des facteurs les plus importants qui influent sur les performances du forage est la vitesse
d’avancement de l’outil, qui est définie comme le taux de pénétration de l’outil (ROP). Ainsi,
l’amélioration du taux de pénétration va conduire à un succès économique incontestable.
L’objet de ce mémoire est l’optimisation des paramètres de forage dans un puits horizontal.
Dans un premier temps, on commence par définir le forage horizontal, puis on va présenter les
différentes études et publications qui définissent et traitent le sujet de l’optimisation et les
paramètres concernés.
1
Introduction générale
On passe par la suite au développement de notre modèle d’optimisation, puis son utilisation
dans un logiciel qu’on développera pour cette raison. On termine par une partie calcul puis une
conclusion générale.
2
Chapitre 1 : Overview sur le forage horizontal
1.1. Introduction
Dans la plupart du temps, les puits sont forés avec les formes suivantes :
• Verticale (figure 1),
• type S (figure 2),
• type J (figure 3),
• horizontale (figure 4).
3
Chapitre 1 : Overview sur le forage horizontal
pieds (6 à 10 m) dans le réservoir, tandis qu’un puits horizontalforé dans le même réservoir peut
donner une idée sur les réserves de 2000 à 3000 pieds (600 à 1000 m).
1.2. Avantagesduforagehorizontal
-leforagehorizontalpermetledéveloppementdechampsquin'auraientpuêtre exploités
commercialement autrement.
-dansla pluparts deréservoirs,leforagehorizontalpermetd'augmenterlaproduction et
d'améliorerletauxderécupération,ceciparunmeilleurdrainageetenretardantl'arrivé d'eau. Et voici
quelques exemples de l’utilisationduforage horizontal :
1.2.1. Réservoirsfracturés
Les fractures de ces réservoirs étant sub-verticales, une conséquence directe est que le
meilleurmoyend'enintercepterleplusgrandnombreestdeforer unpuitshorizontal perpendiculairement à
leur direction principale (figure 5).
Figure5:Réservoirs fracturrés
Danslaplupartdesréservoirsmulti–couches,unpuitshorizontalpeutremplacerplusieurs puits
verticaux ou déviés.
4
Chapitre 1 : Overview sur le forage horizontal
La figure (6) illustre un tel cas de réservoir compartimenté ou un seul puits horizontal
remplacesixpuitsverticauxetdeplusaméliorelaproductionenretardantleconing(figure7).
Figure6:Reservoirsmulti–couches.
Figure7:Waterconing
: Water conin
Leforagehorizontaldansunréservoiràbasseperméabilitéestunealternativeàlafracturation de ce
réservoir.
Le drain horizontal se comporte comme une fracture, avec plusieurs avantages :
- ilestplusfacileetpluséconomiquedeforerunlongdrainplutôtqued'essayerdecréerune
fractureéquivalente.
-la direction est parfaitement contrôlée, ce qui n'est pas possible avec la fracturation.
5
Chapitre 1 : Overview sur le forage horizontal
Laproductiondesablesnonconsolidésprésentedesérieuxproblèmespourlimiterlaquantité de sable
entrant dans le puits.
UndrainhorizontalCetteproductionforédansuntelréservoirpermetderéduirelavitesseàlaparoieten
conséquence la production de sable, qui peut aller jusqu'à être totalement éliminée.
Uneautreconséquenceestquecespuitspeuventquelquefoisêtreproduitssansmiseenplace de
coûteusescrépines calibrées.
Beaucoupderéservoirssontproduitsgrâceàunaquifèreactifouparinjectionartificielle.La
productiondiminiertrèsrapidement si le niveau d'eau remonte trop vite dans le puits.
1.3. Inconvénientsduforagehorizontal
6
Chapitre 1 : Overview sur le forage horizontal
1.3.1. Coûtsadditionnels
Il est évident qu'un forage horizontal a un coût plus élevé qu'un forage vertical ou peu dévié.
Les coûts additionnels sont dusà deux facteurs principaux :
Un puits puits foré dans un réservopir on shore peu profond sera comparativement plus cher
qu’un puits foré offshore dans un réservoir profond dont le surcoût pourrait ne pas dépasser 20 à
30%.
1.3.2. Risquesopératoires
Parrapportauxpuitsverticauxoupeudéviéslespuitshorizontauxreprésentent,aucoursde leur
réalisation, un certain nombrede risques supplémentaires.
• Le comportement des formations : Le comportement des formations peut être très différent
lorsqu'elles sont forées à haute inclinaison, en particulier leur stabilité.
7
Chapitre 1 : Overview sur le forage horizontal
• L'évaluation du potentiel de production
1.4. Typesdepuitshorizontaux
- lespuitsàlongrayonutilisantdesrayonsdecourbureexcédant200mètres(gradients<3 deg/10m),
- lespuitsàmoyenrayonutilisantdesrayonsdecourburecomprisentre80et200mètres(3<deg/10m< 7),
- lespuitsàrayonscourtouultracourtavecdesrayonsdecourburepouvantatteindre5à6 mètres,
seulement réalisable avec des équipements spécifiques.
- à build-up simple avec une seule montée eninclinaison de 0 à 90 degrés (figure 8),
- àbuild-upmultiple(deuxsectionsdemontéeouplus,séparéespardessectionsdroitesdites
"tangentielles")
Figure8: Puitshorizontalàunsimplebuild-up
8
Chapitre 1 : Overview sur le forage horizontal
Les puits à long déplacement ont une longue section forée à haute inclinaison jusqu'au toit du
réservoir,suivied'undrainhorizontal.Lalongueurtotalepeutatteindreplusieurskilomètres, le record
étant à ce jour de 12 kmavec BP sur le champ de WitchFarm.
Figures9: puitshorizontauxàdifférentsdéplacements
Et il y a aussi :
Lespuitsenré-entréesontuneautreséquelleduforagehorizontaletconsistentunrepartir
d'unpuitsdéjàexistant,verticalou dévié, qui est alors abandonné etbouchéafindeforerun puits latéral
se terminant le plus souvent à l'horizontale.
Lespuitsmultilatérauxsontundéveloppementdu foragehorizontal.Ilsconsistentàforerune
ouplusieursbranchesàpartird'undrainprimairequipeutêtre vertical,déviéouluimême horizontal (figure
10).
9
Chapitre 1 : Overview sur le forage horizontal
Figure10:Puitsmultilatéralenarrêtedepoisson
1.5.Conclusion
10
Chapitre2:Paramètresde forage
2.1. Introduction
Lesparamètresdeforagesontdesdonnéesprimordiales pourplanifieretréaliserunforage.
Plusieurstentativesthéoriques etpratiquesontétéeffectuées pourcapteretcontrôlerces
données.Lesétudesthéoriquesconsistentàétablirdesrelationsrationnalesouempiriques
entreslesdifférentsparamètres afindepouvoirlesoptimiserparexemple,tandisqueles
appareilsdemesureontlerôle de mesuré avec précisionlesparamètresdanslesdifférentesphases
deforage.
2.2. Paramètresinfluençantsurlesperformancesdel’outildeforage
Figure11:Formationnonconsolidée Figure12:Formationfracturée
11
Chapitre2:Paramètresde forage
Figure13:Formationabrasive
Laduretédesrochesestliéeàl'efforteffectifd’emprisonnementsurlaroche.Enoutre,la
différenceentrelapressiondepuitsetlapressiondeporeestcouranteentantquel’effort
efficaced’emprisonnement. Àmesurequel'effortd'emprisonnementaugmente,l'effortetla
contraintederupturedelarocheaugmentent.L'augmentation delacontraintederupture
conduitàl’augmentationdutravailexigépourramenerlarocheàlarupture.Parconséquent,
quandladifférenceentrelapressiondepuits,etlapressiondeporeaugmentent,leROP diminue[2].
Letypedel’outildeforagechoisietsescaractéristiquesdeconceptionontuneinfluence
significativesurleROP.Nombre,taille,formede dents,laconceptiondesduses,etc. sont des
caractéristiquesd'outilde foragequi affectentleROPetlesperformancesdel’outil[2,5,6,7].
L’étatdel’outil(figures14et15),notammentl’étatdel’usureaunegrandeinfluencesur
l’efficacitéduforagecarl’augmentationdel’usurevadiminuerleROPetlesperformances
del’outil[2,5].
12
Chapitre2:Paramètresde forage
Figure14:OutiluséFigure15:Outilneuf
lefluidedeforageaffectentfortementleROP.Ladensité,lespropriétés
rhéologiquesdel’écoulement,caractéristiquesdefiltrat,distributiondelataille etcontenues de
solidessontcertainespropriétésquiontuneinfluenceélevéesurlesperformancesde l'outil de
forage[2,5,8].Dansle tableausuivant,nousconstatonsl’influencede la naturede la boue sur
lecoefficientdefrottement:
Organigramme:L’Coefficientdefrottementcorrespondantàchaquetypedeboue[Schlumberg
er]
13
Chapitre2:Paramètresde forage
Typiquement,unecourbedeROPenfonctiondeWOB,obtenueexpéri
ent,unecourbedeROPenfonctiondeWOB,obtenueexpérimentale
entalementavectous
lesautresparamètresaurontlafor
ètresaurontlaformecaractéristiquereprésentéedanslafigure(16
16)[5].
AucunROPsignificatifn'estobtenujusqu
estobtenujusqu'àcequeleseuilWOBsoitappliqué.Puis,letaux
ppliqué.Puis,letaux
depénétrationaugmenterapidem
mentavecl'augmentationdeWOBpourdesvaleurs
entationdeWOBpourdesvaleursmodérées
deWOB.Etàdesvaleursplusélevées
sélevéesdeWOB,onobserveseulement
ent delégères
améliorationsdeROP.LepointcorrespondantàROP
éliorationsdeROP.LepointcorrespondantàROPmaximumestappeléflounderpoint.
mestappeléflounderpoint. Dans
lesschistes,le WOBplusélevéqueleflounderpointdi
BplusélevéqueleflounderpointdiminueleROP [6, 9].
Figure16:ROPenfoncti
:ROPenfonctionduWOB(comportementnormal)
Lorsquetoutesautresvaria
ariablesdeforagesonttenuesconstantes,
lesdeforagesonttenuesconstantes, leROPhabituelle
leROPhabituellement
augmente
ente avecdesbassesvaleursdeRPM.AdesvaleursplusélevéesdeRPM,leROP di
diminue[5].
CettefaiblerépondedeROP estdueau mauvaisnettoyagedu puits[5].
LechoixduWOBetduRPMappropriés
BetduRPMappropriés estforte
estfortementinfluencéparlestypesderoches
encéparlestypesderoches
traversées.Parexemple,onforegénérale
ple,onforegénéralement lesterrainstendresavecunWO
OBfaibleetun RPM
élevé[10], tandisquelesterrainsdures sontmieuxforés
sont euxforés avecun WOB élevéetun RPM faible[10].
2.2.6.3. Débit
14
Chapitre2:Paramètresde forage
L’augmentationdudébitestconsidéréecommel’undesparamètressignificatifsinfluençant
surleROP.Letransportdesdéblaisetlenettoyagesontliésétroitement audébitutilisé
pendantleforage(figure17).
Figure17:Déblaisànettoyer
Laplupartdesoutilsdeforagetendentàforerpluslentpendant quelacoursedel’outilde
forageprogresseenraisondel'usurededent.Cetteusureesttrèsimportantedanscertains
terrainsdursetabrasifs,d’oùladiminutiondutauxdepénétrationdecettepartieduforage[5].
2.3. RelationentrelesdonnéesdelasurfaceetlesdonnéesdeFonddeTrou
15
Chapitre2:Paramètresde forage
AfindemesurerleWOBàlasurface,uncalibrage doitêtreeffectuéavantleforageen
ramenantlavaleurinitialedelachargeàzéro. Ilfautaussiinitialiserlavaleur delacharge au
crochet.Ainsi,pendantleforage,ladifférenceentrelavaleurdelachargeaucrochetinitiale
etactuellereprésente leWOB.Maisenréalité,lepoidsappliquésurl’outilestgénéralement
différentauWOB,danscecas,ils’agitduDWOB.LeDWOBesttoujoursinférieureau
WOB.Ceciestduaufrottement etaudragquiagissentsurlagarniture.Lecoupleest
généralementmesuréenutilisantun courantélectriquefourniàlatablederotation.
L'instrumentation conventionnelledutauxdepénétrationnefournitpasunemesurecorrecte
indiquantleprogrèsdel’outil[11].Elle mesuresimplementleprogrèsdumouvementdehaut
enbasdel'extrémitésupérieure delagarnituredeforage[11].Plusieursphénomènes
apparaissentlorsdel’application dupoidstelqueleflambage,danscecas,onnepeutpas comptersur
unelecturefidèledu ROP.
2.4. InfluencedesparamètresdeforagesurleROP
LetransfertdupoidsestunparamètreimportantquipeutdiminuerleROP.Danslespuits
déviés,lesvaleursau-dessusde20000lbdepoidset25000ft-lbdetorque,pourrontêtre
perdusparledragetlesfrottements. Lafigure(18)montreunexempled’unefaible transmission
depoidsdelasurfaceaufonddupuitsdueprobablementàl’accumulationdes cuttings.
16
Chapitre2:Paramètresde forage
Figure18:Exempled’unefaibletransmissiondepoids
Les moteurs de boue de fond utilisés actuellement sont pratiquement tous du type
déplacementpositif.Si∆Pestlapressiondifférentielleàtraverslemoteurdueaupassage
entrelerotoretlestator.LeRPM diminuescommesuit:
RPM= 4.01.0×ΔP
Pour un moteurusé :
RPM= 4.02.0×ΔP
Letorque(T)augmenteavec∆P commesuit:
T=2.0×ΔP
Lapuissancetendvers le maximumsuivantcetteformule:
Moteurneuf:MHP=8.0×ΔP2.0×ΔP 2
17
Chapitre2:Paramètresde forage
Figure19:Caractéristiquesd’unmoteurdeboueàdéplacement
positifdanslesdeuxétats:neufpuisusé
∗ ∗ ∗
FORS= ∗
∗
18
Chapitre2:Paramètresde forage
Cetteéquationdonnelavaleurdeladuretédelaformationquiestindépendante del’étatde
l’usuredel’outiletd’autresparamètres mesurablesdeforage.Lafigure(20)estune schématisation
dutorqueadimensionnel (TD)enfonctiondeladuretéapparentedela formation(FORS).
Afindecompareretcombinerlesdonnéesdedifférentsintervallesetpuitsforés,leROPdoit
êtrenormalisé.Lanormalisationestdéveloppéeenutilisantle modèlesuivant:
=
²
OU :
ROPactuel= taux de pénétration en fb par heurs .
ROPn= taux de pnétration normalisé.
WOB=poids sur l’outil en klbf.
Db= diamétre de l’outil en in .
N= vitesse de rotation en rpm.
2.5. Conclusion
La compréhension des relations entre les différents paramètres de forage, permet d’établir
des concepts théoriques et pratiquespourmieuxmaîtriser l’opérationdeforage.Ces méthodes sont
très utiles pour gagner en terme de temps et matériel (outil de forage par exemple).
19
Chapitre3:Modèlesd’optimisationdanslalittérature
3.1. Introduction
Danslemomentoùl'industriepétrolière essayederéduirelescoûtsetaméliorer la
performancedeforage,lescompagniespétrolièrescherchentdeplusenplusà mieuxmaîtriser lesrisques
d’investissementdeleurcapitald'équipement.
Letempsconsacrédansl’avancementduforagereprésenteuneportionimportantedansle coûtde
forage.DansdespuitstypiquesforéesparBPExploration(BPX),le tempsperdudans
l’avancementdel’outilaétéestiméde10%à30%ducoûtdupuits[12].Celasignifiequela
maîtrisedutauxdepénétration (ROP)aitunegrandeinfluencesurlaréductiondescoûtsde
forage.Eneffet,dansunenvironnementnormal,c’estleROPquiauneinfluencedirectesur
l’efficacitéopérationnelledu forage[12].
3.2. OptimisationduROPenutilisantlathéoriedel’énergiespécifique
20
Chapitre3:Modèlesd’optimisationdanslalittérature
Les spécialistes de l’optimisation des performances et de la conception des outils
comprennentl'effetdeconfinerlapressionsurl’augmentationdeladuretéapparentedelarocheetsonimpac
tsignificatifsurlesperformancesdel’outil.Entantquetels,beaucoupont
développédesalgorithmesdeprévisiondesperformancesbaséssurl'UCS,d'autrespropriétés
deroche,pressiondeboue,pressiondepore,propriétés deboue,etainsidesuite.Quelques
fournisseursd’outilsdeforageontappliquédesméthodes pourcalculeretutiliserune
résistanceàlapressionconfinée. Cependant, touteslesdeuxapproches,ontunouplusieurs
inconvénients, d'êtres de propriétés industrielles et ainsi pas largement disponibles,
incorrectes,demanquer delarobustessepourl'applicationglobalesanscalibragelocal,oude
manquerdelasimplicitéintuitivebaséesur lespremiersprincipes.
L’optimisation desperformancesdeforagepasseparlabonnesélectiondel’outilet
l’optimisationdes paramètresopérationnels.Ceprocessus d'optimisationestsouventaccompli parun
testdevalidation.
Lasélectiondel’outildeforage,prévisiondesperformances,etl'optimisationbaséesurle
MEMestprincipalementlimitéeaux spécialisteschezles compagniesde fournisseursd’outils
etquelquesopérateurs.Danslamesureoùdesprévisionsprécisesdesperformances d’outil
sontdéveloppées,ellessontsouventcomplexes, depropriétésindustrielles,oubaséessurdes
corrélationsempiriqueslocales. Enplus,ellespeuventnepasêtreglobalementapplicablesde
lasurfaceàdegrandesprofondeurs, àtraverslagammedestaillesdetrou,etpourtousles
typesprédominantsd’outils.Lesméthodesdeprévisiondesperformancesdel’outilpeuvent
nepasêtreliéesousensiblesauxcapacitésdel’appareildeforage(puissance disponible et
caractéristiquesdelagarniture).
21
Chapitre3:Modèlesd’optimisationdanslalittérature
tiondanslalittérature
3.2.3. L'environnement du fond de trou
Lafigure(21)illustrel'envirronnementgénéraldufonddutroupourunpuits
entgénéraldufonddutroupourunpuits verticaldansune
rocheporeuse/perméable.Laroche,etladuretéassociéeduterrainderecouvre
éable.Laroche,etladuretéassociéeduterrainderecouvrement,ontété
ent,ontété
remplacéesparuntrouremplidefl
plidefluidedeforage. Bienqu'ilpuisseyavoirdesexceptions,la
ilpuisseyavoirdesexceptions,la
pressionexercée parlefluidedeforage, esten
estengénéralplusgrande
usgrande quelapression
quelapressiondepore,et
inférieureàlapressiondesterrainsderecouvre
insderecouvrement.Danscettecondition
ent.Danscettecondition deforage,laroche
s’étendlégèrement
ent aufonddutrouenraisondelaréductiond effort(lapressiondufluidede
aufonddutrouenraisondelaréductiond'effort(lapressiondufluidede
forageestmoinsquelapression
ression exercéeparlesterrainsderecouvrement).
exe
Figure21:Environnementduf
nementdufonddutroud’unpuitsvertical dansunerocheporeuse/per
eporeuse/perméable
Silarocheestporeuse,ilyaurauneréductionimmédiatedepressiondeporedelaroche
édiatedepressiondeporedelaroche expansée.Si la
rocheestperméable,la réductiondepressionde
tiondepressionde poreacommeconséquencele mouvementdu
mouv fluidedu
champlointain(réservoir)vers
lointain(réservoir)vers larégionexpansée.
Letauxetledegréauxquelslefluidedeporecouledanslarégionexpanséeetégalisela
pressiondeporedelarégionexpa
pressiondeporedelarégionexpansée àcelleduchamplointain(
plointain(pressionderéservoir),est
unefonctiondeplusieursfacteurs.Cesfacteursincluentletauxde
rs.Cesfacteursincluentletauxdel’altérationdelaroche
l’altérationdelaroche
(corrélatifautauxdepénétrationn)etlaperméabilitérelativedelarocheaufluidedep
idedepore
(supposantquelevolumederéservoirestrelati
ederéservoirestrelativementgrandcomparéàlaprofondeurdelarégionforée).En
paréàlaprofondeurdelarégionforée).En
mêmetemps,silapressiondufluidedeforagee
ssiondufluidedeforageestsupérieure àlapression
Pourlarocheimperméabletellequeleschiste,ilyaunmouvementinsignifiantdefluidedu
poreoufiltratd’invasiondanslaprofondeurdelazoneforée.Parconséquent,lapression
instantanéedeporedanslaprofondeurdelazoneforéeestunefonctionduchangement de contraintes,
despropriétésderoche,etdespropriétésinsitudeliquidesdepore (principalementcompressibilité).
Leparamètredel’Esestunemesureutilepourprévoirlapuissancerequise(letorquede
l’outiletlerpm)pouruntypeparticulier d’outilpourforeràunROPdonnédansuntypede rochedonné.
Lathéoried'Esn'estpasnouvelle, elleaétéemployéepourl'évaluationrapidedes
performancesdel’outilpendantdesannées.L’équation(1)montrelarelationdel’énergie
spécifiquedeTeale[13] dérivéepour leforagerotarydans lesconditionsatmosphériques.
∗ ∗ ∗
Es= + (1)
B B∗
Pessier[14]avalidél'équation(1)pourunforagesouslapressionhydrostatique. Puisquela
majoritédedonnéesréellesestsousformedemesuresensurfacetellesque:poidssurl’outil
(WOB),rpm(N),etletauxdepénétration(ROP).Uncoefficientspécifiquedefrottementd’un
= (2)
∗
!
EFFM= ∗ $%% (&)
"#
"**ᴍ 0∗ ² -4%∗5
⎠
SubstituantEsentermesd'efficacité etcouplemécaniquesenfonctionduWOBetrésolvant
24
Chapitre3:Modèlesd’optimisationdanslalittérature
Ilexisteuneméthodelargementpratiquéeetadmiseen mécaniquedesrochespourcalculerla
résistance à la pression confinée de la roche. Elle est comme suit :
CCS=UCS+DP+2DP*sinFA/(1sinFA)(6)
CCS_DP=UCS+DP+2DP*sinFA/(1sinFA)(7)
25
Chapitre3:Modèlesd’optimisationdanslalittérature
Pour le cas de la roche imperméable et d'un puits vertical, une relation décrite par Skempton
[15]aétéappliquéecommeunmoyendecalculer la pression de pore dans la roche expansée dans la
profondeur de la zone forée.Ceci a comme conséquence ce qui suit :
= CCS_Sksi Phie≤0.50
perméableestessentiellementégalàladiffére
entégalàladifférence
ce entre la pression de boue et la pression in situ de
pore, excepté des différences mineuresprovoquées
mineuresprovoquées par le profil inférieur de trou.
Lafigure(22)illustrela diff
fférenceentrelapressiondeboueetlapressiondeporechangée
érenceentrelapressiondeboueetlapressiondeporechangée (dans la
profondeur
ondeur expansée de la zone
zon forée)
ée) pour un ensemble donné de conditions pour la roche
imperméable.
Figure22:Pressiondifférentiellesaufonddutroupouruneroche
essiondifférentiellesaufonddutroupourunerocheimperméable,
méable,comparaisonentrela
méthodesimplifiéede
edeSkempton.etlesrésultatsreportéesparWarrenetSmith
rWarrenetSmith[16]
27
Chapitre3:Modèlesd’optimisationdanslalittérature
tiondanslalittérature
Figure2
Figure23:ROPenfonctiondeCCSpourun micro–outil
de1.25indanslelaboratoire
Pourlesessaisspécifiquesre
esreprésentésdansletableau(2),unelignede
tésdansletableau(2),unelignede tendance
tendancedelaloide
puissanceestassortieetlaformulerésultantee
ulerésultanteest indiquée dans le tableau (2). La formule de
ROPreprésentéedansletableau(
leau(2)estspécifiqueaumicro–outildelaboratoire1.25inet
outildelaboratoire1.25inet aux paramètres
para
de forage (poids surr l’outil, le rpm,
r et le débit).
Tableau2:Valeurscalculééesdelapressiondifférentielle,effortdecompression
ressionconfinée,
ettauxdepénétration
28
Chapitre3:Modèlesd’optimisationdanslalittérature
tiondanslalittérature
Depth=[Feet]
PP =Pore Pressure,[ppg] 0B=Over
erburden,[ppg]
ECD = EquivalentCirculati
tingDensity,[ppg] UCS =UnconfinedCompressiveStrength,[psi]
ressiveStrength,[psi] FA
=FrictionAngle,[degrees]
CCS_DP = ConfinedCompressiveStr
ressiveStrength,[psi],basedonDP
CCS_Sk= ConfinedCompressiveStre
ressiveStrength,[psi], basedon DP_Sk
CCS_ECD =ConfinedCompressiveStren
ressiveStrength,[psi], basedon DP_ECD
P_ECD ROP_DP =
=Rateof
Penetration,[ft/hr],basedon CCS_DP
DP
ROP_Sk= RateofPenetration,[ft/hr],ba
hr],basedonCCS_Sk
ROP_ECD =RateofPenetration,[ft//hr],basedonCCS_ECD
DP =(ECD - porepressure),[psi]
DP_Sk=[ECD- {PP-(OB-ECD)/3}],[p
ECD)/3}],[psi] DP_ECD=ECD Pressure,[psi]
29
Chapitre3:Modèlesd’optimisationdanslalittérature
tiondanslalittérature
Figure
Figure24:LeROPenfonctiondeladensitédeboue
Figure
Figure25:LeROPenfonctiondelapressiondepore
Lafigure(26)illustreuneapplic
neapplicationréelle.Cettefiguremontrelelogdelamécaniquede
écaniquede roche et les
performances
ances du foraged’un puits dans le golf du Mexique.
30
Chapitre3:Modèlesd’optimisationdanslalittérature
tiondanslalittérature
Figure26:Logdela
:Logdela mécaniquedesrochesetdesperformancesde
formancesdeforage
d’unpuitsdanslegolfduMexique
UnexempledelafaçondontleCoefficientSpécifiqueduFrotte
façondontleCoefficientSpécifiqueduFrottementenSlidingd’un
entenSlidingd’un Outil de PDC
avec plus de sept lames
es a été déterminé
déter est illustré dans la figure (27).
31
Chapitre3:Modèlesd’optimisationdanslalittérature
tiondanslalittérature
Figure27:Coefficientspécifiquedufr
:Coefficientspécifiquedufrottementensliding
d’unoutilP
d’unoutilPDCàplusdeseptlames
Leséchantillonsderochesprovenantdugrèsdever
Leséchantillonsderochesprovenantdugrèsdevergerdecrabe,schistedeCatoo
ecrabe,schistedeCatoosa,etdu marbre de
Carthage ont été employés
ployés pour les essais
e multiples avec un outil
til de PDC avec plus de sept lames.
Tous les essais ont employé
ployé un poids de boue de 9.5 ppg. Les valeurs correspondantes de CCS à
6.000 psi de pression de trou inférieurétaientde 18.500 psi pour le
schistedeCatoosa,de36.226psipourle
atoosa,de36.226psipourlemarbre
arbre deCarthage,etde66.000psipourleverger de crabe.
Lacorrélations'estdéveloppéepourcalculer
veloppéepourcalculerμpourunoutildePDCavecplusdeseptla
pourunoutildePDCavecplusdeseptlames, dérivée de la
μ=0.9402*EXP (8E06*CCS
CCS) (13)
Les mêmesprocéduresetessais
etessais completsdesimulateurontétéexécutéspourdéter
completsdesimulateuronté éexécutéspourdéterminerla relation
de μen
en fonction de l’effort de compression
co confinée pour tous les types d’outils.
32
Chapitre3:Modèlesd’optimisationdanslalittérature
tiondanslalittérature
Lafigure(28)illustrelesrapp
)illustrelesrapportsdesefficacitésminimumetmaximumpourl’
mumpourl’outilPDC avec plus de
sept lames.
FFM minimum)etl'efficacitémaximum(maximumEFF
L'efficacitéminimum(EFF aximumEFFM)pourun outil
PDC avec plus de sept lames
es sont calculées
calculé en utilisant les équationss (14 et 15)
1 :
EFFMmax=0.0008*
0.0008*CCS+8.834(14 et 15)
EFFMmax=0.00
0.001*CCS+13.804
Figure28
28:Efficacitémécaniqueminimumet
inimumet maximum
d’unoutilPDCàplusdeseptlames
L’efficacitémécaniquenom
minale(Nom EFFM)est uneefficacitémoyennedérivéedes
efficacitésminimumetmaximum
m.L'équation(16)indiqueNomEFFM pourlesoutilsPDC avec plus de
sept lames.
NomEFFM=0.00095*CCS+10.319(16)
33
Chapitre3:Modèlesd’optimisationdanslalittérature
tiondanslalittérature
Ceuxquicomprennentlesperfor
Ceuxquicomprennentlesperformancesdel’outilsaventquel'efficacitédel’outil
efficacitédel’outildeforage
estaffectéeparlepoidsdelaboue.D’aprèslesessais,lavaleurde pourl’outildePDCest
estaffectéeparlepoidsdelaboue.D’aprèslesessais,lavaleurdeμpourl’outildePDCest
réduiteapproximativementà49%enaug
ativementà49%enaugmentant
entant lepoidsdebouedu9.5ppgà16.5ppg.
lepoidsdebouedu9.5ppgà16.5ppg.En
conséquence,lavaleurdeμdoitêtrecorrigéesilepoidsdeboueestdifférentduppg9.5[18].
doitêtrecorrigéesilepoidsdeboueestdifférentduppg9.5[18]. La figure (29)
(
Figure29:Correctiondufacteurμd’un
d’unPDC
enfonctiondupoidsdelaboue
μ=[0.9402*EXP (8E06*CCS
)]*[0.8876*Ln(MudWeight)+2.998
+ ](17)
34
Chapitre3:Modèlesd’optimisationdanslalittérature
tiondanslalittérature
Figure30:Correctiondufacteurμd’unoutil
d’unoutilPDC
enfonctiondupoidsdelaboue
Leséquations(18et19)montrentlescorrél
ontrentlescorrélationsréviséespourlesefficacitésmécaniques
efficacitésmécaniques Min et la
Max pour un outil PDCavec plus de sept lames
la et pourn'importequel
portequel poids de la boue.
MaxEFFM=(0.0008*CCS+8.834
0.0008* ) *[1.0144*Ln(MudWeight )+3.283
3.2836](18)
3.2.10 .Facteur
Facteur de correction pour les outils PDC
DC dû à la taille de la dent :
Plusieurs essais ont été réalisés pour déterminer l’effet des tailles de dents des outils PDC
sur le ROP (figure 31). Les résultats d'essais ont montré que le coefficient de frottement d’unoutil
en sliding (µ)) est diminué ou augmenté par 1.77% quand la taille de dent est diminuée ou
augmentée chaque millimètre au-dessus
au ou au-dessous de 19 mm.
35
Chapitre3:Modèlesd’optimisationdanslalittérature
tiondanslalittérature
Figure31:Effetdelatailledeladentd’unoutilPDCsurlecoeffic
ffetdelatailledeladentd’unoutilPDCsurlecoeffic
fetdelatailledeladentd’unoutilPDCsurlecoefficientμ
μ=0.0177*CutterSize+0.6637
0. (20)
Encombinanttouslesfacteursde
nttouslesfacteursdecorrection,lacorrélationfinalepourμpourl’outilPDC
pourl’outilPDC avec plus
μ=[0.9402*EXP(8E06*CCS)]*[0.8876*Ln(MudWeight
MudWeigh )+2.998]*
[0.0177*CutterSize+0.6637](21)
36
Chapitre3:Modèlesd’optimisationdanslalittérature
Dans un scénario d'optimisation d’outil en temps réel, le ROP prévu et les valeurs d'énergie
Es peuvent être employés pour évaluer les performances de l’outil. Ceci peut être accompli si les
propriétés de roche sont connues, par corrélation ou directement mesurées et calculées à partir
des données de LWD. Les performances et les conditions peuvent être évalués en comparant le
Es réelle à Es prévue, aussi bien que par la comparaison du ROP réel au ROP prévu. L'analyse
des performances des outils en utilisant le Es prévue en temps réel et les valeurs réelles de Es
peuvent être également employées pour détecter et corriger des problèmes de forage, tels que la
vibration de l’outil et les chocs sur l’outil, analyse de la rupture, et les outils inactifs.
3.3. Conclusions
Plusieurs méthodes ont été développées pour optimiser les paramètres de forage. Bien
entendu, le détail de la plupart de ces méthodes reste confidentiel, car elle restent la
propriété industrielle de la société qui a investi sur cette recherche. Néanmoins, les éléments de
base de chaque recherche sont publiés dans les revues internationales.
37
Chapitre 4 : Conception et structures du modèle à optimiser
d’optimisation
4.1. Introduction
Dans le chapitre précédent, nous avons présenté trois méthodes d’optimisation des paramètres
de forage. Pour notre étude, nous avons choisi la méthode d’optimisation du ROP par la
théorie de l’énergie spécifique. Cette méthode utilise directement les paramètres mesurés en surface
dans la cabine Mud logging.
Dans la partie verticale, on peut utiliser directement la formule (5), car l’auteur [14] a
considéré les conditions dans un puits vertical.
. ∗ ∗
= (5)
−
ᴍ∗
Avec :
ROP : Taux de pénétration.
= (2)
∗
38
Chapitre 4 : Conception et structures du modèle à optimiser
d’optimisation
Pour éviter le calcul de CCS et EFFM, on va les déterminer à partir des premières données
disponibles lors du forage. A partir de l’équation (5), on a :
. ∗ ∗
=% + *∗ (33)
ᴍ ∗& ' )
On pose : ∝=
ᴍ
Donc:
. ∗ ∗
∝= % + *∗ (34)
∗& ' )
Les premiers mètres forés peuvent nous données les valeurs suivantes : µ, N, ROP et le WOB.
Dans ce cas, on peut calculer facilement ∝. La valeur de ∝ sera utilisée par la suite pour
calculer le ROP optimal, en utilisant l’équation (5). A son tour, la valeur de ROP optimal sera
En effet, les paramètres de forage sont limités par des intervalles donnés par les constructeurs.
Par exemple, le WOB est limité par la capacité de l’outil de forage. Ainsi, on obtient le
modèle d’optimisation suivant :
. ∗ ∗
• Fonction objectife : = (/)
ᴍ∗
.
NMIN≤N≤Nmax
39
Chapitre 4 : Conception et structures du modèle à optimiser
d’optimisation
L’organigramme illustré dans la figure (32), montre le mécanisme de calcul pour obtenir le ROP
optimal.
Valeurs ROP,WOB
N ,TOR
N,TOR
Calcul α
Calcul ROP
Si ROP=ROPn+1 Si ROP≠ROPn+1
Calcul WOB,N
optimal
ROP réel
WOB,N optimal
Figure 32 : Processus simplifié de calcul pour l’obtention du ROP optimal, cas d’un forage vertical
40
Chapitre 4 : Conception et structures du modèle à’optimisation
Cependant la détermination des paramètres de forage sera moins coûteuses et plus rapide mais
l’information obtenue est moins fiable que celle obtenue expérimentalement.
Dans ce chapitre, nous présentons le programme informatique élaboré pour l’optimisation des
paramètres de forage. Ce programme se compose de deux parties principales :
41
Chapitre 4 : Conception et structures du modèle à’optimisation
NV=Nmin:PN:Nmax;
m(1)=(36*TOR)/(D*WOBO);
alpha(1)=(((13.33*m(1)*NO)+(1/AB))/(D*ROPO))*WOBO;
ROP(1)=(13.33*m(1)*NO)/(((alpha(1)/WOBO)-(1/AB))*D);
j=2;PX=2;epsi=1;
while PX>epsi
m(j)=(36*TOR)/(D*WOBV(j-1));
alpha(j)=(((13.33*m(j)*NV(j-1))/(1/AB))+(D*ROP(j-1)))*WOBV(j-1);
ROP(j)=(13.33*m(j)*NV(j-1))/((alpha(j)/WOBV(j-1))-(1/AB)*D);
epsi=ROP(j)*0.05;
PX=ROP(j)-ROP(j-1);
j=j+1;
end ;
ROPREAL=ROP(j);
wobOPT=WOB(j-1);
Nopt=N(j-1);
input('les valeurs ROPRAL=',ROPREAL);
input('les valeyrs WOBopt=',WOBopt,'les valeyrs Nopt=',Nopt);
Les Tableau suivants présentes les données de calcul des paramètres de forage
Bit Time
Bit Duse TFA Depth In Depth Out ROP WOB RPM LPM SPP
Type
N°/ Diam 2 (m) (m) (m/h) (t) (tr/mn) (l/mn) (psi)
(in ) (h)
01 / 26’’ GTXC3G1 6*18 00 491 36.37 13.5 15-27 110 3900 1900
Bit Time
Bit Duse TFA Depth In Depth Out ROP WOB RPM LPM SPP
Type
N°/ Diam 2 (m) (m) (m/h) (t) (tr/mn) (l/mn) (psi)
(in ) (h)
42
Chapitre 4 : Conception et structures du modèle à’optimisation
MKS-
05/ 12¼ 7*16/32 3080 3086 5.04 1.19 1-5 115 1600 2900
57SDG2
06/ 12¼ HCM408 8*15/32 3086 3284 198 2.51 4-8 115 1620 2900
HCM408
07/ 12¼ 8*15/32 3284 3301 12.54 1.36 8-12 130 1700 3100
7004676
Bit Tim
Bit Duse TFA Depth In Depth Out ROP WOB RPM LPM SPP
Type
N°/ Diam 2 (m) (m) (m/h) (t) (tr/mn) (l/mn) (psi)
(in ) (h)
8 / 8⅜ SE3841I 4*16 3301 3431 46.07 2.82 6-10 200-210 1600 2600
1
" = 36 23∗456
43
Chapitre 4 : Conception et structures du modèle à’optimisation
µ0=36*(362.6433)/(0.31115*15*103) µ0=2.80
78.88∗µ∗9 7
∝= + ∗ WOB
2∗:;< =3
78.88∗µ∗9
ABC DDE K
23 .
FGGᴍ∗HIJ LJ
ES=(WOB)/(AB +120*π*N*T)/(AB*ROP)
ES1=(28660.09)/(290.1606)+(120*3.14*160*2623)/(290.1606*44.901) ES1=12242.51486 psi
Esmin=(WOBmin)/(AB +120*π*Nmin*T)/(AB*ROPmin)
Esmin1=(33069.34 /290.1606) +((120* π*110*2623) / (290.1606*44.29134))
Esmin1=8577.758212psi
On présentera dans cette partie les valeurs des paramètres de forage optimisé par la méthode
d’énergie spécifique.
Après la lecture des données concernant les paramètres de forage et de la géométrie de puits
(well bore), ainsi que les paramètres initial, on fait appel au programme OPTIFOR pour le calcul
des paramètres optimaux.
Le tableau (7) suivant représenté les paramètres de forage optimisées par la méthode
d’énergie spécifique dans le cas verticale de 3 phase :
44
Chapitre 4 : Conception et structures du modèle à’optimisation
Tableau (7) : les paramètres de forage optimisées par la méthode d’énergie spécifique dans le cas
verticale
Partie horizontal :
RPM = 4.0 −1.0 × ∆P ∆P= 3bar
RPM=4.0-1.0× 3 RPM=1
Mouteur neuf :
MHP=8× ∆P -2 × ∆P2=8 ×3- 2 ×32 MHP=6
ABCactuel
ABCn = = 44.90125/(13/12.24996)² ∗ (160/60)
TBU Y
²
VW 60
ROPn=14.9510 ft/hr
On présentera dans cette partie les valeurs des paramètres de forage optimisé par la méthode
d’énergie spécifique.
Après la lecture des données concernant les paramètres de forage et de la géométrie de puits
(well bore), ainsi que les paramètres initial, on fait appel au programme OPTIFOR HORIZONTAL
pour le calcul des paramètres optimaux.
Le tableau (8) suivant représenté les paramètres de forage optimisées par la méthode
45
Chapitre 4 : Conception et structures du modèle à’optimisation
4.5 Conclusion
Le développement du code de calculs, a nécessité le travail sur deux parties. La première est
liée à la fonction objectife et ses contraintes, et la seconde concerne le modèle mathématique à
appliquer.
Pour la partie verticale du forage, il n’y avait pas de problème puisque le modèle
d’optimisation peut être appliqué directement. En revanche, la partie horizontale, a nécessité des
adaptations sur le modèle précèdent.
46
Conclusion générale
Conclusion générale
Pour consolider la position de l’ENTP dans le marché du forage algérien, et faire des pas surs
dans le marché international, il faudra travailler sur plusieurs fronts. L’un de ces fronts, est la
maîtrise et la minimisation du temps alloué au forage.
L’un des défis du forage en Algérie, est la vérification des plafonds de production dressés par
l’état. L’objectif tracé pour l’année 2012 est d’atteindre 2millions barils par jours, et 2.5 millions
barils par jour en 2016.
Par notre modeste travail, nous essayons de contribuer à optimiser le temps de forage, et
ce, par l’optimisation des paramètres de forage, en construisant un logiciel destiné à ce type de
calculs.
49
Références bibliographiques
Références Bibliographiques
50
Références bibliographiques
51
Références bibliographiques
52