Skip to content

Build biomass potential

build_biomass_potential_xing(biomass_potentials_path)

build potential from Xing et al. https://doi.org/10.1038/s41467-021-23282-x

Parameters:

Name Type Description Default
biomass_potentials_path PathLike

the path to the Xing SI data (xlsx).

 required
Source code in workflow/scripts/build_biomass_potential.py 
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
def build_biomass_potential_xing(biomass_potentials_path: PathLike):
    """build potential from Xing et al. https://doi.org/10.1038/s41467-021-23282-x

    Args:
        biomass_potentials_path (PathLike, optional): the path to the Xing SI data (xlsx).
    """

    df = read_xing_si_data(biomass_potentials_path)

    # select only relevant part of potential
    df = df[df.columns[df.columns.str.contains("Agricultural residues burnt as waste")]].sum(axis=1)

    # convert t biomass yr-1 to MWh, heat content is from paper reference 92
    heat_content = 19  # GJ (t biomass−1)
    heat_content *= 1000 / 3600  # GJ/t -> MWh
    df = df * heat_content

    return df

estimate_co2_intensity_xing()

Estimate the biomass Co2 intensity from the Xing paper

Returns:

Name Type Description
float float

the biomass co2 intensity in t/MWhth
Source code in workflow/scripts/build_biomass_potential.py 
12
13
14
15
16
17
18
19
20
21
22
23
24
25
def estimate_co2_intensity_xing() -> float:
    """Estimate the biomass Co2 intensity from the Xing paper

    Returns:
        float: the biomass co2 intensity in t/MWhth
    """

    biomass_potential_tot = 3.04  # Gt
    embodied_co2_tot = 5.24  # Gt
    heat_content = 19 * 1000 / 3600  # GJ/t -> MWh_th/t
    unit_co2 = embodied_co2_tot / biomass_potential_tot  # t CO2/t biomass
    co2_intens = unit_co2 / heat_content  # t CO2/MWh_th

    return co2_intens

read_xing_si_data(biomass_potentials_path)

read and prepare the xing SI data

Parameters:

Name Type Description Default
biomass_potentials_path PathLike

the path to the Xing SI data (xlsx).

 required
Source code in workflow/scripts/build_biomass_potential.py 
28
29
30
31
32
33
34
35
36
37
38
39
40
41
def read_xing_si_data(biomass_potentials_path: PathLike):
    """read and prepare the xing SI data

    Args:
        biomass_potentials_path (PathLike): the path to the Xing SI data (xlsx).
    """
    # data is indexed by province and county
    df = pd.read_excel(biomass_potentials_path, sheet_name="supplementary data 1")
    df = df.groupby("Province name").sum()

    df = df.rename(index={"Inner-Monglia": "InnerMongolia", "Anhui ": "Anhui"})
    df = df.add_suffix(" biomass")

    return df