FR 2808909 A1 20011116 FR 0006113 A 20000511 Domaine technique L'invention se rapporte au domaine technique du Risk Management Industriel et plus spécifiquement - aux processus statistiques/probabilistes de prédiction de la loi de probabilité de la distribution statistique des variations d'un indicateur de comportement futur et de l'optimisation des paramètres de rendement d'une multitude de systèmes industriels locaux distants soumis à un environnement global commun aléatoire évolutif, - du type par simulation, c'est-à-dire fonctionnant par échantillonnage préalable de scénarios probables des états futurs des variables d'environnement. L'invention se rapporte spécifiquement au cas - où une multitude de sites industriels décentralisés ont un comportement incertain du fait qu'ils sont soumis à un environnement global commun aléatoire évolutif (que l'on peut dans la majorité des cas qualifier d'environnement de risque), - et où cet environnement aléatoire évolutif est néanmoins quantifiable, c'est-à-dire qu'il peut être majoritairement décrit numériquement par des variables d'état communes (au sens défini ci-après) à l'ensemble des sites industriels. Ces variables d'état sont en effet: - soit, communes à l'ensemble des sites industriels, - soit, lorsqu'elles sont en apparence spécifiques à un site industriel, quasi-commune du fait que leur comportement est dicté par des relations complexes avec d'autres variables d'état du modèle (non spécifiques au dit site), de telle sorte que leur loi de probabilité ne peut être décrite indépendamment de celle des autres variables d'état communes ou spécifiques d'autres sites. Par la suite, et notamment dans la description qui va suivre ainsi que dans le texte des revendications, on désignera (par extension) du terme de variables d'état communes un ensemble de variables d'état de l'un des deux types décrits ci-dessus. Par la suite, on utilisera également indifféremment le terme d'environnement aléatoire ou d'environnement de risque. Par site industriel on entend, selon l'invention : tout lieu d'activité humaine ayant pour objet la production d'une valeur ajoutée (flux de production) - ou exprimable directement ou indirectement en terme de valeur ajoutée -,que celle-ci soit liée à l'industrie manufacturière, énergétique, agroalimentaire, financière, de transport, de santé publique ou d'exploitation des matières premières, etc. Par la suite on utilisera fréquemment comme exemple d'application de l'invention, et comme domaine d'analyse de l'art antérieur, le domaine financier. Car c'est dans ce domaine que l'application des processus de prédiction probabiliste de comportements aléatoires est le plus fréquent du fait de l'adoption de réglementations l'y contraignant (directives CAD du Conseil des Communautés Européennes, recommandations du Comité de Bâle , recommandation du Groupe des 30 ). Mais les enseignements de l'invention peuvent s'appliquer dans presque tous les domaines industriels. Si l'on prend comme exemple le domaine financier, on entend par risque quantifiable (environnement aléatoire ) les risques de marché, de crédit, de contrepartie, de liquidité... Par opposition, on peut considérer comme non quantifiables les risques légaux, de fraudes, médiatiques, L'invention concerne spécifiquement les cas - où le comportement des variables d'état de l'environnement est incertain ou incontrôlable, en sorte que celles-ci constituent des facteurs de risque aléatoires qui peuvent causer une variation et notamment une baisse d'efficacité (ou un risque sur le comportement de chacun) des sites industriels décentralisés, pouvant se manifester notamment par des pertes (d'argent, de vie humaine, de production,...) ; - et où l'efficacité de chaque système local peut en outre être décrite par un indicateur numérique de rendement selon un modèle numérique systémique local (généralement du type paramétré) préalablement établi, dépendant pour chacun des systèmes, et à chaque instant, de tout ou partie des variables d'état et de paramètres locaux d'action permettant de régler le niveau de fonctionnement du système industriel local. L'invention concerne les domaines industriels, aussi bien manufacturiers que financiers, publics ou privés, où il convient - d'estimer statistiquement pour chaque site industriel, afin de le circonscrire, un indicateur numérique de rendement aléatoire qui peut être de type soit industriel manufacturier (production, nuisance), soit financier (profit/perte) ou autre, et qui est lié à l'exercice d'une production de valeur ajoutée (production manufacturière, agricole ou financière...), - et d'optimiser l'indicateur numérique de rendement (ou de risque) de chaque site industriel, en réglant les paramètres locaux d'action de ce site industriel, - ceci, en tenant compte notamment de critères réglementaires imposant des limites à cet indicateur de rendement (ou de risque), Description de l'art antérieur La deuxième moitié du XXè" siècle a vu la prise de conscience des effets pervers du risque et ) par voie de conséquence de son coût et la naissance d'une nouvelle science : le Risk Management . Le Risk Management est la science - qui étudie les facteurs quantifiables ou non qui menacent le rendement ou le comportement d'un système industriel ou privé, - et qui apporte des solutions stratégiques (organisationnelles et de procédure) et tactiques (méthodes quantitatives de gestion ... ) afin de combattre ces facteurs. L'art antérieur connaît des méthodes de Risk Management de systèmes industriels. Ainsi le brevet US 6021402 décrit une méthode de gestion des risques d'un réseau de systèmes industriels de production électrique. Le brevet US 5930762 décrit un système de contrôle de risque d'un dispositif dont le comportement est lié â des paramètres multiples d'environnement. L'art antérieur connaît principalement cinq méthodes différentes de mesure et/ou prise de décision relative au risque des systèmes aléatoires - les processus à arbre décisionnel, - les processus statistiques/probabilistes, (ex : Variance-covariance, Historique, Monte- Carlo), - les méthodes variationnelles (ex : La valeur de 1 point de base (PVO I )), - les méthodes par scénarios extrêmes (ex : Stress testing , Scénario catastrophe ), - les méthodes d'intelligence artificielle. a) Les processus à arbre décisionnel sont utilisés dans les procédés de prise de décision séquentielle. Un diagramme qui ressemble aux branches d'un arbre montre l'ensemble des événements possibles et l'enchaînement des options de décision qui sont le résultat de décisions précédentes. Cette méthode à arbre décisionnel peut être utilisée uniquement dans les cas très simples où les événements anticipés et la probabilité de chaque événement sont déjà connus. Les calculs numériques mis en oeuvre par cette méthode sont simples et peuvent d'ailleurs souvent être implémentés sans calculateur. Mais elle est insuffisante pour décrire la structure de risque de systèmes aléatoires complexes. b) Dans de nombreux cas, les événements anticipés dépendent de différentes variables d'état dont la valeur ne peut pas être connue exactement. Ce type de problèmes est habituellement analysé par un processus statistique/probabiliste. On appellera par la suite processus statistique/probabiliste d'estimation du comportement futur d'un système soumis à un environnement aléatoire pouvant être décrit à l'aide de variables d'état, une méthode consistant - à effectuer une estimation statistique de la loi de probabilité de la distribution statistique des variations d'un indicateur numérique de rendement de ce système, - à partir d'un modèle de distribution de probabilité jointe des variables d'état de l'environnement (facteurs de risque) préalablement établie. Ceci est en opposition avec le processus du type à arbre décisionnel (ou toute autre mesure du risque) dont les résultats sont des valeurs uniques. c) Les méthodes variationnelles quantifient la conséquence sur le comportement d'un système de la variation d'un ou plusieurs paramètres d'environnement (variables d'état). Ainsi dans le domaine financier, la méthode dite valeur de 1 point de base (PVOI) quantifie la sensibilité d'un portefeuille de titres financiers à un mouvement parallèle et unilatéral de la courbe des taux d'intérêt vers le haut ou vers le bas pour un décalage de 0,01% (un point de base ). On notera que dans ce cas la probabilité d'occurrence et l'amplitude du mouvement ne sont pas appréhendés. Ceci est insuffisant pour décrire la réalité du comportement de risque d'un système aléatoire complexe. d) Les méthodes par scénarios extrêmes visent à estimer le résultat sur un système de l'occurrence de scénarios particuliers. Ainsi, dans le domaine financier, la méthode stress testing permet de faire varier les taux d'intérêt à plusieurs endroits de la courbe, par exemple : concomitamment à la hausse pour une des échéances et à la baisse pour une autre. On utilise cette méthode à des fins de simulation de ce qui se passerait suite à des scénarios catastrophes (prévision de ce qui est sensé n'arriver que très rarement). On notera que, dans ce cas, les aspects temps et probabilité n'interviennent pas. Ceci est également insuffisant pour décrire la réalité du comportement de risque d'un système aléatoire complexe. e) Dans les cas différents où l'on souhaite non pas estimer un risque, mais plutôt prédire une évolution aléatoire telle que par exemple le prix du pétrole ou la variation de prix d'une action, trop de facteurs entrent en jeu qui sont eux-mêmes trop incertains et difficiles à modéliser avec des outils mathématiques. Pour aborder ces problèmes de prédiction, on utilise des méthodes dites d'intelligence artificielle à réseau de neurones et des systèmes experts . Ces techniques d'intelligence artificielle ont la capacité de raisonner sur la base d'informations floues ou incomplètes. Récemment, ces technologies ont été utilisées dans le domaine de la prédiction géologique de production et celui des prédictions boursières. Mais elles sont inadaptées pour des applications générales d'analyse de risque de systèmes industriels, eu égard à des critères (notamment réglementaires) communs. Ainsi le brevet US 5521813 décrit un système de prédiction des conséquences des variations climatiques sur un système productif, utilisant une méthode à régression multiple et réseau de neurones. La présente invention se rapporte spécifiquement au domaine technologique des processus de statistiques/probabilistes qui visent à effectuer une estimation statistique de la loi de probabilité de la distribution des variations d'un indicateur numérique de rendement de ce système, à partir d'un modèle de distribution de probabilité jointe des variables d'état de l'environnement. Le domaine de l'art antérieur le plus courant dans lequel les processus statistiques/probabilistes d'analyse de risque sont utilisés est celui de la finance et plus spécialement celui des risques de marché. Dans cet environnement global de l'industrie financière, le risque de marché est le risque auquel est confronté l'investisseur en raison de sa méconnaissance de l'évolution future des variables d'état communes fondamentales du marché que sont : les cours des titres, les taux d'intérêts, les taux de change.... Ces variables de marché, également appelées facteurs de risque, déterminent le prix des titres, des actifs contingents, des portefeuilles, ...etc. Les risques associés à différents produits financiers peuvent être évalués par quelques nombres caractéristiques : duration et convexité pour les obligations, bêta et performance espérée pour les actions, valeur-temps, volatilité, delta, gamma, thêta, véga (dérivée par rapport à la volatilité), etc... pour les options. Mais ces différentes approches du risque dites par sensibilité - sont spécifiques à un type de produit financier, - s'appuient sur des hypothèses de distribution spécifiques, - ne sont pas faciles à interpréter, - sous-estiment souvent le risque global de perte sur un portefeuille qui est pourtant celui qui importe pour l'investisseur, - n'évaluent pas le niveau de probabilité de survenance de ces risques, - et ne prennent pas l'horizon temporel en considération. Il est préférable d'avoir une approche unifiée probabiliste et temporelle du risque. Ainsi, l'art antérieur utilise la méthode dite Value at Risk (VaR) comme indicateur de risque d'un système financier soumis à un environnement aléatoire de marché. Pour un actif financier (ou un portefeuille financier), la VaR est un montant tel que la perte encourue sur cet actif (ou ce portefeuille) durant un certain intervalle de temps ne dépassera VaR qu'avec une probabilité de niveau (faible) donné. Le calcul de la VaR nécessite la connaissance de la loi de probabilité de la perte subie (indicateur numérique de rendement). Le problème fondamental des processus statistiques/probabilistes d'analyse de risque (et notamment des processus de calcul de la VaR) est donc de prévoir la distribution de loi de probabilité de la distribution statistique des variations d'un indicateur numérique de rendement (perte dans le cas de la VaR). Les processus statistiques/probabilistes d'estimation statistique de loi de probabilité de la distribution des variations d'un indicateur numérique de rendement d'un système sont principalement de deux types - Ceux, que nous appellerons par la suite par simulation , # qui procèdent par génération d'une multitude d'échantillons statistiques multi-variés (de n-uplets) de description de valeurs probables conjointes des variables d'état (paramètres de risques), # qui aboutissent à la description de l'histogramme point par point de la structure de la loi de probabilité de la distribution statistique des variations de l'indicateur numérique de rendement du système, (c'est-à-dire dans son intégralité), # et qui permettent notamment de calculer la VaR de manière générale, sans hypothèse préalable sur la loi de probabilité de cette distribution statistique des variations de l'indicateur numérique de rendement ; - Par opposition à ceux, que nous appellerons par la suite analytiques (et que l'on pourrait également qualifier de paramétriques, normatifs ou synthétiques,...), # qui ne fournissent que quelques paramètres (en nombre fixé) descriptifs de la loi de probabilité de la distribution statistique des variations de l'indicateur numérique de rendement du système, # qui imposent une forme particulière au modèle d'évaluation de l'indicateur numérique de rendement du système en fonction des variables d'état (paramètres de risques), (notamment une hypothèse de linéarité, quadraticité, etc.), # et qui aboutissent à des hypothèses limitatives de description paramétrique spécifique de la loi de probabilité de la distribution statistique des variations de l'indicateur numérique de rendement. Ainsi, l'art antérieur connaît principalement quatre méthodes classiques d'estimation de la VaR (et plus généralement d'estimation statistique de loi de probabilité de la distribution statistique des variations d'un indicateur numérique de rendement d'un système soumis à un environnement aléatoire) - deux méthodes statistiques/probabilistes analytiques # la technique delta normal ou technique de la matrice des variance-covariance estimées , qui suppose le rendement linéaire ; # la technique delta-gamma , qui tient compte d'un terme quadratique dans le calcul du rendement du système ; - et deux méthodes de type statistique/probabiliste par simulation # la technique de simulation de Monte-Carlo ; # la technique de l'analyse historique. L'invention se rapporte exclusivement aux domaines techniques des procédés de type statistiques / probabilistes par simulation . Pour le calcul de la VaR d'un portefeuille financier, les entrées sont classiquement - les modèles mathématiques d'évaluation de prix des différents actifs, ex # le modèle de prix par combinaison linéaire de variation du prix d'un actif en fonction des variations des facteurs de risque (tel que le modèle à indice simple de Sharpe pour les actions, qui postule que la variation du prix d'une action est une fonction du premier degré de la variation d'un indice général du marché), # les modèles de prix non linéaire, tel que le modèle en développement de Taylor indépendant du temps, # et les modèles de prix fonction du temps, notamment du type modèle aléatoire stationnaire, tel que le modèle de Black et Sholes pour l'évaluation des options ; - et les historiques, c'est-à-dire les observations des prix des différents actifs pour un certain nombre de périodes passées. La méthode delta normal (ou linéaire) dite de la matrice des variances-covariances estimées est celle qui a été proposée et développée par JP Morgan, avec son système RiskMetrix T#'1. Elle consiste en trois étapes - le relevé des facteurs de risque élémentaires en lesquels les actifs financiers du portefeuille peuvent être décomposés, - la répartition des flux financiers associés à ces facteurs de risque élémentaires en flux agrégés correspondant à des échéances standards, - le calcul effectif de la VaR. Cette méthode fait les hypothèses théoriques spécifiques limitatives suivantes - les facteurs de risque liés à chaque actif obéissent à une loi normale, - la performance des actifs dépend linéairement des facteurs de risque, - et par voie de conséquence la performance globale du portefeuille dépend linéairement des facteurs de risque. Le calcul de la VaR s'effectue alors à partir de la matrice M dite de variance-covariance , de dimension nxn dans le cas de n facteurs de risque, et constituée à partir des nx(rt+t)/2 variances et covariances des variations des facteurs de risque, calculées sur la base d'un historique passé. Si l'on appelle A;(`-' ") les sensibilités du portefeuille sur les différents facteurs de risque (X = X;('' -la VaR est donnée par la formule où z est le quantile de la distribution normale réduite pour la valeur de l'intervalle de confiance de la VaR retenue. En sorte que, selon la méthode de la matrice de variance- covariance, on calcule la structure du risque, dans un environnement à il facteurs de risques, en générant un ensemble de n' paramètres fixes de la matrice de variance-covariance, descriptifs des variances et covariances des variations des facteurs de risque. C'est en ce sens que la méthode de la matrice des variances-covariances est une méthode statistique/probabiliste analytique. L'hypothèse de linéarité retenue nécessite pour être mise en oeuvre des hypothèses strictes sur les sensibilités du portefeuille, la rendant inapplicable pour des portefeuilles incluant des risques dits gamma ou de convexité . La méthode dite delta-gamma est une forme plus élaborée de la méthode de la matrice des variance-covariance. Elle abandonne l'hypothèse de linéarité au profit d'une hypothèse quadratique. Elle incorpore un second ordre de sensibilité dans les analyses. Elle nécessite la mise en oeuvre de calculs mathématiques sophistiqués. Malgré sa supériorité vis-à-vis de la méthode delta normal , la méthode delta-gamma est inadaptée dans le cas général, notamment pour évaluer le risque de portefeuilles incluant des options arrivant à maturité ou des dérivés exotiques tels que des options knock-out . Au contraire, selon la méthode de simulation de Monte-Carlo, prônée notamment par Bankers Trust, avec son système RaRoc 2020 Tm: - On fait le choix d'une famille de modèles d'évaluation pour les différents facteurs de risque ou les actifs du portefeuille, justifiés à partir d'observations historiques. - A partir des historiques des actifs en présence et des facteurs de risque, on estime la distribution des différents facteurs de risque, ainsi que des paramètres qui lui sont associés : moyennes et matrice de variance-covariance..., # ainsi que les relations liant les prix des actifs et les facteurs de risque. - Pour l'ensemble des différents facteurs de risque et sur la base des modèles probabilistes retenus et des paramètres de distribution déterminés, on émet grâce à un générateur pseudo-aléatoire, un grand nombre (quelques centaines à quelques milliers) d'échantillons pseudo-aléatoires extraits d'une distribution estimée des variations probables futures des facteurs de risque. - On introduit les valeurs des échantillons de chacune de ces simulations dans les relations de comportement des actifs en fonction des facteurs communs de risque. - On reconstitue pour chacun des échantillons de simulation la valeur correspondante du portefeuille. - On construit l'histogramme global de loi de probabilité de la distribution statistique des variations de valeur du portefeuille, sans hypothèse préalable sur cette dernière. C'est en ce sens que la méthode de simulation de Monte-Carlo est une méthode statistique/probabiliste par simulation. La méthode d'analyse historique, prônée notamment par Chase Manhattan, avec les systèmes Charisma Tm et Risk$ Tm est une simplification de la méthode Monte-Carlo. Selon cette méthode d'analyse historique, on émet la multitude d'échantillons des variations probables futures des facteurs de risque directement à partir de l'historique du marché. A titre d'exemple, la VaR historique à un jour d'un portefeuille peut être estimée à partir de cent échantillons de scénarios de la manière suivante - On collecte les valeurs reconstituées des facteurs de risque à partir de l'historique du marché sur les cent derniers jours de cotation. - Pour chaque jour, on détermine le pourcentage de variation de chacun des facteurs de risque. Cela permet de simuler cent échantillons multi-valués de valeurs probables de chacun des facteurs de risque. - On introduit les valeurs des échantillons correspondant à chacune de ces observations historiques dans les relations de comportement des actifs en fonction des facteurs communs de risque. - On reconstitue pour chacun des échantillons historiques la valeur correspondante du portefeuille. - On construit l'histogramme global de loi de probabilité de la distribution statistique des variations de valeur du portefeuille, sans hypothèse préalable sur cette dernière. C'est en ce sens que la méthode par analyse historique est également une méthode statistique/probabiliste par simulation. Le principe de l'analyse de risque utilisant une méthode statistique/probabiliste par simulation est connu de l'art antérieur et utilisé dans de nombreux domaines aussi bien financiers qu'industriels. La présente invention qui comporte (en combinaison) la mise en aeuvre d'une méthode de ce type ne revendique nullement ce moyen pris seul. On peut citer le brevet FR2769703 qui décrit un procédé pour déterminer selon la méthode Monte-Carlo l'incertitude de mesure d'une machine â mesurer en coordonnées. Le brevet EP0837408 décrit un procédé de simulation de pulvérisation par méthode Monte-Carlo. Le brevet W09905685 décrit un système automatique de surveillance industrielle par méthode Monte-Carlo. Le brevet W09859322 décrit un système de génération des valeurs de pixels d'une image par méthode Monte-Carlo. Le brevet W09852132 décrit un système d'affectation de ressources organisationnelles et d'évaluation des coûts et prix des ressources par méthode Monte-Carlo. Le brevet W09850835 décrit un contrôleur programmable pour superviser les activités d'équipements industriels par méthode Monte-Carlo. Le brevet W09836364 décrit un processus d'évaluation d'intégrale par méthode Monte-Carlo. Le brevet US5940810 décrit une méthode d'estimation de la valeur d'actifs financiers complexes par méthode Monte- Carlo. On constate que l'on peut répartir les méthodes statistiques/probabilistes d'estimation de loi de probabilité de la distribution statistique des variations d'un indicateur numérique de rendement d'un système soumis à un environnement aléatoire (et notamment les méthodes d'estimation de la VaR) en deux groupes distincts: - les méthodes analytiques selon lesquelles on génère un nombre fixé et petit de paramètres descriptifs du comportement statistique des paramètres de risque, et aboutissant à une description paramétrique de la loi de probabilité de la distribution statistique des variations de l'indicateur numérique de rendement d'un système, (c'est le cas notamment de la méthode de la matrice des variances-covariances de RiskMetrics TM delta normal ou delta-gamma ) ; - et les méthodes par simulation qui procèdent par échantillonnage de scénarios probables de variables d'environnement et selon lesquelles on génère une multitude (dont le nombre peut être choisi) de tirages répétés d'échantillons de valeurs futures probables plausibles des variables d'environnement, aboutissant à une description globale absolue de l'histogramme de loi de probabilité de la distribution statistique des variations d'un indicateur numérique de rendement d'un système (c'est le cas notamment de la méthode de Monte-Carlo et de la méthode d'analyse historique). L'invention se rapporte exclusivement aux procédés d'analyse de comportement - de type statistique/probabiliste par simulation , - effectuant un échantillonnage (tirage) de scénarios probables multi-valués des variables d'environnement . L'art antérieur sait, au niveau d'un site unique, mettre en oeuvre une méthode de type statistique/probabiliste par simulation , et notamment effectuer l'ensemble des acquisitions de données, calculs statistiques, etc... en vue de déterminer la loi de probabilité de la distribution statistique des variations d'un indicateur numérique de rendement d'un système lorsqu'il est soumis à un environnement aléatoire. L'art antérieur sait également calculer la VaR du portefeuille d'une institution financière seule, par méthode . Il est par ailleurs connu par l'art antérieur d'émettre répétitivement des valeurs numériques d'un site centralisé en direction d'une multitude de sites distants. C'est notamment le cas des crans de données financières de type Reuters ou Bloomberg. Un site central est relié à des sources de données variables et aléatoires d'environnements (cours d'actifs financiers, etc...). Il concentre ces données et les émet répétitivement de manière décentralisée en direction d'une multitude de sites distants (écrans d'opérateurs financiers). Mais on constatera que selon l'art antérieur, les données télétransmises sont des données réelles observables (valeurs de titres) ou sont quasi-observables, c'est-à-dire reconstituées par calcul (valeurs d'indices, taux de rendement internes...) et destinées à être utilisées comme telles. L'art antérieur ne connait pas le principe de la télétransmission systématique et répétitive d'échantillons probables multi-valués de valeurs futures possibles de variables d'environnement en vue de reconstituer un histogramme. Les méthodes de l'art antérieur telles qu'exposées ci-dessus ne permettent pas d'effectuer de manière satisfaisante la prédiction probabiliste du comportement futur et/ou l'optimisation des paramètres de rendement d'une multitude de systèmes industriels locaux distants soumis à un environnement global commun. Notamment elles ne permettent pas de calculer concomitamment dans de bonnes conditions la VaR d'une multitude de sites financiers distincts gérant des portefeuilles de titres. En effet D'une part, comme on l'a vu plus haut, seuls les processus statistiques/probabilistes qui visent à effectuer une estimation statistique de loi de probabilité de la distribution statistique des variations d'un indicateur numérique de comportement, donnent une mesure satisfaisante fiable et fidèle des risques. D'autre part, parmi les processus statistiques/probabilistes - A) La méthode analytique, et notamment la méthode de la matrice de variance-covariance estimée, est simple et nécessite peu de calcul. Elle peut d'ailleurs être mise en oeuvre conjointement au niveau de plusieurs sites distants. Ainsi, selon une variante de mise en oeuvre de l'art antérieur, le système RiskMetrics Tm fournit à partir d'un site centralisé, par Internet, et en direction d'une multitude de sites distants qui utilisent sa méthode, une émission répétitives d'un nombre fixé de paramètres de la matrice de variance-covariance. En sorte que le calcul global de la VaR est hiérarchisé entre le site central et les sites distants. Par contre les valeurs transmises ne sont pas des échantillons probables multi- valués de variables d'environnement mais de simples paramètres (voir plus haut). Le défaut de cette méthode est qu'elle nécessite des hypothèses de base très strictes (normalité conditionnelle, linéarité, stationnarité) qui ne sont pas satisfaites pour une grande partie des instruments financiers. En sorte que le calcul de la VaR selon la méthode de la matrice de variance-covariance produit des erreurs dans le cas de portefeuilles comprenant des instruments ayant un gamma , une convexité (voire des inversions de convexité ou des discontinuités) tels que les options, les notes structurées, les mortgage-backed securities ... Plus généralement, les VaR analytiques, y compris les plus perfectionnées, résument le comportement d'une structure de risque à l'émission d'un nombre restreint de paramètres, ses grecques (delta, gamma, véga, thêta, ... ). Ceci est en général insuffisant pour la description d'environnements ayant des facteurs de risque complexes. - B) La méthode classique par simulation, et notamment la technique de simulation de Monte-Carlo est considérée par la majorité des spécialistes comme la plus pure , la plus globale, la plus tournée vers l'avenir dans sa conception. Elle est incontestablement la plus juste théoriquement, la plus précise et elle s'applique à la totalité des instruments financiers. Cependant elle nécessite des choix complexes de modèles d'évaluation ainsi que des processus stochastiques de génération d'échantillons aléatoires correctement choisis. Ceci est très complexe étant donné la difficulté de l'évaluation de certains instruments et leurs constantes évolutions. Elle mobilise des équipes de spécialistes nombreuses pour être correctement implémentée sur un site. En outre elle nécessite pour être précise la génération d'un nombre très important de tirages pseudo-aléatoires d'échantillons. Ceci rend nécessaire la mise à disposition de systèmes informatiques de grande taille et coûteux. En résumé, selon la méthode Monte-Carlo de l'art antérieur, chaque site doit disposer d'un système informatique coûteux et d'une large et coûteuse équipe de mise en oeuvre. En sorte qu'elle est peu accessible à la majorité des institutions financières ou sites industriels qui pourtant souhaiteraient et auraient intérêt à la mettre en couvre pour réduire leur risque. - C) L'inconvénient majeur de la méthode de simulation par analyse historique est sa très grande sensibilité à la qualité des données. En effet, une ou quelques observations sortant de l'ordinaire (qu'elles soient vraiment exceptionnelles ou qu'une erreur soit à leur origine) influencent fortement la valeur de la VaR pendant une longue période (égale à la durée des historiques). On a dit que la théorie des valeurs extrêmes pouvait pallier cet inconvénient. Malheureusement, la masse des calculs que son application entraîne est prohibitive. En outre, les VaR historiques sont imprécises et doivent être complétées par des stress tests adéquats. Par conséquent, selon la méthode historique de l'art antérieur, chaque site doit disposer d'un système informatique coûteux et d'une large équipe de mise en couvre. En sorte qu'elle est également peu accessible à la majorité des institutions financières ou sites industriels qui pourtant souhaiteraient et auraient intérêt à la mettre en couvre pour réduire leur risque. Il est connu par l'art antérieur de diffuser un scénario le plus probable de variables d'environnement en direction de sites distants (ex : téléviseurs). C'est le cas notamment des prévisions de variables météorologiques (vents, températures ... ) en plusieurs points. Mais ceci ne permet pas d'avoir une vision probabiliste globale (loi de probabilité de la distribution statistique des variations d'un indicateur numérique de comportement) mais limitée à un petit nombre de paramètres au niveau de chacun des sites. En outre, les procédures de contrôle des risques requièrent une uniformisation des mesures de risque, donc une unité de méthode dans leur évaluation. En sorte que l'analyse de l'art antérieur fait ressortir le besoin d'une technique de prédiction probabiliste du comportement futur (et/ou d'optimisation des paramètres de rendement) d'une multitude de systèmes industriels locaux distants (notamment d'une multitude de sites financiers) soumis à un environnement global commun, donnant une mesure satisfaisante fiable et fidèle des risques, de type statistique/probabiliste par simulation c'est-à-dire fonctionnant par échantillonnage de scénarios probables de variables d'environnement, en utilisant globalement le minimum - de moyens humains, - et de moyens en matériel de calcul informatique. Résumé de l'invention Le Procédé statistique/probabiliste de simulation stochastique selon l'invention vise à effectuer de manière hiérarchisée mais conjointe (c'est-à-dire mutualisée) la prédiction probabiliste du comportement futur et/ou l'optimisation des paramètres de rendement d'une multitude de systèmes industriels locaux distants (SJ - 1 ..-")). Il s'applique au cas - où ces systèmes industriels locaux (Su) sont très distants les uns des autres, mais sont soumis à un environnement global commun (E) évolutif, pouvant être décrit numériquement par p variables d'état (X;(t)"-' ...P'), constituant les composantes d'un vecteur d'état (X(t) X((I)(1 - '"-Pl), - où les p variables d'état (X;(t)('- '...Pl), sont communes à la multitude des n (ii étant très grand 1) systèmes industriels locaux distants (S,('= '.."l ) décentralisés, - et où l'efficacité de chaque système local (S,,Iu= ' -.n) peut en outre être décrite par un indicateur numérique de rendement (lu""- ' ... nl) @ selon un modèle numérique systémique local paramétré (Du) préalablement établi (lu = Du(u I -"l (1u, X(t))), dépendant pour chacun des systèmes (Su) et â chaque instant (t) # de tout ou partie des variables d'état (X;(t)(' --n') du vecteur d'état (X(t)), # et d'un vecteur de q paramètres (réglables) locaux d'action (a,,= ..., @iuq)). Le procédé de simulation stochastique selon l'invention procède selon une méthode classique de type statistique/probabiliste par simulation qui consiste tout d'abord â effectuer un échantillonnage de scénarios probables des variables de l'environnement , puis à effectuer une estimation statistique de la loi de probabilité de la distribution des variations d'un indicateur numérique de rendement de chacun des systèmes industriels locaux (S") à partir de cet échantillonnage de description de la distribution de probabilité jointe des variables d'état de l'environnement préalablement établi. Mais le procédé selon l'invention vise en outre - à centraliser (mutualiser) et à mettre en commun au sein d'un site central les moyens lourds (en matériel informatique et en personnel) nécessaires à la génération de cet échantillonnage de scénarios probables des variables de l'environnement, - et à ne réaliser, au niveau de la multitude des systèmes industriels locaux, que des calculs simples et peu onéreux en matériel informatique et en personnel spécialisé. Dans sa forme la plus générale, le procédé de simulation stochastique selon l'invention est remarquable par la combinaison des deux moyens suivants: I a D'une part, au niveau d'un site centralisé émetteur (U), très distant (de plusieurs kilomètres, voire de plusieurs centaines de kilomètres) d'une majorité des st systèmes industriels locaux distants (S"(,= i...n>)@ 1 on génère répétitivement de manière centralisée à des instants (t,) successifs, une swap série de simulations probabilistes (W;,j@ ('=1...p,j_1...s,k=l...n) multi-valuée de description des variables d'état (X;(t)@'= @...p@) de l'environnement (E). Cette snap série de simulations probabilistes est, selon l'invention, constituée # soit d'une série de pseudo-états probables futurs des variables (X;(t)(' I -P)) de l'environnement (E), # soit d'une série de pseudo-variations probables futures de ces variables (Axi(t)(i = 1 ...nl), (ces séries étant simulées à un (ou plusieurs an, avec ni >_ 1) horizons de temps (tk(k 1-ni) > t0) ; 1 et on télétransmet répétitivement ces swap séries de simulations probabilistes (Wl,jk) multi-valuée de description des pseudo-états ou pseudo- variations probables futur(e)s des variables de l'environnement, ce de manière répartie décentralisée # du site émetteur (U), # en direction de l'ensemble de la multitude des systèmes locaux distants (Su). o D'autre part, de manière indépendante, au niveau de chacun des systèmes locaux (Su), majoritairement très distants du site centralisé émetteur (U), répétitivement on collecte lesdites smp séries de simulation probabiliste (W;jk) multi-valuées de description des pseudo-états ou pseudo-variations probables futur(e)s, 1 et on effectue un traitement numérique statistique local spécifique sur lesdites séries de simulation probabiliste (W,,@) reçues pour prédire le comportement futur probabiliste des systèmes industriels locaux (S"). Eventuellement, on règle le niveau du vecteur des paramètres locaux d'action (#,u) ...,@.uN) de chaque système local (Su) à une valeur optimale calculée en fonction des critères préétablis, afin notamment d'anticiper le comportement probable futur et l'efficacité de chaque système (S") soumis à l'environnement (E) au cours du ou des horizon(s) de temps (tk) retenu(s). Figures Sur les dessins - La figure 1 : décrit une méthode classique recommandée par l'invention pour la saisie de mesures observables de l'état d'un environnement et leur transmission par un fournisseur de ces données observables à un site émetteur de données probabilistes selon l'invention. - La figure 2 : décrit les principales étapes du processus de génération d'un échantillonnage de scénarios probables de variables d'état d'environnement par un site émetteur de données probabilistes selon l'invention procédant par méthode de Monte Carlo. - La figure 3 : décrit les principales étapes du processus selon l'invention de traitement de séries d'échantillonnages de scénarios probables de variables d'état d'environnement reçus par un groupe de trois systèmes industriels locaux distants fonctionnant selon l'invention, l'un des trois systèmes étant décrit en détail. - La figure 4 : décrit globalement les principales étapes du procédé selon l'invention de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, pour la prédiction probabiliste du comportement futur et/ou l'optimisation des paramètres de rendement d'une multitude de systèmes industriels locaux distants. - La figure 5 : décrit globalement les principales étapes d'une variante du procédé selon l'invention de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios variationnels probables, avec téléadaptation des paramètres numériques de calcul du rendement de chaque système industriel distant. - Les figures 6, 7 et 8: décrivent trois méthodes pour effectuer selon l'invention la transmission de smp séries probabilistes du site centralisé émetteur vers la multitude des sites industriels locaux distants et trois méthodes pour filtrer sélectivement le contenu des smp séries probabilistes transmises. - La figure 9 : décrit une application particulière du procédé selon l'invention dans le domaine de l'industrie financière. - La figure 10 : décrit une variante de mise en oeuvre de l'invention pour effectuer une optimisation probabiliste groupée de certains groupes de systèmes industriels locaux. Description de l'invention Les figures 1, 2 et 3 illustrent l'organisation générale des principales étapes du procédé par simulation stochastique de génération répétitive centralisée d'échantillons statistiques et de télétransmission de scénarios probables, pour la prédiction probabiliste du comportement futur et/ou l'optimisation des paramètres de rendement d'une multitude de systèmes industriels locaux distants (Su). Dans l'exemple retenu, les systèmes industriels locaux distants (S") sont constitués par des centrales électriques (S,,...,S",.... S") placés dans un même environnement (E) et formant, dans ce cas particulier, un réseau de distribution électrique d'un état ou d'une région. On voit en référence aux figures 1 et 4, l'environnement (E) d'usines de production électrique, connectée â un même réseau, influencées par des phénomènes aléatoires que sont le climat (1) agissant directement sur les performances des usines (rendement des centrales, niveau des barrages), les caractéristiques de la demande (quantité, prix) (2) affectant les revenus de la centrale électrique, la disponibilité du réseau (3) influençant le régime de fonctionnement des usines, et les caractéristiques des énergies primaires (4) influençant les coûts d'approvisionnement. Ces phénomènes aléatoires liés au comportement de l'environnement (B) du réseau sont donnés à titre d'exemple. Mais en comprendra que dans la pratique les phénomènes aléatoires à prendre en compte sont beaucoup plus nombreux et peuvent être différents. Même lorsqu'ils sont locaux (comme les phénomènes climatiques), ces phénomènes aléatoires ne peuvent être considérés indépendamment des phénomènes communs à toutes les usines (comme le prix du pétrole, le niveau de la demande sur le réseau). Ainsi, la température à proximité d'une usine agit non seulement sur le rendement de l'usine (donc sur son coût marginal), mais interagit avec le niveau de la demande et le prix de l'électricité. On a de fortes chances de voir évoluer cette température locale de manière corrélée avec la température à proximité des autres usines, et donc de voir le coût marginal de ces dernières varier dans le même sens et interagir par ailleurs avec le niveau de la demande. Rechercher à optimiser la marge d'une usine en faisant une analyse statistique des paramètres locaux de climat indépendamment de l'analyse statistique sur les paramètres globaux conduirait à des erreurs manifestes d'optimisation. Car paramètres locaux et paramètres globaux interagissent dans le cadre d'un même environnement global commun. En référence aux figures 1, 4 et 5, on constate qu'un centre fournisseur de données (Fo) est chargé de réunir des données observables permettant de décrire l'état de l'environnement global (E) du réseau électrique. Dans ce cas particulier de l'industrie électrique, le centre fournisseur de données (Fo) peut faire partie de la compagnie gérant l'organisation du réseau électrique. Mais dans d'autres applications, notamment financières, ce centre fournisseur de données (Fo) est une organisation indépendante (telle que les compagnies Reuters, Bloomberg) fournissant commercialement des données observables à une multitude de clients. On a installé en différents points de l'environnement (E), une série de r capteurs (y,,("- 1...') ). Ils constituent des observateurs de l'état d'environnement (E). Le fournisseur de données (Fo) collecte à distance, à chaque instant (t), les r mesures de l'état de l'environnement (eth(t)" - 1 ...r)) en provenance des capteurs. Ces r mesures constituent, à chaque instant d'échantillonnage (tj, les composantes d'un vecteur temporel de mesure d'environnement (M(tf)) (uh(tk#)l" - 1...r). Le fournisseur de données (Fo) est équipé d'un émetteur (Ao) de données observables. Celui-ci émet sur un réseau (Ro) d'émission de données observables - soit directement les mesures recueillies par les capteurs (variables observables), - soit les résultats de calculs algébriques déterministes (tels que moyennes, ... ) opérés par le fournisseur de données (Fo) sur ces mesures (variables quasi -observables). Par la suite, on assimilera variables observables et quasi-observables que l'on appellera globalement variables observables . Dans le cas représenté figure 1, le réseau (Ro) d'émission de données observables est schématisé par une antenne parabolique d'émission (5) transmettant les données (mesures) observables (M(tk)) (p,,(tk#)l" --r@) de l'environnement (E) à un satellite (6). Un site centralisé émetteur (U) selon l'invention est décrit dans sa forme générale en référence à la figure 2. Il est généralement constitué d'un système informatique (7) comprenant principalement - une interface de modélisation (S2) permettant à un opérateur mathématicien de faire le choix des différents modèles et variables, - un calculateur (C) de simulation pour effectuer les calculs numériques et statistiques, - une mémoire de masse (H) centrale de données probabilistes, - un récepteur (Bo) de mesures observables, - et un émetteur (Ap) de séries de simulations probabilistes. Ce site centralisé émetteur (U) est très distant (de plusieurs kilomètres, voire de plusieurs centaines de kilomètres) d'une majorité des n systèmes industriels locaux distants (S"("- 1...n) ), les centrales électriques (S,,.... S",.... Sn) qui sont réparties sur tout le territoire de l'état et apparaîtront figures 3 à 8. Le récepteur (Bo) de mesures observables est schématisé par une antenne parabolique de réception de données observables (8) qui reçoit en provenance du satellite (6), en flux continu périodique, les r mesures d'environnement (M(tk)) (Ph(tk)(h='@..r émises par le fournisseur de données (Fo) via le réseau (Ro) d'émission de données observables. L'interface de modélisation (S2) est préférentiellement située sur le site centralisé émetteur (U). Dans le cas plus général, elle est située au sein d'un site expert (U'), qui peut être distinct du site émetteur (U) mais en tout cas distant (de plusieurs kilomètres, voire de plusieurs centaines de kilomètres) d'une majorité des n systèmes industriels locaux distants (Su'u-@-..R). Au niveau de l'interface de modélisation (S2) ou du site expert (U'), on fait le choix - des p variables (X;li='...Pl) de l'environnement (E), également appelées facteurs de risque, - des horizons de temps (tk(k-'--"' des différentes simulations, - du modèle numérique descriptif d'environnement (D), définissant l'ensemble de p variables d'état (X;(t)''-'..-P) _ (@(,u@(t),...,@Cr(t))) en fonction des mesures observables (ph(t)(h= 1...r) (avec généralement p S r mais suffisamment grand, en sorte que l'ensemble des variables d'environnement (Xi(t)"-'...Pl) décrivent physiquement, à chaque instant d'échantillonnage (tk), la globalité de l'état de l'environnement (E)), - et d'une loi statistique paramétrée (modèle) (fa(X,,...,XP; t; Atk)) du comportement joint des facteurs de risque, au temps (t) et par tranches de temps (dtk(k =' "' , des variables d'état (X;(t)"-'...Pl), décrite par sa densité de distribution jointe. La loi statistique (fa) paramétrée (décrite numériquement) est généralement choisie par un mathématicien spécialiste des modèles statistiques. Elle est justifiée à partir d'observations historiques - il peut s'agir notamment de lois normales, log-normales, etc., ou plus énéralement d'une distribution statistique jointe des p variables d'environnement (X;(tk) (' ..P') ou de leurs variations temporelles (AX;(tk) (' - '...Pl), validée par des tests statistiques classiques adéquats (moyennes, variances, corrélations, sauts, comportements extrêmes, etc.). Au fur et à mesure de la réception des données observables d'environnement que l'on reçoit sur demande ou en continu, on met à jour et on enregistre dans la mémoire de masse (H) l'historique périodique des valeurs des mesures (Ph(t)lhLe calculateur (C) calcule, à chaque période d'échantillonnage (tk) successive, à l'aide du modèle numérique descriptif d'environnement (D) et du vecteur de mesures (M(t)), les p variables d'état (X;(t)'`-'P@) _ (@(,u,(t),..."ur(t))). Préférentiellement, le calculateur (C) calcule et met également à jour répétitivement l'historique des variations des p variables d'état (OX;(t)('='...P) qu'il enregistre dans la mémoire de masse (H). On entend par variations des variables d'état, soit la différence, soit le ratio de deux valeurs d'une variable à des instants successifs, ou plus généralement toute valeur numérique décrivant les variations (absolues ou relatives) de la variable entre deux instants successifs. Le calculateur (C) est muni d'un système simulateur (9) de scénarios, d'un système calibreur (10) de la loi statistique paramétrée (fa) pour qu'elle soit conforme aux observations historiques mémorisées dans la mémoire (H), et d'un générateur aléatoire (11) permettant de construire répétitivement des smp séries de simulations probabilistes (W;,,@) multi-valuées de description des pseudo-états ou pseudo-variations probables futur(e)s des variables d'état (Xt(tfi ,...Pl). Diverses méthodes de fonctionnement du simulateur (9) de scénarios seront décrites plus loin. Certaines génèrent des séries de simulations des variables d'état (X;(t)('-'"), et d'autres des séries de simulation des variations des p variables d'état (OX;(t)l'=' ,@Pl). En sorte qu'au niveau du site centralisé émetteur (U), à partir de la densité de distribution jointe de la loi statistique paramétrée (fa (Xi,..., XP; t; Atk)) retenue du comportement joint dans le temps (t) des variables d'état (Xt), on génère numériquement répétitivement à des instants t. successifs, une snip série de simulation probabiliste (W; ,j@ l'-'".p'j-' ..-'s' k-'"'") multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s de l'environnement (E). Selon un premier mode de mise en aeuvre, cette snip série de simulation probabiliste est une snip série probabiliste (X;,j@) de pseudo-états formée de sxrn p-uplets d'échantillons multi-valués de pseudo-états futurs (X;.j{tk)l'='".p'j='".S, k='"'"')), comprenant sxnixp valeurs numériques possibles des pseudo-valeurs futures probables des p variables d'état (X,(t)(' et constituées - pour chacun des m (ni >_ 1) horizons de temps (tk"=''n')) de la simulation, - de s (avec s très grand 1) pseudo-états futurs possibles (Xjk@= 1..s)) i (X;,j(tk)('-''"p'j-''.'S) du vecteur des p variables d'état (X;(tk)"-'"'p) . Selon un second mode de mise en oeuvre décrit figure 5, cette snip série probabiliste variationnelle (AX;,@k) est formée de sxm p-uplets d'échantillons multi-valués de pseudo- variations (AX;j{tk)(f-'...p, j=l...s, k=I...n:)) futures à partir de la valeur courante des variables d'état (X;(to)"-'".p), comprenant sxnixp valeurs numériques de pseudo-variations futures probables des p variables d'état (X;(tk)"-'"'p) et constituées - pour chacun des m (ni > 1) horizons de temps (tk(k-'".n')) de la simulation, - de s (avec s très grand 1) pseudo-variations numériques futures possibles (AXj,k@ _' "S)) (OX;,j(tk)"='"'p'j=''') du vecteur des p variables d'état (Xi(tk)"='..'p)). Comme cela apparaît en référence aux figures 2, 3, 4 et 5, grâce à l'émetteur (Ap) de série de simulations probabilistes équipant le site centralisé émetteur (M, répétitivement, à des instants 1o successifs, on diffuse ladite snip série probabiliste (W;,jk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s, comprenant les sxinxp pseudo-valeurs numériques constituée par : (cas 1) ladite snip série probabiliste (X;,jk) des pseudo-valeurs des variables d'état (Xi,j(Ik)"-'".p'j-'----` k=' '.' "') ), et/ou (cas 2) ladite snip série probabiliste variationnelle (AX;,jk) des pseudo-variations (OX;,j(tk)('-'".p,j='"-'s, k-1...nr)). Comme on le constate en référence aux figures 2, 3, 4 et 5, cette émission de simulations probabilistes a lieu de manière décentralisée du site émetteur (n en direction de l'ensemble de la multitude des systèmes locaux distants (Su) via le réseau (Rp) d'émission de données probabilistes. Dans le cas représenté, l'émetteur (Ap) de séries de simulations probabilistes est schématisé par une antenne parabolique (12) d'émission de données probabilistes transmettant successivement les données probabilistes, à savoir les sine séries probabilistes (W;,jA) au satellite (6). C'est-à-dire que les réseaux (Ro) et (RI sont confondus. On comprendra, que l'on pourrait tout aussi bien utiliser un autre satellite pour l'émission de données probabilistes que celui utilisé par la réception de données observables, ou un autre type de transmission de données, notamment via le réseau Internet. En référence à la figure 3, on voit l'organisation générale, selon l'invention, d'un ensemble de systèmes industriels locaux (S"). Il s'agit, dans l'exemple retenu lié à l'industrie électrique, d'un ensemble de centrales électriques (S,,...,S",.... Sn) dont le comportement industriel (production de KWs énergétiques) est influencé par l'environnement (E) décrit à la figure 1. Chaque système industriel local (S"("= 1...n)) comprend un système informatique local (13, 14, 15) sensiblement équivalent. L'architecture d'un système informatique local (15) du système industriel (S") est décrite figure 3. Celui-ci comprend principalement - une interface locale de modélisation (S2,u) pour configurer le site (S"( =''")) en fonction de différents variables et modèles, - un calculateur local (Cu) pour effectuer les calculs numériques et statistiques, - une mémoire de masse locale (H" ), - un récepteur local (Bp") de données probabilistes, conçu et réglé de manière à recevoir le flux de données probabilistes en provenance du satellite (6). On constate que chaque système industriel local (S") est équipé d'une antenne parabolique locale (16,17,l8) de réception de données probabilistes et d'un récepteur local de données probabilistes (Bp"). De manière indépendante, au sein de chacun des systèmes industriels locaux (S"), répétitivement on collecte, grâce à l'antenne parabolique (18) et au récepteur local (Bp"), la succession des snap séries probabilistes (W,jk) multi-valuées de description des pseudo-états ou pseudo-variations probables futur(e)s à sxnrxp valeurs numériques (W,jk X;-j' ou OX;jk) transmises par le satellite (6). Par l'intermédiaire de l'interface locale (S2"), on choisit m' (rn' ni) horizons de temps, retenus parmi les m dits horizons de temps (tk(k-'..."')). Eventuellement (dans le cas (2) décrit ci-dessus et représenté figure 5) où le site industriel local (S") reçoit du satellite (6) une sntp série probabiliste variationnelle (AX;,,k), on calcule localement ladite snzp série probabiliste des pseudo-états (X;ik) constituée de pseudo-échantilIons d'états futurs (X;,j(tk)1'-''..n,'-@`k-1.... correspondante (au moins pour les ni" horizons de temps (tk) retenus). Pour chacun des m' horizons de temps (tk) retenus, on applique numériquement le modèle systémique local (D") préalablement fixé afin de déterminer dans chacun des s scénarios de simulation reçus (j 1, .... s) l'indicateur de rendement (I"jk D"(X", (X;,j(tk)) - pour au moins une valeur particulière (k,") = (A,,,,...,du vecteur (,'#") des paramètres locaux d'action (que l'on pourra être amené à faire varier), - et pour chacun des s p-uplets simulés de variables d'état (X;j{tk)"-'...n), correspondant audit horizon de temps (tk) retenu de la série probabiliste de pseudo-état (X;,j(tk)#i-'...s . On détermine ainsi s pseudo-valeurs simulées correspondantes de l'indicateur de rendement (I",j(tk,@.")) D"(X",(X;,j{tkf'-1--`j-1...St). On regroupe et on ordonne les s valeurs correspondantes de l'indicateur de rendement (lu, j(tk,@,u)@ - `@)> selon le niveau de rendement (1",j). On calcule et on mémorise dans la mémoire de masse locale (H") l'histogramme (1'(I")) ainsi que la distribution de probabilité (en nombre ou pourcentage) ('h"(1")) du niveau atteint de l'indicateur de rendement local (I") dans la configuration selon le modèle local (D"). On en déduit, par un calcul statistique, la valeur d'au moins un paramètre statistique de réglage choisi (V(k")) (tel que : moyenne, écart type, Value at Risk, utilité, etc...) caractéristique de la distribution de probabilité (`P(?,,",I)) des niveaux de rendement (I"(tk,@,"), en fonction du vecteur des paramètres locaux d'action (#,") = (;,",, ...,@."q). Eventuellement, on effectue localement, au sein de chaque système local (S"), une optimisation numérique (selon des critères prédéfinis) du (ou des) paramètre(s) statistique(s) de réglage choisi (V(a,")). Puis on règle le niveau du vecteur des paramètres locaux d'action (@-") ...,/ 1,"G) de chaque système local (S") à une valeur optimale calculée en fonction des critères préétablis, afin notamment d'anticiper le comportement probable futur et l'efficacité de chaque système (S") soumis à l'environnement (E) au cours du ou des horizon(s) de temps (tk) retenu(s). Un mode de mise en oeuvre particulière de l'invention décrit figures 3 et 5 vise â rendre adaptatif pour chaque système industriel local (S"("-"--" la modélisation de l'indicateur de rendement correspondant à chaque paramètre local d'action Rh). A cet effet, au niveau du site expert (U'), on associe à chaque paramètre local d'action (;"h(h-'...q)) un code générique (Lh) (ne dépendant pas du système local (S")). On établit pour chaque code (Lh), et on actualise régulièrement, un évaluateur de rendement (11h), c'est-à-dire un programme (matérialisé sous forme d'un fichier numérisé) permettant de calculer, en fonction du p-uplet des variables d'état (X) (X;('- 1...n et du paramètre local d'action (Â"h), la composante correspondante du rendement (1"h) = (IIh(A"h,X)), de telle sorte que, pour chaque système local (S"), ledit modèle systémique local (D") consiste simplement à additionner les composantes (1"h) de l'indicateur de rendement correspondant à chaque paramètre local d'action (A"h) Le site expert (U') est relié au site centralisé émetteur par un réseau (Rl et lui transmet les définitions numériques adaptatives des composantes de rendement (I"hjk) = (IZh(Auh,Xjk)). Au niveau du site centralisé émetteur (U), d'une part, on génère numériquement et on émet, répétitivement à des instants (to) successifs, à partir du modèle statistique retenu (fa(Xl,...,Xp; t; At)) du comportement joint dans le temps (t) des variables d'état (X;('='...P , la sntp série probabiliste (W;Jk) multi-valuée de description des pseudo-états ou pseudo- variations probables futur(e)s. D'autre part, à l'apparition de chaque nouveau code (Lh) ou de chaque modification de l'évaluateur de rendement (IIh), on émet de manière centralisée la paire (Lh,IIh) sous forme d'un fichier numérisé. Au niveau de chaque système industriel local distant (Su) on calcule la composante de rendement (luhjk) (11h(Âuh,Xjk)) - pour chacun des in ' (m `S m) horizons de temps (tk) de la simulation retenus à l'échelon local (S"), - pour chacun des j = 1,...,s (avec s très grand 1) pseudo-états numériques possibles (Xjk (i='-s)) reçus, ou reconstitués à partir de la série des pseudo-variations (OX;,jk), - et pour chacun des h = 1, .... q paramètres locaux d'action (Ah). t Puis, par sommation, selon la formule ci-dessus (*) de l'indicateur de rendement (I"jk), on calcule son histogramme, puis le paramètre statistique retenu (V(@")). Enfin on applique le vecteurs de paramètres locaux d'action (@,") qui optimise la valeur du paramètre statistique (V). Un cas particulier du mode de mise en oeuvre ci-dessus consiste en ce que les évaluateurs de rendement (11h) soient proportionnels au paramètre local d'action (Âuh) c'est-à-dire que IIh(@uh,X) = Âuh 7rh(X ) où (irh) est appelé évaluateur de rendement marginal . Dans ce cas, le modèle systémique local (Du) revient à faire un produit scalaire entre le vecteur de paramètres locaux d'action (@,") et le vecteur de rendements marginaux (n(X)) (ir,,(X)(' Du(Xu,X) (X, I n(X)) (**) Le calcul, au niveau du site local (S" ), du paramètre statistique (V(%.")) s'en trouve simplifié et son optimisation accélérée par l'enregistrement préalable, dans la mémoire de masse (H") du calculateur local (C") des séries de rendements marginaux (Irh(Xjk)(h -'"'9'j l ...s, k =1 ... > . Une première variante particulière de mise en oeuvre de l'invention est décrite également en référence à la figure 4. Elle concerne le cas où, au niveau d'un site centralisé émetteur (U), on génère numériquement répétitivement à des instants to successifs, à partir du modèle statistique retenu (fa(Xl,...,Xn; t; dt)) du comportement joint dans le temps (t) des variables d'état (X;), une série probabiliste variationnelle (fX;Jk) formée de sxnt p-uplets d'échantillons multi-valués de pseudo-variations (AX,,,{tJ'-'...p,'=l...s,k=1..."'@) d'états futurs à partir de l'état courant des variables d'état (X;(to)('-'---n . Cette série de pseudo-variations comprend sxntxp valeurs numériques de pseudo-variations futures probables des p variables d'état (X;(tk)"='...P'). Elle est constituée pour chacun des m (ni >_ 1) horizons de temps (tktk=' de la simulation, de s (avec s très grand 1) pseudo-variations numériques possibles (OX;,j(tk)I# ='...s') des p variables d'état (X;(tk)('-'...n). On diffuse de manière répétitive cette sntp série probabiliste variationnelle (OXi,jk) des pseudo-variations d'états (OX;,j(tk)('-'...P. j=I...s, k=I...m)), ce de manière décentralisée du site émetteur (U) en direction de l'ensemble de la multitude des systèmes locaux distants (Su). Ceci est effectué par l'émetteur (Ap) de séries de simulations probabilistes grâce â l'antenne parabolique (12) d'émission de données probabilistes et via le satellite (6). De manière indépendante, au niveau de chacun des systèmes locaux (S"), le récepteur local (Bpu) de données probabilistes collecte de manière répétitive via l'antenne parabolique locale (16,17,18) de réception de données probabilistes, à chaque instant (ta), et selon une fréquence déterminée, voire à chaque modification de l'une des variables (Xi), la valeur courante des variables d'état (Xi(to)('='...p). Il collecte également ladite srnp série probabiliste variationnelle (AXiJk) des variations d'états (OXij(Otk)('='...p' j=I...s, k=I...m) en fonction de l'incrément de temps Atk , émise par le site central (U). Ces données sont mémorisées à l'intérieur de la mémoire de masse locale (H"). Le calculateur local (C") évalue localement ladite smp série probabiliste (Xi,jk) de pseudo-états (Xi,j(tk)(i=I...p, j=l...s, k=1...`)) aux instants futurs (tk) (t, + Atk), correspondant à la série variationnelle par calcul numérique incrémental du type Xij@tk) Xi(to) (1 + AXi,j(tk)). On peut utiliser toute autre formule (addition, exponentielle, etc...) faisant intervenir ladite snip série probabiliste variationnelle (AXi,jk) émise par le site centralisé émetteur (U) et la valeur de l'état courant des variables d'état (Xi(t)('='...p). Une méthode particulière de mise en oeuvre de cette première variante est celle décrite plus haut en référence à la figure 4. Au niveau du site centralisé émetteur (U), on génère numériquement répétitivement â des instants to successifs, une silip série probabiliste variationnelle (AXi,jk) telle que décrite ci-dessus. On collecte de manière répétitive et on mémorise la valeur courante des variables d'état (Xi(to)('=',**p)). On diffuse de manière répétitive et décentralisée du site émetteur (U) en direction de l'ensemble de la multitude des systèmes locaux distants (Su), à la fois - la siiip série probabiliste variationnelle (OXi,jk) des pseudo-variations d'état (Arij(tk)(i-I...p, j=l...s, k=I...ni)), - et la valeur courante des variables d'état (Xi(to)('='-- p). Une seconde variante de combinaison particulière de mise en oeuvre de l'invention est décrite en référence aux figures 1 et 4. Comme cela apparaît figure 1, préalablement, et en différents points de l'environnement (E), on installe une série de r capteurs (y,,(h=1 ...r)), observateurs de l'état de l'environnement (E), émettant à chaque instant (t) un ensemble de r mesures de l'état de l'environnement (ph(t)(h='...')). Dans le cas décrit figures 1 et 4, le site émetteur (U) s'est abonné à la fourniture de données d'environnement fournies par un fournisseur (Fo) de données observables qui lui-même opère le réseau de capteurs (yh(h=' ...')). Au niveau du site concepteur (U'), très distant d'une majorité des n systèmes industriels locaux distants (Su(u= I...n)), on définit préalablement un modèle numérique descriptif d'environnement définissant l'ensemble de p variables d'état (Xi(t)('='...p)) (D(ftl(t),.... p,(t))) en fonction des mesures (jch(t)(h='...`)) (avec généralement p 5 r), en sorte que l'ensemble des variables d'environnement(Xi(t)('='...p)) décrivent physiquement, à chaque instant d'échantillonnage (tk), la globalité de l'état de l'environnement (E). Mais on comprendra que l'on peut également effectuer cette opération directement au niveau du site émetteur (U), ou pour le compte de celui-ci, notamment lorsque le site concepteur (U') est situé au sein du site émetteur (U). Au niveau du site centralisé émetteur (U), on collecte à distance et on mémorise les r mesures d'environnement (Ph(t)(h=' -r)), ces r mesures constituant à chaque instant d'échantillonnage (tk) les composantes d'un vecteur temporel de mesure d'environnement (1VI(tk)) (tlh(tk)(h= I....), puis on calcule, à chaque période d'échantillonnage (tk) successive, à l'aide du modèle numérique descriptif d'environnement (D) et du vecteur de mesures (M(t)), les p variables d'état (Xi(t)(i='..'p) (e(y,(t),...,p,(t))). L'invention recommande deux méthodes particulières de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, pour la prédiction probabiliste du comportement futur et/ou l'optimisation des paramètres de rendement d'une multitude de systèmes industriels locaux distants (S" ). La première méthode proposée par l'invention de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables consiste, en combinaison, à utiliser la méthode connue dite de Monte-Carlo. Selon cette première méthode, qui est décrite en référence à la figure 2, au niveau du site centralisé émetteur (U) on maintient et on met à jour un historique des valeurs des p variables d'état (X;(t)('-"@@P de l'environnement global commun (E) au cours du temps. On génère numériquement répétitivement à des instants to successifs, à partir - du modèle statistique retenu (fa(X,,...,Xp; t; At)) du comportement joint dans le temps (t) des variables d'état (X;), - et de cet historique mémorisé des valeurs des p variables d'état (X;(t)('- 1...P), la smp série probabiliste (W,jk) multi-valuée de description des pseudo-états ou pseudo- variations probables futur(e)s de l'environnement (E), selon la méthode dite de Monte- Carlo . C'est-à-dire que l'on procède par tirages pseudo-aléatoires des pseudo-valeurs (WJk) grâce à un générateur aléatoire (11) de nombres situé sur le site centralisé émetteur (U), voire sur le site concepteur (U'). On peut utiliser un processus - soit pseudo-aléatoire de nature ergodique, - soit à discrépance faible (suites de Sobol, Hamersley, etc.), - soit toute autre méthode produisant des nombres asymptotiquement répartis selon une distribution spécifiée (par exemple uniforme sur un intervalle). Et on diffuse ladite smp série probabiliste (W;jk) multi-valuée de description des pseudo- états ou pseudo-variations probables futur(e)s de manière décentralisée du site émetteur (U) en direction de l'ensemble de la multitude des systèmes locaux distants (S" ). Cette méthode de mise en oeuvre est celle qui est recommandée par l'invention. La seconde méthode proposée par l'invention de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables consiste, en combinaison, à utiliser la méthode connue d'analyse historique. C'est une simplification de la méthode Monte-Carlo décrite ci-dessus. Selon cette méthode d'analyse historique, on émet la multitude d'échantillons des variations probables futures des facteurs de risque directement à partir de l'historique mémorisé. Au niveau du site centralisé émetteur (U), on maintient et on met à jour un historique des valeurs des p variables d'état (X;(t)"-' ...P') de l'environnement global commun (E) au cours du temps. On génère numériquement répétitivement à des instants t, successifs, à partir - du modèle statistique retenu (fx(X,,...,XP; t; At)) du comportement joint dans le temps (t) des variables d'état (X;), - et de cet historique mémorisé des valeurs des p variables d'état (X;(t)"- P'), la smp série probabiliste (W,ik) multi-valuée de description des pseudo-états ou pseudo- variations probables futur(e)s de l'environnement (E). Selon la méthode dite historique , on reconstruit l'estimation d'états futurs probables des p variables d'état (X;(t)"- ' ...P') à partir de l'état courant (X;(to)"-''.-Pt) et de variations historiques passées (AX;(tl"-'...P) pour un ensemble de dates passées (t" sntp série probabiliste (WJk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s, de manière décentralisée du site émetteur (U) en direction de l'ensemble de la multitude des systèmes locaux distants (S" ). Une troisième variante de mise en oeuvre de l'invention consiste à introduire des biais dans les simulations stochastiques émises. Ces biais visent à prendre en compte des situations extrêmes qui n'arrivent que rarement et/ou constituent des scénarios catastrophes . Selon cette méthode, on biaise la génération des snip séries probabilistes (W;,jk) multi-valuées de description des pseudo-états ou pseudo-variations probables futur(e)s émises par le site émetteur (U) en direction des multiples systèmes industriels locaux distants (S"), ce par rapport à la statistique de l'historique des variables d'états ou des variations de variables d'état. Selon cette variante de l'invention, au niveau du site centralisé émetteur (U), on maintient et on met â jour un historique des valeurs des p variables d'état (Xi(t)i`- ' ...P') de l'environnement global commun (E) au cours du temps. On génère numériquement répétitivement à des instants to successifs la snip série probabiliste (W,à multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s de l'environnement (E), à partir - du modèle statistique retenu (fa(X,,...,XP; t; At)) du comportement joint dans le temps (t) des variables d'état (X;), - et de cet historique mémorisé des valeurs des p variables d'état (X;(t)"- i- P)). Mais on introduit volontairement dans la snip série probabiliste (W,ik) des biais de comportement statistique joint des variables d'état (Xi(t)('-i...P) vis-à-vis de leur comportement statistique historique. Les séries biaisées émises reproduisent ainsi des scénarios catastrophes . Ceci permet de tenir compte de l'impact desdits scénarios catastrophes sur les systèmes distants (S"). Et on diffuse ladite snip série probabiliste (W,j@) multi-valuée biaisée de description des pseudo-états ou pseudo-variations probables futur(e)s de manière décentralisée du site émetteur (U) en direction de l'ensemble de la multitude des systèmes locaux distants (S"). Selon une première sous-variante de la troisième variante décrite ci-dessus, l'invention recommande d'adapter sélectivement la teneur des snip séries probabilistes (W,ik) multi- valuées biaisées de description des pseudo-états ou pseudo-variations probables futur(e)s, émises par le site émetteur (U) en direction de l'ensemble de la multitude des systèmes locaux distants (S"), ce en fonction de requêtes effectuées sélectivement par les différents systèmes distants (Su). A cet effet, certains des systèmes distants (S") envoient au système centralisé émetteur (U) des requêtes sous forme de paramètres numériques de scénarios catastrophes selon des biais spécifiques définis par eux. On génère au niveau du site centralisé (U), en fonction de ces requêtes, éventuellement en fusionnant celles qui sont identiques, des sflip séries probabilistes (W,lk) multi-valuées de description des pseudo-états ou pseudo- variations probables futur(e)s, biaisées en fonction des paramètres spécifiés requis par les systèmes distants. Le site émetteur (U) envoie sélectivement vers chaque système local distant (S" ), ladite snip série probabiliste (W;,,@) multi-valuée biaisée de description des pseudo- états ou pseudo-variations probables futur(e)s qui correspond à sa requête. Cet envoi est effectué de manière électronique (via Internet, un réseau local, le téléphone, ou tout autre technique permettant la transmission de données numérisées). Selon une deuxième sous-variante de la troisième variante décrite ci-dessus, l'invention recommande que les systèmes distants (S") altèrent eux-même tout ou partie des snip séries probabilistes (W,j@) qu'ils reçoivent du site émetteur (U) pour prendre en compte des situation possibles extrêmes qui les concernent. Selon cette sous-variante, certains des systèmes distants (Su) altèrent eux-mêmes tout ou partie des snip séries probabilistes (W,J@) multi-valuées de description des pseudo-états ou pseudo-variations probables futur(e)s. Ceci permet de construire localement des séries biaisées reproduisant des scénarios catastrophes dont les paramètres sont : soit spécifiés au niveau du système émetteur central (E) et diffusés en direction des systèmes locaux (S" ), soit définis au niveau des systèmes locaux (S" ). Une quatrième variante de mise en oeuvre préférée de l'invention consiste en ce qu'au niveau du site centralisé émetteur (M, on génère numériquement la sinp série probabiliste (Wjk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s de l'environnement (E), de manière périodique, selon une fréquence d'échantillonnage p). Selon une cinquième variante de mise en oeuvre, l'invention recommande trois méthodes pour effectuer la transmission des smp séries probabilistes (W,j) multi-valuées de description des pseudo-états ou pseudo-variations probables futur(e)s, du site centralisé émetteur (0 vers la multitude des sites industriels locaux distants (Su(" = 1...n) ) . Une première sous-variante d'émission de données probabilistes est décrite figure 6. Au niveau du site centralisé émetteur (U), on diffuse ladite snlp série probabiliste (W,ik) multi- valuée de description des pseudo-états ou pseudo-variations probables futur(e)s, de manière t décentralisée du site émetteur (M en direction de l'ensemble de la multitude des systèmes locaux distants (S") par émission sous une forme broadcast , notamment par émission radio. De manière indépendante, chacun des systèmes locaux (S"), majoritairement très distants du site centralisé émetteur (U), collecte répétitivement ladite slllp série probabiliste (Wijk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s à sxnlxp valeurs numériques (Wjk X;,lk ou AX;Jk), par réception quasiment en continu de cette émission broadcast en provenance du site centralisé émetteur (U). Une deuxième sous-variante d'émission de données probabilistes est décrite figure 7. Au niveau du site centralisé émetteur (0, on diffuse ladite snlp série probabiliste (W;ik) multi- valuée de description des pseudo-états ou pseudo-variations probables futur(e)s, de manière décentralisée du site émetteur (U) en direction de l'ensemble de la multitude des systèmes locaux distants (S") par émission sous une forme dite PP aval , c'est-à-dire par émission sur un réseau (R) -notamment de type Internet - du site centralisé émetteur (U) à l'adresse spécifique de chacun des systèmes industriels locaux distants (S,("="," . A cet effet, on constate que le site émetteur (U) est équipé d'un modem central (20) et que chaque système industriel local distant (S"("- 1-n) est équipé d'un modem local (21) par lesquels transitent les données probabilistes. De manière indépendante, chacun des systèmes locaux (S"), majoritairement très distants du site centralisé émetteur (M, collecte répétitivement ladite slip série probabiliste (W,jk) multi-valuée de description des pseudo-états ou pseudo- variations probables futur(e)s à sxnlxp valeurs numériques (Wjk X;jk ou AX;,ik), par enregistrement successifs dans la mémoire de masse locale (H") du calculateur local (C"). Une troisième sous-variante d'émission de données probabilistes est décrite figure 8. Au niveau du site centralisé émetteur (n, on diffuse ladite slip série probabiliste (W,jk) multi- valuée de description des pseudo-états ou pseudo-variations probables futur(e)s, de manière décentralisée du site émetteur (U) en direction de l'ensemble de la multitude des systèmes locaux distants (S"), par mise à disposition de ladite slip série probabiliste (Wjk) multi- valuée sur un serveur (Z) de données probabilistes. Celui-ci est disponible sur un réseau (R) d'émission (notamment de type Internet), ce par accès PP de chacun des systèmes industriels locaux distants (S"t" - I ...n . De manière indépendante, chacun des systèmes locaux (S") collecte répétitivement ladite smp série probabiliste (W;,jk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s, par accès sous une forme dite PP amont . C'est-à-dire qu'ils accèdent au serveur (Z) du site centralisé émetteur (U) via le réseau (R) par requête à son adresse spécifique (notamment de type Internet). Puis ils acquièrent les smp séries probabilistes par transfert (notamment de type FTP ) via le réseau (R) et par enregistrement dans la mémoire de masse locale (H") de leur calculateur local (C"). Selon une sixième variante de mise en oeuvre, l'invention recommande trois méthodes pour organiser des slllp séries probabilistes (Wl,jk) transmises par le site centralisé émetteur (U) vers la multitude des systèmes industriels locaux distants (S"(" - '-- '), ceci afin de compresser le volume des données probabilistes transmises et d'accélérer leur prise en compte par la multitude des systèmes industriels locaux distants (Su(" -1 ...n)). Selon une première sous-variante d'organisation de l'émission des données probabilistes, au niveau du site centralisé émetteur (U), on génère numériquement répétitivement à des instants to successifs, à partir du modèle statistique retenu (fjXi,...,Xn; t; At)) du comportement joint dans le temps (t) des variables d'état (X;(`-1 ...P'), une snip série probabiliste (WJk) multi- valuée de description des pseudo-états ou pseudo-variations probables futur(e)s de l'environnement (E). Celle-ci est constituée de s (avec s très grand 1) pseudo-valeurs numériques des p composantes (W,j(tk)('-1 ... P' j=1...5) pour chacun des iii (n1 > 1) horizons de temps (tk(k= i...Pn) ), sous la forme d'un smp fichier de sxmxp données alphanumériques successives . Mais, d'une part, on subdivise ce smp fichier en ni sous-fichiers temporels successifs de niveau 1 regroupant les pseudo-valeurs d'un même horizon de temps (tk(k-1."n` . D'autre part, chacun des sous-fichiers temporels est constitué d'une succession de s sous-fichiers de tirage de niveau 2 constitués chacun d'un p-uplet de pseudo-valeurs. Selon une deuxième sous-variante d'organisation de l'émission des données probabilistes au niveau du site centralisé émetteur (U), on génère numériquement répétitivement à des instants to successifs une snip série probabiliste (W;,-jk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s de l'environnement (E) en séparant les ni sous- fichiers temporels successifs de niveau 1 regroupant les pseudo-valeurs d'un même horizon de temps (tk(k-1..."')) par un indicateur d'horizon de temps (Ilk(k=L.."') correspondant. Et, au niveau de chacun des systèmes locaux (S" ), on distingue les uns des autres les différents dits sous-fichiers temporels d'une même smp série probabiliste (W,ik) multi-valuée collectée en identifiant leurs indicateurs spécifiques d'horizon de temps (17k(k= 1..."'@). Selon une troisième sous-variante d'organisation de l'émission des données probabilistes au niveau du site centralisé émetteur (0, on génère numériquement répétitivement à des instants ta successifs les sinp séries probabilistes (W;.J#, ) multi-valuées de description des pseudo- états ou pseudo-variations probables futur(e)s de l'environnement (E), en séparant les s sous-fichiers de tirage successifs de niveau 2 (d'un même fichier temporel) regroupant chacun un p-uplet de pseudo-valeurs par un indicateur de tirage (zj'#= `-`) correspondant. Et, au niveau de chacun des systèmes locaux (S"), on distingue les uns des autres les différents sous-fichiers de tirage successifs de niveau 2 (d'un même fichier temporel) regroupant chacun un p-uplet de pseudo-valeurs d'une même sinp série probabiliste (W,jk) multi- valuée collectée, en identifiant leur indicateur spécifique de tirage (zjU- 1...s . Selon une septième variante de mise en oeuvre, l'invention recommande trois méthodes pour filtrer sélectivement le contenu des smp séries probabilistes (WJk) transmises par le site centralisé émetteur (U) vers la multitude des systèmes industriels locaux distants (S"(" =1 ...n)). Ceci vise à adapter leur contenu aux besoins et aux particularités des systèmes industriels locaux distants (Su(" 1...n)) ainsi qu'aux liens hiérarchiques, industriels ou commerciaux entre le site centralisé émetteur (U) et la multitude des systèmes industriels locaux distants (Su(u = 1...n)). Une première sous-variante de filtration de l'émission des données probabilistes est décrite figure 6. Au niveau du site centralisé émetteur (U), on diffuse dans sa globalité ladite snip série probabiliste (Wjk) multi-valuée de manière décentralisée du site émetteur (U) en direction de l'ensemble de la multitude des systèmes locaux distants (S") par émission sous une forme broadcast . On peut notamment procéder par émission radio, accessible à tous les systèmes locaux (S"). Et, de manière indépendante, au niveau de chacun des systèmes locaux (Su), on installe et on utilise un filtre sélectif préréglé (22) de mise à disposition d'une fraction (W ijk) prédéfinie de ladite snip série probabiliste (W,ik) collectée par réception en provenance du site centralisé émetteur (U). Les caractéristiques de filtrage des données par ce filtre sélectif peuvent notamment être acquises par un abonnement spécifique pris par chaque système local (Si,) auprès du site émetteur (U). Elles peuvent être modifiées par téléaction à travers le réseau reliant le site centralisé émetteur (U) aux systèmes locaux (Su). Une deuxième sous-variante de filtration de l'émission des données probabilistes est décrite figure 7. Au niveau du site centralisé émetteur (U), on diffuse de manière sélective (selon une règle d'attribution prédéfinie) diverses fractions différentes ((W,jk)g(g-i ,,'' prédéfinies (notamment acquises par abonnement spécifique) de ladite snip série probabiliste (Wi,ik). Cette diffusion sélective est effectuée de manière décentralisée du site émetteur (U) en direction de plusieurs groupes (Gglg-""*" attributifs de la multitude des systèmes locaux distants (S"). Chaque groupe (G8) recoit une fraction particulière ((Wi,@k)g). Ceci s'effectue par émission sous une forme dite PP aval , c'est-à-dire par émission sur un réseau (R) (notamment de type Internet) du site centralisé émetteur (U) à l'adresse spécifique de chacun des systèmes industriels locaux distants (Si, lu = i..." . De manière indépendante, chacun des groupes (G8'8-""" attributifs de systèmes locaux (S") collecte répétitivement ladite fraction prédéfinie ((Wijk)g) de la sinp série probabiliste (W ,,i k). Ils enregistrent successivement cette fraction prédéfinie ((Wi,,@)b) dans la mémoire de masse locale (Hu) de leur calculateur local (C,"). Une troisième sous-variante de filtration de l'émission des données probabilistes est décrite figure 8. Au niveau du site centralisé émetteur (U), on diffuse ladite snip série probabiliste (Wijk) de manière décentralisée du site émetteur (U) en direction de l'ensemble de la multitude des systèmes locaux distants (S") par mise à disposition de ladite snip série probabiliste (W,jk) sur un serveur (Z) disponible sur un réseau (R) (notamment de type Internet) auquel chacun des systèmes industriels locaux distants (Sue" - '- --") peut accéder par accès PP amont . Mais on met de manière sélective à disposition de différents groures (Gg'g- '---"') de systèmes locaux (Su) sur le serveur (Z) diverses fractions limitatives ((Wi,j' )g) prédéfinies (notamment acquise par un abonnement spécifique) de ladite siiip série probabiliste (Wi,.jk). Ceci peut notamment s'opérer par l'intermédiaire d'un système d' identificateurs et de mots de passe (23) transmis par chaque système local (S"@" --"@) lors de son accès au serveur (Z) . De manière indépendante, chacun des systèmes locaux (Su) collecte répétitivement ladite fraction ((Wi.k)g) prédéfinie attribuée de ladite snip série probabiliste (Wjk) par accès sous une forme dite PP aval , c'est-à-dire par accès via le réseau (R) du site centralisé émetteur (U) à son adresse spécifique. Ils acquièrent par transfert ladite fraction ((W,j@)g) prédéfinie attribuée de ladite snip série probabiliste (W,lk) multi- valuée qui leur est attribuée. Ils enregistrent successivement cette fraction prédéfinie ((Wi,ik)g) dans la mémoire de masse locale (Hu) de leur calculateur local (C"). Selon une huitième variante de mise en oeuvre, décrite en référence à la figure 10, l'invention recommande d'effectuer une optimisation probabiliste groupée de certains groupes de systèmes industriels locaux distants (Su(" -- " . Ceci se rapporte notamment au cas où les systèmes industriels d'un même groupe d'optimisation ont des objectifs industriels (de production, d'économie, etc.) communs qu'il convient d'optimiser globalement. On choisit préalablement n' (n'5 n) systèmes industriels locaux, retenus parmi les n dits systèmes industriels locaux (S"). On choisit également préalablement ni' (ni ` horizons de temps, retenus parmi les m dits horizons de temps (tk(k - 1 ... "`)), identiques pour tous les systèmes locaux. Au niveau de chacun des n' systèmes locaux (Su) retenus, et pour chacun des in' horizons de temps (tk) retenus, on applique numériquement le modèle systémique local (Du) préalablement fixé afin de déterminer dans chacun des s scénarios (j = 1,...,s) le niveau de l'indicateur de rendement (I"jk D"(Â,", (Xi,j(tk)) - pour au moins une valeur particulière (;,") du vecteur des paramètres locaux d'action, - et pour chacun des s p-uplets simulés de variables d'état (X;,j(tk)" _', -"n'), correspondant audit horizon de temps (tk) retenu, de la série probabiliste de pseudo-états (Xj(tk)"= '...S>) (Xi j(tk)cr='...n,j=l...s . On détermine les s pseudo-valeurs simulées de l'indicateur de rendement (lu,j(tk,Xu)) (Du(;#u,Xj(tk))Ii ''..sl) correspondantes. On émet la série de simulations locales constituée des s pseudo-valeurs simulées de l'indicateur de rendement (luj(tk,@u)) (Du(Xu,Xj(tk))'I '...s') correspondantes du système local (S"), de manière centralisante en direction d'un site centralisateur (B). Au niveau du site centralisateur (B), on collecte les sxn' pseudo-valeurs simulées d'indicateurs de rendements (lu,j(tk,@,"))@=l...s) en provenance de chacun des n' systèmes locaux (S") et pour chacune des s simulations, on calcule, pour chaque simulation (j 1,...,s), un indicateur de rendement global (Jj) de l'ensemble des systèmes locaux à partir d'un modèle d'évaluation global (Jj)=(x(lu,j(Ik,@,")@"='..-n:j=i...s>). On regroupe et on ordonne les s valeurs correspondantes de l'indicateur de rendement selon son niveau. On calcule et on mémorise l'histogramme (n ainsi que la distribution de probabilité (en nombre ou pourcentage) ('E(J)) du niveau atteint de l'indicateur de rendement global (J) dans la configuration choisie des vecteurs de paramètres locaux d'action (A) = (@,"@"= 1...n@), selon les modèles locaux (Du) et le modèle d'évaluation global (x) retenus. On en déduit, par un calcul statistique, la valeur d'au moins un paramètre statistique de réglage global choisi (V(A)) (tel que : moyenne, écart type, Value-at-Risk, utilité, etc...) caractéristique de la distribution de probabilité (x(?@,,...,@,n ; J)) des niveaux de rendement globaux (J(tk ;#,,,...,k" )) en fonction des valeurs du multi-vecteur (#j = (J 1, .... Lino des vecteurs de paramètres locaux d'action (glu). Eventuellement on effectue globalement une optimisation numérique (selon des critères prédéfinis) du paramètre statistique de réglage global choisi (V('_])) par rapport au multi- vecteur (_) = (J 1,...,n). Et on transmet à chaque système local (Su) le niveau correspondant de son vecteur optimisé des paramètres locaux d'action (@,"), afin notamment d'anticiper globalement le comportement probable futur et l'efficacité de l'ensemble des système locaux (S") soumis à l'environnement (E) au cours du (ou des) horizon(s) de temps (tk) retenu(s). On décrit en référence à la figure 9 une application du procédé selon l'invention dans le domaine de l'industrie financière. Le procédé de l'invention effectue une simulation stochastique répétitive hiérarchisée par échantillonnage de pseudo-états futurs probables de variables de marché, pour la prédiction probabiliste du comportement futur et/ou l'optimisation des paramètres de rendement d'une multitude de n (n étant très grand 1) d'opérateurs financiers locaux (S"@ = @...n)), L'invention s'applique dans le cas où les opérateurs financiers locaux (S"(" = I...n@) sont décentralisés et très distants les uns des autres. Chaque opérateur financier local (S"("= '...n>) possède un portefeuille (P"( -"* n constitué d'actifs financiers (6h@h='N@), dont la valeur (1"("- 1...n dépend des mouvements du marché des capitaux (E). La valeur (lu@"=1...n>) d'un portefeuille (P@Iu-'...n)J)eut être calculée par combinaison d'un vecteur de variables de marché (X(t)) (X;(t)l' - --nl) et, pour chaque opérateur financier local distant (S"), d'un vecteur de quantités (a,") = (A"h(h ' ..-N)) de chacun des actifs financiers (9h(h='...9I) détenus dans son portefeuille (Pü). Les variables de marché (X;(t)('='-1) sont communes aux différents opérateurs financiers locaux (S"(u= 1-..n , Elles varient au cours du temps (t). Il est possible pour chaque opérateur financier local distant (S") de reconstituer avec une précision acceptable la valeur de tous les actifs financiers pouvant être présents dans son portefeuille (P"), selon un modèle local d'évaluation (D"), à partir du vecteur des variables de marché (X(t)) = (X;(t)" -""#P . Les variables de marché peuvent représenter des prix d'actifs ou de commodités, des taux d'intérêt ou des coefficients de courbe de taux, des volatilités ou des coefficients de surfaces de volatilité, des indicateurs de défaut d'un émetteur ou d'une contrepartie, ou encore tout autre nombre apportant une information pertinente sur l'état des marchés de capitaux. Par exemple la valeur (1"I" - 1 ... ") d'un portefeuille (pJu - ... R)) peut être calculée selon la combinaison linéaire où irh(t,X(t),6h) représente le prix de l'actif(6h) à l'instant (t) selon le modèle local (Du) lorsque le marché est dans l'état (X(t)), ou une autre formule pouvant faire apparaître une non-linéarité par rapport aux quantités (Â"h), pour tenir compte notamment du risque de liquidité. Selon une première sous-variante d'application du procédé selon l'invention au domaine de l'industrie financière, au niveau d'un site expert (U'), on fait préalablement le choix - des p variables de marché (Xi" = ' ... P)) décrivant l'ensemble du marché (E), - et d'un modèle statistique du comportement joint au temps (t) et par tranches de temps (Atk(k=]-,-"')) des variables de marché (Xi(t)"-'...P)), décrit par sa densité de distribution jointe (fa(Xl,...,XP; t; Atk)). Le modèle statistique du comportement joint au temps (t) et par tranches de temps (Atk - l...m)) des variables de marché est justifié à partir d'observations historiques - il peut s'agir notamment de lois normales, log-normales, etc., ou plus généralement d'une distribution statistique jointe des p variables d'environnement (Xi(tk)I' - 1 ... )) ou de leurs variations temporelles (AXi(tj' - 1---P)). Il est validé par des tests statistiques classiques adéquats (moyennes, variances, corrélations, sauts, comportements extrêmes, etc.). Au niveau d'un site centralisé émetteur (U), on génère numériquement, répétitivement à des instants (to) successifs, à partir du modèle statistique retenu (fa(Xl,...,XP; t; Atk)) du comportement joint dans le temps (t) des variables de marché (Xi -1 ...P)), une snip série probabiliste (Wi,lik) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s du marché (E). Selon un premier mode de mise en aeuvre, cette snip série probabiliste (Wijk) est une snip série probabiliste (Xi,ik) de pseudo-états formée de sxni p-uplets d'échantillons multi-valués de pseudo-états futurs (Xi,!(tk)I`=1...P'j-l...s,k-1-..n)@ comprenant sxrrixp valeurs numériques possibles des p variables de marché (Xi(t)('-1 ."P)). Elle est constituée - pour chacun des ni (m >_ 1) horizons de temps (tk(k-1..."')) de la simulation, - de s (avec s très grand 1) pseudo-états futurs possibles des p variables de marché (Xi(tk)Ir 1...P)) . Selon un second mode de mise en oeuvre, cette snip série probabiliste (Wiik) est une snip série probabiliste variationnelle (AXi,ik) formée de sxni p-uplets d'échantillons multi-valués de pseudo-variations (AXi,j(tk)('-1"'P'J-l ...s, k=l...m) futures à partir de l'état courant des variables de marché (Xi(to)('-1 ...P)), comprenant sxmxp valeurs numériques de pseudo-variations futures probables des p variables d'état (Xi(tk)('-1 '..P). Elle est constituée - pour chacun des m (m >_ 1) horizons de temps (tk(k= l..@n')) de la simulation, - de s (avec s très grand 1) pseudo-variations numériques possibles (AXi,j(tk)(i-1 ...s)) des p variables d'état (Xi(tk)(i-1--.P)). Le procédé de l'invention consiste alors, au niveau du site centralisé émetteur (U) très distant de la majorité des n opérateurs financiers locaux distants (S"1u-1-..n)), répétitivement à des instants (ta) en ce qu'on diffuse ladite snip série probabiliste (Wi,ik) multi-valuée de description des pseudo-états ou pseudo variations probables futur(e)s comprenant sxmxp pseudo-valeurs numériques constituée par: - (cas 1) ladite smp série probabiliste (Xi,jk) des pseudo-valeurs des variables de marché (XtrJ(tk)(i=l...p, j=I...s, _l ...m))@ - (cas 2) et/ou ladite snap série probabiliste variationnelle (AXiik) des pseudo-variations (AXi.j(tk)(i-l...p, j=l_s, k=I...nt)), Cette diffusion de ladite sntp série probabiliste (Wi,jk) est effectuée de manière décentralisée du site émetteur (U) en direction de l'ensemble de la multitude des opérateurs locaux distants (S") . *De manière indépendante, au niveau de chacun des opérateurs locaux (Su), majoritairement très distants du site centralisé émetteur (U), on choisit m' (in `:5 in) horizons de temps, retenus parmi les ni dits horizons de temps (tk(k-1...n`)). Eventuellement, dans le cas (2) où l'on reçoit uniquement ladite smp série probabiliste variationnelle (4Xi,jk) d'échantillons variationnels (AXi,,{tk)), on calcule localement ladite smp série probabiliste (Xi,jk) des pseudo-états constituée de pseudo-échantillons d'états futurs (Xi,j(tk)(`-l...p,'=l...`,k=1..."')) correspondante, (au moins pour les m' horizons de temps (tk) retenus). Puis, pour chacun des l' horizons de temps (tk) retenus, on calcule numériquement, selon le modèle local d'évaluation (D") préalablement fixé dans chacun des s scénarios (j la valeur globale (I"jk) (Du(@,", (Xi,j(tk)) du portefeuille (P"). On effectue ce calcul : pour au moins une valeur particulière du vecteur (Â,")= (/@ul-..,;l,"y) des quantités d'actifs détenus et pour chacun des s p-uplets simulés de variables d'état (Xi,j(tk)), correspondant audit horizon de temps (tk) retenu, de la série probabiliste de pseudo-état (Xi,j(tk)J=1 ...5)). Et on détermine les s pseudo-valeurs simulées du portefeuille (Iu,j(tk,2l.u)) (Du(k",Xj(tk))v-""s)) correspondantes. On regroupe et on ordonne ces s valeurs du portefeuille (P"). Puis on calcule et on mémorise l'histogramme (T'") ainsi que la distribution de probabilité (en nombre ou pourcentage) ('f"(1" )) de la valeur (I") du portefeuille local (P") dans la configuration (#,J, selon le modèle local (D"). On en déduit, par un calcul statistique, la valeur d'au moins un paramètre statistique de réglage choisi (V(k")) (tel que : moyenne, écart type, Value at Risk, utilité, etc...) caractéristique de la distribution de probabilité (`Y(;,",I)) des valeurs (I"(tk,(2#-")) en fonction des variations du vecteur (k.") = (#"1, ...,I."q) décrivant le portefeuille (P"). Eventuellement, on effectue localement, au niveau de chaque opérateur local (S"), une optimisation numérique (selon des critères prédéfinis) du (ou des) paramètre(s) statistique(s) de réglage choisi (V(7@")), et on ajuste la composition du portefeuille (@.") _ (@,"1, ...,@,"y) de chaque opérateur local (S") à une valeur optimale calculée en fonction des critères préétablis, afin notamment d'anticiper le comportement probable futur et le rendement de chaque portefeuille (P") soumis au marché (E) au cours du ou des horizon(s) de temps (tk) retenu(s). Selon une deuxième sous-variante d'application du procédé selon l'invention au domaine de l'industrie financière, pour évaluer la valeur (I"jk) du portefeuille (P") d'un opérateur financier (S") à la date future (tk) selon le modèle (D") et dans l'état du marché simulé (Xj(tk)) correspondant au scénario j, on prend en compte dans le modèle (D") les états de marché simulés, dans le même scénario j, aux dates ultérieures (ti(i=k+l...m)), selon le principe suivant - on définit pour chaque actif financier ((9h)("=1 @)) une quantité maximale (V,,) qu'il est possible de traiter (à l'achat ou à la vente) dans une journée, - si aucune des quantités (A,",,) d'actifs (6,,) détenus dans les portefeuilles (P") n'excède la valeur limite (V,,), alors on applique la formule (***) pour évaluer, à la date (tk) et dans la simulation (j), la valeur (1",;k) du portefeuille (Pu) - si une des quantités (1.üh) d'actifs (Oh) détenus dans un des portefeuilles (P") excède en valeur absolue la valeur limite (Vh), alors on remplace la quantité (Auh) par (4), affectée de son signe, à la date (tx) et on reporte le reliquat (Auh Vh) sur la date suivante (tk+1) en utilisant le vecteur de variables de marché (Xj(tk+,)) ; - si cette quantité est encore supérieure à la limite (Vh) en valeur absolue, on la fixe à la valeur maximale (vh), avec le même signe, et on reporte à nouveau le reliquat (Âuh 2Vh) sur la date suivante (tk+2) en utilisant le vecteur de variables de marché (Xj(tk.+2)), et ainsi de suite. Cette disposition particulière permet de tenir compte du risque de liquidité imposant des contraintes de borne supérieure à la quantité de chaque actif (B,,) que l'on est en droit de valoriser selon le modèle local (Du) à une date donnée (tk), et entrainant un report éventuel de la position virtuelle résultante dans chaque simulation sur des dates ultérieures. Applications industrielles de l'invention L'invention trouve une application industrielle dans tous les domaines tels que l'industrie manufacturière, énergétique, agroalimentaire, financière, de transport, de santé publique ou d'exploitation des matières premières, ... soumis à un environnement global commun et où il convient de mesurer et/ou de contrôler l'impact sur la production des évolutions possibles des variables d'environnement globales et locales en une multitude de sites industriels éloignés. A titre d'exemple, on peut avantageusement appliquer les enseignements de l'invention pour l'optimisation de l'arrosage de cultures, et pour le calcul des temps de transport routier ou la gestion d'une flotte automobile. Ainsi, une première application de l'invention concerne l'agriculture. Pour cette application, les systèmes industriels locaux distants (Su@u= @...n@) sont les exploitations agricoles de différentes parcelles cultivées. Les variables d'état (X;(t)"='...P) de l'environnement sont la température de l'air, la pluviométrie, l'hygrométrie, la sécheresse du sol, etc. en différents points géographiques. Les paramètres locaux d'action (a,") (A",--,Auq ) sont principalement les caractéristiques des parcelles et des plantations. Le but de l'invention est de limiter la surconsommation d'eau en déclenchant les différents systèmes d'arrosage automatiques uniquement lorsque la probabilité que les cultures aient à souffrir de sécheresse dans les heures ou les jours à venir est supérieure à un seuil d'acceptance défini par l'exploitant de chaque parcelle. Ainsi, les systèmes d'arrosage automatique ne se déclencheront pas si la probabilité d'une averse dans les heures à venir est forte. Inversement, il pourra par exemple se déclencher automatiquement un soir si la probabilité d'une forte sécheresse le lendemain est élevée. La mutualisation entre les exploitants de la structure de génération d'échantillons de variables météorologiques et donc des frais liés au site centralisé émetteur probabiliste (M rend la mise en place de ce système de déclenchement d'arrosage automatique moins coûteuse. Une seconde application de l'invention concerne l'optimisation de trajet routier. L'objectif est d'estimer la loi de probabilité du temps d'atteinte de la destination d'un automobiliste ou d'un transporteur compte tenu de l'ensemble des facteurs aléatoires extérieurs. Pour un conducteur qui doit se rendre en un point à une heure donnée, seule l'optimisation d'un critère probabiliste complexe (probabilité d'atteindre le point fixé à l'heure donnée) peut lui permettre de réellement optimiser l'itinéraire à choisir. En effet, statistiquement, il peut arriver qu'un tronçon de route soit généralement fluide, mais que lorsqu'un ralentissement survient, celui-ci reste encombré longtemps : la vision du temps de trajet moyen (espéré) est de ce fait trop réductrice. Selon l'invention, chacun des véhicules d'un réseau routier constitue alors un système industriel local distant (Su). Ceux-ci sont soumis aux facteurs aléatoires (X;(t)"= ' -..P) qui sont par exemple : la fluidité de la circulation sur les différents tronçons routiers, la situation météorologique, la survenue d'une panne ou d'un accident... Les paramètres locaux d'action (#,") = (A,,,, ...,@.uq ) de chaque véhicule sont la situation sur le réseau routier, l'itinéraire â venir, la vitesse du véhicule, le type de véhicule et de conduite, etc. Au niveau des différents systèmes industriels locaux, l'interface (S2") de saisie des paramètres (@,ul, ...,@."9 ) peut être reliée aux différents instruments de bord du véhicule (notamment à un système de navigation par satellite de type GPS et des capteurs du véhicule) et le récepteur (Rp") de données probabilistes peut être lié à un terminal Internet embarqué. Une évolution de cette seconde application concerne la gestion d'une flotte automobile (transport en commun, transport de fret, loueur de voiture, etc.). Si chaque automobile ou système roulant est équipé d'un système de localisation (par exemple de type GPS) et d'un émetteur, il est possible de connaître la répartition de toute la flotte à un instant passé ou présent. Il est aussi possible de connaître la situation la plus probable (l'espérance) de la position des véhicules à une date future donnée. En centralisant au niveau d'un site centralisateur (B) les pseudo-valeurs simulées des indicateurs de rendement (lu,@{t@,@,u)), l'invention permet au gestionnaire du parc d'affiner sa vision sur la probabilité des évènements à venir. Ainsi, le gestionnaire du parc peut contrôler automatiquement la probabilité de la survenue d'évènements peu fréquents mais lourds de conséquences (probabilité que tous les camions d'un transporteur de fret tombent en panne simultanément, probabilité pour que tous les véhicules d'un loueur de voiture se retrouvent concentrés dans une région). Une troisième application particulièrement avantageuse de l'invention concerne la maîtrise des paramètres d'unités industrielles fonctionnant sur un même réseau et soumises à un environnement global aléatoire commun. A titre d'exemple, on peut avantageusement appliquer les enseignements de l'invention au domaine de la production d'électricité, dans un environnement dérégulé, où le prix de l'électricité est fixé par confrontation de l'offre de n producteurs (S"( - "*-" et de la demande d'un nombre élevé de consommateurs sur un marché central (le pool énergétique ). La rencontre de l'offre et de la demande peut se faire sur la base d'un planificateur central, indiquant à chaque producteur quels doivent être ses paramètres de fonctionnement. Une telle approche pose un double problème - Le planificateur central doit prendre en compte une masse considérable d'informations, non seulement des informations communes (le niveau de la demande, l'état du réseau, les prix d'énergie primaire), des informations locales interagissant fortement avec l'environnement global (le climat), mais également de multiples informations spécifiques à chaque producteur et indépendantes de l'environnement global (l'état de marche et les performances de telle turbine, les problèmes de limite de rejet polluant dans telle rivière, etc ... ). Selon les méthodes classiques, ce problème d'optimisation devient insoluble dans ces conditions, sauf à faire des hypothèses grossièrement simplificatrices sur les contraintes locales. - Chaque producteur local doit s'efforcer d'être capable de répondre à la requête du planificateur central. Ignorant par définition le comportement du modèle d'optimisation central, il se trouve confronté à un aléa de type administratif (quelle requête va adresser le planificateur central ?) en sus d'un environnement global aléatoire. Cet aléa administratif est d'autant plus difficile â modéliser qu'il ne prend en compte qu'imparfaitement les contraintes opérationnelles locales. Le risque opérationnel devient alors difficile à gérer. Dans ces conditions, l'optimum global selon les méthodes classiques devient illusoire, car il repose sur une série d'états opérationnels locaux loin d'un réel optimum. Cette approche, prédominante dans le domaine de la production d'électricité, a été analysée par les demandeurs. Il ressort de cette étude qu'elle conduit en pratique à des situations de surcapacités industrielles et énergétiques significatives. L'alternative classique pour délocaliser le problème d'optimisation est d'établir un marché de rencontre d'offre et de demande, avec un prix d'équilibre. La théorie économique montre que la recherche d'optimums locaux correspond, sous certaines hypothèses très générales de convexité, à un optimum global. En matière d'électricité, le problème est techniquement ardu, car l'électricité est un bien non stockable. Il y a donc autant de marchés que de périodes de l'année. Le mouvement général récent de déréglementation des marchés de l'électricité conduit aujourd'hui à plusieurs modèles d'organisation de marché, qu'il serait fastidieux de décrire ici. Mais ils se rapprochent tous plus ou moins du modèle ci-dessous. Le pool énergétique est segmenté en (in) sous-marchés, correspondant à des horizons i temps de préavis distincts (Atk(k - 1...` . A titre d'exemple, on distingue un marché très court terme (quelques minutes) destiné à maintenir la tension sur le réseau, un marché court terme (1 heure, balancing market ) pour compenser les aléas de l'offre et de la demande, un marché normal à 24 heures, etc...Chaque producteur (Sj"= '...") (usine de production électrique) doit s'engager, pour chaque sous-marché, sur le respect d'une courbe de coût de production (Ck"(w)), donnant un coût de production pour chaque niveau de capacité (w) offert au réseau (ce coût pouvant être infini, au-delà d'un niveau maximum de capacité). Le pool énergétique classe sur chacun des sous-marchés les offres de capacité par coût croissant et les confronte avec la demande sur ce sous-marché : le prix (Pk(k - .."')) du kilowattheure sur le sous-marché résulte de cette confrontation. Les producteurs dont le coût de production est supérieur au prix (Pk) ne sont pas appelés par le gestionnaire du réseau. Si un producteur, suite à une indisponibilité accidentelle, constate, un temps (Atk) avant la livraison, son incapacité à respecter son engagement sur le sous-marché 1 (avec Atk il doit racheter les quantités correspondantes sur le sous-marché (k) (ou tout autre sous-marché ayant un préavis plus court). Dans ce cadre d'application de l'invention, l'indicateur de rendement du producteur est sa marge brute. Les paramètres de commandes d'un producteur (Su(u= 1...n) sont notamment - l'allocation de sa capacité disponible (Cap") entre les in sous-marchés et un n1+1-ième sous-marché spécifique qu'est l'entretien programmé de la centrale. - le choix de l'énergie primaire à utiliser, lorsqu'il a cette faculté (gaz / charbon, fuel / charbon par exemple). Le choix des paramètres est soumis à des contraintes opérationnelles locales tel que le temps et le coût minimum de démarrage d'une unité ou de changement de type de combustible, les minima et maxima techniques, les limites de rejets polluants dans l'air et dans l'eau, etc. Le problème d'optimisation local est complexe, dans la mesure où il fait appel à la prévision sur la base de variables aléatoires relevant non seulement d'un environnement spécifique à l'installation (disponibilité des installations de production), mais aussi d'un environnement global commun à tous les producteurs, comprenant - le climat, et notamment la température pour une centrale thermique (affectant le rendement thermodynamique de la centrale), la pluviométrie pour une centrale hydraulique (affectant débit et réserves)... - les prix sur les m sous-marchés d'électricité, - les prix des combustibles utilisés par le producteur, etc. Cet environnement global est bien commun à l'ensemble des producteurs : les prix sont les mêmes pour tous et les variables climatiques sont non seulement fortement corrélées entre les différents sites mais interagissent de surcroît avec les prix (on conçoit par exemple qu'une journée très froide prévue à 24 heures se traduise par une envolée du prix de l'électricité sur le marché à 24 heures). Le problème d'optimisation global par un planificateur central est remplacé par n problèmes d'optimisation locaux, afin par exemple de maximiser l'espérance de marge brute et de garantir à x% la marge brute minimale pour supporter la charge de la dette (typiquement un calcul de Value-at-Risk), ceci impliquant des calculs complexes. Selon l'art antérieur, les solutions possibles sont les trois suivantes - Chaque producteur contracte, après mise en concurrence, avec un organisme ( asset manager ) disposant des outils d'analyse de l'environnement global, qui dictera le choix des paramètres en contrepartie d'une garantie d'une marge brute prédéterminée ( tolling agreement ). Cette solution permet d'éliminer par mise en concurrence une éventuelle inefficacité d'un planificateur central, mais pose fondamentalement le même type de problèmes : les accords de ce type ne peuvent jamais parfaitement prendre en compte les contraintes opérationnelles locales. Elle conduit par ailleurs â communiquer à l'asset manager des informations confidentielles sur les contraintes de fonctionnement du site. - Chaque producteur effectue son optimisation sur la base d'une modélisation analytique de l'environnement global et de son impact sur son équation de rendement. Ceci se fait au prix de fortes hypothèses simplificatrices qui, lorsqu'elles ne sont pas respectées, conduisent à des situations de crise (pertes très élevées pour le producteur, déficit de production globale, etc ... ). - Chaque producteur simule tous les états possibles de l'environnement (global et spécifique pour chacun des horizons de temps) correspondant aux sous-marchés d'électricité. Il calcule sa marge brute pour chaque état et chaque trajectoire de paramètres, et sélectionne le vecteur de paramètres optimisant le couple Espérance / VaR sur sa marge brute. La simulation des états possibles de l'environnement global implique la mise en oeuvre de calculs complexes, mobilisant des moyens informatiques et des compétences coûteux, à l'évidence dupliqués entre les différents producteurs. La mise en eeuvre des enseignements de l'invention permet à chaque producteur d'électricité de choisir la troisième approche, c'est-à-dire la seule permettant de faire une recherche précise de l'optimum local, tout en mutualisant la partie complexe des calculs. Selon l'invention - Il collecte répétitivement (avec une fréquence correspondant à la fréquence de cotation des prix sur le pool) ladite sntp série probabiliste (W;jk) multi-valuée de description des pseudo-états probables futurs â sxmxp valeurs numériques (W ;.Jl' X;.,@ global) qui lui est fournie par le centralisé émetteur (M. - Pour chaque sfnp série probabiliste , et pour chaque vecteur possible de paramètres de commandes, il tire une fonction aléatoire de défaut reflétant la probabilité de défaillance technique de l'installation. Il en déduit pour chaque état et chaque vecteur de paramètre de commande une marge brute (MB"(W;,lz,#." )), calculée par une fonction arithmétique simple. - Il en déduit par un tri simple classant pour chaque vecteur (#,") les tirages de sa marge brute (MB" ), une fonction VaR"(@,") et une fonction espérance de marge brute Enth"(;,"). Le problème d'optimisation complexe en environnement aléatoire est ramené pour une centrale (Su(u - I...n)) à un problème d'optimisation simple : rechercher le maximum de Ernb"(@.") sous la contrainte VaR"(kn) Vn,ar. Avantages de l'invention L'invention permet de minimiser les moyens informatiques et humains nécessaires à la mise en oeuvre de simulations probabilistes sur un ensemble de sites industriels distants soumis à un environnement aléatoire commun et vue d'aboutir notamment à - une optimisation de la production en respectant des contraintes réglementaires, - une optimisation des coûts de la production en contrôlant la survenue d'évènements défavorables, - un contrôle et une minimisation de la survenue d'évènements rares défavorables. Un autre avantage de l'invention est de simplifier le calcul du risque global d'un groupe de systèmes industriels locaux distants. Un troisième avantage de l'invention est de permettre des calculs de simulation statistique locaux au niveau d'une multitude de sites industriels à partir de tirages aléatoires mutualisés sans qu'aucun des systèmes locaux n'ait à transmettre d'informations au site centralisé émetteur sur sa situation ou le réglage de ses paramètres industriels. Ceci permet d'assurer une meilleure confidentialité industrielle et commerciale. En outre, l'exemple d'application dans le domaine de l'énergie électrique donné ci-dessus permet de constater les avantages suivants pour un des producteurs industriels - Le producteur peut se rapprocher de l'optimum local, en mettant en oeuvre simplement et avec peu de moyens locaux, les techniques d'optimisation les plus sophistiquées. Ceci se traduira par un temps de retour plus rapide sur investissement ou par une baisse des prix au profit des consommateurs. - Il n'a pas besoin d'investir dans de coûteux moyens de modélisation et de simulation de l'environnement global. II peut se concentrer sur la simulation des paramètres spécifiques locaux, indépendants de l'environnement global. - Il n'a pas besoin de communiquer à un organisme central ses paramètres de fonctionnement. Dans un environnement concurrentiel, il peut donc préserver la confidentialité de sa stratégie. L'invention ayant maintenant été décrite, et son intérêt justifié sur des exemples détaillés, les demandeurs s'en réservent l'exclusivité pendant toute la durée du brevet, sans limitation autre que celle des termes des revendications ci-après. REVENDICATIONS 1) Procédé de simulation stochastique répétitive hiérarchisée et mutualisée par échantillonnage de scénarios probables, pour la prédiction probabiliste du comportement futur et/ou l'optimisation des paramètres de rendement d'une multitude de systèmes industriels locaux distants (S"(" - 1 ... n)) o très distants les uns des autres, o soumis à un environnement global commun (E) évolutif, pouvant être décrit numériquement par p variables d'état (X;(t)"-i P), constituant les composantes d'un vecteur d'état (X(t)) (X;(t)(` - '- - -P), o et dont l'efficacité de chaque système local (Si,(" - i...n)) peut en outre être décrite par un indicateur numérique de rendement (I"""- '...n)), selon un modèle numérique systémique local paramétré (D") préalablement établi (I") (D"(?.",X(t))), dépendant pour chacun des systèmes (S") et à chaque instant (t) 1 de toutou partie des variables d'état (X;(t)"- '..-n) du vecteur d'état (X(t)), 1 et d'un vecteur de q paramètres (réglables) locaux d'action (k,") ce procédé étant caractérisé en ce qu'en combinaison u au niveau d'un site centralisé émetteur (n, très distant (de plusieurs kilomètres, voire de plusieurs centaines de kilomètres) d'une majorité des n systèmes industriels locaux distants (S"(" - I ...n)) : 1 on génère répétitivement de manière centralisée à des instants (to) successifs une snip série de simulations probabilistes (W,J#" ('-i..-P''-i...s'k-i..."')) multi-valuée de description des variables (X;(t)(' - "-,P)) de l'environnement (E) constituée a. soit de pseudo-états probables futurs des variables (Xi(t)"= '---P') de l'environnement (E), b. soit de pseudo-variations probables futures de ces variables (X;(t)(`-i...n)) de l'environnement (E), communs à la multitude des n (n étant très grand 1) systèmes industriels locaux distants (S,,("- i"") ) décentralisés, et simulés à un (ou plusieurs iii, avec 1) horizons de temps (tk(k - 1 ... n') > t'), 1 et on télétransmet répétitivement cette snip série de simulations probabilistes (WJk) multi-valuée de description des pseudo-états ou pseudo- variations probables futur(e)s des variables de l'environnement, ce de manière répartie décentralisée # du site émetteur (U), # en direction de l'ensemble de la multitude des systèmes locaux distants (S"); u et, de manière indépendante, au niveau de chacun des systèmes locaux (S"), majoritairement très distants du site centralisé émetteur (U), répétitivement on collecte ladite snip série probabiliste (W,ik) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s, 1 on effectue un traitement numérique statistique local spécifique sur ladite série probabiliste (W;,jk) reçue pour prédire le comportement futur probabiliste des systèmes industriels locaux (S"), et (éventuellement) on règle le niveau du vecteur des paramètres locaux d'action (#.") = (A"i, ...,de chaque système local (S") à une valeur optimale calculée en fonction des critères préétablis, afin notamment d'anticiper le comportement probable futur et l'efficacité de chaque système (S") soumis à l'environnement (E) au cours du ou des horizon(s) de temps (tk) retenu(s). 2) Procédé selon la revendication 1 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables pour la prédiction probabiliste du comportement futur et/ou l'optimisation des paramètres de rendement d'une multitude de systèmes industriels locaux distants (Su), ce procédé étant du type spécifique consistant a au niveau d'un site expert (U'), très distant (de plusieurs kilomètres, voire de plusieurs centaines de kilomètres) d'une majorité des n systèmes industriels locaux distants (Su(u = I ... n # à faire préalablement le choix 1 des p variables d'environnement (X;"='---P)) de l'environnement (E), et d'un modèle statistique du comportement joint au temps (t) et par tranches de temps (Atk(k - '...m)) des variables d'état (X;(t)"='".P), décrit par sa densité de distribution jointe (tXI,.... XP; t; Atk)) # justifié à partir d'observations historiques - il peut s'agir notamment de lois normales, log-normales, etc., ou plus généralement d'une distribution statistique jointe des p variables d'environnement (X;(tk) '' P) ou de leurs variations temporelles (AX;(tk)(`='-..P) # et validé par des tests statistiques classiques adéquats (moyennes, variances, corrélations, sauts, comportements extrêmes, ...), o au niveau d'un site centralisé émetteur (U), très distant (de plusieurs kilomètres, voire de plusieurs centaines de kilomètres) d'une majorité des n systèmes industriels locaux distants (Su(' = 1...n)) 1 â générer numériquement répétitivement à des instants t,, successifs, â partir du modèle statistique retenu (fa(XI,...,XP; t; Atk)) du comportement joint dans le temps (t) des variables d'état (X;), une snip série probabiliste (W;,ik ('=I...P,j=I...s, k=i...in)) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s de l'environnement (E) constituée # cas 1 : d'une snip série probabiliste (X;,jk) de pseudo-états formée de sxni p- uplets d'échantillons multi-valués de pseudo-états futurs Comprenant sxnixp valeurs numériques possibles des pseudo-états futurs probables des p variables d'état (X;(t)"='...P)) et constituées # pour chacun des in (ni >_ 1) horizons de temps (tk(k='...m)) de la simulation, # de s (avec s très grand 1) pseudo-états futurs possibles (X;,,{tk)Ii ... s)) de chacune des p variables d'état (Xt(tk)(' .--P) , # cas 2 : et/ou éventuellement, d'une snip série probabiliste variationnelle (AX;,Jk) formée de sxni p-uplets d'échantillons multi-valués de pseudo-variations (AX;,j(tk)''='...P''='...s' k='...") futures à partir de la valeur courante des variables d'état (X;(to)"='---P)), comprenant sxirixp valeurs numériques de pseudo- variations futures probables des p variables d'état (X;(tk)''=',* P)) et constituées # pour chacun des m (m _> 1) horizons de temps (tk(k='..."')) de la simulation, # de s (avec s très grand 1) pseudo-variations numériques futures possibles (AX ;,j{tk)Ii ... s) de chacune des p variables d'état (X;(tk)'' ... P), le procédé étant caractérisé en ce qu'en combinaison u au niveau du site centralisé émetteur (U) très distant de la majorité des n systèmes industriels locaux distants (Su(" = I...n) ), répétitivement à des instants (to) on diffuse ladite snip série probabiliste (W;j@) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s, comprenant sxnixp pseudo- valeurs numériques constituée par a. (cas 1) ladite snip série probabiliste (X;,,@) des pseudo-valeurs des variables d'état (Xi,j(tk)"='...P.1=I...s, k=I...ni)), b. (cas 2) et/ou ladite snip série probabiliste variationnelle (dX;,ik) des pseudo- variations (AX;,j(tk)(`='...P,j=i...s,k=1...n), ce de manière décentralisée 0 du site émetteur (U), # en direction de l'ensemble de la multitude des systèmes locaux distants (S") ; u et, de manière indépendante, au niveau de chacun des systèmes locaux (S"), majoritairement très distants du site centralisé émetteur (U), répétitivement on collecte ladite snip série probabiliste (W;,jk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s à sxitixp valeurs numériques (Wi,jk Xr.jk ou OXi,jk), on choisit m'(m'5 m) horizons de temps, retenus parmi les m dits horizons de temps (tk(k = I ..."t)) , éventuellement, dans le cas (2) où l'on reçoit uniquement ladite smp série probabiliste variationnelle (OXi,jk) d'échantillons variationnels (AX;j(tk)), on calcule localement ladite snip série probabiliste (X;,jk) des pseudo-états constituée de pseudo-échantillons d'états futurs (X;,j(tk)"-'...P'j_'..-` k-'..."`) correspondante (au moins pour les m' horizons de temps (tk) retenus), pour chacun des i' horizons de temps (tk) retenus # on applique numériquement le modèle systémique local (D") préalablement fixé afin de déterminer dans chacun des s scénarios j = 1,...,s l'indicateur de rendement Iujk Du(ku,Xj(tk)) # pour au moins une valeur particulière (Xu)Q (@.",,...,@,"y) du vecteur (a,") des paramètres locaux d'action, # et pour chacun des s p-uplets simulés de variables d'état (X;,j{tj'-'...n), correspondant audit horizon de temps (tk) retenu, de la série probabiliste de pseudo-état (X;, j(tk)U = 1 ... S)), et on détermine les s pseudo-valeurs simulées de l'indicateur de rendement (lu,j(tk,@u)) _ correspondantes, # on regroupe et on ordonne les s valeurs correspondantes de l'indicateur de rendement selon le niveau de rendement (1,j), # on calcule et on mémorise l'histogramme (I'") ainsi que la distribution de probabilité (en nombre ou pourcentage) (LF"(I")) du niveau atteint de l'indicateur de rendement local (I") dans la configuration (@."), selon le modèle local (Du), # on en déduit par un calcul statistique la valeur d'au moins un paramètre statistique de réglage choisi (V(Â.u)) (tel que : moyenne, écart type, Value-at- Risk, utilité, etc...) caractéristique # de la distribution de probabilité ('F(î,",I)) des niveaux de rendement (I"(tk,;,u), # en fonction des variations du vecteur des paramètres locaux d'action (éventuellement) on effectue localement, au sein de chaque système local (S" ), une optimisation numérique (selon des critères prédéfinis) du (ou des) paramètre(s) statistique(s) de réglage choisi (V(,,")), et on règle le niveau du vecteur des paramètres locaux d'action (#.") = (/1u@,...,/luG) de chaque système local (S") à une valeur optimale calculée en fonction des critères préétablis, afin notamment d'anticiper le comportement probable futur et l'efficacité de chaque système (S") soumis à l'environnement (E) au cours du ou des horizon(s) de temps (tk) retenu(s). 3) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison u au niveau du site expert (U'), on associe â chaque paramètre local d'action (1,"h@h=' _`') un code générique (L,,) (ne dépendant pas du système local (S")), 1 on établit pour chaque code (Lh), et on actualise régulièrement, un évaluateur de rendement (IIh), c'est-à-dire un programme (matérialisé sous forme d'un fichier numérisé) permettant de calculer, en fonction du p-uplet des variables d'état (X) (X;(`- '...P') et du paramètre local d'action (,,uh), la composante correspondante du rendement (luh) = (Ilh(l#.uh,X)), 1 de telle sorte que, pour chaque système local (Su), ledit modèle systémique local (Du) consiste simplement à additionner les composantes (luh) de l'indicateur de rendement correspondant à chaque paramètre local d'action (fui,) o au niveau d'un site centralisé émetteur (U): on génère numériquement et on émet, répétitivement à des instants (to) successifs, à partir du modèle statistique retenu (fa(X,,...,XP; t; At)) du comportement joint dans le temps (t) des variables d'état (X;), la snip série probabiliste (Wi,jk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s, et on émet, à l'apparition de chaque nouveau code (Lh) ou de chaque modification de l'évaluateur de rendement (IIh), la paire (Lh,IIh) sous forme d'un fichier numérisé, o au niveau de chaque système industriel local distant (Su) on calcule # pour chacun des ni' (ni ` horizons de temps (tk) de la simulation retenus à l'échelon local (Su), # pour chacun des = 1, .... s (avec s très grand 1) pseudo-états numériques possibles (Xjk(j-@`t) reçus, ou reconstitués à partir de la série des pseudo- variations (OX,jk), # pour chacun des h = 1,...,q paramètres locaux d'action (@,",,), # la composante de rendement (luijk) (IIh(Auh,Xjk)), # puis, par sommation, selon la formule (*), l'indicateur de rendement (I"jk), # le paramètre statistique retenu (V(ku)), 1 on applique le vecteur de paramètres locaux d'action (ku) qui optimise la valeur du paramètre statistique (V). 4) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison a au niveau du site expert (U'), on associe à chaque paramètre local d'action (Auh(h- 1...`i un code générique (Li,) (ne dépendant pas du système local (Su)), 1 on établit pour chaque code (Lh), et on actualise régulièrement, un évaluateur de rendement marginal (ich), c'est-à-dire un programme (matérialisé sous forme d'un fichier numérisé) permettant de calculer, en fonction du p-uplet des variables d'état (X)=(X;"-@ P')), la composante correspondante du rendement marginal (l,h) _ (nh(X)), de telle sorte que, pour chaque système local (Su), ledit modèle systémique local (Du) consiste simplement à effectuer un produit scalaire entre le vecteur de paramètre locaux d'action (ku) =(Auh(h "*"@) et le vecteur des rendements marginaux (Il a au niveau d'un site centralisé émetteur (U) 1 on génère numériquement et on émet, répétitivement à des instants (to) successifs, à partir du modèle statistique retenu (fa(X,,.... XP; t; 4t)) du comportement joint dans le temps (t) des variables d'état (Xi), la snîp série probabiliste (Wijk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s, et on émet, à l'apparition de chaque nouveau code (Lh) ou de chaque modification de l'évaluateur de rendement (nh), la paire (Lh, irh) sous forme d'un fichier numérisé, 1 et on calcule # pour chacun des m horizons de temps (tk), # pour chacun des j = 1,...,s (avec s très grand 1) pseudo-états numériques possibles (Xjk (j= 1'..s), # pour chacun des h = 1,...,q paramètres locaux d'action (A"h), # la composante de rendement marginal (I hjk) _ (nh(Xjk)), et on diffuse la snîq série probabiliste (W h,jk) multi-valuée de description des pseudo-états (1 hjk) ou pseudo-variations probables futur(e)s (01 hjk) , u au niveau de chaque système industriel local distant (S") on collecte la snîq série probabiliste (W;,j k) multi-valuée de description des pseudo-états (l'yjk) ou pseudo-variations probables futur(e)s (AI 'hjk) puis, par sommation, selon la formule (**), on en déduit la série des indicateurs de rendement (I"jk) et le paramètre statistique retenu (V(k" )), 1 on calcule, puis on applique le vecteur de paramètres locaux d'action (Xu) qui optimise la valeur du paramètre statistique (V). 5) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison a au niveau d'un site centralisé émetteur (M: 1 on génère numériquement répétitivement à des instants (to) successifs, à partir du modèle statistique retenu (fa(XI,...,XP; t; Dt)) du comportement joint dans le temps (t) des variables d'état (Xi('= I ...P), une série probabiliste variationnelle (AXi,jk) formée de sxwt p-uplets d'échantillons multi-valués de pseudo-variations (AXi,j(tk)c'=1...p'j=l...s' k=i..."') d'états futurs à partir de l'état courant des variables d'état (Xi(to)"=I ..-P), comprenant sxnîxp valeurs numériques de pseudo-variations futures probables des p variables d'état (Xi(tk)"- I ...P) et constituées: # pour chacun des ni (ni ? 1) horizons de temps (tk(k-1..."')) de la simulation, # de s (avec s très grand 1) pseudo-variations numériques possibles (AXi, j(tk)" =1 - . -s) des p variables d'état (Xi(tk)(i -P'), i on diffuse de manière répétitive ladite sînp série probabiliste variationnelle (AXi,jk) des pseudo-variations d'états (OXi,j(tk)l'=I...P, j=l...s, k=1 _.."`)), ce de manière décentralisée # du site émetteur (U), # en direction de l'ensemble de la multitude des systèmes locaux distants (S" ), a et, de manière indépendante, au niveau de chacun des systèmes locaux (Su) on collecte et on mémorise de manière répétitive à chaque instant (to), selon une fréquence déterminée, voire à chaque modification de l'une des variables (Xi), la valeur courante des variables d'état (Xi(to)('=l-@-P), on collecte ladite smp série probabiliste variationnelle (4Xiik) des variations d'états (dXi,j(Atk)l'=I ...P,j=l...s, k=1...`)) en fonction de l'incrément de temps Atk , émise par le site central (n, on calcule localement ladite snîp série probabiliste (Xi,jk) de pseudo-états (Xirj(tk)(`-I...P.j-l...s' k=1--) aux instants futurs (tk) (t, + Atk), correspondant à la série variationnelle par calcul numérique incrémental du type Xi,j(tk) =Xi(to) (1 + AXi,j(tk)) ou toute autre formule (addition, exponentielle, etc...) faisant intervenir # ladite smp série probabiliste variationnelle (AXi,jk) émise par le site centralisé émetteur (U), # et la valeur courante des variables d'état (Xi(t)('='...p)). 6) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison u au niveau du site centralisé émetteur (U): 1 on génère numériquement répétitivement à des instants (tj successifs, à partir du modèle statistique retenu (fa(Xl,...,Xp; t; At)) du comportement joint dans le temps (t) des variables d'état (Xi), une snip série probabiliste variationnelle (OXi,jk) formée de sxin p-uplets d'échantillons multi-valués de pseudo-variations (AXt'J(tk)(ï=I...p, j=l...s, k=l...ne)) d'états futurs à partir de l'état courant des variables d'état (Xi(to)('='---p)), comprenant sxmxp valeurs numériques de pseudo-variations futures probables des p variables d'état (Xi(tk)(`='...p)) et constituées # pour chacun des in (ni >_ 1) horizons de temps (tk(k='..."')) de la simulation, # de s (avec s très grand 1) pseudo-variations numériques possibles (AXi,j(tk)"='...5)) des p variables d'état (Xi(tk)(i='---p), 1 on collecte de manière répétitive et on mémorise la valeur de l'état courant des variables d'état (Xi(to)(' ... p), on diffuse de manière répétitive et décentralisée du site émetteur (U) en direction de l'ensemble de la multitude des systèmes locaux distants (Su ), à la fois # la snip série probabiliste variationnelle (OXi,jk) des pseudo-variations (AXi,j(tk)(i-l...p, j=l...s, k=l...nt) ' # et la valeur courante des variables d'état (Xi(to)(` _' p). 7) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables caractérisé en ce qu'en combinaison a préalablement, et en différents points de l'environnement (E), on installe une série de r capteurs (y, (h=' .-.')) observateurs de l'état d'environnement (E), émettant à chaque instant (t), r mesures de l'état de l'environnement (ui,(t)(h ...r)), o au niveau du site concepteur (U'), on définit préalablement un modèle numérique descriptif d'environnement (D), définissant l'ensemble de p variables d'état (Xi(t)('='...p) _ (@(ul(t),...,,ur(t))) en fonction des mesures (ph(t)(''='',-')) (avec généralement p S r), en sorte que l'ensemble des variables d'environnement (Xi(t)(`='...p)) décrivent physiquement avec une précision acceptable, à chaque instant d'échantillonnage (tk), la globalité de l'état de l'environnement (E), o au niveau du site centralisé émetteur (n: 1 on collecte à distance et on mémorise les r mesures d'environnement (ui,(t)(h='...r))ces r mesures constituant à chaque instant d'échantillonnage (tj les composantes d'un vecteur temporel de mesure d'environnement (M(tk)) (uh(tk)(h ... ')), et on calcule, à chaque période d'échantillonnage (tj successive, à l'aide du modèle numérique descriptif d'environnement (D) et du vecteur de mesures (M(t)), les p variables d'état (Xi(t)('='...p) (D(ul(t),...,,ur(t))). 8) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison u au niveau du site centralisé émetteur (U): on maintient et on met â jour un historique des valeurs des p variables d'état (Xi(t)('='...p) de l'environnement global commun (E) au cours du temps. 1 on génère numériquement répétitivement à des instants (tj successifs, à partir # du modèle statistique retenu (fa(X,,...,XP; t; Ot)) du comportement joint dans le temps (t) des variables d'état (X;(' -' ".Pl), # et de cet historique mémorisé des valeurs des p variables d'état (X;(t)('-'...Pl), la smp série probabiliste (W,à multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s de l'environnement (E), 1 on diffuse ladite smp série probabiliste (W,J@) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s de manière décentralisée du site émetteur (U) en direction de l'ensemble de la multitude des systèmes locaux distants (S"). 9) Procédé selon la revendication 8 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison u au niveau du site centralisé émetteur (U) on maintient et on met à jour un historique des valeurs des p variables d'état (X;(t)"-'---Pl) de l'environnement global commun (E) au cours du temps, on génère numériquement répétitivement à des instants (t ) successifs, à partir # du modèle statistique retenu (fa(X,,.... XP; t; At)) du comportement joint dans le temps (t) des variables d'état (X;('-'...Pl), # et de cet historique mémorisé des valeurs des p variables d'état (X;(t)(' - I ..-Pl), la snap série probabiliste (WJ@) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s de l'environnement (E), selon la méthode dite de Monte-Carlo , par tirages pseudo-aléatoires des pseudo-valeurs (WJk) grâce à un générateur de nombres situé sur le site centralisé émetteur (U), voire sur le site concepteur (U' ), reposant sur une méthode # soit pseudo-aléatoire de nature ergodique, # soit à discrépance faible (suites de Sobol, Hamersley, etc.), # soit sur toute autre méthode produisant des nombres asymptotiquement répartis selon une distribution spécifiée (par exemple uniforme sur un intervalle), 1 on diffuse ladite snap série probabiliste (W,ik) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s de manière décentralisée du site émetteur (U) en direction de l'ensemble de la multitude des systèmes locaux distants (S" ). 10) Procédé selon la revendication 8 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison u au niveau du site centralisé émetteur (U): 1 on maintient et on met â jour un historique des valeurs des p variables d'état (X;(t)l' - ' ... P) de l'environnement global commun (E) au cours du temps, 1 on génère numériquement répétitivement à des instants (to) successifs, à partir # du modèle statistique retenu (fa(X,,...,XP; t; At)) du comportement joint dans le temps (t) des variables d'état (X;('-'...Pl), # et de cet historique mémorisé des valeurs des p variables d'état (X;(t)" -' Pl), la smp série probabiliste (Wjk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s de l'environnement (E), selon la méthode dite Historique d'estimation d'états futurs probables des p variables d'état (X;(t)'`-'..-P) reconstruits à partir de l'état courant (X;(to)(`-'---P@) et de variations historiques passées (OX;(t j)"-'...P,'-'---" pour un ensemble de dates passées t.), ce en affectant à chaque état reconstruit une probabilité de réalisation pouvant dépendre de l'âge de la variation historique correspondante, et on diffuse ladite srnp série probabiliste (Wj) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s de manière décentralisée du site émetteur (U) en direction de l'ensemble de la multitude des systèmes locaux distants (Su). 11) Procédé selon la revendication 8 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison o au niveau du site centralisé émetteur (U) 1 on maintient et on met à jour un historique des valeurs des p variables d'état (X;(t)"- '...P) de l'environnement global commun (E) au cours du temps, 1 on génère numériquement répétitivement à des instants (to) successifs la snip série probabiliste (W,ik) multi-valuée de description des pseudo-états ou pseudo- variations probables futur(e)s de l'environnement (E), à partir # du modèle statistique retenu (fa(XI,...,XP; t; At)) du comportement joint dans le temps (t) des variables d'état (X;), # et de cet historique mémorisé des valeurs des p variables d'état (X;(t)('- I P)), # mais en introduisant volontairement dans la snip série probabiliste (Wjk) des biais de comportement statistique joint des variables d'état (X;(t)" - 'vis-à- vis de leur comportement statistique historique, les séries biaisées reproduisant ainsi des scénarios catastrophes , afin de tenir compte de l'impact desdits scénarios catastrophes sur les systèmes distants (Su), 1 on diffuse ladite snip série probabiliste (W;Jk) multi-valuée biaisée de description des pseudo-états ou pseudo-variations probables futur(e)s de manière décentralisée du site émetteur (U) en direction de l'ensemble de la multitude des systèmes locaux distants (Su). 12) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison o certains des systèmes distants (Su) envoient au système central (U) des requêtes sous forme de paramètres de scénarios catastrophes selon des biais spécifiques définis par eux, o au niveau du système central on génère en fonction de ces requêtes, éventuellement en fusionnant celles qui sont identiques, des snip séries probabilistes (Wjk) multi-valuées de description des pseudo-états ou pseudo-variations probables futur(e)s, biaisées en fonction des paramètres spécifiés par les systèmes distants, et on envoie du site émetteur (U) vers chaque système local distant (S") ladite snip série probabiliste (W,J@) multi-valuée biaisée de description des pseudo-états ou pseudo-variations probables futur(e)s qui correspond à sa requête. 13) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison u chaque système distant (Su) altère lui-même tout ou partie des sinp séries probabilistes (W;jk) multi-valuées de description des pseudo-états ou pseudo- variations probables futur(e)s, de façon à en faire des séries biaisées reproduisant des scénarios catastrophes dont les paramètres sont soit spécifiés au niveau du système central (U) et diffusés en direction des systèmes locaux (Su), soit définis au niveau des systèmes locaux (S"). 14) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison o au niveau du site centralisé émetteur (U) on génère numériquement la snip série probabiliste (Wjk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s de l'environnement (E) de manière périodique, selon une fréquence d'échantillonnage p). 15) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison o au niveau du site centralisé émetteur (M: on diffuse ladite snip série probabiliste (W,jk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s de manière décentralisée du site émetteur (U) en direction de l'ensemble de la multitude des systèmes locaux distants (S") par émission sous une forme broadcast , notamment par émission radio, o et, de manière indépendante, au niveau de chacun des systèmes locaux (Su), majoritairement très distants du site centralisé émetteur (U), répétitivement on collecte ladite snip série probabiliste (W,jk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s â sxmxp valeurs numériques (W,jk k ou A & â, par réception quasiment en continu de cette émission broadcast en provenance du site centralisé émetteur (U). 16) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison o au niveau du site centralisé émetteur (U): 1 on diffuse ladite smp série probabiliste (W;,jk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s, de manière décentralisée du site émetteur (U) en direction de l'ensemble de la multitude des systèmes locaux distants (S") par émission sous une forme dite PP aval , c'est-à-dire par émission sur un réseau (R) (notamment de type Internet) du site centralisé émetteur (U) à l'adresse spécifique de chacun des systèmes industriels locaux distants (Sj"-'@. n), o et, de manière indépendante, au niveau de chacun des systèmes locaux (S"), majoritairement très distants du site centralisé émetteur (U), répétitivement 1 on collecte, ladite snip série probabiliste (W,jk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s à sxirixp valeurs numériques (W,jk X;,jk ou A & i'), par enregistrement successifs dans la mémoire (H") de l'ordinateur local (Cu). 17) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison a au niveau du site centralisé émetteur (Ü): 1 on diffuse ladite snip série probabiliste (W,jk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s de manière décentralisée du site émetteur (U) en direction de l'ensemble de la multitude des systèmes locaux distants (S") par mise à disposition de ladite snip série probabiliste (W,jk) multi- valuée sur un serveur (Z) disponible # sur un réseau (R) (notamment de type Internet), )a et, # de par manière accès indépendante, PP par chacun au des niveau systèmes de industriels chacun des locaux systèmes distants locaux (S"(" = i...n) (Su) répétitivement 1 on collecte ladite snip série probabiliste (W,jk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s à sxnixp valeurs numériques (W-j#, Xi,jk ou AX;,jk), par accès sous une forme dite PP amont , c'est-à-dire par accès via le réseau (R) (notamment de type Internet) au site centralisé émetteur (U) à son adresse spécifique, par acquisition par transfert (notamment de type FTP ) de ladite smp série probabiliste (W,ik) multi-valuée et par enregistrement dans la mémoire (Hu) de l'ordinateur local (Q. 18) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison u au niveau du site centralisé émetteur (M, on génère numériquement répétitivement à des instants (to) successifs, à partir du modèle statistique retenu (fa(XI,...,Xp; t; 4t)) du comportement joint dans le temps (t) des variables d'état (X;(' ..n)), une snip série probabiliste (W;,ik) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s de l'environnement (E) constituée de s (avec s très grand 1) pseudo-valeurs numériques des p composantes (S,,j(tk)('-I...n''-I...S) pour chacun des m (m >_ 1) horizons de temps (tk(k-I...m)), sous la forme d'un fichier de sxmxp données alphanumériques successives subdivisé 1 en m sous-fichiers temporels successifs de niveau 1 regroupant les pseudo- valeurs d'un même horizon de temps (tk(k - I...n')), 1 chaque sous-fichier temporel incluant une succession de s sous-fichiers de tirage de niveau 2 constitués chacun d'un p-uplet de pseudo-valeurs. 19) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison a au niveau du site centralisé émetteur (n, on génère numériquement répétitivement à des instants (t,,) successifs, une snip série probabiliste (W,Jk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s de l'environnement (E) en séparant les in sous-fichiers temporels successifs de niveau 1 regroupant les pseudo-valeurs d'un même horizon de temps (tk(k-I..."')) par un indicateur d'horizon de temps (r7k(k - I .-.n')) correspondant, o et, au niveau de chacun des systèmes locaux (S"), on distingue les uns des autres les différents dits sous-fichiers temporels d'une même srnp série probabiliste (W,ik) multi-valuée collectée en identifiant leurs indicateurs spécifiques d'horizon de temps (Ilk(k I...n:)). 20) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison a au niveau du site centralisé émetteur (n, on génère numériquement répétitivement à des instants (tJ successifs, une sinp série probabiliste (W,jk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s de l'environnement (E), en séparant les s sous-fichiers de tirage successifs de niveau 2 (d'un même fichier temporel) regroupant chacun un p-uplet de pseudo-valeurs par un indicateur de tirage (zJ(j - I .-.S correspondant, o et, au niveau de chacun des systèmes locaux (Si,), majoritairement très distants du site centralisé émetteur (U), on distingue les uns des autres les différents sous-fichiers de tirage successifs de niveau 2 (d'un même fichier temporel) regroupant chacun un p- uplet de pseudo-valeurs d'une même snip série probabiliste (Wjk) multi-valuée collectée en identifiant leur indicateur spécifique de tirage (zj'# - I -S). 21) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison u au niveau du site centralisé émetteur (U): on diffuse dans sa globalité ladite snip série probabiliste (Wjk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s de manière décentralisée du site émetteur (U) en direction de l'ensemble de la multitude des systèmes locaux distants (Su) par émission sous une forme broadcast , notamment par émission radio, accessible à tous les systèmes locaux (Su), 1 et, de manière indépendante, au niveau de chacun des systèmes locaux (S"), on installe et on utilise un filtre sélectif préréglé de mise à disposition d'une fraction (W ,J k)' prédéfinie (notamment acquise par un abonnement spécifique) de ladite snip série probabiliste (Wik) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s, collectée par réception en provenance du site centralisé émetteur (M . 22) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison o au niveau du site centralisé émetteur (U), on diffuse de manière sélective (selon une règle d'attribution prédéfinie) 1 diverses fractions différentes ((W;ik)sis-'...h1) prédéfinies (notamment acquises par un abonnement spécifique) de ladite snap série probabiliste (Wjk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s, 1 de manière décentralisée du site émetteur (U) en direction de plusieurs groupes (G,'9='---") attributifs de la multitude des systèmes locaux distants (S"), chaque groupe (G.) recevant une fraction particulière ((W,ik)x), 1 par émission sous une forme dite PP aval , c'est-à-dire par émission sur un réseau (R) (notamment de type Internet) du site centralisé émetteur (U) à l'adresse spécifique de chacun des systèmes industriels locaux distants (S"("-'-'), a et, de manière indépendante, au niveau de chacun des groupes (G,,'s-'---") attributifs de systèmes locaux (S"), majoritairement très distants du site centralisé émetteur (U), répétitivement 1 on collecte ladite fraction prédéfinie ((W,jk)X) de la snip série probabiliste (W,ik) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s par enregistrements successifs dans la mémoire (H") de l'ordinateur local (Cu). 23) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison a au niveau du site centralisé émetteur (U) 1 on diffuse ladite stnp série probabiliste (W;,ik ) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s de manière décentralisée du site émetteur (U) en direction de l'ensemble de la multitude des systèmes locaux distants (S") par mise à disposition de ladite srnp série probabiliste (W;jk) multi- valuée sur un serveur (Z) disponible sur un réseau (R) (notamment de type Internet) par accès PP par chacun des systèmes industriels locaux distants (S"("- 1...n , 1 en mettant de manière sélective à disposition de différents groupes (Gg(g-'...")) de systèmes locaux (S") sur le serveur (Z) diverses fractions limitatives ((Wjk)K) prédéfinies (notamment acquises par un abonnement spécifique) de ladite snip série probabiliste (W,j#, ) multi-valuée de description des pseudo-états probables futurs, ce notamment par l'intermédiaire d'un système d'identificateurs et de mots de passe transmis par chaque système local (Su(" -' - - `), a et, de manière indépendante, au niveau de chacun des systèmes locaux (S"), répétitivement on collecte ladite fraction ((W,jk)s) prédéfinie attribuée de ladite snip série probabiliste (Wjk) multi-valuée de description des pseudo-états ou pseudo- variations probables futur(e)s par accès sous une forme dite PP amont , c'est-à- dire par accès via le réseau (R) (notamment de type Internet) du site centralisé émetteur (U) à son adresse spécifique, acquisition par transfert de ladite fraction ((W;,-jk)x) prédéfinie attribuée de ladite sinp série probabiliste (W,ik) multi- valuée et par enregistrement dans la mémoire (H") de l'ordinateur local (C"). 24) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de scénarios probables, caractérisé en ce qu'en combinaison 0 on choisit préalablement n'(n'S n) systèmes industriels locaux, retenus parmi les il dits systèmes industriels locaux (Su(" - i ...n , a on choisit préalablement m' (ni `S m) horizons de temps, retenus parmi les m dits horizons de temps (tk(k - 1 ... "')), identiques pour tous les systèmes locaux, a au niveau de chacun des n' systèmes locaux (S") retenus, pour chacun des m' horizons de temps (tk) retenus on applique numériquement le modèle systémique local (D") préalablement fixé afin de déterminer dans chacun des s scénarios (j = 1, .... s) l'indicateur de rendement (lujk - Du(ku, (Xi,j(tk)) # pour au moins une valeur particulière (;.u) du vecteur des paramètres locaux d'action, # et pour chacun des s p-uplets simulés de variables d'état (Xij(tk)" -' ... P'), correspondant audit horizon de temps (tk) retenu, de la série probabiliste de pseudo-états (Xj(tk)(j -' ... s'), # et on détermine les s pseudo-valeurs simulées de l'indicateur de rendement (lu.@\tk,@u)) = (Du(Xu,Xj(tk))U-'...s))COrrespondanteS, 1 on émet la série de simulations locales constituée des s pseudo-valeurs simulées de l'indicateur de rendement (I",j{tk,X"))=(D"(a,",Xj(tk))#i '"" correspondantes du système local (Su), de manière centralisante en direction d'un site centralisateur (B), a au niveau du site centralisateur (B), on collecte les sxn' pseudo-valeurs simulées d'indicateurs de rendement (I"@(tk,@,ü)@-@s@) en provenance des chacun des n' systèmes locaux (S") et pour chacune des s simulations, on calcule, pour chaque simulation (j = 1, .... s) un indicateur de rendement global (Jj) de l'ensemble des systèmes locaux â partir d'un modèle d'évaluation global (Jj) _ (D(lu.j(tk,%1,u)@u-l...nI)), on regroupe et on ordonne les s valeurs correspondantes de l'indicateur de rendement selon son niveau, on calcule et on mémorise l'histogramme (I') ainsi que la distribution de probabilité (en nombre ou pourcentage) ('I'(J)) du niveau atteint de l'indicateur de rendement global (J) dans la configuration choisie des vecteurs de paramètres locaux d'action (A) = (@ju-1...n), selon les modèles locaux (D") et le modèle d'évaluation global (D) retenus, 1 on en déduit par un calcul statistique la valeur d'au moins un paramètre statistique de réglage global choisi (V(A)) (tel que : moyenne, écart type, Value at Risk, utilité, etc...) caractéristique # de la distribution de probabilité (x(ki,.... ;.n, J)) des niveaux de rendement globaux # en fonction des valeurs du multi-vecteur (A) des vecteurs de paramètres locaux d'action (@,ü), (éventuellement) on effectue globalement une optimisation numérique (selon des critères prédéfinis) du paramètre statistique de réglage global choisi (V(A)) par rapport au multi-vecteur (A) = (T.j,.... Xn) et on transmet à chaque système local (S") le niveau correspondant de son vecteur optimisé des paramètres locaux d'action (,,"), afin notamment d'anticiper globalement le comportement probable futur et l'efficacité de l'ensemble des système locaux (S") soumis à l'environnement (E) au cours du (ou des) horizon(s) de temps (tk) retenu(s). 25) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de pseudo-états futurs probables de variables de marché, pour la prédiction probabiliste du comportement futur et/ou l'optimisation des paramètres de rendement d'une multitude de ri (n étant très grand 1) d'opérateurs financiers locaux (S"( =1...n)) o décentralisés et très distants les uns des autres, o possédant chacun un portefeuille (P,,('= 1 ... n)) constitué d'actifs financiers (9h(h=1 ...4)), 1 dont la valeur (I,,('= 1 ...n)) dépend des mouvements du marché des capitaux (E), et dont la valeur (In(u =1 ... n)) peut être calculée par combinaison # d'un vecteur de variables de marché (X(t)) = (X;(t)('= '... ') communes variant au cours du temps (t), à partir duquel il possible, pour chaque opérateur financier local distant (Su), de reconstituer avec une précision acceptable la valeur de tous les actifs financiers pouvant être présents dans son portefeuille (P" ), selon un modèle local d'évaluation (D"), (ces variables de marché pouvant représenter des prix d'actifs ou de commodités, des taux d'intérêt ou des coefficients de courbe de taux, des volatilités ou des coefficients de surfaces de volatilité, des indicateurs de défaut d'un émetteur ou d'une contrepartie, ou tout autre nombre apportant une information pertinente sur l'état des marchés de capitaux), # et, pour chaque opérateur financier local distant (S"), d'un vecteur de quantités ( U (Xuh(h 1...y)) de chacun des actifs financiers (9h(h=1 ...9)) détenus dans son portefeuille (P"), par exemple selon la combinaison linéaire où (ir"(t,X(t),6h)) représente le prix de l'actif(9j,) à l'instant (t) selon le modèle local (Du) lorsque le marché est dans l'état X(t), ou une autre formule pouvant faire apparaître une non-linéarité par rapport aux quantités (@",,), pour tenir compte notamment du risque de liquidité; ce procédé étant du type consistant o au niveau d'un site expert (U'), très distant (de plusieurs kilomètres, voire de plusieurs centaines de kilomètres) d'une majorité des n opérateurs financiers locaux distants (Su(u ), à faire préalablement le choix 1 des p variables de marché ((X;)(' - 1, -P) décrivant l'ensemble du marché (E), et d'un modèle statistique du comportement joint au temps (t) et par tranches de temps (Atk(k - 1..-m)) des variables de marché (X;(t)('=1 -P), décrit par sa densité de distribution jointe (fa(X1,.... XP; t; OtA)) # justifié à partir d'observations historiques - il peut s'agir notamment de lois normales, log-normales, etc., ou plus généralement d'une distribution statistique jointe des p variables d'environnement (X;(tk)(' ou de leurs variations temporelles (àX ;(tx)('= 1-..P)), # et validé par des tests statistiques classiques adéquats (moyennes, variances, corrélations, sauts, comportements extrêmes, ... ), a au niveau d'un site centralisé émetteur (U), très distant (de plusieurs kilomètres, voire de plusieurs centaines de kilomètres) d'une majorité des n systèmes industriels locaux distants (Su(, - 1...n)) à générer numériquement répétitivement à des instants (to) successifs, à partir du modèle statistique retenu (fa(X1,...,XP; t; Atk)) du comportement joint dans le temps (t) des variables de marché (Xi(' ...P)), une sntp série probabiliste (Wjz) multi- valuée de description des pseudo-états ou pseudo-variations probables futur(e)s du marché (E) constituée # cas 1 : d'une srnp série probabiliste (X;,jk) de pseudo-états formée de sxna p- uplets d'échantillons multi-valués de pseudo-états futurs (X;,j(tk)l'-l...comprenant sxmxp valeurs numériques possibles des p variables de marché (X;(t)('-1 ...P)) et constituées # pour chacun des m (m >_ 1) horizons de temps (tk(k= 1 ... "')) de la simulation, # de s (avec s très rand 1) pseudo-états futurs possibles des p variables de marché (X;(tk)(' ... P) , # cas 2 : et/ou éventuellement, d'une snip série probabiliste variationnelle (AX ;,j k) formée de sxm p-uplets d'échantillons multi-valués de pseudo-variations (AX`'j(tk)(i=l...p, j=l...s, k=1 ... m)) futures à partir de l'état courant des variables de marché (X;(to)(ï ...P)), comprenant sxmxp valeurs numériques de pseudo- variations futures probables des p variables d'état (X;(tk)(i= 1...P)) et constituées # pour chacun des m (ni > 1) horizons de temps (tk(k=1..."')) de la simulation, # de s (avec s très grand 1) pseudo-variations numériques possibles (OX;,j(tk)(#=1.__s) des p variables de marché (Xl(tk)('= 1...P)), le procédé étant caractérisé en ce qu'en combinaison o au niveau du site centralisé émetteur (M très distant de la majorité des il opérateurs financiers locaux distants (Si,(' l -n)), répétitivement â des instants (t,,) 1 on diffuse ladite smp série probabiliste (W;J) multi-valuée de description des pseudo-états ou pseudo variations probables futur(e)s, comprenant sxrnxp pseudo- valeurs numériques, constituée par a. (cas 1) ladite srnp série probabiliste (X;,jk) des pseudo-valeurs des variables de marché (Xi,j(tk)(i-l...p, j=1...s, k=1 ... m)) b. (cas 2) et/ou ladite snip série probabiliste variationnelle (4X;,jk) des pseudo- variations (AXi.j(tk)('-1...p, j=l...s, k=l...m))' ce de manière décentralisée # du site émetteur (n, # en direction de l'ensemble de la multitude des opérateurs locaux distants (S") ; a et, de manière indépendante, au niveau de chacun des opérateurs locaux (Si,), majoritairement très distants du site centralisé émetteur (n, répétitivement on collecte ladite snip série probabiliste (W,jk) multi-valuée de description des pseudo-états ou pseudo-variations probables futur(e)s à sxmxp valeurs numériques (Wijk Xr.jk ou A & i'), 1 on choisit nt'(m' 1...nr)) éventuellement, dans le cas (2) où l'on reçoit uniquement ladite snip série probabiliste variationnelle (OXj,jk) d'échantillons variationnels (OX;j(tk)), on calcule localement ladite snip série probabiliste (X;,jk) des pseudo-états constituée de pseudo-échantillons d'états futurs (X;j{tk)('=1 ...P''=1...`,k=1..."')) correspondante, (au moins pour les m' horizons de temps (tk) retenus), 1 pour chacun des m' horizons de temps (tk) retenus # on calcule numériquement selon le modèle local d'évaluation (Du) préalablement fixé, dans chacun des s scénarios (j = 1,...,s), la valeur globale (jujk) (Du(?l,u,Xj(tk)) du portefeuille (Pu) # pour au moins une valeur particulière (@u)= (Â.udu vecteur (@u) des quantités d'actifs détenus, # et pour chacun des s p-uplets simulés de variables d'état (X;,j(tk)), correspondant audit horizon de temps (tk.) retenu, de la série probabiliste de pseudo-état (X;,j(tk)(i 1...s)), et on détermine les s pseudo-valeurs simulées du portefeuille (Iu,j(tk,Xu)) Du(#,u,Xj(tJ'I-'...s' coriespondantes, # on regroupe et on ordonne les s valeurs correspondantes (Iu,j(t#,,#,u)Y=i._.s@) du portefeuille (P" ), # on calcule et on mémorise la distribution de probabilité (en nombre ou pourcentage) `Y"(lu) de la valeur (I") du portefeuille local (P") dans la configuration (Xu), selon le modèle local (Du), # on en déduit par un calcul statistique la valeur d'au moins un paramètre statistique de réglage choisi (V(ku)) (tel que : moyenne, écart type, Value-at- Risk, utilité, etc ... ) caractéristique # de la distribution de probabilité `h((@."),1) des valeurs (I"(t@,@,u), # et du vecteur (;.u) = ...,@,"q) décrivant le portefeuille (P"), (éventuellement) on effectue localement, au niveau de chaque opérateur local (Su), une optimisation numérique (selon des critères prédéfinis) du (ou des) paramètre(s) statistique(s) de réglage choisi (V(X")), et on ajuste la composition du portefeuille (a,.) = (Â. 1, .... Â,,,) de chaque opérateur local (S") â une valeur optimale calculée en fonction des critères préétablis, afin notamment d'anticiper le comportement probable futur et le rendement de chaque portefeuille (P") soumis au marché (E) au cours du ou des horizon(s) de temps (tR) retenu(s). 26) Procédé selon la revendication 2 de simulation stochastique répétitive hiérarchisée par échantillonnage de pseudo-états futurs probables de variables de marché, pour la prédiction probabiliste du comportement futur et/ou l'optimisation des paramètres de rendement d'une multitude de n (n étant très grand 1) d'opérateurs financiers locaux (S"t"= @...n)) o décentralisés et très distants les uns des autres, o possédant chacun un portefeuille (P"('= @...n)) ; 1 constitué d'actifs financiers dont la valeur unitaire (I h(" = '-- ` dépend des mouvements du marché des capitaux (E), et est une fonction d'un vecteur de variables de marché (X(t)) (X;(t)"=1...P) communes variant au cours du temps (t), à partir duquel il possible, avec une précision acceptable, de reconstituer la valeur de tous les actifs financiers pouvant être présents dans les portefeuilles (P") des acteurs locaux, selon un modèle commun d'évaluation (D), (ces variables de marché pouvant représenter des prix d'actifs ou de commodités, des taux d'intérêt ou des coefficients de courbe de taux, des volatilités ou des coefficients de surfaces de volatilité, des indicateurs de défaut d'un émetteur ou d'une contrepartie, ou tout autre nombre apportant une information pertinente sur l'état des marchés de capitaux), 1 et dont la valeur (I"In='...n) peut être calculée par produit scalaire entre # le vecteur des valeurs unitaires (l') = (I"' =' ...q') # et, pour chaque opérateur financier local distant (S" ), d'un vecteur de quantités (,u) ('uh(h 1...q)) de chacun des actifs financiers (6,,(h= 1...q) détenus dans son portefeuille (P") Du(ku,X) (X, 1 I@) (*x) où (1 h) = (7rh(t,X(t),6h)) représente le prix de l'actif(91,) à l'instant (t) selon le modèle commun (D) lorsque le marché est dans l'état (X(t)), ce procédé étant du type consistant u au niveau d'un site expert (U'), très distant (de plusieurs kilomètres, voire de plusieurs centaines de kilomètres) d'une majorité des n opérateurs financiers locaux distants (Su(u= 1...n) )@ à faire préalablement le choix des p variables de marché (X;('= 1 ...P) décrivant l'ensemble du marché (E), 1 et d'un modèle statistique du comportement joint au temps (t) et par tranches de temps (Otk(k '.'n) des variables de marché (X;(t)('=1-..n), décrit par sa densité de distribution jointe (fa(XI,.... Xp; t; Atk)) # justifié à partir d'observations historiques - il peut s'agir notamment de lois normales, log-normales, etc., ou plus généralement d'une distribution statistique jointe des p variables d'environnement (X;(tk)(' I@'@) ou de leurs variations temporelles (DX;(tk)(' ...P>), # et validé par des tests statistiques classiques adéquats (moyennes, variances, corrélations, sauts, comportements extrêmes, etc.), u au niveau d'un site centralisé émetteur (U), très distant (de plusieurs kilomètres, voire de plusieurs centaines de kilomètres) d'une majorité des n systèmes industriels locaux distants (S,('= ".")) ; 1 à générer numériquement répétitivement à des instants (ta) successifs, à partir du modèle statistique retenu (fa.(XI,...,Xn; t; Atk)) du comportement joint dans le temps (t) des variables de marché (Xi(' = 1--n), une smp série probabiliste (W;,jk) multi- valuée de description des pseudo-états ou pseudo-variations probables futur(e)s du marché (E), constituée # cas n' l : d'une sntp série probabiliste (X;.j) de pseudo-états formée de sxm p-uplets d'échantillons multi-valués de pseudo-états futurs (X;.j{tk)(i=I...p,;=l...s,k=1...rrt @ comprenant sxntxp valeurs numériques représentant # pour chacun des m (m > 1) horizons de temps (tk(k =1 ... ni) de la simulation, # s (avec s très grand 1) pseudo-états futurs possibles des p variables de marché (X;(tk)(' =1...Pl) , # cas n 2 : d'une smp série probabiliste variationnelle (fX;,jk) de pseudo-états formée de sxm p-uplets d'échantillons multi-valués de pseudo-variations futures (aX;,j(tk)('=I ..-Pj=l_s,k=I...n'l), comprenant sxrnxp valeurs numériques représentant # pour chacun des m (in >_ 1) horizons de temps (tk(k=1...") de la simulation, # s (avec s très grand 1) pseudo-variations futures possibles des p variables de marché (X;(tk)('= "*-n>), 1 à calculer pour chaque vecteur (X,jk) de pseudo-état des variables de marché une srnq série probabiliste (W hjk) multi-valuée de description # cas n 1 : des pseudo-prix probables futurs d'actif (I hljk (h-I...q.j=I...s,k=I...rn)) - D(tk,X,j(tk),eh) # cas n 2 : et/ou des pseudo-variations (fil h,j(tk)(h=I...q, j=l...s, k=l...rrr>) qui en découlent le procédé étant caractérisé en ce qu'en combinaison o au niveau du site centralisé émetteur (U) très distant de la majorité des n opérateurs financiers locaux distants (Si,(' "..n>), répétitivement à des instants (ta) 1 on diffuse, ladite smq série probabiliste (W h,jk) multi-valuée de description des pseudo-états ou pseudo variations probables futur(e)s, comprenant sxnlxq pseudo- valeurs numériques constituée par # (cas 1) ladite srnq série probabiliste (I',,,j) des pseudo-valeurs des prix des actifs (I h,j@ (1=1...q, j=I...s, k=l...m)) # (cas 2) et/ou ladite srnq série probabiliste variationnelle des pseudo- k (h=l...q, j=I...s, k=1 ... m)) variations (AI hj ce de manière décentralisée # du site émetteur (U), # en direction de l'ensemble de la multitude des opérateurs locaux distants (S") ; o et, de manière indépendante, au niveau de chacun des opérateurs locaux (Su), majoritairement très distants du site centralisé émetteur (U), répétitivement 1 on collecte ladite srnq série probabiliste (W 'h jk) multi-valuée de description des pseudo-états ou pséudo-variations probables futur(e)s à sxmxq valeurs numériques ( W . k , k , k h, j = I h, j Ou AI h,j ), on choisit in '(m' horizons de temps, retenus parmi les ni dits horizons de temps (tk(k=1...n,)) 1 éventuellement, dans le cas (2) où l'on reçoit uniquement ladite smq série probabiliste variationnelle (DI j,,jk) d'échantillons variationnels (@I h,j{tk)), on calcule localement ladite smq série probabiliste (1 h,jk) des pseudo-états constituée de pseudo-échantillons d'états futurs (I h,j(tk)(h-1...y, j=l...ç, k=1...n) correspondante, (au moins pour les m' horizons de temps (tk) retenus), 1 pour chacun des ni " horizons de temps (1k) retenus # on calcule numériquement selon le modèle local d'évaluation (Du) préalablement fixé dans chacun des s scénarios (j = 1,...,s) la valeur globale (Iujk) lek) du portefeuille (Pu) # pour au moins une valeur particulière du vecteur (Xu)= (/l,u,,...,@uy) des quantités d'actifs détenus, # et pour chacun des s q-uplets simulés de prix d'actifs (I h,,{tk)), correspondant audit horizon de temps (tk) retenu, de la série probabiliste de pseudo-états (I h'J(tk)U I...s)), et on détermine les s pseudo-valeurs simulées du portefeuille (I",j (tk,ku)) (@u I I@k)@=1 _.S)) correspondantes, # pour la valeur (#.u)= (Au,,...,du vecteur des paramètres locaux d'action retenus, # on regroupe et on ordonne les s valeurs correspondantes (Iu,j(tk,%,u)@- @`)) du portefeuille (Pu), # on calcule et on mémorise l'histogramme (ru) ainsi que la distribution de probabilité (en nombre ou pourcentage) ('Fu(lu)) de la valeur (lu) du portefeuille local (Pu) dans la configuration (@,u), selon le modèle local (Du), # on en déduit par un calcul statistique la valeur d'au moins un paramètre statistique de réglage choisi (V(k")) (tel que : moyenne, écart type, Value-at- Risk, utilité, etc...) caractéristique # de la distribution de probabilité @F, ,,",lu)) des valeurs (lu(tk,@.u), # en fonction du vecteur (;Q _ (@u,, ...,puy) décrivant le portefeuille (Pu), 1 (éventuellement) on effectue localement, au niveau de chaque opérateur local (Su), une optimisation numérique (selon des critères prédéfinis) du (ou des) paramètre(s) statistique(s) de réglage choisi (V(#.u)), et on ajuste la composition du portefeuille (a,u) = (Au,, .... Au.) de chaque opérateur local (Su) à une valeur optimale calculée en fonction des critères préétablis, afin notamment d'anticiper le comportement probable futur et le rendement de chaque portefeuille (Pu) soumis au marché (E) au cours du ou des horizon(s) de temps (tk) retenu(s). 27) Procédé selon la revendication 25 de simulation stochastique répétitive hiérarchisée par échantillonnage de pseudo-états futurs probables de variables de marché, pour la prédiction probabiliste du comportement futur et/ou l'optimisation des paramètres de rendement d'une multitude de n (n étant très grand 1) d'opérateurs financiers locaux, ce procédé étant caractérisé en ce que u pour évaluer la valeur (Iujk) du portefeuille (Pu) d'un opérateur financier (Su) à la date future (tk) selon le modèle (Du) et dans l'état du marché simulé (Xj(tk)) correspondant au scénario (j), on prend en compte dans le modèle (Du) les états de marché simulés dans le même scénario (j), aux dates ultérieures (tl =k+1...m>), selon le principe suivant on définit pour chaque actif financier (9h(h-1-q) une quantité maximale (vh) qu'il est possible de traiter (à l'achat ou à la vente) dans une journée, 1 si aucune des quantités (Â"h) d'actifs (6h) détenus dans les portefeuilles (P") n'excède en valeur absolue la valeur limite (vh), alors on applique la formule (***) pour évaluer, à la date J et dans la simulation (j) la valeur (I",,) du portefeuille (P") si une des quantités (A"h) d'actifs (9h) détenus dans un des portefeuilles (P") excède en valeur absolue la valeur limite (vh), alors on remplace la quantité (;"h) par (vh), affectée de son signe, à la date (tk) et on reporte le reliquat (@."h vh) sur la date suivante (tk+,) en utilisant le vecteur de variables de marché (XJ(tk+,)) ; si cette quantité était encore en valeur absolue supérieure à la limite (vh), on la fixe à la valeur maximale (vh), avec le même signe, et on reporte à nouveau le reliquat (Auh 2vh) sur la date suivante (tk+2) en utilisant le vecteur de variables de marché (Xj(tk+2)), et ainsi de suite... ; ceci, notamment pour tenir compte du risque de liquidité imposant des contraintes de borne supérieure à la quantité de chaque actif (9h) que l'on est en droit de valoriser selon le modèle local (D") à une date donnée (tk), et entraînant un report éventuel de la position virtuelle résultante dans chaque simulation sur des dates ultérieures.