@open3cl/engine 1.8.42 → 1.8.44

package/13.2_generateur_combustion.js

@@ -45,6 +45,31 @@ function evaluateFormula(formulaOrValue, pn, E, F) {
   return evaluateTvFormula(formulaOrValue, pn, E, F);
 }
 
+/**
+ * Puissance nominale maximale (W) retenue dans la formule Qp0 = Pn × (E + F × log Pn) / 100 des
+ * chaudières gaz standard et basse température installées après 2015.
+ */
+export const PN_MAX_QP0_CHAUDIERE_GAZ = 400000;
+
+/**
+ * Puissance (W) à utiliser dans la formule qp0 de la ligne de table.
+ *
+ * Pour les chaudières gaz après 2015 (qp0 = Pn × (E + F × log Pn) / 100), le moteur de référence
+ * (Tribu, `Calcul_generateur.cs` Rpn_chaudieregaz, `Math.Min(400.0, Pn)`) plafonne Pn à 400 kW :
+ * au-delà, la formule n'est plus représentative (elle décroît et devient négative vers 1 500 kW).
+ * Le plafond ne s'applique pas aux chaudières fioul (Rpn_chaudierefioul : Pn non plafonné).
+ * @see : Methode_de_calcul_3CL_DPE_2021-338.pdf - §13.2.2
+ *
+ * @param row {object} ligne de la table generateur_combustion
+ * @param pnCollectif {number} puissance nominale du générateur (W)
+ * @returns {number}
+ */
+export function puissanceFormuleQp0(row, pnCollectif) {
+  const formuleEF = /E\s*\+\s*F/.test(row.qp0_perc || '');
+  const chaudiereGaz = /chaudi[eè]re gaz/i.test(row.type_generateur || '');
+  return formuleEF && chaudiereGaz ? Math.min(pnCollectif, PN_MAX_QP0_CHAUDIERE_GAZ) : pnCollectif;
+}
+
 /**
  * Caractéristiques forfaitaires (rpn, rpint, qp0) d'un générateur à combustion.
  *
@@ -78,7 +103,7 @@ export function caracteristiquesGenerateurCombustion(row, pn, ratio, E, F, type,
 
   if (![4, 5].includes(methodeSaisie)) {
     if (row.qp0_perc) {
-      const qp0_calc = evaluateFormula(row.qp0_perc, pn / ratio, E, F);
+      const qp0_calc = evaluateFormula(row.qp0_perc, puissanceFormuleQp0(row, pn / ratio), E, F);
       // Certaines chaudières ont un qp0 en % de pn, d'autres ont des valeurs constantes
       if (row.qp0_perc.includes('Pn')) {
         caracteristiques.qp0 = qp0_calc * 1000 * ratio;

package/13.2_generateur_combustion.spec.js

@@ -64,7 +64,9 @@ const {
   tv_generateur_combustion,
   findGenerateurChMixteJumeau,
   updateGenerateurCombustion,
-  caracteristiquesGenerateurCombustion
+  caracteristiquesGenerateurCombustion,
+  puissanceFormuleQp0,
+  PN_MAX_QP0_CHAUDIERE_GAZ
 } = await import('./13.2_generateur_combustion.js');
 const { tv, tvColumnLines } = await import('./utils.js');
 const { conventionPuissanceNominale } = await import('./13.2_generateur_combustion_pn.js');
@@ -1099,3 +1101,72 @@ describe('caracteristiquesGenerateurCombustion - selon la méthode de saisie', (
     ).toBeCloseTo(50, 9);
   });
 });
+
+/**
+ * Plafond de Pn à 400 kW dans Qp0 = Pn × (E + F × log Pn) / 100 des chaudières gaz après 2015
+ * (moteur de référence Tribu : Calcul_generateur.cs, Rpn_chaudieregaz, Math.Min(400.0, Pn)).
+ * @see : Methode_de_calcul_3CL_DPE_2021-338.pdf - §13.2.2
+ */
+describe('puissanceFormuleQp0 - plafond de Pn à 400 kW pour les chaudières gaz', () => {
+  const ROW_GAZ = {
+    type_generateur: 'Chaudière gaz basse température après 2015',
+    qp0_perc: 'Pn*(E+F*log10(Pn))/100'
+  };
+  const ROW_FIOUL = {
+    type_generateur: 'Chaudière fioul standard après 2015',
+    qp0_perc: 'Pn*(E+F*log10(Pn))/100'
+  };
+
+  test('chaudière gaz de 500 kW : puissance plafonnée à 400 kW', () => {
+    expect(PN_MAX_QP0_CHAUDIERE_GAZ).toBe(400000);
+    expect(puissanceFormuleQp0(ROW_GAZ, 500000)).toBe(400000);
+  });
+
+  test('chaudière gaz de 400 kW ou moins : puissance inchangée', () => {
+    expect(puissanceFormuleQp0(ROW_GAZ, 400000)).toBe(400000);
+    expect(puissanceFormuleQp0(ROW_GAZ, 24000)).toBe(24000);
+  });
+
+  test('chaudière fioul : pas de plafond', () => {
+    expect(puissanceFormuleQp0(ROW_FIOUL, 500000)).toBe(500000);
+  });
+
+  test('chaudière gaz avec qp0 forfaitaire (sans E, F) : pas de plafond', () => {
+    expect(
+      puissanceFormuleQp0(
+        { type_generateur: 'Chaudière gaz standard 2001-2015', qp0_perc: '1%' },
+        5e5
+      )
+    ).toBe(500000);
+  });
+
+  test('ligne sans libellé ni formule qp0 : pas de plafond', () => {
+    expect(puissanceFormuleQp0({}, 500000)).toBe(500000);
+  });
+
+  test('qp0 d’une chaudière gaz BT de 500 kW avec ventilateur (E = 1.75, F = -0.55) : Pn = 400 kW', () => {
+    const { qp0 } = caracteristiquesGenerateurCombustion(ROW_GAZ, 500000, 1, 1.75, -0.55, 'ch', 1);
+    // 400 × (1.75 − 0.55 × log10(400)) / 100 kW = Qp0 Tribu IC4-0-21 -- valeur de référence de régression
+    expect(qp0).toBeCloseTo(1275.468019078, 6);
+  });
+
+  test('qp0 d’une chaudière fioul de 500 kW : formule appliquée sur Pn = 500 kW', () => {
+    const { qp0 } = caracteristiquesGenerateurCombustion(
+      ROW_FIOUL,
+      500000,
+      1,
+      1.75,
+      -0.55,
+      'ch',
+      1
+    );
+    // 500 × (1.75 − 0.55 × log10(500)) / 100 kW -- valeur de référence de régression
+    expect(qp0).toBeCloseTo(500 * (1.75 - 0.55 * Math.log10(500)) * 10, 6);
+  });
+
+  test('chaudière gaz virtualisée : plafond appliqué sur Pn collectif puis qp0 ramené au logement', () => {
+    const { qp0 } = caracteristiquesGenerateurCombustion(ROW_GAZ, 50000, 0.1, 1.75, -0.55, 'ch', 1);
+    // Pn collectif = 500 kW plafonné à 400 kW, puis × ratio 0.1
+    expect(qp0).toBeCloseTo(127.5468019078, 6);
+  });
+});

