Brevet FR3038091A1 (20161230)
REVENDICATIONS1. Procédé de détermination d'une pluralité de premières lois reliant la perméabilité à la porosité au sein d'un réservoir souterrain (1), ledit procédé comprend : une étape d'obtention (E100) d'une première pluralité de points de mesure pour ledit réservoir (1), chaque point de mesure comprenant une donnée de porosité et une première donnée de perméabilité ; une étape de définition (E200) d'une famille de lois reliant la porosité à au moins une perméabilité ; une première étape de dénombrement (E500), pour chaque loi de ladite famille, des points de mesure de la pluralité de points reproduits par ladite loi ; et - une étape de sélection(E700) d'une pluralité de premières lois dans ladite famille sur la base d'au moins le résultat de ladite étape de dénombrement (E500).\n<CLM>\1\tProcédé selon la revendication 1, dans lequel ladite pluralité de premières lois sélectionnée est choisie parmi les lois de ladite famille reproduisant au moins un nombre minimal de points de mesure de ladite pluralité de points de mesure.\n<CLM>\1\tProcédé selon l'une des revendications 1 à 2, dans lequel un résultat de dénombrement est pondéré de manière à être plus important si des points de mesures reproduits par ladite loi sont répartis le long d'un voisinage de la courbe représentant ladite loi.\n<CLM>\1\tProcédé selon l'une des revendications 1 à 3, dans lequel lesdites lois liant la porosité à la perméabilité sont paramétrées par au moins deux paramètres.\n<CLM>\1\tProcédé selon l'une des revendications 1 à 4, dans lequel dans lequel les lois de ladite famille de loi sont des lois log-log définies par deux paramètres A et B et ladite étape de dénombrement (E500) comprend : - une étape de transformation (E300) de ladite pluralité de points de mesure par application de la fonction logarithme aux dites premières données de perméabilité et auxdites données de porosité ; - une étape de représentation (E300) de ladite pluralité de points transformés sous forme d'une image d'intensité ; et 3 03 809 1 41 ledit dénombrement comprenant un comptage du nombre de points de ladite pluralité de points transformés représentée dans ladite image d'intensité et dont la distance à une droite ayant pour équation Y=A.X+B est inférieure à un seuil. 5\n<CLM>\1\tProcédé selon l'une des revendications 1 à 4, dans lequel dans lequel les lois de ladite famille de loi sont des lois semi-log définies par deux paramètres A et B et ladite étape de dénombrement comprend : une étape de transformation (E300) de ladite pluralité de points de mesure par application de la fonction logarithme aux dites première données de perméabilité ; 10 une étape de représentation (E300) de ladite pluralité de points transformés sous la forme d'une image d'intensité ; et ledit dénombrement comprenant un comptage du nombre de points de ladite pluralité de points transformés représentée dans ladite image d'intensité et dont la distance à une droite ayant pour équation Y=A.X+B est inférieure à un seuil. 15\n<CLM>\1\tProcédé selon l'une des revendications 5 à 6 comprenant en outre une étape de lissage (E350) de ladite image d'intensité préalablement audit dénombrement.\n<CLM>\1\tProcédé selon l'une des revendications 1 à 7, dans lequel l'obtention desdites 20 données de porosité et desdites premières données de perméabilité comprend une analyse de carottes d'échantillons issues dudit réservoir ou une analyse de mesures de diagraphie, et ladite étape d'obtention (E100) comprend en outre une étape (E150) d'ajout de points de mesure supplémentaires dans ladite première pluralité de points, lesdits points de mesure supplémentaires ajoutés étant sélectionnés parmi ladite 25 première pluralité de points de mesure à partir de l'analyse de deuxièmes données de perméabilité obtenues à partir d'au moins un test de formation effectué au sein dudit réservoir.\n<CLM>\1\tProcédé de détermination selon la revendication 8 dans lequel ladite étape 30 d'ajout (E150) comprend en outre : une étape (F400) d'obtention d'un histogramme réel (Dist2) du logarithme desdites premières données de perméabilité obtenues par analyse desdites carottes ou par analyse desdites mesures de diagraphie ; une étape (F300, F450) d'obtention d'un histogramme théorique (Dist1) du logarithme desdites deuxièmes données de perméabilité obtenues à partir dudit au 3 03 80 9 1 42 moins un test de formation, lesdits intervalles dudit histogramme théorique étant égales auxdits intervalles dudit histogramme réel ; une étape (F500) d'obtention d'un histogramme agrégé de probabilité obtenu par le calcul du produit entre lesdits histogrammes théorique et réelle ; 5 une étape (F600) d'obtention pour au moins un intervalle de perméabilité pour lequel ledit histogramme agrégé est non nul d'un ensemble de points de mesure supplémentaires, le nombre de points de mesure supplémentaires dudit ensemble étant fonction de la valeur dudit histogramme agrégé évaluée pour ledit au moins un intervalle de perméabilité, lesdits points de mesure supplémentaires étant 10 sélectionnés aléatoirement parmi les points de mesure de ladite pluralité pour lesquels ladite première donnée de perméabilité correspond audit au moins un intervalle de perméabilité.\n<CLM>\1\tProcédé selon l'une des revendications 8 ou 9, dans lequel ladite première 15 donnée de perméabilité et ladite deuxième donnée de perméabilité sont des perméabilités horizontales.\n<CLM>\1\tProcédé selon l'une des revendications 8 ou 9, dans lequel ladite première donnée de perméabilité et ladite deuxième donnée de perméabilité sont des 20 perméabilités verticales.\n<CLM>\1\tProcédé selon la revendication 10, dans lequel ladite famille de lois est paramétrée par une pluralité de paramètres et le procédé comprend en outre : une étape d'obtention d'une deuxième pluralité de points de mesure, lesdits points de 25 mesure comprenant une desdites données de porosité et une troisième donnée de perméabilité verticale ; une deuxième étape de dénombrement, pour chaque loi de ladite famille, des points de ladite deuxième pluralité de points reproduits par ladite loi ; une étape de représentation des résultats desdites première et deuxième étapes de 30 dénombrement sous forme d'un premier et d'un deuxième signal d'intensité dépendant de ladite pluralité de paramètres ; une étape d'estimation d'un vecteur de translation par analyse d'une corrélation entre ledit premier signal et ledit deuxième signal d'intensité ; une étape de translation dudit vecteur de translation dudit deuxième signal 35 d'intensité ;ladite étape de sélection prend en compte au moins l'analyse dudit premier signal d'intensité et dudit deuxième signale d'intensité, ledit procédé 3 03 809 1 43 comprenant en outre la détermination d'une pluralité de deuxièmes lois reliant la perméabilité verticale à la porosité, ladite pluralité de deuxième loi étant obtenu à partir de ladite pluralité de premières lois par translation desdits paramètres par ledit vecteur. 5\n<CLM>\1\tProcédé selon la revendication 11, dans lequel ladite famille de lois est paramétrée par une pluralité de paramètres et le procédé comprend en outre : une étape d'obtention d'une deuxième pluralité de points de mesure, lesdits points de mesure comprenant une desdites données de porosité et une troisième donnée de 10 perméabilité horizontale ; une deuxième étape de dénombrement, pour chaque loi de ladite famille, des points de ladite deuxième pluralité de points reproduits par ladite loi ; une étape de représentation des résultats desdites première et deuxième étapes de dénombrement sous forme d'un premier et d'un deuxième signal d'intensité 15 dépendant de ladite pluralité de paramètres ; une étape d'estimation d'un vecteur de translation par analyse d'une corrélation entre ledit premier signal et ledit deuxième signal d'intensité ; une étape de translation dudit vecteur de translation dudit deuxième signal d'intensité ; 20 ladite étape de sélection prend en compte au moins l'analyse dudit premier signal d'intensité et dudit deuxième signal d'intensité , ledit procédé comprenant en outre la détermination d'une pluralité de deuxièmes lois reliant la perméabilité horizontale à la porosité, ladite pluralité de deuxième loi étant obtenu à partir de ladite pluralité de premières lois par translation desdits paramètres par ledit vecteur. 25\n<CLM>\1\tProcédé selon la revendication 12 ou la revendication 13 comprenant en outre une étape de normalisation dudit premier et dudit deuxième signal d'intensité préalablement à ladite étape d'estimation dudit vecteur de translation. 30\n<CLM>\1\tProcédé d'estimation d'au moins une perméabilité moyenne pour un ensemble de puits d'un réservoir souterrain, ledit procédé comprend : une étape (H200) d'obtention d'une distribution de données de porosité pour ledit ensemble de puits ; une étape (H300) d'obtention d'une pluralité de premières lois reliant la porosité à la 35 perméabilité pour ledit ensemble selon un procédé de détermination selon l'une des revendications 1 à 14 ; 3 03 809 1 44 une étape (H400) d'obtention d'une loi de probabilité approximant ladite distribution sur la base d'au moins ladite pluralité de premières lois ; et une étape (H500) d'estimation de ladite au moins une perméabilité moyenne pour ledit ensemble de puits à partir d'au moins ladite loi normale asymétrique et de ladite 5 pluralité de premières lois.\n<CLM>\1\tProcédé selon la revendication 15, dans lequel ladite étape (H400) d'obtention d'une loi de probabilité s'effectue en minimisant une fonction objectif prenant en compte au moins un terme parmi les trois termes suivants : 10 un premier terme favorisant les lois de probabilité approchant au mieux ladite distribution de données de porosité ; un deuxième terme favorisant les lois de probabilité dont la valeur moyenne approche au mieux la valeur moyenne de ladite distribution de données de porosité ; et 15 un troisième terme favorisant les lois de probabilité qui minimisent pour chaque loi sélectionnée de ladite première pluralité de lois la somme des différences entre la valeur de la perméabilité moyenne calculée après application de la loi sélectionnée à ladite distribution de données de porosité et la valeur de la perméabilité moyenne calculée après application de la loi sélectionnée à ladite loi de probabilité. 20\n<CLM>\1\tProcédé selon la revendication 16, dans lequel : ladite pluralité de premières lois est définie par la relation log(Kh) = f(cp ) + Bi) où Kh est une perméabilité horizontale, Ai et BI sont deux paramètres réels appartenant à une région déterminée de l'espace défini par les paramètres A et B et 25 où f est la fonction identité ou la fonction log ; et - ladite fonction objectif est une combinaison linéaire ccEi + (1 - a)[(1 13-1E_3 + IRE ] - + où a et 13 sont deux coefficients positifs et inférieurs à un et où : 2 ledit premier terme El est égal à Er=i (Fi - LP((pi , S1, Sg)) ; ledit deuxième terme E2 est égal à (wi - U(Si, Sg))2 ; et 2 30 ledit troisième terme E3 est égal à el (Rffli, Bi) - Bi, m, S1, ..., Sg)) avec log(Khi (Ai, Bi, (pi )) = f((pi) + Bi) et iZh(Ai, Bi) = Er,.1 Khi (Ai, Bi, (pi ). Fi et où : N est le cardinal de ladite région déterminée de l'espace défini par les paramètres A et .B, n est le nombre d'intervalles de ladite distribution de données de porosité, 3 0 3 8 0 9 1 45 Fi est la fréquence d'occurrence associée à chacun des intervalles de porosité (pi , (TI, est la valeur associée auxdits intervalles de porosité, et Kh(A1,131, S1, Sg) est la perméabilité moyenne horizontale calculée après 5 application de la loi 1og(Kh ) = (Ai. f(cp ) + Bi) à la loi de probabilité LP (g) ; S9) dépendant des paramètres S1, ..., S9 ; et ladite au moins une perméabilité moyenne est une perméabilité moyenne horizontale Khs donnée par la formule : Khs (Ai, Bi) = exp (Ai. f((pi) + Bi) . LPs((pj ; S1, ..., Se), où LPs est la loi de 10 probabilité minimisant la fonction objectif.\n<CLM>\1\tProcédé selon la revendication 16, dans lequel : ladite pluralité de premières lois est définie par la relation log(Kv) = (Avi. f(cp ) + Bvi) où Kv est une perméabilité verticale, Avi et Bvi sont deux paramètres réels 15 appartenant à une région déterminée de l'espace défini par les paramètres A et B et où f est la fonction identité ou la fonction log ; et ladite fonction objectif est une combinaison linéaire aE1 + (1 - oc)[(1 - 0)E3 + 8E2], où a et 13 sont deux coefficients positifs et inférieurs à un et où : 2 - ledit premier terme El est égal à Er,' (F1 - LP((pi , S1, ..., Se)) ; 20 - ledit deuxième terme E2 est égal à ((T), - LP(S1, Sg))2 ; et 2 - ledit troisième terme E3 est égal à l (1v(Avi, Bvi) - Kv(Avi, Bvi, ... Se)) E" avec log(Kvi(Avi, Bvi, (pi )) = (Ai. f((pi) + Bi) et 1/1-7/(Avi, Bvi) = Er=i Fi/ Kvi(Avi, Bvi, (pi ) et où : N est le cardinal de ladite région déterminée de l'espace défini par les paramètres 25 Av et Bv, n est le nombre d'intervalles de ladite distribution de données de porosité, F1 est la fréquence d'occurrence associée à chacun des intervalles de porosité (pi , Tp-; est la valeur associée auxdits intervalles de porosité, et 30 1-(Avt, Bvi, S1, , Sg) est la perméabilité moyenne verticale calculée après application de la loi log(Kv ) = (Avi. f(cp ) + Bvi) à la loi de probabilité LP (cps ; S1, ..., S9) dépendant des paramètres S1, ..., S9 ; et ladite au moins une perméabilité moyenne est une perméabilité moyenne verticale Kvs donnée par la formule : 3 03 809 1 46 1/Kvs(Avi, Bvi) = LPs(cpi ; Si, ..., Sg )/exp(Avi. f(cpj) + Bvi), où LPs est la loi de probabilité minimisant ladite fonction objectif.\n<CLM>\1\tProcédé selon la revendication 15 comprenant en outre une étape 5 d'obtention d'une pluralité de troisième lois reliant la porosité à la perméabilité verticale sur la base d'au moins le résultat dudit deuxième dénombrement et dans lequel : - ladite pluralité de premières lois relie la porosité à la perméabilité horizontale, et - ladite étape d'obtention d'une loi normale de probabilité s'effectue en minimisant une fonction objectif prenant en compte au moins un terme parmi les trois termes 10 suivants : un premier terme favorisant les lois de probabilité approchant au mieux ladite distribution de données de porosité ; un deuxième terme favorisant les lois de probabilité dont la valeur moyenne approche au mieux la valeur moyenne de ladite distribution de données de porosité ; 15 - un troisième terme favorisant les lois de probabilité qui minimisent pour chaque loi sélectionnée de ladite première pluralité de lois et pour chaque loi sélectionnée de ladite troisième pluralité de lois la somme des différences entre la valeur de la perméabilité moyenne horizontale totale calculée après application des lois sélectionnées à ladite distribution de données de porosité et la valeur de la 20 perméabilité moyenne horizontale totale calculée après application des lois sélectionnées à ladite loi de probabilité ; et ladite au moins une perméabilité moyenne est une perméabilité moyenne horizontale totale. 25\n<CLM>\1\tProcédé selon la revendication 19, dans lequel : ladite pluralité de premières lois est définie par la relation log(Kh) = (A1. f(9 ) + Bi) où, Kh est une perméabilité horizontale, Ai et Bi sont deux paramètres réels appartenant à une région déterminée de l'espace défini par les paramètres A et B et où f est la fonction identité ou la fonction log ; 30 ladite pluralité de troisième lois est définie par la relation log(Kv) = (Av.. f(cp ) + Bv,',) où, Kv est une perméabilité verticale, Avm et Bv, sont deux paramètres réels appartenant à une région déterminée de l'espace défini par les paramètres Av et Bv ; et ladite fonction objectif est une combinaison linéaire aEi + (1 - a)[(1 - (3)E3 + 13E21, 35 où a et sont deux coefficients positifs et inférieurs à un et où : 3038091 47 2 - ledit premier terme El est égal à Er.,i (Fi - LP(q); , S1, ..., Sg)) ; - ledit deuxième terme E2 est égal à W., - LP(Si, ...,Sg))2 ; et - ledit troisième terme E3 est égal à s 2 U.1.1E14-1 (Rht(Ai, B1, Avm, Bvm) - Kht(Ab Bi, Avm, Bvm, S1, , Sg ) avec 5 Kht(Ab Avm, Bvm) = Ch. Kh(Ab Bi) + (1 - Ch). Kv(Avm, Bvm), log(Kvi(Avm, Bvm, (pi )) = (AVm. f((pi) + Bvm) , 1/Kv(Avm, BVm) = Er =1 Fj/Kvi (Avm, BVm, (Pi )1 log(Khi(Ai, Bi, (pi )) = (Ab f(cpj) + Bi), Ih(Jeki, Bi) = E 1 Fi. Khi (Ai, Bi, (pi ) et où : 10 Ch est un coefficient positif compris entre 0 et 1, N est le cardinal de ladite région déterminée de l'espace défini par les paramètres A et B, M est le cardinal de ladite région déterminée de l'espace défini par les paramètres Av et Bv, 15 n est le nombre d'intervalles de ladite distribution de données de porosité, Fi est la fréquence d'occurrence associée à chacun des intervalles de porosité (pi , - est la valeur associée auxdits intervalles de porosité, et Kht(Ab Bi, Avm, Bvm, S1, ..., Sg) est la perméabilité moyenne horizontale calculée 20 après application des lois log(Kh ) = (Ai. f(cp ) + B1) et log(Kv ) = (Avm. f(cp ) + Bvhn) à la loi de probabilité LP (q:e ; S1, S9) dépendant des paramètres S1, Sg ; et ladite au moins une perméabilité moyenne est une perméabilité moyenne horizontale totale donnée par l'une des formules : 25 Khts(Ai, Bi, AVm, Bvm) = Z=1 Kht] (Ab Bi, Avm, BVm, (pi ). LPs((Pi ; S1, , Sg ), où LPs la loi de probabilité minimisant ladite fonction objectif.\n<CLM>\1\tProcédé selon la revendication 15 comprenant en outre une étape d'obtention d'une pluralité de troisièmes lois reliant la porosité à la perméabilité verticale 30 sur la base d'au moins le résultat dudit deuxième dénombrement et dans lequel : - ladite pluralité de premières lois relie la porosité à la perméabilité horizontale, et - ladite étape d'obtention d'une loi normale de probabilité s'effectue en minimisant une fonction objectif prenant en compte au moins un terme parmi les trois termes suivants : 3 0 3 8 0 9 1 48 un premier terme favorisant les lois de probabilité approchant au mieux ladite distribution de données de porosité ; un deuxième terme favorisant les lois de probabilité dont la valeur moyenne approche au mieux la valeur moyenne de ladite distribution de données de porosité ; 5 un troisième terme favorisant les lois de probabilité qui minimisent pour chaque loi sélectionnée de ladite première pluralité de lois et pour chaque loi sélectionnée de la pluralité de troisième lois la somme des différences entre la valeur de la perméabilité moyenne verticale totale calculée après application des lois sélectionnées à ladite distribution de données de porosité et la valeur de la perméabilité moyenne verticale 10 totale calculée après application des lois sélectionnées à ladite loi de probabilité ; et ladite au moins une perméabilité moyenne est une perméabilité moyenne verticale totale.\n<CLM>\1\tProcédé selon la revendication 21, dans lequel : 15 ladite pluralité de premières lois est définie par la relation log(Kh) = f(cp ) + Bi) où, Kh est une perméabilité horizontale, Ai et B1 sont deux paramètres réels appartenant à une région déterminée de l'espace défini par les paramètres A et B et où f est la fonction identité ou la fonction log ; ladite pluralité de troisièmes lois est définie par la relation log(Kv) = (Avm. f(q) ) + 20 Bvm) où, K,' est une perméabilité verticale, Avn, et Bvn, sont deux paramètres réels appartenant à une région déterminée de l'espace défini par les paramètres Av et Bv ; et ladite fonction objectif est une combinaison linéaire aEi + (1 - a)[(1- (3)E3 + [3E2], où a et 13 sont deux coefficients positifs et inférieurs à un et où : \ 2 25 - ledit premier terme El est égal à Er=i (Fi - LP(4); , S1, ...,Sg)) ; - ledit deuxième terme E2 est égal à (-cp - LP(Si, Sg))2 ; et - ledit troisième terme E3 est égal à , 2 (Rvt(Ai, Bi, Avm, Bvm) - Kvt(A1, Bi, Avm, Bvm, S1, ..., Sg )) avec Kvt(A1, Bi, Avm, Bvm) = Cv. Kh(At, Bi) + (1 - Cy). Kv(Avm, Bvm), 30 log(Kvi(Avm, BVm, (pi )) = (Avm. Kg)) + Bvm) , 1/Kv(Avm, BVm) = Er,,i Fi/Kvi (Avm, BVm, (pi ), log(Khi (Ai, Bi, (pi )) = (Ai. f(cpi) + Bi), kh(Ai, Bi) = Er.,1 Fj. Khi (Ai, Bb (p1 ) et où : Cv est un coefficient positif compris entre 0 et 1, 3 03 80 9 1 49 N est le cardinal de ladite région déterminée de l'espace défini par les paramètres A et B, M est le cardinal de ladite région déterminée de l'espace défini par les paramètres Av et Bv, 5 n est le nombre d'intervalles de ladite distribution de données de porosité, Fi est la fréquence d'occurrence associée à chacun des intervalles de porosité (pi , est la valeur associée auxdits intervalles de porosité, et Kvt(Ai, Bi, Avm, Bvm, S1, , Sg) est la perméabilité moyenne verticale calculée 10 après application des lois log(Kh ) = (Ai. f(cp ) + Bi) et log(Kv ) = (Avm. f((p ) + Bvim) à la loi de probabilité LP (q) ; S1, ..., S9) dépendant des paramètres S1, Sg ; et ladite au moins une perméabilité moyenne est une perméabilité moyenne verticale totale donnée par l'une des formules : 15 1/Kvts (Ai, Bi, Avm, Bvm) 1/Kvti (Ai, Bi, Avm, Bvm, (pi ). LPs(cpj ; S1, ..., Sg ), où LPs la loi de probabilité minimisant ladite fonction objectif.\n<CLM>\1\tProcédé selon la revendication 15, dans lequel ladite étape d'obtention d'une loi de probabilité s'effectue sur la base d'au moins ladite pluralité de deuxièmes lois, ledit 20 procédé comprend en outre une étape d'estimation d'au moins une perméabilité moyenne verticale totale sur la base d'au moins ladite loi de probabilité et de ladite pluralité de deuxièmes lois.\n<CLM>\1\tProcédé selon la revendication 23, dans lequel ladite étape d'obtention d'une 25 loi de probabilité s'effectue en minimisant une fonction objectif prenant en compte au moins un terme parmi les trois termes suivants : un premier terme favorisant les lois de probabilité approchant au mieux ladite distribution de données de porosité ; un deuxième terme favorisant les lois de probabilité dont la valeur moyenne 30 approche au mieux la valeur moyenne de ladite distribution de données de porosité ; et un troisième terme favorisant les lois de probabilité qui minimisent pour chaque loi sélectionnée de ladite première pluralité de lois : o la somme des différences entre la valeur de la perméabilité moyenne 35 horizontale totale calculée après application de la loi sélectionnée à ladite distribution de données de porosité et la valeur de la perméabilité moyenne 3038091 50 horizontale totale calculée après application de la loi sélectionnée à ladite loi de probabilité, ainsi que o la somme des différences entre la valeur de la perméabilité moyenne verticale totale calculée après application de la loi sélectionnée à ladite 5 distribution de données de porosité et la valeur de la perméabilité moyenne verticale totale calculée après application de la loi sélectionnée à ladite loi de probabilité.\n<CLM>\1\tProcédé selon la revendication 24, dans lequel : 10 ladite pluralité de premières lois est définie par la relation log(%) --= f(q) ) + Bi) où Kh est une perméabilité horizontale, A et Bi sont deux paramètres réels appartenant à une région déterminée de l'espace défini par les paramètres A et B et où f est la fonction identité ou la fonction log ; ladite pluralité de deuxièmes lois est définie par la relation log(Kv) = ((A1 + 15 dA). f(q) ) + Bi + dB) où Kv est une perméabilité verticale et dA et dB sont deux paramètres réels ; et ladite fonction objectif est une combinaison linéaire aE1 + (1 - a)[(1- B)E3 + 13E2], où a et (3 sont deux coefficients positifs et inférieurs à un et où : \ 2 ledit premier terme El est égal à Er=i (Fi- LP(cp; , S1, Sg)) ; 20 ledit deuxième terme E2 est égal à (iT - LP(51, Sg))2 ; et s 2 ledit troisième terme E3 est égal à ri% (Kht(Ai, Bi) - Kht(Ai, Bi, m, Sl, , Sg)) + 2 (Kvt(Ai, BD - Kvt(Ai, Bi, m, S1, ..., Sg)) avec Kht(Ai, = Ch. Kh(A1, + (1 - Ch). Kv(Ai, Kvt(Ai, = Cv. Kh(A4, + - Cv). Kv(Ai, 25 log(Kvi(A, Bi)) = ((A1 + dA). f(cpi) + Bi + dB) et 1/Kv(Ab = Ej 1 Fj/Kvi (Ai, BD, log(Khi (Ai, Bi)) = f(cpD + B1) et k«Ai, Bi) = E 1 Fj. Khj (Ai, , et où : Ch et C sont des coefficients positifs compris entre 0 et 1, N est le cardinal de ladite région déterminée de l'espace défini par les paramètres A et B, 30 n est le nombre d'intervalles de ladite distribution de données de porosité, Fi est la fréquence d'occurrence associée à chacun des intervalles de porosité (pi , est la valeur associée auxdits intervalles de porosité, et 3038091 51 Kht(Ai, B1, m, Si, Si/S2) est la perméabilité moyenne horizontale totale calculée après application des lois log(Kh ) = (Ai. f(q) ) + Bi) et log(Kv ) = ((Ai+dA). f((p ) + Bi + dB) à la loi de probabilité LP (q) ; S1, S9) dépendant des paramètres S1, Sg ; 5 Kvt(Ai, Bi, m, Si, Si /S2) est la perméabilité moyenne verticale totale calculée après application de la loi log(Kv ) = ((Ai + dA). f(cp ) + Bi + dB) à la loi de probabilité LP (q:1; S1, ..., S9) dépendant des paramètres S1, ..., Sg ; ladite au moins une perméabilité moyenne est une perméabilité moyenne horizontale totale donnée par l'une des formules : 10 Khts (Ai, Bi) = Er,1 Khtj (Ai, Bi). LPs @pi , S1, , Sg) où LPs la loi de probabilité minimisant ladite fonction objectif ; et ladite au moins une perméabilité moyenne verticale totale est donnée par l'une des formules : 1/Kvts (Ai, Bi) = Er=, 1/Kvti (Ai, Bi). LPs (cpi , Si, ..., Sg)). 15\n<CLM>\1\tProcédé selon l'une des revendications 15 à 25 dans lequel ladite loi de probabilité est une loi normale ou une combinaison linéaire de lois normales.\n<CLM>\1\tProcédé selon l'une des revendications 15 à 25 dans lequel ladite loi de 20 probabilité est une loi normale asymétrique.\n<CLM>\1\tProcédé selon la revendication 25, dans lequel le coefficient Ch est supérieur à 0.75 et inférieur à\n<CLM>\1\t25\n<CLM>\1\tProcédé selon la revendication 22 ou 25, dans lequel le coefficient C., est supérieur à 0 et inférieur à 0.25.\n<CLM>\1\tProcédé de calcul d'une perméabilité moyenne en un emplacement d'un réservoir souterrain, ledit procédé comprend : 30 une étape (3100) de sélection d'un ensemble de puits dudit réservoir, ledit ensemble comprenant au moins un puits ; une étape (J200) de détermination d'une pluralité de premières lois reliant la perméabilité à la porosité pour ledit ensemble de puits selon l'une des revendications 1 à14 35 pour chacun desdits puits dudit ensemble : 3 03 80 9 1 52 o une étape (3300) d'obtention d'une distribution de données de porosité pour le puits ; o une étape (3300) d'obtention d'une loi de probabilité approximant ladite distribution sur la base d'au moins ladite pluralité de premières lois ; et 5 une étape (3400) de calcul d'une loi de probabilité en ledit emplacement à partir desdites lois de probabilité obtenues pour chacun des dits puits ; et une étape (3500) de calcul de ladite perméabilité moyenne en ledit emplacement à partir d'au moins ladite loi de probabilité en ledit emplacement et d'au moins ladite pluralité de premières lois. 10\n<CLM>\1\tProcédé de calcul selon la revendication 30 comprenant en outre une étape (J600) de calcul d'une porosité moyenne en ledit emplacement à partir d'au moins ladite loi de probabilité en ledit emplacement 15\n<CLM>\1\tProgramme d'ordinateur comportant des instructions pour l'exécution des étapes du procédé de détermination d'au moins une première loi reliant la perméabilité à la porosité selon l'une quelconque des revendications 1 à 14 lorsque ledit programme est exécuté par ordinateur. 20\n<CLM>\1\tSupport d'enregistrement lisible par un ordinateur sur lequel est enregistré un programme d'ordinateur comprenant des instructions pour l'exécution des étapes du procédé de détermination d'au moins une première loi reliant la perméabilité à la porosité selon l'une quelconque des revendications 1 à\n<CLM>\1\t25\n<CLM>\1\tProgramme d'ordinateur comportant des instructions pour l'exécution des étapes du procédé d'estimation d'une perméabilité moyenne selon l'une quelconque des revendications 15 à 29 lorsque ledit programme est exécuté par ordinateur.\n<CLM>\1\tSupport d'enregistrement lisible par un ordinateur sur lequel est enregistré 30 un programme d'ordinateur comprenant des instructions pour l'exécution des étapes du procédé d'estimation d'une perméabilité moyenne selon l'une quelconque des revendications 15 à\n<CLM>\1\t36. Programme d'ordinateur comportant des instructions pour l'exécution des 35 étapes du procédé de calcul d'une perméabilité moyenne selon l'une quelconque des revendications 30 à 31 lorsque ledit programme est exécuté par ordinateur. 303 809 1 53\n<CLM>\1\tSupport d'enregistrement lisible par un ordinateur sur lequel est enregistré un programme d'ordinateur comprenant des instructions pour l'exécution des étapes du procédé de calcul d'une perméabilité moyenne selon l'une quelconque des revendications 5 30 à\n<CLM>\1\t38. Dispositif de détermination (3) d'une pluralité de premières lois reliant la perméabilité à la porosité au sein d'un réservoir souterrain, ledit dispositif 3 comprenant : un module d'obtention (3B1) d'une première pluralité de points de mesure pour ledit 10 réservoir, chaque point de mesure comprenant une donnée de porosité et une première donnée de perméabilité ; un module de définition (3B2) d'une famille de lois reliant la porosité à au moins une perméabilité ; un premier module de dénombrement (3B3), pour chaque loi de ladite famille, des 15 points de mesure de la pluralité de points reproduits par ladite loi ; et un module de sélection (3B4) d'une pluralité de premières lois dans ladite famille sur la base d'au moins le résultat du dénombrement effectué par ledit premier module de dénombrement (3B3). 20\n<CLM>\1\tDispositif d'estimation (4) d'au moins une perméabilité moyenne pour un ensemble de puits d'un réservoir souterrain, ledit dispositif comprend : un module d'obtention (4B1) d'une distribution de données de porosité pour ledit ensemble de puits ; un dispositif de détermination (4B2) d'une pluralité de premières lois reliant la 25 porosité à la perméabilité pour ledit ensemble de puits selon la revendication 38 ; un module d'obtention (4B3) d'une loi de probabilité approximant ladite distribution sur la base d'au moins ladite pluralité de premières lois ; et - un module d'estimation (4B4) de la perméabilité moyenne pour ledit ensemble de puits à partir d'au moins ladite loi normale asymétrique et de ladite pluralité de 30 premières lois.\n<CLM>\1\tDispositif de calcul (5) d'une perméabilité moyenne en un emplacement d'un réservoir souterrain, ledit dispositif comprend : - un module de sélection (5B1) d'un ensemble de puits dudit réservoir, ledit ensemble 35 comprenant au moins un puits ; 3038091 54 un dispositif de détermination (5B2) d'une pluralité de premières lois reliant la perméabilité à la porosité pour ledit ensemble de puits selon la revendication 38; un module d'obtention (5B3) d'une distribution de données de porosité pour chacun des puits dudit ensemble ; 5 un module d'obtention (5B4) d'une loi de probabilité approximant ladite distribution pour chacun desdits puits dudit ensemble sur la base d'au moins ladite pluralité de premières lois ; et un module de calcul (5B5) d'une loi de probabilité en ledit emplacement à partir desdites lois de probabilité obtenues pour chacun des dits puits 10 un module de calcul (5B6) de ladite perméabilité moyenne en ledit emplacement à partir d'au moins ladite loi de probabilité en ledit emplacement et d'au moins ladite pluralité de premières lois.