FR 2849211 A1 20040625 FR 0216456 A 20021220 La présente invention a pour objet une méthode pour réaliser un modèle pour simuler le remplissage sédimentaire de bassins sur des échelles de temps et d'espace importants. La méthode selon l'invention concerne plus particulièrement la réalisation d'un 10 modèle stratigraphique numérique permettant la simulation en deux ou trois dimensions ( 2 D ou 3 D) du remplissage multi-lithologique d'un bassin, dans le but de simuler la réponse stratigraphique du système sédimentaire à des variations au cours du temps de l'eustatisme, de la subsidence, des apports sédimentaires et des paramètres physiques régissant le transport des sédiments dans le bassin. ART ANTERIEUR Les progrès récents en géologie, qui ont donné naissance au cours de ces vingt dernières années à la stratigraphie sismique, puis à la stratigraphie génétique, ont profondément modifié la compréhension de l'histoire du remplissage sédimentaire des bassins sur de grandes échelles de temps et d'espace en montrant l'influence primordiale de 20 trois paramètres principaux: l'eustatisme, la tectonique et le flux sédimentaire. De nombreux modèles et notamment des modèles numériques déterministes, ont été élaborés afin de pouvoir appréhender les implications géométriques et lithologiques de ces nouvelles approches. Ces modèles numériques simulent le transport et la sédimentation (ou l'érosion) des 25 sédiments dans le bassin en se basant sur une description plus ou moins complexe de la nature, à partir d'une estimation de l'eustatisme, de la subsidence et de l'apport sédimentaire On rappelle que l'on désigne par eustatisme les variations de la surface des océans enregistrées simultanément sur l'ensemble de la surface terrestre, et par subsidence le déplacement absolu du fond d'un bassin sédimentaire par rapport à un niveau repère fixe au cours du temps Ces modèles permettent ainsi de tester l'influence de différents concepts (relation entre le climat et l'eustatisme, etc) sur l'agencement des unités sédimentaires De plus, appliqués sur des cas réels, ils permettent de tester la cohérence des paramètres 5 introduits dans le modèle tels que les variations eustatiques et tectoniques et de conforter l'interprétation géologique du bassin étudié. On attend d'un modèle qu'il puisse définir les grandes tendances de variations de faciès (variation du rapport sable/argile, teneur en carbonate, etc) au sein des unités génétiques Pour ce faire, il est nécessaire de concevoir une modélisation multi10 lithologique Le modèle doit pouvoir simuler le transport et la sédimentation de différentes lithologies siliciclastiques (sable, argile, etc) et carbonatées (récif, boue pélagique carbonatée, bioclastes, etc) afin que ces tendances de variations de faciès soient un résultat du modèle, indépendant d'a priori géologique tel que plaine alluviale sableuse et domaine marin argileux Dans ce modèle idéal, les sédiments siliciclastiques sont introduits dans le 15 bassin par les fleuves et les courants marins, tandis que les carbonates sont produits en domaine marin en prenant en compte la bathymétrie, la turbidité de l'au, et l'action des vagues Après avoir été introduits au niveau des frontières du bassin, ou produits dans le bassin, les sédiments sont transportés puis sédimentés. Il existe trois grandes familles de modèle déterministe qui régissent le transport des 20 sédiments: les modèles particulaires basés sur la résolution du mouvement des particules (calcul de l'écoulement de l'eau puis relation entre flux d'eau et flux de sédiments); les modèles diffusifs basés sur une équation de diffusion dans laquelle la définition du coefficient de diffusion est plus ou moins raffinée (prise en compte de plusieurs 25 lithologies, du flux d'eau, de l'environnement de dépôt, etc); et les modèles géométriques basés sur une définition du profil géométrique des environnements de dépôt (longueur et pente de la plaine alluviale, etc) ou sur une définition géométrique du taux de sédimentation (décroissance exponentielle en domaine marin, etc). Les modèles particulaires utilisent une description soignée des processus sédimentaires et sont donc aussi chaotiques que la nature Ils permettent essentiellement d'obtenir des simulations à l'échelle du réservoir (longueur de l'ordre de 0,5 à 50 km et durée de l'ordre de 5 à 500 ka) Les modèles diffusifs et géométriques fournissent tous deux des approximations plus ou moins grossières des processus naturels Ils sont beaucoup plus stables que la réalité dont ils rendent compte; mais ils ne fournissent qu'une 5 estimation lissée de la nature Ces modèles sont applicables préférentiellement à l'échelle du bassin (longueur de l'ordre de 10 à 1000 km et durée de l'ordre de 0,1 à 100 Ma) Les modèles géométriques, basés sur une approximation de la morphologie des bassins, sont essentiellement applicables sur des cas simples en 2 D dans lesquels la subsidence, nature des sédiments et le climat ne perturbent pas trop la définition du profil d'équilibre Plus 10 généraux, les modèles diffusifs, fondés sur une approximation de la physique des sédiments, sont applicables en 3 D et peuvent traiter le transport et la sédimentation de multiples lithologies. Différents modèles diffusifs connus sont décrits par exemple par: Kenyon et Turcotte, 1985, Morphology of a delta prograding by bulk sediment 15 transport in Geol Soc Amer Bull, 96, 1457-1465. Begin, Z B, 1988, Application of a diffusion erosion model to alluvial channels which degrade due to base-level lowering Earth Surface Processes and Landforms, 13, 487-500. Rivenaes, J C, 1988, Application of a duallithology, depth dependent diffusion 20 equation in stratigraphic simulation Basin Research, 4, 133-146. Un avantage des modèles diffusifs est de permettre un retour vers les concepts géologiques en quantifiant certaines relations telles que la durée variable des phases de progradation et de rétrogradation des unités génétiques, ou l'évolution de la sablosité en fonction soit de la bathymétrie, soit de la tendance progradante ou rétrogradante du littoral. Les modèles antérieurs ont essentiellement pour but une expérimentation numérique de concepts théoriques le long de profils en 2 D. La méthode concerne la réalisation d'un modèle déterministe numérique de type diffusif permettant la simulation en 2 D ou 3 D du remplissage des bassins sédimentaires par des sédiments siliclastiques et carbonatés. Par le brevet FR 2 744 224 du demandeur, on connaît une méthode pour simuler le remplissage d'un bassin sédimentaire A partir de données connues sur l'architecture d'un bassin et de données de mesures: diagraphies de puits, campagnes sismiques, etc, on constitue un ensemble de données d'entrée portant sur l'espace disponible créé par 5 subsidence et eustatisme, sur l'apport et la production de sédiments fluviatiles ou marins, et sur des paramètres physiques de transport tels que des coefficients de diffusion des différentes lithologies Cet ensemble de données est appliqué à un modèle numérique Les résultats que l'on peut en déduire sur la géométrie et les lithologies des unités sédimentaires, sont confrontés avec les données de mesure et, de proche en proche, par 10 inversion, on affine les données d'entrée. LA METHODE SELON L'INVENTION La méthode de modélisation selon l'invention permet de constituer un modèle déterministe de type diffusif (en 2 D ou 3 D) qui simule le remplissage multi-lithologique d'un bassin sédimentaire Elle comporte, à partir de données de terrain connues relatives à 15 l'architecture du bassin et à de données de mesure telles que des données de diagraphie de puits ou des données sismiques, la constitution d'un ensemble de données d'entrée portant sur un espace disponible par subsidence et eustatisme, sur l'apport de sédiments fluviatiles ou marins et leur transport, et sur des paramètres physiques tels que des coefficients de diffusion des différentes lithologies, par un processus itératif incluant le maillage du bassin 20 en mailles de dimensions régulières; une modélisation suivant un schéma explicite en volumes finis avec des pas de temps constants, de façon à simuler le flux de chaque lithologie déposée sur chaque maille; une comparaison des résultats de la simulation aux données de terrain et une modification par inversion et de proche en proche des données d'entrée. La méthode est caractérisée principalement en ce que à chaque pas de temps et pour chaque lithologie: on modélise par des équations déterministes respectivement: ( 1) la moyenne des processus de transport agissant de manière permanente et continue sur le long terme (supérieures à quelques années); ( 2) la moyenne des processus de transport agissant de manière permanente et continue sur le court terme (de l'ordre de quelques jours à quelques mois) ; ( 3) la moyenne des processus de transport agissant de manière catastrophique sur le très court terme (de l'ordre de quelques heures à quelques jours); et on détermine la stratigraphie résultante du bassin en tenant compte de la conservation des masses. La méthode comporte par exemple: une modélisation des flux à long terme de sédiments transportés par l'eau en tenant compte d'un modèle de répartition des flux d'eau sur le bassin et de la capacité de transport des sédiments par de tels flux; une modélisation des flux à court terme de sédiments transportés par l'eau en utilisant 10 un modèle exponentiel pour estimer la vitesse de l'eau; une modélisation des flux catastrophiques transportés par l'eau dans des zones instables repérées en appliquant un critère de pente critique. Le modèle diffusif obtenu par la méthode permet d'appréhender en 3 D le remplissage multi-lithologique d'un bassin Afin de l'appliquer dans un contexte industriel, 15 une méthodologie d'inversion des paramètres est utilisée Une analyse des données disponibles, telles que des données de puits ou de sismiques, permet d'estimer un premier jeu de paramètres (eustatisme, subsidence, apports, etc) requis par le modèle Le modèle diffusif est alors utilisé pour obtenir une première simulation du bassin Cette simulation est comparée aux données Les paramètres sont alors modifiés afin de réduire les écarts 20 entre données et simulation Cette boucle d'inversion (simulation comparaison modification) est poursuivie et permet un affinement progressif de la quantification de chaque paramètre A la fin de cette boucle d'inversion, l'utilisateur du modèle dispose ainsi d'une simulation calée sur les données disponibles, mais également d'un ensemble de paramètres physiques améliorant sa connaissance du bassin étudié. La modélisation étant effectuée de façon détaillée selon les types de transport, on obtient une distribution plus précise des sédiments à l'intérieur de chaque géométrie et aussi une quantification plus fine de chacun d'eux. Les principales applications du modèle diffusif et de sa boucle d'inversion selon l'invention sont: la quantification des paramètres physiques et géologiques (eustasie, subsidence, localisation des sources de sédiments, paramètres de transport, etc) pour différents modèles géologiques, ce qui permet: de tester la cohérence de différents modèles géologiques; de déterminer le meilleur schéma de corrélation entre différents puits; et de déterminer la meilleure interprétation de données sismiques; la reconstitution de l'évolution au cours du temps de la stratigraphie d'un bassin sédimentaire, ce qui permet en particulier: d'estimer l'extension et la connectivité des différentes unités sédimentaires; ou d'estimer la distribution régionale des faciès sédimentaires. En exploration, le modèle ainsi formé permet d'identifier les régions potentiellement les plus intéressantes, grâce à la localisation et à la caractérisation des couches sédimentaires contenant les roches mères, les drains, les réservoirs et/ou les couvertures Les ingénieurs de réservoir peuvent de plus, à partir du modèle, déterminer les 15 propriétés pétrophysiques d'un réservoir étudié (porosité moyenne et perméabilité) et sélectionner les zones de production les plus prometteuses. PRESENTATION SUCCINCTE DES FIGURES Les caractéristiques et avantages de la méthode selon l'invention, apparaîtront plus clairement à la lecture de la description ci-après d'un exemple non limitatif de réalisation, 20 en se référant aux dessins annexés o: la figure 1 montre un organigramme de la méthode permettant de modéliser par itérations successives, la stratigraphie d'un bassin; la figure 2 illustre le mode de transport à long terme (permanents) des sédiments induits par les fleuves, les courants marins, ou l'action lente de la gravité (reptation, 25 etc); il est simulé à l'aide d'un modèle diffusif, contraint par le flux d'eau et une disponibilité réduite; dans le domaine marin, la fraction fine des sédiments (composante argileuse), est d'abord transportée par suspension, puis remaniée par les processus à long terme agissant sur le fonds marin. la figure 3 illustre le mode de transport à court terme des sédiments lors de phénomènes épisodiques mais intenses, tels que des pluies importantes induisant des crues fluviales et des courants hyperpycnaux en domaine marin; ces processus sont modélisés à l'aide d'un modèle diffusif, contraint par le flux d'eau, la disponibilité 5 réduite et la vitesse de l'eau; cette vitesse de l'eau étant calculé à l'aide d'un modèle exponentiel; et la figure 4 illustre le mode de transport en masse ou catastrophique des sédiments suite à des éboulements de falaise ou de pente; dans un premier temps, on détermine les zones instables à l'aide d'un critère sur la pente du sol; tout sédiment situé au dessus 10 d'une pente critique est considéré comme instable; on le détache du sol, et on le déplace à l'aide de l'équation de Newton (Accélération = Gravité Friction) , tant que sa vitesse est supérieure à une vitesse critique. DESCRIPTION DETAILLEE DE LA METHODE La méthode de modélisation selon l'invention permet d'obtenir une simulation 15 stratigraphique déterministe sur de grandes durées (de l'ordre de la dizaine de milliers à quelques millions d'années), et sur de grandes distances (quelques dizaines à centaines de kilomètres), en 3 D. On va considérer ici que les sédiments présents dans le bassin étudiés peuvent être définis à l'aide de lithologies de base Une lithologie correspond à une classe 20 granulométrique et une nature de grains On distingue ainsi par exemple des sables grossiers, des argiles fines, des carbonates moyens, etc. Comme le montre l'organigramme de la Fig 1, la modélisation stratigraphique du bassin est obtenu itérativement du passé jusqu'au présent (ou au moins jusqu'à un temps plus récent), en une série de pas de temps constant, de l'ordre de la dizaine à la centaine de 25 milliers d'années Au cours de chaque pas de temps, le modèle stratigraphique calcule la déformation du bassin et les apports de chaque lithologie, puis quantifie le transport de chacune des lithologies L'application de la conservation de la masse permet d'obtenir une nouvelle géographie, et termine ainsi le pas de temps. Dans une première étape, on définit l'espace disponible aux sédiments Ensuite, on 30 introduit les sédiments, et on simule la production sédimentaire dans le bassin Puis on répartit ces sédiments en accord avec des lois de macro transport qui seront définis dans la suite de la description en limitant l'érosion à la couche d'altération. I) La création de l'espace disponible L'espace disponible pour le remplissage sédimentaire du bassin, est la somme de 5 l'eustatisme et de la subsidence On a choisi de le définir point par point, soit à partir de courbes eustatiques et de cartes de subsidence, soit directement à partir de cartes d'accommodation, sans utiliser donc de modèles physiques reliant l'eustatisme au climat, ou la subsidence à des processus tectoniques, isostatiques ou thermiques. Bien que ne modifiant pas l'espace disponible, la compaction influe sur le 10 remplissage sédimentaire en modifiant l'épaisseur des couches sédimentaires Pour rendre compte de la compaction mécanique, on a choisi de relier directement la porosité des sédiments à l'enfouissement maximal atteint par eux, par une relation de forme exponentielle qui permet d'obtenir une bonne approximation de la compaction telle que définie par X Parasnis, D S, 1960, The compaction of sediments and its bearing on some geophysical problems Journal of the Royal Astronomical Society, v 3, n 1, pp 1-28. Pour tenir compte du cas de sédiments multi-lithologiques, on a considéré que chaque lithologie a une porosité indépendante des autres, ce qui revient à assimiler les couches sédimentaires composées d'un mélange de plusieurs lithologies telles que le sable 20 et l'argile à une superposition de multiples sous-couches composées de lithologies pures. La porosité liée à chaque lithologie est ainsi traitée individuellement en fonction de l'enfouissement maximal atteint par la couche sédimentaire. pour chaquelithologiei = (I)i {Xr,i + (i Fo,i i l', porosité de la lithologie i (en %) Dj porosité résiduelle de la lithologie i (en %) avec (Do porosité de dépôt de la lithologie i (en %) z enfouissement maximale atteint par la couche sédimentaire étudiée (en m) z enfouissement de référence de la lithologie i (en m) Cette définition d'une porosité liée aux lithologies permet en particulier de simuler 25 le transport individuel de chaque lithologie et d'en déduire les conséquences sur la porosité des couches sédimentaires. II) L'introduction et la production de sédiments L'espace disponible ayant été créé, la deuxième étape de formation du modèle recherché consiste à introduire des sédiments dans le bassin, ou à les produire à l'intérieur du domaine marin. a) L'introduction de sédiments dans le domaine d'étude Les flux de sédiments aux frontières du domaine étudié représentent géologiquement l'apport sédimentaire responsable du remplissage du bassin Ils peuvent être perçus physiquement en termes de conditions limites du flux de sédiments. On fixe ces conditions limites soit en imposant la valeur exacte du flux en un 10 secteur de la frontière, soit encore en imposant au flux de sédiments d'avoir une évolution continue Le premier cas représente une zone d'apport imposée par des conditions externes au bassin tels que le débouché d'un fleuve drainant un bassin versant extérieur au bassin. Le deuxième cas représente une zone libre le long de laquelle le flux de sédiments est régi par des paramètres physiques internes au bassin tels que les caractéristiques des vagues. 15 Lors de la simulation d'un bassin, il est possible de combiner ces deux types de conditions limites, en distinguant par exemple une zone continentale o le flux est imposé par les apports fluviaux externes et une zone marine o le flux est défini par les lois de transport interne. b) La production sédimentaire à l'intérieur du domaine d'étude Les sédiments peuvent également être produits à l'intérieur du bassin, et plus particulièrement dans le cas de sédiments carbonatés Pour cette phase de la construction du modèle selon l'invention, on a choisi de simuler cette production à l'aide d'une formule empirique reliant le taux de production en chaque point du bassin à la bathymétrie et aux flux de sédiments telle que définie par exemple par: Lawrence et al, 1990, Stratigraphic simulation of sedimentary basins: concepts and calibration in Amer Assoc Petrol Geol Bull, 74, 3, 273-295. pour chaque lithologiei Pl = Pj i Fj j, taux de production de la lithologie i (en m/s) avec c Po J taux de production maximal de la lithologie i (en m/s) Fl OF si b 5 O lBi influence de la bathymétrie (adimensionnelle) influence de la bathymétrie Bi si Ob bi r b bathymétrie (en m) lbt b' seuil bathymétrique (en m) le-b-' si b>bi lA coefficientd'atténuation (en m) Fji sensibilité de la lithologie i au flux de la lithologie j I ifunedlaltooij jtsi Q n S J Q flux de la lithologie j (en m 2 Is) inluencedela holoie i F =I Fjj = e-jj(QJ-ji) si Qj >Sji Sj, seuildedébutd'inhibition (enm 2/s) lt. yt pr 1 coefficient de sensibilité (en m-2 s) Afin de prendre en compte l'érosion des récifs lorsque ceux-ci sont émergés, le processus d'altération mécanique de ces récifs a été simulé en considérant que tout sédiment récifal situé dans la couche d'altération en domaine continental est transformé en 5 sédiment "bioclastique" Les lithologies récifales et bioclastiques sont transportées dans le bassin à l'instar des sédiments siliciclastiques. I) Le transport des sédiments a) Le transport sur le long terme Les processus à long terme correspondent au transport permanent des sédiments, 10 induits par les fleuves, les courants marins (figure 2), ou l'action lente de la gravité (reptation, etc) Dans le cadre de la méthode, nous considérons que ces processus peuvent être modélisés par un équation de diffusion, contrainte par le flux d'eau et par une disponibilité réduite. Nous commençons par définir la valeur du flux d'eau Qeau en chaque point M du 15 bassin Ce calcul se fait de l'amont vers l'aval, en supposant que toute l'eau arrivant en un point donné du bassin est redistribué vers les voisins en aval de ce point, au prorata de la pente; le voisin ayant la pente la plus forte recevra ainsi plus d'eau que les autres. Dans un deuxième temps, nous calculons la capacité de transport long terme en tout point M du bassin Cette capacité de transport est définie par un modèle diffusif contraint 20 par le flux d'eau, c'est-à-dire que la capacité de transport de chaque lithologie est définie par la formule suivante: Qmax, litho = (Kgravité + K fluviale Qeau)grad z avec: Qmax, litho, la capacité de transport de la lithologie " litho " en M; Qeau, le flux d'eau s'écoulant à la surface du sol en M; Kgravité, le coefficient de diffusion lié aux processus gravitaires permanents (reptation du sol, etc); Kfluviale, le coefficient de diffusion lié au transport fluvial et océanique permanent; z, l'altitude du sol en M. Pour définir le flux réel Qlitho de chaque lithologie transporté au niveau du point M, nous parcourons le bassin d'amont en aval En chaque point M de ce parcours, nous connaissons ainsi le flux réel de chaque lithologie fournit par l'amont (étant donné que 10 nous avons précédemment calculé ce qui se passait en amont du point M). Nous considérons de plus que localement, l'érosion ne peut excéder une valeur maximale E, définie par la teneur en lithologie dans le sol (un récif carbonaté s'érodera ainsi moins vite qu'une argile indurée, qui s'érodera moins vite qu'un grès, etc), et par la bathymétrie (l'érosion étant différente en milieu aérien et sous-marin). Nous calculons le flux réel " sortant " en aval du point M en considérant que ce flux ne peut pas dépasser la capacité de transport maximale en aval de M, et ne peut pas dépasser la disponibilité maximale en sédiment, définie comme étant la somme du flux " entrant " en amont et de l'érosion maximale Nous avons ainsi: Qlitho = min (Qlitho amont + Elitho; Qmax,litho aval) Nous supposons de plus que le flux reste diffusif, c'est-à-dire qu'il suit la plus grande pente Nous pouvons définir ainsi un coefficient pondérateur A,, compris entre O et 1, tel que: Qlitho = Qmax, litho = (Kgravité + Kfluviale Qeau)gra z b) Le transport sur le court terme Les processus court termes correspondent (figure 3) au transport des sédiments lors de phénomènes épisodiques mais intenses, tels que les pluies importantes induisant des crues fluviales et des courants hyperpycnaux en domaine marin Dans le cadre de cette invention, nous considérons que ces processus peuvent être modélisés par un équation de diffusion améliorée, contrainte par le flux d'eau, par une disponibilité réduite, et par la vitesse de l'écoulement de l'eau. Un tel transport épisodique de sédiment a depuis fort longtemps été très étudié pour des problèmes d'ingénierie fluviale, côtière et marine Le but de ces études est 5 généralement de déterminer la capacité de transport d'un courant d'eau, afin de savoir si les installations humaines, telles que des maisons, des ponts, des ports ou des plate-formes marines, résisteront à un tel courant On admet généralement qu'un cours d'eau peut être caractérisé par trois grandeurs principales: sa vitesse u, sa profondeur h et sa largeur w. On admet de plus généralement que les équations de Saint-Venant, correspondant à la 10 formulation pour un fluide de l'équation de Newton (somme des forces = accélération), permet de déterminer ces caractéristiques en tout point d'un écoulement Les équations de Saint Venant se présentent en 2 D sous la forme simplifiée suivante: db h h auh conservation du volume d'eau d = + = O dt dt ax I d uh auh au 2 h 2 quantité de mouvement et bilan des forces d u + = g h S-Cu dt at 3 x o S représente la pente de la surface du sol, et C le coefficient de friction de l'eau contre 15 le sol. La résolution de ces équations sur de très grandes durées (dizaines à centaines de milliers d'années), et sur de très grandes distances (dizaines à centaines de kilomètres), posent des problèmes numériques énormes pour une application pratique Le principal problème est le temps de calcul nécessaire pour résoudre ces équations Il peut se chiffrer 20 en jours pour un pas de temps Une simulation stratigraphique durerait ainsi des semaines. Pour lever cet obstacle rédhibitoire, la méthode proposée permet de simplifier grandement le calcul de la vitesse de l'eau. A cet effet, nous considérons dans un premier temps qu'un régime d'équilibre existe virtuellement Ce régime d'équilibre est défini par la solution permanente uniforme 25 des équations ci-dessus Nous pouvons ainsi définir une vitesse d'équilibre Ue, une hauteur d'équilibre He, et une distance d'équilibre Le, par: Ue =(Qeau gs) C He_=Q eau =C Qeau 2 I ue g S Le= Qeau He 3 Ue C 3 C Nous considérons dans un second temps que la vitesse d'un cours d'eau essaye en tout point de rejoindre sa vitesse d'équilibre Cette tendance d'évolution est supposée suivre un modèle exponentiel: a U(M) Ue i 34 Le o t la distance le long de l'écoulement. Dans le cas o la vitesse et la distance d'équilibre restent sensiblement uniforme entre un point de référence et le point M, la relation ci-dessus aboutit au modèle exponentiel suivant: U(M) = Ue +(Uo -Ue) e-/Le o Uo représente la vitesse au niveau d'un point de référence de l'écoulement. Nous considérons en outre que la vitesse de l'eau agit comme un pondérateur P(U) du modèle de vitesse, suivant la relation suivante: Qmax, litho (Kgravité + Kfiuviale Qeau P(U))grad z _ c,litho avec P(U) = U o S U >$U 2 c,litho -U 2 3/2 et P(U)=O si U Nous appliquons finalement la règle de disponibilité donnée plus haut au sujet du transport long terme pour contraindre et quantifier le flux réel de sédiment transporté par un courant épisodique d'eau. c) Le transport catastrophique Les processus catastrophiques correspondent au transport en masse de sédiments suite à des éboulements de falaise ou de pente Dans un premier temps, nous déterminons les zones instables en appliquant un critère de pente critique Cette pente critique varie localement et dépend de la nature lithologique des sédiments dans le sol. Toute zone instable est ensuite déplacée en masse, en appliquantl'équation de 10 Newton: l'accélération de l'éboulement est égale à la somme des forces s'appliquant sur l'éboulement (la gravité et la friction) Cette équation nous permet de mouvoir chaque éboulement à travers le bassin Nous considérons que l'éboulement se fige dès que sa vitesse décroît en dessous d'un seuil critique. d) Résolution numérique La réalisation du modèle selon l'invention comporte la résolution numérique des équations de transport basée sur une discrétisation spatiale du bassin étudié, et une discrétisation temporelle de la formation. Le bassin étudié est décomposé en mailles carrées de taille constante, tandis que le remplissage de ce maillage est simulé au niveau d'une succession de temps de calcul, 20 séparés par un pas de temps constant La largeur des mailles est par exemple de l'ordre de 1 à 10 kilomètres en fonction du bassin simulé, tandis que le pas de temps est de l'ordre de 000 ans. Les équations de transport sont ensuite résolues à l'aide d'un schéma numérique explicite o les paramètres physiques tels que le flux ou l'altitude au temps (n+ 1) sont 25 exprimés en fonction de paramètres mesurées au temps (n) et dont les valeurs sont ainsi connues. Dans le cas monolithologique en deux dimensions, cette résolution s'écrit ainsi de la manière suivante. ooh h+n ( 1) Q =-K A en chaquemaillei Qin+l K dx Ah _ a Q hln+') h(n) Q(n+l) Q&sn+l) ( 2) Q > en chaque maille i + dt dt dx d(n)_, K(n, ,hln)_ -n soit h In+l) = h(n) + dt K h, -(n) hn i, h(?) altitude de la maille i et au temps N (en m) K>n) coefficient de diffusion au niveau de la maille i (en m 21 s) avec Q) flux de sédiments entre les mailles i et i+l (en m 2 /s) d pas de temps (en s) dx largeur des mailles (en m) Cette écriture assez simple devient cependant plus complexe lors du passage en multi-lithologie et en 3 D, et lors de la prise en compte de la restriction du transport des sédiments à la couche altérée Dans un tel cas, on écrit que: h max, soit =-_V, max, soit h (n+l) h (n) V 1 dt La résolution explicite par volumes finis est la méthode de calcul la plus rapide à mettre en oeuvre Elle fournit de plus des résultats très précis si on reste dans le domaine de stabilité de l'algorithme, ce qui nécessite des pas de temps de calcul internes très petits, de 10 l'ordre du siècle au millénaire. L'application du modèle aux données d'entrée formées à partir des données de champ, permet de simuler le transport des sédiments sur tout le bassin étudié La validité de la simulation est ensuite testée en confrontant les résultats fournis par le modèle aux données collectées sur le terrain et principalement aux épaisseurs de sédiments et aux 15 facies observés sur les diagraphies de puits En cas de désaccord, on cherche par inversion le jeu de paramètres d'entrée du modèle: espace disponible et quantité de sédiments transportés, tel que l'écart entre les résultats obtenus par ces paramètres et les contraintes imposées soit minimale, tel que schématisé ci-après modèle direct =ô S = M( p) fonction écart = E(p) = Ils cll p ensemble des paramètres d'entrée du modèle M ensemble des équations régissant le modèle S ensemble des résultats en sortie du modèle avec C ensemble des contraintes géologiques E fonction d 'écart entre les résultats et les contraintes lIx mesure de l'écart entre les résultats et les contraintes inversion = trouver p tel que = E(j) soit le minimum absolu de la fonction E Vp, E(p) E Dans le cadre de la présente méthode, l'inversion a ainsi pour but de quantifier en particulier les valeurs de l'accommodation, des apports sédimentaires et des coefficients de diffusion afin d'obtenir une simulation dont les géométries des corps sédimentaires, les 5 épaisseurs et les faciès mesurés au droit des puits de calage soient le plus près possible des contraintes géologiques. On utilise une méthode d'inversion de type essais et erreurs On définit un jeu de paramètres initial On calcule la solution du modèle liée à ces paramètres, et on modifie la valeur des paramètres en fonction de l'écart entre la solution et les contraintes géologiques. 10 On continue ce procédé jusqu'au moment o l'écart devient suffisamment faible. L'optimum de concordance étant atteint, la simulation conduit à des données quantitatives sur la géométrie du bassin et sur les lithologies des unités sédimentaires Elle permet aussi de vérifier la cohérence du schéma de corrélation des puits. IV) Validation industrielle du modèle obtenu: La validité du modèle ainsi constitué a été testée sur des cas théoriques simples à deux ou trois dimensions On a supposé constants par exemple la vitesse de subsidence, le taux d'apport sédimentaire et la diffusivité des sédiments En définissant l'eustatisme suivant un ou deux ordres de cyclicités, on a ainsi pu montrer à l'aide de ces cas simples mais réalistes que le modèle diffusif permet de restituer les implications des concepts de la 20 stratigraphie génétique, telles que le partitionnement volumétrique des sédiments au sein d'une unité génétique et la distorsion des unités génétiques La cohérence géologique des résultats de ces simulations permet de valider empiriquement l'utilisation de l'équation de diffusion multi-lithologique sur de grandes échelles de temps et d'espace. Le modèle a été testé aussi sur des cas réels actuels et anciens, dans différents contextes sédimentaires, dans le cadre de thèses et de projets en collaboration avec des compagnies pétrolières Les tests ont notamment été effectués sur: ( 1) des systèmes siliciclastiques (sable, argile, ) dans des environnements deltaque à marin peu profond (bathymétrie Ces applications ont permis de montrer que le modèle diffusif présenté dans le cadre de cette invention permet d'obtenir des simulations très précises, avec un écart moyen des épaisseurs cumulées des sédiments de l'ordre de 5 mètres à 25 kilomètres des 15 puits de contraintes (pour une formation ayant une épaisseur moyenne de l'ordre de 100 mètres), et un écart entre la position des rivages simulés et observés inférieur à 10 kilomètres (la taille des mailles de calcul étant de 10 kilomètres). Du fait de cette restitution géométrique et faciologique, et de l'accès à la physique à grandes échelles des processus sédimentaires, le modèle diffusif selon l'invention permet 20 de nouveau un retour sur les bases de données géologiques en affinant par exemple la valeur des bathymétries et des flux de sédiments, en restituant l'évolution de l'accommodation au cours du temps, et en confirmant ou en infirmant le choix d'un schéma de corrélation. REVENDICATIONS 1) Méthode pour constituer un modèle déterministe de type diffusif permettant de simuler le remplissage multi-lithologique d'un bassin sédimentaire, comportant, à partir de 5 données de terrain connues relatives à l'architecture du bassin et à des données de mesure telles que des données de diagraphie de puits ou des données sismiques, la constitution d'un ensemble de données d'entrée portant sur un espace disponible par subsidence et eustatisme, sur l'apport de sédiments fluviatiles ou marins et leur transport, et sur des paramètres physiques tels que des coefficients de diffusion des différentes lithologies, par 10 un processus itératif incluant le maillage du bassin en mailles de dimensions régulières, une modélisation suivant un schéma explicite en volumes finis avec des pas de temps constants, de façon à simuler le flux de chaque lithologie déposée sur chaque maille, une comparaison des résultats de la simulation aux données de terrain et la modification par inversion et de proche en proche des données d'entrée, caractérisée en ce que, à chaque pas 15 de temps et pour chaque lithologie: on modélise par des équations déterministes respectivement: a) la moyenne des processus de transport agissant de manière permanente et continue sur le long terme; b) la moyenne des processus de transport agissant de manière permanente et 20 continue sur le court terme; c) la moyenne des processus de transport agissant de manière catastrophique sur le très court terme; et on détermine la stratigraphie résultante du bassin en tenant compte de la conservation des masses. 2) Méthode selon la revendication 1, caractérisée en ce que, à chaque pas de temps et pour chaque lithologie: on modélise les flux à long terme de sédiments transportés par l'eau en tenant compte d'un modèle de répartition des flux d'eau sur le bassin et de la capacité de transport des sédiments par de tels flux; on modélise les flux à court terme de sédiments transportés par l'eau en utilisant un modèle exponentiel pour estimer la vitesse de l'eau; on modélise les flux catastrophiques de sédiments transportés par l'eau dans des zones instables repérées en appliquant un critère de pente critique. 3) Méthode selon la revendication 1 ou 2, caractérisée en ce que l'on modélise les flux à long terme par la relation Qlitho = Qmax, litho = X (Kgravité + Kfluviale Qeau) grad z, O Qmax, litho est la capacité de transport de la lithologie "litho" en tout point du bassin; Qeau, le flux d'eau s'écoulant à la surface du sol en ce point; Kgravité, le coefficient de diffusion lié aux processus gravitaires permanents; Kfluviale, le coefficient de diffusion lié au transport fluvial et océanique permanent; z, l'altitude du sol en ce point; et A, un coefficient de pondération. 4) Méthode selon la revendication 1 ou 2, caractérisée en ce que l'on modélise les 15 flux à court terme de sédiment transporté par l'eau en appliquant pour la vitesse de l'eau le modèle exponentiel: U(M) = Ue +(Uo -Ue) e /Le O Uo représente la vitesse au niveau d'un point de référence de l'écoulement, Ue, une vitesse d'équilibre de l'eau, et , une distance le long de l'écoulement, les flux étant 20 estimés par la relation: Q max, litho = (Kgravité + K fluviale Qeau P(U))grad z avec P(U) ( clitho 32 si U> U avec P(U) = clitho si U > Uclitho et P(U) = O si U