FR 2989200 A1 20131011 FR 1201057 A 20120410 La présente invention concerne le domaine technique de l'industrie pétrolière, et plus particulièrement l'optimisation de placement des puits, ou zones de drainage, dans un réservoir pétrolier pour produire des hydrocarbures. Pour produire un gisement pétrolier, les ingénieurs réservoirs cherchent à définir les meilleurs positions et trajectoires des puits à forer correspondant aux meilleures zones de drainages. Ces puits sont généralement des puits producteurs qui permettent de faire remonter les hydrocarbures depuis le réservoir jusqu'à la surface, ou parfois des puits injecteurs, par exemple d'eau ou de gaz, pour améliorer la récupération des hydrocarbures contenus dans le réservoir. Le forage des puits est une décision importante dans la vie d'un gisement pétrolier, vu que forer un puits est une action couteuse à la fois en argent et en temps. Pour cette raison, les ingénieurs réservoir utilisent généralement deux outils. - un premier outil appelé "modèle de réservoir" qui constitue une image représentative de la structure du réservoir ainsi que de son comportement. Cette image ou représentation est construite à partir des données collectés sur le champ pétrolier considéré (carottes, diagraphies, sismiques,...) réalisées par des géologues, des géophysiciens, des pétrophysiciens et un ordinateur. Cependant, en pratique, les géologues, les géophysiciens et les pétrophysiciens ne possèdent pas suffisamment des données pour construire un modèle qui correspond parfaitement au sous-sol. Donc, l'approche usuelle bien connue par les spécialistes est de construire plusieurs modèles de réservoir représentant des propriétés, par exemple pétrophysiques, différentes. Chaque modèle de réservoir est considéré comme une image possible de la structure du réservoir et de son comportement. L'incertitude sur la connaissance précise du réservoir est donc prise en compte en produisant plusieurs images possibles du réservoir. - un deuxième outil appelé "simulateur d'écoulement". Un simulateur d'écoulement est un logiciel qui a pour objectif de modéliser les écoulements des fluides au sein du réservoir pétrolier représenté par un modèle de réservoir. Par exemple, ce logiciel permet d'estimer pour une configuration de puits donnée et pour un intervalle de temps donné, les quantités produites d'eau, d'huile et de gaz. PumaFlow® (IFP Energies nouvelles, France) est un exemple de simulateur d'écoulement. On connait des techniques permettant d'optimiser l'emplacement des zones de drainage en utilisant plusieurs modèles de réservoir possibles et un simulateur d'écoulement. On connait par exemple les méthodes suivantes: - Guyaguler, B. and Home, R.N. 2001. Uncertainty assessment of well placement optimization. In: SPE annual technical conference and exhibition. SPE 71625 - Ozdogan, U. and Home, R.N. 2006. Optimization of Well Placement Under TimeDependent Uncertainty. SPE Res Eval & Eng 9 (2): 135-145. SPE-90091-PA - Alhuthali, A.H., Datta-Gupta, A., Yuen, B. and Fontanilla, J.P. 2010. Optimizing smart well controls under geologic uncertainty. Journal of Petroleum Science and Engineering 73 (1-2): 107-121. - Schulze-Riegert, R., Bagheri, M., Krosche, M., Kuck, N. and Ma, D. 2011. Multiple- objective optimization applied to well path design under geological uncertainty. In SPE reservoir simulation symposium. SPE 141712 Cependant, ces approches représentent un processus long et fastidieux nécessitant un nombre très élevé de simulations d'écoulement au moyen du simulateur d'écoulement. En effet, pour chaque configuration des zones de drainage considérée, ces méthodes effectuent une simulation d'écoulement pour chaque modèle de réservoir possible, pour les combiner après en utilisant un critère de qualité, classiquement, défini comme étant un paramètre appelé « Net Present Value », ou NPV. Le NPV est la différence des flux de trésorerie généré par l'investissement correspondant au placement des zones de drainage. Pour répondre aux besoins des ingénieurs et des spécialistes en charge de l'exploitation du réservoir pétrolier, il est essentiel de réduire le nombre de simulations d'écoulements pour optimiser le placement des zones de drainage sur plusieurs modèles de réservoir possibles. Pour réduire le nombre de simulations d'écoulement, on connaît la méthode décrite dans le document suivant : - Wang, H., Echeverria Ciaurri, D., Durlofsky, L.J. and Cominelli, A. 2011. Optimal well placement under uncertainty using a Retrospective Optimization framework. In: SPE Reservoir Simulation Symposium. SPE 141950-MS Cette méthode propose d'utiliser un nombre limité (choisi par l'ingénieur ou l'utilisateur de la méthode) de modèles de réservoir au début de l'optimisation, puis d'utiliser de plus en plus (d'une façon choisi par l'ingénieur ou l'utilisateur de l'approche) de modèles en fonction de l'évolution de l'optimisation. En pratique, il est difficile pour l'ingénieur de choisir ce nombre de modèles à utiliser à chaque étape ou itération de la méthode. De plus, au cours des dernières générations, on considère pour chaque configuration de zones de drainage, un nombre élevé (21 ou 104 dans l'exemple proposé) de modèles de réservoir possibles, et donc de faire un nombre élevé de simulations de réservoir pour chaque configuration des zones de drainage. En pratique, il est très difficile aussi d'effectuer ce nombre élevé de simulations d'écoulement. Ainsi, l'objectif de l'invention concerne un procédé pour optimiser le placement des puits à partir de plusieurs modèles de réservoir possibles, dans laquelle le nombre de simulations d'écoulement est réduit. Le procédé se base sur l'utilisation des résultats de simulations d'écoulement effectués tout au long de la procédure d'optimisation, et ainsi défini un nouveau critère de qualité qui approche le NPV mais nécessitant un nombre réduit (voir une seule) de simulations d'écoulement. De façon générale, l'invention concerne un procédé d'exploitation d'un réservoir souterrain, à partir d'un ensemble d'images du réservoir, dans lequel on définit un nombre de puits à forer au sein du réservoir. On détermine des emplacements des puits en réalisant les étapes suivantes : i) on génère aléatoirement des configurations de puits, chaque configuration correspondant à des emplacements possibles des puits ; ii) on optimise lesdites configurations en calculant un premier critère de qualité de l'exploitation du réservoir au moyen d'un simulateur d'écoulement appliqué à un nombre restreint d'image du réservoir et au moyen de résultats de simulations de configuration obtenus préalablement pour des configurations proches, par rapport à une distance, de la configuration à évaluer, ladite distance étant fonction des coordonnées des puits ; iii) on calcule pour lesdites configurations optimisées, un second critère de qualité de l'exploitation du réservoir au moyen d'un simulateur d'écoulement appliqué à chaque image du réservoir, et l'on choisit l'emplacement des puits correspondant à la configuration optimisée ayant le meilleur second critère de qualité; iv) on fore les puits aux emplacements choisis. Selon l'invention, les résultats de simulations de configuration obtenus préalablement peuvent être obtenus au cours de l'optimisation, ou en calculant le second critère de qualité. Selon l'invention, les premier et second critères de qualité peuvent correspondre à une valeur présente nette NPV caractéristique de toutes les images du réservoir. Le second critère C peut alors être calculé par la formule suivante : C = - 1 NR I avec : le NPV issu d'une simulation associée à l'image de réservoir i ; NR : le nombre d'images du réservoir. Le premier critère Q peut alors être calculé par la formule suivante : piCi j=1 i=1 Avec : NS : nombre de configurations déjà simulées NM : nombre d'images de réservoir considérées, NM NR , NR étant le nombre d'images du réservoir. : résultat de simulation correspondant à chacune des Ns configurations. : résultat de simulation correspondant à la configuration en cours, pour chacune des NM images de réservoir considérées. : poids affecté à chaque Ci 23i poids affecté à chaque é; . Selon l'invention, on peut réitérer le procédé en faisant varier le nombre de puits, et on peut optimiser les configurations au moyen d'un algorithme de type CMA-ES. Enfin, les puits à forer peuvent comporter des drains multi-branches. D'autres caractéristiques et avantages du procédé selon l'invention, apparaîtront à la lecture de la description ci-après d'exemples non limitatifs de réalisations, en se référant aux figures annexées et décrites ci-après. Présentation sommaire des figures La figure 1 montre le maillage du réservoir utilisé. Elle comporte une vue en perspective du réservoir et une vue de dessus, et présente l'élévation verticale du réservoir. La figure 2 montre l'évolution de la meilleure valeur de critère de qualité approché trouvé en utilisant la méthode selon l'invention. L'axe des abscisses présente le nombre de simulations de réservoir utilisées. L'axe des ordonnées présente la valeur du critère de qualité approché. La figure 3 montre l'évolution de la meilleure valeur de vrai critère de qualité trouvé en utilisant la méthode selon l'invention. L'axe des abscisses présente le nombre de simulations de réservoir utilisées. L'axe des ordonnées présente la valeur du vrai critère de qualité. La figure 4 montre l'évolution de la meilleure valeur du critère de qualité approché trouvé en utilisant la méthode selon l'invention. Les points correspondent à chaque évaluation du vrai critère de qualité. L'axe des abscisses présente le nombre de simulations de réservoir utilisées. L'axe des ordonnées présente la valeur du critère de qualité (approché ou vrai). La figure 5 montre l'évolution de la meilleure valeur de NPV trouvé en utilisant la méthode selon l'invention et la compare avec celle trouvée en utilisant l'approche conventionnelle. L'axe des abscisses présente le nombre de simulations de réservoir utilisées. L'axe des ordonnées présente la valeur du NPV. La figure 6 montre une vue de dessus du maillage du réservoir. Elle montre aussi les emplacements des puits à forer obtenus en utilisant la méthode selon l'approche conventionnelle: le producteur, noté P1, et l'injecteur, noté 11. La figure 7 montre une vue de dessus du maillage du réservoir. Elle montre aussi les emplacements des puits à forer obtenus en utilisant la méthode selon l'invention: le producteur, noté P1, et l'injecteur, noté 11. Description du procédé L'objet de l'invention concerne une technique de décision de placement des puits dans un réservoir, en prenant en compte l'incertitude sur la connaissance du réservoir avec un nombre réduit de simulations de réservoir. On appelle « puits » toute zone de drainage. La méthode utilise une optimisation d'un critère de qualité approché permettant d'obtenir un résultat satisfaisant du vrai critère de qualité (moyenne des NPV ou des productions des zones de drainage sur tous les modèles de réservoir). La réduction du nombre de simulations de réservoir permet l'application industrielle par l'ingénieur réservoir. A partir d'un ensemble d'images du réservoir (modèles de réservoir), le procédé d'exploitation d'un réservoir souterrain selon l'invention comporte les étapes principales suivantes : 1. on définit un nombre de puits à forer au sein du réservoir 2. on génère aléatoirement des configurations de puits 3. on optimise lesdites configurations en calculant un premier critère de qualité (critère approché) 4. on calcule pour les configurations optimisées, un second critère de qualité (vrai critère) 5. on fore les puits aux emplacements choisis. On commence par définir un ensemble de modèles de réservoir (d'images du réservoir) qui correspondent chacun à une grille définie par les propriétés du réservoir comme la porosité, la perméabilité... Ces modèles sont construits à partir de mesures collectées du réservoir à l'aide d'une méthode géostatistique. On choisit également un simulateur d'écoulement. On peut utiliser par exemple le logiciel PumaFlowTM (IFP Energies nouvelles, France). 1. Définition du nombre de puits à forer au sein du réservoir On définit dans cette étape le nombre de puits à forer au sein du réservoir. Ces puits définissent des zones de drainage. On définit dans cette étape le type de chaque puits : producteur ou injecteur. Pour les puits injecteurs, on définit aussi le type du fluide à injecter (eau, gaz,...). 2. Génération aléatoire des configurations de puits Une configuration de puits correspond à l'emplacement, au sein du réservoir, de l'ensemble des puits à forer dans le réservoir. On définit un paramètre À, dit taille de population, caractérisant le nombre de configurations de puits utilisées à chaque génération. La population initiale (première itération) est générée aléatoirement. La configuration de puits est représentée par la localisation de chaque puits au sein du réservoir. Il peut s'agir des coordonnées des points définissants les extrémités, au sein du réservoir, des puits et de ses drains latéraux. La dimension du problème est ainsi égale au nombre de variables utilisées pour représenter la configuration de puits. 3. Optimisation des configurations en calculant un critère de qualité approché 3.1 Choix d'un critère de qualité approché de l'exploitation du réservoir La qualité de l'exploitation de réservoir est classiquement, défini comme étant un paramètre appelé "Net Present Value", ou NPV. Le NPV est la différence des flux de trésorerie généré par l'investissement correspondant au placement des zones de drainage. Dans cette étape, on définit un critère de qualité de l'exploitation du réservoir remplaçant le critère de qualité vrai de l'étape 4.1. La qualité approchée de l'exploitation du réservoir est calculée en utilisant un nombre réduit de nouvelles simulations d'écoulement et un nombre de résultats de simulations déjà existantes avec d'autres configurations de zones de drainages qui sont proches de la configuration considérée. Ainsi, on calcule un premier critère de qualité de l'exploitation du réservoir au moyen d'un simulateur d'écoulement appliqué à un nombre restreint d'image du réservoir et au moyen de résultats de simulations de configuration obtenus préalablement pour des configurations proches, par rapport à une distance, de la configuration à évaluer, ladite distance étant fonction des coordonnées des puits. Les résultats de simulations de configuration obtenus préalablement sont obtenus au cours de l'optimisation ou obtenus en calculant un second critère de qualité (voir étape 4.1). En effet, pour une configuration de zones de drainage donnée, on effectue les étapes suivantes : 3.1.1 on sélectionne un nombre NN ("Number of neighbor points") de configurations de zone de drainage les plus proches, au sens d'une distance, de la configuration de zones de drainage considéré. Ces NN configurations sont sélectionnées parmi les configurations de zones de drainage déjà simulées avec un simulateur d'écoulement. On note par Ns le nombre de configurations de zones de drainage déjà simulé pour lesquelles on dispose les résultats de simulation d'écoulement. NN doit être inférieur ou égal à Ns. On note par c, les résultats de simulation correspondants, à chaque configuration de zones de drainage sélectionnée. Typiquement NN est égale à min(2*NR, Ns); 3.1.2 on sélectionne aléatoirement un nombre NM de modèles de réservoir, sur chacune on effectue une simulation d'écoulement de la configuration de zone de drainage considéré. On note par Ci les résultats de simulation correspondants, à la configuration de zone de drainage considérée, à chaque modèle de réservoir sélectionné. Typiquement, NM est égale à 1. Ceci signifie qu'on effectue une seule simulation d'écoulement; 3.1.3 on choisit une fonction d'attribution de poids aux différents critères de qualités, i.e., on associe une valeur de poids à chaque résultat de simulation ê, ou C. Les valeurs de poids sont normalisées entre 0 et 1. Typiquement, les valeurs des poids correspondants aux Ci sont plus importantes que celles des poids correspondants aux C; . Typiquement aussi, les valeurs des poids correspondants aux C; sont égales. Typiquement enfin, les valeurs des poids correspondants aux ê, , doivent être décroissante en fonction de la distance de la configuration sélectionnée de la configuration considérée. On note dans ce qui suit par pi les NM valeurs de poids correspondants à chaque C; et par p; les NS valeurs de poids correspondants à chaque . En particulier, un moyen de définir les NM valeurs de poids pi et les NS valeurs de p, est de définir une distance maximale dite de sélection égal à h (par exemple égale à 4000). Ensuite, pour chaque /37 , on note par di la distance entre la configuration sélectionnée et la configuration considérée. On peut définir donc pi et p, de la façon suivante: On peut également utiliser une autre forme de qualité comme par exemple utiliser la différence entre la moyenne des NPV sur tous les modèles de réservoir et l'écart type des NPV sur tous les modèles réservoir. D'une façon générale, le critère de qualité de l'exploitation de réservoir, noté C, est évalué comme suit: C avec Ci le résultat de simulation associé au modèle de réservoir i. Typiquement, C correspond au NPV, et chaque C, correspond à un NPV évalué avec un modèle de réservoir donné. Le calcul d'une valeur de Ci nécessite une simulation d'écoulement. D'une façon générale, le calcul du critère de qualité de l'exploitation du réservoir, permettant de mesurer la qualité d'une configuration de zones de drainages donnée, nécessite donc un nombre de simulations d'écoulements égal au nombre de modèles de réservoir considéré qui est égal à NR. 4.2 Choix des meilleurs emplacements de chaque puits Dans cette étape on commence par choisir un nombre d'emplacements de chaque zone de drainage candidats à être parmi les meilleurs emplacements parcourus pendant l'étape précédente. Un moyen est de choisir le meilleur emplacement obtenu à chaque itération. Un autre moyen peut être par exemple de choisir seulement le meilleur emplacement obtenu à chaque génération s'il représente le meilleur par rapport à tous les points parcourus (avant cette génération) durant toute l'étape précédente. Une fois que l'on a choisi des emplacements candidats à être les meilleurs emplacements, on effectue pour chaque emplacement choisi des simulations d'écoulements avec tous les modèles de réservoir, et on calcule le critère de qualité comme défini à l'étape 4.1. Puis, on choisit un des emplacements défini à l'étape précédente. Généralement, il s'agit de prendre la configuration des zones de drainage qui offre le meilleur critère de qualité et qu'il est techniquement possible de forer. 5. Forage des puits aux emplacements choisis. Dans cette étape, on fore les puits définis dans l'étape précédente. Exemple d'application La méthode de placement des puits dans un réservoir incertain selon l'invention peut être utilisée pour placer de nouveaux puits ou ajouter des nouveaux latéraux à des puits déjà existants. Elle peut être appliquée sur des nouveaux champs ou des champs matures (contenant des puits déjà forés). L'incertitude du réservoir est définie moyennant un nombre de modèles de réservoir équiprobables. Ici, nous présentons en particulier un exemple d'application pour le placement d'un puits producteur et un puits injecteur dans un nouveau champ avec 20 modèles de réservoir possibles. Cette méthode est appliquée sur un réservoir synthétique. Le réservoir est de taille 3420m x 5040m x 90m. Le maillage est cartésien avec 19 mailles dans la direction x, 28 mailles dans la direction y et 5 mailles dans la direction z. La taille des mailles est de 180m x 180m x 18m. La figure 1 montre le maillage du réservoir utilisé. Elle comporte une vue en perspective du réservoir et une vue de dessus, et présente l'élévation verticale du réservoir. Le champ considéré ne contient aucun puits déjà foré. On propose de trouver le meilleur placement pour chacun des deux puits avec uniquement un tronc principal (sans latéraux), avec la méthode selon l'invention. Le vrai critère de qualité C est choisi comme étant la valeur moyenne de NPV sur les 20 modèles de réservoir: 20f Y C = E 1 r 20 APRy x C" + Qn,g X Cn,g + Qn,wa x C no,nt - Cd \ n=1 ii int où Qn,p est la production du champ considéré en phase p dans la période n sur le modèle de réservoir m, Cn,p est le gain ou la perte associé à la production de la phase p pendant la période n, la phase p peut représenter l'huile, le gaz ou l'eau qui sont notés respectivement par o, g, wa. APR est le taux d'intérêt annuel. Y est le nombre de périodes considérées, Cd est le cout nécessaire pour le forage et la complétion pour les puits considérés, Cd est approximé comme suit : Cd = k=0 avec lw la longueur du puits, dw le diamètre du puits, A une constante, N le nombre total de puits. Les productions Qn,p de la phase p pour la période n sur le modèle de réservoir m sont obtenus moyennant le simulateur de réservoir. Dans cet exemple on utilise le simulateur de réservoir PumaFlowTM (IFP, France). Dans cet exemple, N = 2. Les constantes utilisées pour définir la fonction objectif sont représentées dans le tableau suivant : Constante Valeur Cn,, 60 $ / baril Cn,g -4 $ /baril Cn,w 0 APR 0.2 A 1000 dw 0.1 m Les pressions limites sont imposées au fond de ces deux puits. On utilise l'algorithme d'optimisation CMA-ES pour optimiser les 12 paramètres représentant les 2 puits (6 paramètres par puits correspondant aux coordonnées des extrémités de chaque puits). On choisit une taille de population égale à 40 individus. On définit un critère de qualité approché ê . Pour calculer ê pour une configuration de zone de drainage donnée, on commence par sélectionner un nombre NS de configurations de zone de drainage déjà simulées. Ce nombre est égal au maximum 40, et correspond au nombre de configurations de zone de drainage avec une distance au maximum égale à 4000. La distance entre deux configurations de puits q1 et q2 est calculée comme suit: d(91,92) = -'(9'1- 42)T C (91- 42) V.71,q2 e 9? - On note par i = 1,...,Ns les valeurs de NPV des NS configurations de zones de drainage sélectionnées et par d, les distance entre chacun des NS configurations de zones de drainage sélectionnées et la configuration de zone de drainage à évaluer. On effectue après une seule simulation de réservoir pour la configuration de zone de drainage à évaluer en utilisant un modèle de réservoir sélectionné aléatoirement, et on note C1 la valeur de NPV qui correspond à cette simulation. Le critère de qualité approché est ainsi calculé comme étant: = ± avec : pl et /-3, calculé et normalisé comme suit: Pi ÷- 1, L'Y \z On effectue 20 simulations de réservoir pour chaque configuration de zone de drainage sélectionnée, i.e., pour chaque modèle de réservoir. Ainsi, on effectue 600 (= 20 x 30) simulations de réservoir. Les résultats des 30 évaluations de vrai critère de qualité pour les configurations de zones de drainage sélectionnées sont reportés à la figure 3. Cette figure montre l'évolution de la meilleure valeur de vrai critère de qualité trouvé en utilisant la méthode selon l'invention. L'axe des abscisses présente le nombre de simulations de réservoir utilisées. L'axe des ordonnées présente la valeur du vrai critère de qualité. Chaque point de la courbe correspond à 20 simulations de réservoir effectuées. Les résultats des 30 évaluations de vrai critère de qualité sont montrées dans la figure 4, ensemble avec les résultats des évaluations du critère approché. Cette figure montre l'évolution de la meilleure valeur du critère de qualité approché trouvé en utilisant la méthode selon l'invention. Les points correspondent à chaque évaluation du vrai critère de qualité. L'axe des abscisses présente le nombre de simulations de réservoir utilisées. L'axe des ordonnées présente la valeur du critère de qualité (approché ou vrai). La méthode selon l'invention permet donc d'avoir une meilleure valeur du vrai critère de qualité (moyenne de NPV sur tous les modèles de réservoir) égale à 9.66E9. Avec cette méthode, nous pouvons obtenir une meilleure configuration de zone drainage avec un nombre techniquement abordable de simulations de réservoir. On montre aussi que la méthode selon l'invention est presque 16 fois plus rapide (16 fois de moins en nombre de simulation) que la méthode conventionnelle. La figure 5 montre l'évolution de la meilleure valeur de NPV trouvé en utilisant la méthode selon l'invention et la compare avec celle trouvée en utilisant l'approche conventionnelle. L'axe des abscisses présente le nombre de simulations de réservoir utilisées. L'axe des ordonnées présente la valeur du NPV. Dans cet exemple d'application, la méthode conventionnelle est basée sur une optimisation basée sur la méthode CMA-ES du vrai critère de qualité. La figure 6 montre les positions des configurations de zones de drainage obtenues en utilisant la méthode selon la méthode conventionnelle. La figure 7 montre les positions des configurations de zones de drainage obtenues en utilisant la méthode selon l'invention. Les positions de zones de drainage sont proches et les valeurs de critères de qualités sont équivalentes, mais nous avons besoin de beaucoup moins de temps pour déterminer la zone de drainage selon l'invention.20 REVENDICATIONS1. Procédé d'exploitation d'un réservoir souterrain, à partir d'un ensemble d'images du réservoir, dans lequel on définit un nombre de puits à forer au sein du réservoir, caractérisé en ce que l'on détermine des emplacements desdits puits en réalisant les étapes suivantes : v) on génère aléatoirement des configurations de puits, chaque configuration correspondant à des emplacements possibles des puits ; vi) on optimise lesdites configurations en calculant un premier critère de qualité de l'exploitation du réservoir au moyen d'un simulateur d'écoulement appliqué à un nombre restreint d'image du réservoir et au moyen de résultats de simulations de configuration obtenus préalablement pour des configurations proches, par rapport à une distance, de la configuration à évaluer, ladite distance étant fonction des coordonnées des puits ; vii) on calcule pour lesdites configurations optimisées, un second critère de qualité de l'exploitation du réservoir au moyen d'un simulateur d'écoulement appliqué à chaque image du réservoir, et l'on choisit l'emplacement des puits correspondant à la configuration optimisée ayant le meilleur second critère de qualité ; viii) on fore les puits aux emplacements choisis. 2. Procédé selon la revendication 1, dans lequel les résultats de simulations de configuration obtenus préalablement sont obtenus au cours de l'optimisation. 3. Procédé selon la revendication 1, dans lequel les résultats de simulations de configuration obtenus préalablement sont obtenus en calculant ledit second critère de qualité. 4. Procédé selon l'une des revendications précédentes, dans lequel les premier et second critères de qualité correspondent à une valeur présente nette NPV caractéristique de toutes les images du réservoir.17 5. Procédé selon la revendication 4, dans lequel le second critère C est calculé par la formule suivante : avec : le NPV issu d'une simulation associée à l'image de réservoir i ; NR : le nombre d'images du réservoir. 6. Procédé selon l'une des revendications 4 et 5, dans lequel le premier critère Q est calculé par la formule suivante : 10 Avec : Ns : nombre de configurations déjà simulées NM : nombre d'images de réservoir considérées, NM 5 NR , NR étant le nombre d'images du réservoir. 15 C résultat de simulation correspondant à chacune des Als configurations. Ci résultat de simulation correspondant à la configuration en cours, pour chacune des NM images de réservoir considérées. pi : poids affecté à chaque Ci : poids affecté à chaque ê, . 20 7. Procédé selon l'une des revendications précédentes, dans laquelle on réitère le procédé en faisant varier le nombre de puits. 8. Procédé selon l'une des revendications précédentes, dans lequel on optimise lesdites configurations au moyen d'un algorithme de type CMA-ES. 25 9. Procédé selon l'une des revendications précédentes, dans lequel les puits à forer comportent des drains multi-branches.