Bureau-du-Forestier-en-chef / MoSiR

Python package for the carbon fate of wood products
3 stars 0 forks source link

Mathématique de dégradation #2

Closed Landry-G closed 3 months ago

Landry-G commented 10 months ago

Implémenter l'option de pouvoir choisir l'équation mathématique de dégradation pour les nodes concernées :

Mettre d'une part le choix entre exponentielle et Chi-square. La valeur pourrait être flaggée dans la valeur du temps de demi-vie dans le json, par exemple : Dépotoir [c35] = 35 ans, avec l'équation Chi-square

Par la suite, vérifier la possibilité de transformer un string en équation mathématique pour donner l'option à l'utilisateur d'écrire sa propre formule : https://stackoverflow.com/questions/30775453/converting-a-string-into-equation-and-resolve-it

BFECLucasMoreau commented 10 months ago

Articles sur les différentes courbes de dégradation :

1 - Marland_2010 (Pose la question initiale) Gregg Marland, Eric Marland & Kirk Stellar (2010) " A Distributed Approach to Accounting to Carbon in Wood Products" Mitigation and Adaptation Strategies for Global Change Volume 15 issue 1 pp. 71-91 Version of Record Available From (www.springer.com)

2- Cherubini_2012 (Utilisation Chi-square) GCB Bioenergy (2012) 4, 784–798, https://doi.org/10.1111/j.1757-1707.2011.01156.x

3- Bates_2017 (Déterminer les valeurs des variables d'une distribution Gamma) => Exponentielle + Gamma standard + K=2 + Chi-square http://dx.doi.org/10.1016/j.jfe.2017.02.004

BFECLucasMoreau commented 10 months ago

Utilisation d'une équation Gamma unique basée sur Bates_2017 ==> L'utilisateur ne doit donc pas rentrer les valeurs de sa courbe. Besoin d'un package python qui permet d'utiliser la distribution Gamma. Les variables d'une distribution Gamma sont k ("shape ofthe function") et ϴ ("scales the graph horizontally and vertically")

L'équation va varier selon la valeurs des variables déterminée préalablement et manuellement grâce au fichier Excel R:\07_Outil_moyen_methode\02_Innovation\01_Carbone\ProduitsDuBois\Courbes_degradation\Calcul_Degradation

L'enjeu est de faire correspondre le stock restant à 50% de l'initial à l'année voulue selon la forme de l'équation choisie (Bates_2017 = Gamma standard, Exponentielle, Chi-square, K=2).

La question de flagger quelle équation correspond à quelle node reste identique !

BFECLucasMoreau commented 10 months ago

Capture_ExpDegradation

BFECLucasMoreau commented 10 months ago

Package Python possible : https://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.gamma.html

BFECLucasMoreau commented 10 months ago

Nomenclature dans les nodes Miro à changer :

Remplacer Half-life[35] par gamma(k=X,theta=Y,halflife=35) dans toutes les nodes. (garder le halflife = Z, dans le cas ou on irait plus tard vers un fit de courbe avec une inconnue - Voir issue #9 ) Plus besoin de flagger les nodes, chaque node aura son paramétrage à elle!

L'utilisateur doit déterminer ses variables k et theta avec le fichier Excel (voir réponse un peu plus haut)