package/16.2_production_enr.js

@@ -62,8 +62,18 @@ export class ProductionENR {
    * @param zc_id {string}
    * @param coeff_ep_override {number?}
    * @param ratioPv {number?} proratisation au logement d'une installation d'immeuble (cf. ratioProrataPvLogement)
+   * @param facteursAc {object?} si fourni, reçoit par poste la part non autoconsommée (cf. getFacteursNonAutoconsommes)
    */
-  calculateEnr(productionElecEnr, conso, Sh, th, zc_id, coeff_ep_override, ratioPv = 1) {
+  calculateEnr(
+    productionElecEnr,
+    conso,
+    Sh,
+    th,
+    zc_id,
+    coeff_ep_override,
+    ratioPv = 1,
+    facteursAc = null
+  ) {
     const productionElectricite = {
       conso_elec_ac: 0,
       production_pv: 0,
@@ -92,6 +102,10 @@ export class ProductionENR {
         ratioPv
       );
 
+      if (facteursAc) {
+        Object.assign(facteursAc, this.getFacteursNonAutoconsommes(productionElectricite, conso));
+      }
+
       // Mise à jour des consommations d'énergie finale en minorant l'énergie consommée par l'énergie autoconsommée par le poste
       this.updateEfConso(productionElectricite, conso, Sh);
 
@@ -213,6 +227,55 @@ export class ProductionENR {
     productionElectricite.conso_elec_ac_auxiliaire = consoAcAuxiliaires;
   }
 
+  /**
+   * Part non autoconsommée (1 − Celec_ac_poste / Celec_poste) de la consommation électrique de chaque
+   * poste, à appliquer aux consommations électriques avant calcul des émissions GES et du coût.
+   *
+   * Le moteur de référence calcule les émissions et les coûts sur les consommations électriques
+   * minorées de l'autoconsommation PV (Tribu `Calcul_batiment.cs` l.1188 carbone_total :
+   * (Cch_electricite − Celec_ac_ch) × 0,079 + … ; l.1318 calc_cout.Calcul(…, num19 − num60, …)).
+   * Doit être appelée avant updateEfConso (consommations encore brutes).
+   *
+   * @param productionElectricite autoconsommation par poste (calculateConsoElecAc)
+   * @param conso consommations brutes (ef_conso, sortie_par_energie_collection)
+   * @returns {Object<string, number>} facteur par poste : ch, ecs, fr, eclairage, ventilation,
+   * auxiliaire_generation_ch, auxiliaire_generation_ecs, auxiliaire_distribution_ch,
+   * auxiliaire_distribution_ecs
+   */
+  getFacteursNonAutoconsommes(productionElectricite, conso) {
+    const consoElec = conso.sortie_par_energie_collection.sortie_par_energie.find(
+      (sortie) => sortie.enum_type_energie_id === '1'
+    );
+    const ef = conso.ef_conso;
+    const facteur = (ac, brut) => (brut > 0 ? 1 - ac / brut : 1);
+    return {
+      ch: facteur(productionElectricite.conso_elec_ac_ch, consoElec?.conso_ch),
+      ecs: facteur(productionElectricite.conso_elec_ac_ecs, consoElec?.conso_ecs),
+      fr: facteur(productionElectricite.conso_elec_ac_fr, ef.conso_fr),
+      eclairage: facteur(productionElectricite.conso_elec_ac_eclairage, ef.conso_eclairage),
+      ventilation: facteur(
+        productionElectricite.conso_elec_ac_ventilation,
+        ef.conso_auxiliaire_ventilation
+      ),
+      auxiliaire_generation_ch: facteur(
+        productionElectricite.conso_elec_ac_auxiliaire_generation_ch,
+        ef.conso_auxiliaire_generation_ch
+      ),
+      auxiliaire_generation_ecs: facteur(
+        productionElectricite.conso_elec_ac_auxiliaire_generation_ecs,
+        ef.conso_auxiliaire_generation_ecs
+      ),
+      auxiliaire_distribution_ch: facteur(
+        productionElectricite.conso_elec_ac_auxiliaire_distribution_ch,
+        ef.conso_auxiliaire_distribution_ch
+      ),
+      auxiliaire_distribution_ecs: facteur(
+        productionElectricite.conso_elec_ac_auxiliaire_distribution_ecs,
+        ef.conso_auxiliaire_distribution_ecs
+      )
+    };
+  }
+
   /**
    * Mise à jour des consommations ef en minorant l'énergie finale consommée par l'énergie autoconsommée par chaque enveloppe
    * @param productionElectricite

package/16.2_production_enr.spec.js

@@ -72,6 +72,108 @@ describe('production ENR unit tests', () => {
     }
   );
 
+  test('getFacteursNonAutoconsommes : part non autoconsommée de chaque poste électrique', () => {
+    const productionElectricite = {
+      conso_elec_ac_ch: 250,
+      conso_elec_ac_ecs: 0,
+      conso_elec_ac_fr: 0,
+      conso_elec_ac_eclairage: 50,
+      conso_elec_ac_ventilation: 30,
+      conso_elec_ac_auxiliaire_generation_ch: 0,
+      conso_elec_ac_auxiliaire_generation_ecs: 0,
+      conso_elec_ac_auxiliaire_distribution_ch: 10,
+      conso_elec_ac_auxiliaire_distribution_ecs: 0
+    };
+    const conso = {
+      ef_conso: {
+        conso_fr: 0,
+        conso_eclairage: 100,
+        conso_auxiliaire_ventilation: 120,
+        conso_auxiliaire_generation_ch: 0,
+        conso_auxiliaire_generation_ecs: 0,
+        conso_auxiliaire_distribution_ch: 40,
+        conso_auxiliaire_distribution_ecs: 0
+      },
+      sortie_par_energie_collection: {
+        sortie_par_energie: [
+          { enum_type_energie_id: '2', conso_ch: 5000, conso_ecs: 2000 },
+          { enum_type_energie_id: '1', conso_ch: 1000, conso_ecs: 0 }
+        ]
+      }
+    };
+    expect(productionENR.getFacteursNonAutoconsommes(productionElectricite, conso)).toStrictEqual({
+      ch: 0.75,
+      ecs: 1,
+      fr: 1,
+      eclairage: 0.5,
+      ventilation: 0.75,
+      auxiliaire_generation_ch: 1,
+      auxiliaire_generation_ecs: 1,
+      auxiliaire_distribution_ch: 0.75,
+      auxiliaire_distribution_ecs: 1
+    });
+  });
+
+  test('calculateEnr renseigne les facteurs avant minoration des consommations', () => {
+    const conso = {
+      ef_conso: {
+        conso_ch: 1000,
+        conso_ecs: 0,
+        conso_fr: 0,
+        conso_eclairage: 100,
+        conso_auxiliaire_ventilation: 0,
+        conso_auxiliaire_generation_ch: 0,
+        conso_auxiliaire_generation_ecs: 0,
+        conso_auxiliaire_distribution_ch: 0,
+        conso_auxiliaire_distribution_ecs: 0,
+        conso_totale_auxiliaire: 0,
+        conso_5_usages: 1100,
+        conso_5_usages_m2: 11
+      },
+      ep_conso: {
+        ep_conso_ch: 1900,
+        ep_conso_ecs: 0,
+        ep_conso_fr: 0,
+        ep_conso_eclairage: 190,
+        ep_conso_totale_auxiliaire: 0,
+        ep_conso_5_usages: 2090,
+        ep_conso_5_usages_m2: 20.9
+      },
+      sortie_par_energie_collection: {
+        sortie_par_energie: [
+          { enum_type_energie_id: '1', conso_ch: 1000, conso_ecs: 0, conso_5_usages: 1100 }
+        ]
+      }
+    };
+    const facteursAc = {};
+    productionENR.calculateEnr(
+      {
+        donnee_entree: { presence_production_pv: 1 },
+        panneaux_pv_collection: {
+          panneaux_pv: [
+            {
+              enum_orientation_pv_id: '3',
+              enum_inclinaison_pv_id: '2',
+              surface_totale_capteurs: 10
+            }
+          ]
+        }
+      },
+      conso,
+      100,
+      'maison',
+      1,
+      undefined,
+      1,
+      facteursAc
+    );
+    expect(facteursAc.ch).toBeGreaterThan(0);
+    expect(facteursAc.ch).toBeLessThan(1);
+    expect(facteursAc.eclairage).toBeGreaterThanOrEqual(0);
+    expect(facteursAc.eclairage).toBeLessThan(1);
+    expect(facteursAc.ecs).toBe(1);
+  });
+
   describe('proratisation au logement des installations PV d’immeuble', () => {
     const panneau = (extra = {}) => ({
       panneaux_pv_collection: {

package/README.md

@@ -370,26 +370,26 @@ Ces résultats sont **générés automatiquement** à la fin de `npm run test:co
 
 <!-- CORPUS:START -->
 
-> **Version `1.8.41`** · branche `main` · généré le 2026-10-10
+> **Version `1.8.43`** · branche `main` · généré le 2026-10-10
 > Seuil de tolérance **5%**
 
 <table>
 <tr>
 <td align="center"><strong>9</strong><br/><sub>corpus</sub></td>
 <td align="center"><strong>89 984</strong><br/><sub>DPE analysés</sub></td>
-<td align="center"><strong>61 956</strong><br/><sub>DPE conformes</sub></td>
-<td align="center"><strong>68,85 %</strong><br/><sub>réussite globale</sub></td>
+<td align="center"><strong>61 969</strong><br/><sub>DPE conformes</sub></td>
+<td align="center"><strong>68,87 %</strong><br/><sub>réussite globale</sub></td>
 </tr>
 </table>
 
 |     | Corpus                                                                                                                  |    Réussite |                        |  DPE conformes |
 | :-: | :---------------------------------------------------------------------------------------------------------------------- | ----------: | :--------------------- | -------------: |
-| 🔴  | **Généraliste**<br/><sub>`corpus_dpe.csv`</sub>                                                                         | **46,18 %** | `█████████░░░░░░░░░░░` |  4 618 / 9 999 |
-| 🟢  | **Appartement · chauffage individuel (2025)**<br/><sub>`dpe_appartement_individuel_chauffage_individuel_2025.csv`</sub> | **92,54 %** | `███████████████████░` | 9 254 / 10 000 |
+| 🔴  | **Généraliste**<br/><sub>`corpus_dpe.csv`</sub>                                                                         | **46,24 %** | `█████████░░░░░░░░░░░` |  4 624 / 9 999 |
+| 🟢  | **Appartement · chauffage individuel (2025)**<br/><sub>`dpe_appartement_individuel_chauffage_individuel_2025.csv`</sub> | **92,55 %** | `███████████████████░` | 9 255 / 10 000 |
 | 🟢  | **Logement individuel (2025)**<br/><sub>`dpe_logement_individuel_2025.csv`</sub>                                        | **89,55 %** | `██████████████████░░` |  8 948 / 9 992 |
-| 🟢  | **Maison individuelle (2025)**<br/><sub>`dpe_maison_individuelle_2025.csv`</sub>                                        | **88,61 %** | `██████████████████░░` |  8 860 / 9 999 |
+| 🟢  | **Maison individuelle (2025)**<br/><sub>`dpe_maison_individuelle_2025.csv`</sub>                                        | **88,66 %** | `██████████████████░░` |  8 865 / 9 999 |
 | 🟡  | **Immeuble · chauffage individuel**<br/><sub>`dpe_immeuble_chauffage_individuel.csv`</sub>                              | **84,22 %** | `█████████████████░░░` |  8 420 / 9 998 |
-| 🟡  | **Appartement · chauffage collectif (2025)**<br/><sub>`dpe_appartement_individuel_chauffage_collectif_2025.csv`</sub>   | **71,60 %** | `██████████████░░░░░░` | 7 160 / 10 000 |
+| 🟡  | **Appartement · chauffage collectif (2025)**<br/><sub>`dpe_appartement_individuel_chauffage_collectif_2025.csv`</sub>   | **71,61 %** | `██████████████░░░░░░` | 7 161 / 10 000 |
 | 🟡  | **Immeuble · chauffage collectif**<br/><sub>`dpe_immeuble_chauffage_collectif.csv`</sub>                                | **66,19 %** | `█████████████░░░░░░░` |  6 618 / 9 999 |
 | 🔴  | **Immeuble · chauffage mixte**<br/><sub>`dpe_immeuble_chauffage_mixte.csv`</sub>                                        | **53,74 %** | `███████████░░░░░░░░░` |  5 372 / 9 997 |
 | 🔴  | **Individuel généré depuis l'immeuble (2026)**<br/><sub>`dpe_individuel_a_partir_dpe_immeuble_2026.csv`</sub>           | **27,06 %** | `█████░░░░░░░░░░░░░░░` | 2 706 / 10 000 |
@@ -403,24 +403,24 @@ Ces résultats sont **générés automatiquement** à la fin de `npm run test:co
 
 <table>
 <tr>
-<td align="center"><strong>5,03 ms</strong><br/><sub>moyenne</sub></td>
-<td align="center"><strong>4,37 ms</strong><br/><sub>médiane</sub></td>
+<td align="center"><strong>4,86 ms</strong><br/><sub>moyenne</sub></td>
+<td align="center"><strong>4,24 ms</strong><br/><sub>médiane</sub></td>
 <td align="center"><strong>0,12 ms</strong><br/><sub>min</sub></td>
-<td align="center"><strong>579,2 ms</strong><br/><sub>max</sub></td>
+<td align="center"><strong>549,1 ms</strong><br/><sub>max</sub></td>
 </tr>
 </table>
 
 | Corpus                                         | Moyenne | Médiane |     Min |      Max |     p95 |     p99 |
 | :--------------------------------------------- | ------: | ------: | ------: | -------: | ------: | ------: |
-| **Généraliste**                                | 4,29 ms | 3,83 ms | 0,23 ms |  32,6 ms | 7,90 ms | 10,9 ms |
-| **Appartement · chauffage individuel (2025)**  | 3,69 ms | 3,37 ms | 1,32 ms |  28,5 ms | 6,33 ms | 8,22 ms |
-| **Logement individuel (2025)**                 | 5,06 ms | 4,44 ms | 1,56 ms |  66,4 ms | 9,43 ms | 13,7 ms |
-| **Maison individuelle (2025)**                 | 5,88 ms | 5,24 ms | 1,90 ms |  49,5 ms | 10,6 ms | 14,2 ms |
-| **Immeuble · chauffage individuel**            | 6,22 ms | 5,27 ms | 0,12 ms |  84,0 ms | 12,6 ms | 19,7 ms |
-| **Appartement · chauffage collectif (2025)**   | 4,18 ms | 3,77 ms | 1,37 ms |  36,4 ms | 6,96 ms | 8,88 ms |
-| **Immeuble · chauffage collectif**             | 5,43 ms | 4,57 ms | 0,14 ms | 114,1 ms | 9,98 ms | 17,3 ms |
-| **Immeuble · chauffage mixte**                 | 4,88 ms | 3,99 ms | 0,20 ms | 579,2 ms | 9,79 ms | 16,2 ms |
-| **Individuel généré depuis l'immeuble (2026)** | 5,64 ms | 4,81 ms | 1,59 ms |  57,5 ms | 10,6 ms | 16,2 ms |
+| **Généraliste**                                | 4,32 ms | 3,87 ms | 0,21 ms |  28,9 ms | 8,02 ms | 10,8 ms |
+| **Appartement · chauffage individuel (2025)**  | 3,79 ms | 3,46 ms | 1,50 ms |  34,3 ms | 6,55 ms | 8,57 ms |
+| **Logement individuel (2025)**                 | 4,75 ms | 4,21 ms | 1,61 ms |  34,4 ms | 8,90 ms | 12,2 ms |
+| **Maison individuelle (2025)**                 | 5,40 ms | 4,87 ms | 1,76 ms |  36,9 ms | 9,53 ms | 12,7 ms |
+| **Immeuble · chauffage individuel**            | 6,22 ms | 5,32 ms | 0,14 ms | 115,4 ms | 12,5 ms | 19,3 ms |
+| **Appartement · chauffage collectif (2025)**   | 4,08 ms | 3,65 ms | 1,48 ms |  34,4 ms | 6,76 ms | 8,58 ms |
+| **Immeuble · chauffage collectif**             | 5,30 ms | 4,49 ms | 0,12 ms | 107,4 ms | 9,83 ms | 17,6 ms |
+| **Immeuble · chauffage mixte**                 | 4,78 ms | 3,93 ms | 0,22 ms | 549,1 ms | 9,87 ms | 15,9 ms |
+| **Individuel généré depuis l'immeuble (2026)** | 5,07 ms | 4,39 ms | 1,65 ms |  42,9 ms | 9,57 ms | 14,4 ms |
 
 <sub>La moyenne est tirée vers le haut par les DPE collectifs, dont le coût atteint plusieurs
 dizaines de fois la médiane : c’est la médiane qui décrit le cas courant, et p95/p99 le cas

package/conso.js

@@ -376,7 +376,8 @@ export default function calc_conso(
   dateDpe,
   coeffEp,
   nbLogements = 1,
-  ratioSurfaceCollectif = 1
+  ratioSurfaceCollectif = 1,
+  facteursAc = {}
 ) {
   const gen_ch = ch.reduce((acc, ch) => {
     const generateur_chauffage = ch.generateur_chauffage_collection.generateur_chauffage;
@@ -399,9 +400,6 @@ export default function calc_conso(
 
   ecs = Array.isArray(ecs) ? ecs : [];
 
-  const conso_aux_distribution_ecs_total = ecs.reduce((acc, ecs) => {
-    return acc + (ecs.donnee_intermediaire.conso_auxiliaire_distribution_ecs || 0);
-  }, 0);
   const gen_ecs = ecs.reduce((acc, ecs) => {
     const generateur_ecs = ecs.generateur_ecs_collection.generateur_ecs;
     if (prorataECS === 1) {
@@ -416,6 +414,14 @@ export default function calc_conso(
     return acc.concat(generateur_ecs);
   }, []);
 
+  /**
+   * Coût : le moteur de référence ne minore de l'autoconsommation PV que le chauffage, l'ECS et les
+   * auxiliaires de génération / distribution ; l'éclairage, le froid et la ventilation sont tarifés
+   * bruts (Tribu `Calcul_batiment.cs` l.1318 : num31 = Ceclairage1, conso_elec_*_fr,
+   * conso_elec_*_vent passés sans déduction à calc_cout.Calcul).
+   */
+  const facteursAcCout = { ...facteursAc, eclairage: 1, fr: 1, ventilation: 1 };
+
   /**
    * Barème des coûts : barème 2021 (par poste) avant le 01/07/2024 ; ensuite barème mis à jour
    * avec tranche déterminée sur la consommation totale de l'énergie (voir coefCoutProrata).
@@ -456,7 +462,8 @@ export default function calc_conso(
         masque,
         prorataECS,
         prorataChauffage,
-        ecs
+        ecs,
+        facteursAcCout
       );
       const coll = calc_conso_pond(
         0,
@@ -469,7 +476,8 @@ export default function calc_conso(
         masque,
         prorataECS,
         prorataChauffage,
-        []
+        [],
+        facteursAcCout
       );
       // Appoint individuel des convecteurs bi-jonction : tarifé en individuel (voir coutBijonction)
       const appointIndividuel = gen_ch
@@ -502,8 +510,7 @@ export default function calc_conso(
       null,
       prorataECS,
       prorataChauffage,
-      ecs,
-      conso_aux_distribution_ecs_total
+      ecs
     ),
     ep_conso: calc_conso_pond(
       Sh,
@@ -516,8 +523,7 @@ export default function calc_conso(
       coeffEp,
       prorataECS,
       prorataChauffage,
-      ecs,
-      conso_aux_distribution_ecs_total
+      ecs
     ),
     emission_ges: calc_conso_pond(
       Sh,
@@ -531,7 +537,7 @@ export default function calc_conso(
       prorataECS,
       prorataChauffage,
       ecs,
-      conso_aux_distribution_ecs_total
+      facteursAc
     ),
     cout: calc_conso_pond(
       Sh,
@@ -545,7 +551,7 @@ export default function calc_conso(
       prorataECS,
       prorataChauffage,
       ecs,
-      conso_aux_distribution_ecs_total
+      facteursAcCout
     )
   };
   ret.ep_conso.classe_bilan_dpe = classe_bilan_dpe(
@@ -670,12 +676,14 @@ export function classe_emission_ges(emission_ges_5_usages_m2, zc_id, ca_id, Sh)
  * @param coef {number}
  * @param prorataECS {number}
  * @param prefix {string}
+ * @param facteurAc {number} part non autoconsommée (PV) de l'ECS électrique (1 sans PV)
  * @returns {number}
  */
-function getEcsConso(gen_ecs, field, coef, prorataECS, prefix) {
+function getEcsConso(gen_ecs, field, coef, prorataECS, prefix, facteurAc = 1) {
   return gen_ecs.reduce((acc, gen_ecs) => {
-    const conso = gen_ecs.donnee_intermediaire[field];
     const typeEnergie = getTypeEnergie(gen_ecs.donnee_entree, 'ecs');
+    const conso =
+      gen_ecs.donnee_intermediaire[field] * (typeEnergie === 'électricité ecs' ? facteurAc : 1);
 
     let coeffConsoEcs = { ...coefPourGenerateur(coef, gen_ecs) };
 
@@ -701,13 +709,14 @@ function getEcsConso(gen_ecs, field, coef, prorataECS, prefix) {
  * @param field {string} conso_auxiliaire_generation_ecs[_depensier]
  * @param coef {object}
  * @param prorataECS {number}
+ * @param facteurAc {number} part non autoconsommée (PV) des auxiliaires de génération ECS
  * @returns {number}
  */
-export function getAuxGenerationEcs(ecs, field, coef, prorataECS) {
+export function getAuxGenerationEcs(ecs, field, coef, prorataECS, facteurAc = 1) {
   return ecs.reduce((acc, inst) => {
     const rdim = prorataECS === 1 ? Number(inst.donnee_entree.rdim) || 1 : 1;
     return inst.generateur_ecs_collection.generateur_ecs.reduce((acc2, gen) => {
-      const conso = (gen.donnee_intermediaire[field] || 0) * rdim;
+      const conso = (gen.donnee_intermediaire[field] || 0) * rdim * facteurAc;
       return acc2 + getConso(coef, 'électricité auxiliaire', conso);
     }, acc);
   }, 0);
@@ -720,12 +729,14 @@ export function getAuxGenerationEcs(ecs, field, coef, prorataECS) {
  * @param coef {number}
  * @param prorataChauffage {number}
  * @param prefix {string}
+ * @param facteurAc {number} part non autoconsommée (PV) du chauffage électrique (1 sans PV)
  * @returns {number}
  */
-function getChauffageConso(gen_ch, field, coef, prorataChauffage, prefix) {
+function getChauffageConso(gen_ch, field, coef, prorataChauffage, prefix, facteurAc = 1) {
   return gen_ch.reduce((acc, gen_ch) => {
-    const conso = gen_ch.donnee_intermediaire[field];
     const typeEnergie = getTypeEnergie(gen_ch.donnee_entree, 'ch');
+    const conso =
+      gen_ch.donnee_intermediaire[field] * (typeEnergie === 'électricité ch' ? facteurAc : 1);
 
     // La clé de répartition n'est utilisée que dans le cadre des chauffages collectifs
     const repartition =
@@ -757,11 +768,13 @@ function calc_conso_pond(
   coef,
   prorataECS,
   prorataChauffage,
-  ecs_installations
+  ecs_installations,
+  facteursAc = {}
 ) {
+  const fAc = (poste) => facteursAc[poste] ?? 1;
   const ret = {};
   ret.auxiliaire_ventilation = vt_list.reduce((acc, vt) => {
-    let conso = vt.donnee_intermediaire.conso_auxiliaire_ventilation || 0;
+    let conso = (vt.donnee_intermediaire.conso_auxiliaire_ventilation || 0) * fAc('ventilation');
 
     if (vt.donnee_entree.cle_repartition_ventilation) {
       conso *= vt.donnee_entree.cle_repartition_ventilation;
@@ -771,11 +784,13 @@ function calc_conso_pond(
   }, 0);
 
   const conso_eclairage = calc_conso_eclairage(zc_id) * Sh;
-  ret.eclairage = getConso(coef, 'électricité éclairage', conso_eclairage);
+  ret.eclairage = getConso(coef, 'électricité éclairage', conso_eclairage * fAc('eclairage'));
 
   // aux ch
   ret.auxiliaire_generation_ch = gen_ch.reduce((acc, gen_ch) => {
-    const conso = gen_ch.donnee_intermediaire.conso_auxiliaire_generation_ch || 0;
+    const conso =
+      (gen_ch.donnee_intermediaire.conso_auxiliaire_generation_ch || 0) *
+      fAc('auxiliaire_generation_ch');
     return acc + getConso(coef, 'électricité auxiliaire', conso);
   }, 0);
 
@@ -785,7 +800,9 @@ function calc_conso_pond(
   }, 0);
 
   ret.auxiliaire_distribution_ch = gen_ch.reduce((acc, gen_ch) => {
-    let conso = gen_ch.donnee_intermediaire.conso_auxiliaire_distribution_ch || 0;
+    let conso =
+      (gen_ch.donnee_intermediaire.conso_auxiliaire_distribution_ch || 0) *
+      fAc('auxiliaire_distribution_ch');
 
     /**
      * enum_methode_calcul_conso_id
@@ -801,7 +818,7 @@ function calc_conso_pond(
     return acc + getConso(coef, 'électricité auxiliaire', conso);
   }, 0);
 
-  ret.ch = getChauffageConso(gen_ch, 'conso_ch', coef, prorataChauffage, prefix);
+  ret.ch = getChauffageConso(gen_ch, 'conso_ch', coef, prorataChauffage, prefix, fAc('ch'));
 
   ret.ch_depensier = getChauffageConso(
     gen_ch,
@@ -815,7 +832,8 @@ function calc_conso_pond(
     ecs_installations,
     'conso_auxiliaire_generation_ecs',
     coef,
-    prorataECS
+    prorataECS,
+    fAc('auxiliaire_generation_ecs')
   );
 
   ret.auxiliaire_generation_ecs_depensier = getAuxGenerationEcs(
@@ -826,17 +844,20 @@ function calc_conso_pond(
   );
 
   ret.auxiliaire_distribution_ecs = ecs_installations.reduce((acc, inst) => {
-    const conso = inst.donnee_intermediaire.conso_auxiliaire_distribution_ecs || 0;
+    const conso =
+      (inst.donnee_intermediaire.conso_auxiliaire_distribution_ecs || 0) *
+      fAc('auxiliaire_distribution_ecs');
     return acc + getConso(coef, 'électricité auxiliaire', conso);
   }, 0);
 
-  ret.ecs = getEcsConso(gen_ecs, 'conso_ecs', coef, prorataECS, prefix);
+  ret.ecs = getEcsConso(gen_ecs, 'conso_ecs', coef, prorataECS, prefix, fAc('ecs'));
 
   ret.ecs_depensier = getEcsConso(gen_ecs, 'conso_ecs_depensier', coef, prorataECS, prefix);
 
   ret.fr = fr_list.reduce((acc, fr) => {
-    const conso = fr.donnee_intermediaire.conso_fr;
     const typeEnergie = getTypeEnergie(fr.donnee_entree, 'fr');
+    const conso =
+      fr.donnee_intermediaire.conso_fr * (typeEnergie === 'électricité fr' ? fAc('fr') : 1);
 
     return acc + getConso(coef, typeEnergie, conso);
   }, 0);

package/conso.spec.js

@@ -238,6 +238,121 @@ describe('calc_conso - agrégation des consommations', () => {
     expect(res.ef_conso.conso_eclairage).toBe(100);
   });
 
+  test('autoconsommation PV : GES de tous les postes et coût hors éclairage/froid/ventilation minorés (Tribu Calcul_batiment.cs l.1188 / l.1318)', () => {
+    const ch = [
+      installCh({ cle_repartition_ch: 1 }, [
+        genCh('1', { conso_ch: 1000, conso_ch_depensier: 1200 })
+      ])
+    ];
+    // 25 % du chauffage et 50 % de l'éclairage autoconsommés
+    const facteursAc = { ch: 0.75, eclairage: 0.5 };
+
+    const brut = calc_conso(100, 1, 1, [], ch, [], [], 1, 1, DATE_DPE, coef_ep);
+    const net = calc_conso(100, 1, 1, [], ch, [], [], 1, 1, DATE_DPE, coef_ep, 1, 1, facteursAc);
+
+    // EF / EP inchangées (la minoration EF/EP reste faite par ProductionENR.updateEfConso)
+    expect(net.ef_conso).toStrictEqual(brut.ef_conso);
+    expect(net.ep_conso.ep_conso_5_usages).toBeCloseTo(brut.ep_conso.ep_conso_5_usages, 9);
+    // GES : 750 kWh × 0,079 ; éclairage 100 × 0,5 × coef
+    expect(net.emission_ges.emission_ges_ch).toBeCloseTo(750 * 0.079, 9);
+    expect(net.emission_ges.emission_ges_eclairage).toBeCloseTo(
+      brut.emission_ges.emission_ges_eclairage * 0.5,
+      9
+    );
+    // Coût : chauffage tarifé sur 750 kWh au lieu de 1000, éclairage tarifé brut
+    expect(net.cout.cout_ch).toBeLessThan(brut.cout.cout_ch);
+    expect(net.cout.cout_ch).toBeCloseTo(
+      calc_conso(
+        100,
+        1,
+        1,
+        [],
+        [
+          installCh({ cle_repartition_ch: 1 }, [
+            genCh('1', { conso_ch: 750, conso_ch_depensier: 900 })
+          ])
+        ],
+        [],
+        [],
+        1,
+        1,
+        DATE_DPE,
+        coef_ep
+      ).cout.cout_ch,
+      9
+    );
+    expect(net.cout.cout_eclairage).toBeCloseTo(brut.cout.cout_eclairage, 9);
+  });
+
+  test('autoconsommation PV : froid électrique minoré, froid non électrique non minoré', () => {
+    const frElec = [
+      {
+        donnee_entree: { enum_type_energie_id: '1' },
+        donnee_intermediaire: { conso_fr: 200, conso_fr_depensier: 250 }
+      }
+    ];
+    const frGaz = [
+      {
+        donnee_entree: { enum_type_energie_id: '2' },
+        donnee_intermediaire: { conso_fr: 200, conso_fr_depensier: 250 }
+      }
+    ];
+    const facteursAc = { fr: 0.5 };
+
+    const brutElec = calc_conso(100, 1, 1, [], [], [], frElec, 1, 1, DATE_DPE, coef_ep);
+    const netElec = calc_conso(
+      100,
+      1,
+      1,
+      [],
+      [],
+      [],
+      frElec,
+      1,
+      1,
+      DATE_DPE,
+      coef_ep,
+      1,
+      1,
+      facteursAc
+    );
+    expect(netElec.emission_ges.emission_ges_fr).toBeCloseTo(
+      brutElec.emission_ges.emission_ges_fr * 0.5,
+      9
+    );
+
+    const brutGaz = calc_conso(100, 1, 1, [], [], [], frGaz, 1, 1, DATE_DPE, coef_ep);
+    const netGaz = calc_conso(
+      100,
+      1,
+      1,
+      [],
+      [],
+      [],
+      frGaz,
+      1,
+      1,
+      DATE_DPE,
+      coef_ep,
+      1,
+      1,
+      facteursAc
+    );
+    expect(netGaz.emission_ges.emission_ges_fr).toBeGreaterThan(0);
+    expect(netGaz.emission_ges.emission_ges_fr).toBe(brutGaz.emission_ges.emission_ges_fr);
+  });
+
+  test('sans facteurs d’autoconsommation : résultats identiques', () => {
+    const ch = [
+      installCh({ cle_repartition_ch: 1 }, [
+        genCh('1', { conso_ch: 1000, conso_ch_depensier: 1200 })
+      ])
+    ];
+    const brut = calc_conso(100, 1, 1, [], ch, [], [], 1, 1, DATE_DPE, coef_ep);
+    const net = calc_conso(100, 1, 1, [], ch, [], [], 1, 1, DATE_DPE, coef_ep, 1, 1, {});
+    expect(net).toStrictEqual(brut);
+  });
+
   test('installation collective (méthode 2) : clé de répartition appliquée à conso et distribution', () => {
     const ch = [
       installCh(

package/engine.js

@@ -33,7 +33,7 @@ import getFicheTechnique from './ficheTechnique.js';
 import { ProductionENR, ratioProrataPvLogement } from './16.2_production_enr.js';
 import DpeSanitizerService from './dpe-sanitizer.service.js';
 
-const LIB_VERSION = '1.8.42';
+const LIB_VERSION = '1.8.44';
 
 const dpeSanitizerService = new DpeSanitizerService();
 
@@ -617,6 +617,7 @@ export function calcul_3cl(inputDpe, options) {
     cg.surface_habitable_immeuble
   );
 
+  const facteursAc = {};
   const production_electricite = productionENR.calculateEnr(
     dpe.logement.production_elec_enr,
     conso,
@@ -624,9 +625,35 @@ export function calcul_3cl(inputDpe, options) {
     th,
     zc_id,
     undefined,
-    ratioPv
+    ratioPv,
+    facteursAc
   );
 
+  /**
+   * Émissions GES et coût calculés sur les consommations électriques minorées de l'autoconsommation
+   * PV (Tribu `Calcul_batiment.cs` l.1188 carbone_total, l.1318 calc_cout.Calcul).
+   */
+  if (Object.keys(facteursAc).length) {
+    const consoNette = calc_conso(
+      Sh,
+      zc_id,
+      ca_id,
+      vt_list,
+      instal_ch,
+      ecs,
+      clim,
+      prorataECS,
+      prorataChauffage,
+      dateDpe,
+      coef_ep,
+      th === 'immeuble' ? Nb_lgt : 1,
+      ratioSurfaceCoutCollectif,
+      facteursAc
+    );
+    conso.emission_ges = consoNette.emission_ges;
+    conso.cout = consoNette.cout;
+  }
+
   const conso1_7 = calc_conso(
     Sh,
     zc_id,

package/engine.spec.js

@@ -444,7 +444,8 @@ describe('calcul_3cl - run intégral (maison)', () => {
       'maison',
       'zc1',
       undefined,
-      0.25
+      0.25,
+      {}
     ]);
     expect(calculateEnr.mock.calls[1][5]).toBe(1.7);
     expect(calculateEnr.mock.calls[2][5]).toBe(2.3);
@@ -494,6 +495,29 @@ describe('calcul_3cl - run intégral (maison)', () => {
     expect(calc_clim).toHaveBeenCalledTimes(1);
   });
 
+  test('autoconsommation PV : GES et coût recalculés sur les consommations nettes (Tribu Calcul_batiment.cs l.1188 / l.1318)', () => {
+    const dpe = makeDpe();
+    const facteurs = { ch: 0.8 };
+    vi.mocked(calculateEnr).mockImplementationOnce((enr, conso, Sh, th, zc, coef, ratio, f) => {
+      Object.assign(f, facteurs);
+      return { pv: 1 };
+    });
+    const consoNette = {
+      ...consoResult,
+      emission_ges: { emission_ges_5_usages: 1 },
+      cout: { cout_5_usages: 2 }
+    };
+    vi.mocked(calc_conso).mockReturnValueOnce(consoResult).mockReturnValueOnce(consoNette);
+
+    calcul_3cl(dpe, { sanitize: false });
+
+    expect(calc_conso).toHaveBeenCalledTimes(4);
+    expect(calc_conso.mock.calls[1][13]).toEqual(facteurs);
+    expect(calc_conso.mock.calls[1][10]).toEqual(calc_conso.mock.calls[0][10]);
+    expect(dpe.logement.sortie.emission_ges).toEqual({ emission_ges_5_usages: 1 });
+    expect(dpe.logement.sortie.cout).toEqual({ cout_5_usages: 2 });
+  });
+
   test('bâtiment matériaux anciens => ilpa = 1 ; incohérence inertie => erreur loggée', () => {
     const dpe = makeDpe();
     dpe.logement.meteo.batiment_materiaux_anciens = 1;
@@ -579,6 +603,21 @@ describe('calcul_3cl - surfaces immeuble / appartement et prorata', () => {
     expect(conso_aux_distribution_ecs.mock.calls[0][3]).toBe(200);
   });
 
+  test('immeuble avec autoconsommation PV : conso nette calculée avec le nombre de logements', () => {
+    enums.methode_application_dpe_log = { 1: 'dpe immeuble collectif' };
+    const dpe = makeDpe();
+    vi.mocked(calculateEnr).mockImplementationOnce((enr, conso, Sh, th, zc, coef, ratio, f) => {
+      Object.assign(f, { ch: 0.8 });
+      return { pv: 1 };
+    });
+
+    calcul_3cl(dpe, { sanitize: false });
+
+    // Nombre de logements de l'immeuble (nombre_appartement = 4) pour le recalcul net.
+    expect(calc_conso.mock.calls[1][11]).toBe(4);
+    expect(calc_conso.mock.calls[1][13]).toEqual({ ch: 0.8 });
+  });
+
   test('immeuble : systèmes ECS non tous individuels sans surfaces => repli sur la division par deux', () => {
     enums.methode_application_dpe_log = { 1: 'dpe immeuble collectif' };
     const dpe = makeDpe();

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@open3cl/engine",
-  "version": "1.8.42",
+  "version": "1.8.44",
   "description": "Open Source 3CL-DPE engine",
   "main": "index.js",
   "directories": {