Skip to content

Instantly share code, notes, and snippets.

@cocomice
Last active October 25, 2023 09:20
Show Gist options
  • Select an option

  • Save cocomice/d96709b107ce03ad8cb41526a2e89fce to your computer and use it in GitHub Desktop.

Select an option

Save cocomice/d96709b107ce03ad8cb41526a2e89fce to your computer and use it in GitHub Desktop.
import numpy as np
def calc_urban_wasteWat(N, Fc, eta):
"""
年城镇生活污水原位产生量计算
Args:
N (float): 人口数量,万人
Fc (float): 人均用水量系数,L/(人·d)
eta (float): 折污系数
Returns:
Gc (float): 城镇生活污水产生量,万吨
"""
Gc = N * Fc * eta * 365 / 1000
return Gc
def calc_urban_wasteWat_discharged(Gc, G_rc):
"""
年城镇生活污水排放量计算
Args:
Gc (float): 生活污水产生量,万吨
G_rc (float): 生活污水处理量,万吨
Returns:
Gc (float): 城镇生活污水排放量,万吨
"""
G_pc = Gc - G_rc
return G_pc
def calc_rural_wasteWat(N, Fc):
"""
年农村生活污水原位产生量计算
Args:
N (float): 人口数量,万人
Fc (float): 人均用水量系数,L/(人·d)
Returns:
Gn (float): 农村生活污水产生量,万吨
"""
Gn = N * Fc * 365 / 1000
return Gn
def calc_rural_wasteWat_discharged(G_cn, theta, eta_n):
"""
年农村生活污水排放量计算
Args:
G_cn (float): 农村生活污水产生量,万吨
theta (float): 对生活污水进行处理的行政村比例
eta_n (float): 污染物综合去除率
Returns:
G_pn (float): 农村生活污水排放量,万吨
"""
G_pn = G_cn * (1 - theta * eta_n)
return G_pn
def calc_cropLand_waste_discharged(A_g, e_gj, A_y, e_y, q_j, q_o):
"""
针对j类污染物的年农田面源污染物产生量计算
Args:
A_g (float): 农田面积,万亩
e_gj (float): 农田面源污染物流失系数, kg/(万亩·a)
A_y (float): 园地种植面积,万亩
e_y (float): 园地污染物流失系数, kg/(万亩·a)
q_j (float): 调查年度用于种植业的含氮化肥(含磷化肥)单位面积使用量,公斤/公顷
q_o (float): 用于种植业的含氮化肥(含磷化肥)单位面积建议使用量,公斤/公顷。含氮化肥用量指氮肥和含氮复合肥的折纯用量;含磷化肥用量指磷肥和含磷复合肥的折纯用量。
Returns:
G_gj (float): 种植业第j项污染物排放流失量 ,吨
j 指代污染物类型,如氨氮,总磷,总氮,化学需氧量
"""
G_gj = (A_g * e_gj + A_y * e_y) * q_j / q_o / 1000
return G_gj
def calc_livestock_waste_discharged(q_large, f_large, q_small, f_small):
"""
针对i类畜禽和j类污染物的畜禽养殖业污染物排放量
Args:
q_large (float): 第i类畜禽规模化养殖场的存/出栏量, 头/羽
f_large (float): 第i类畜禽规模化畜禽养殖第j项污染物产污系数(排污系数), kg/(头/羽)
q_small (float): 第i类畜禽养殖户养殖场的存/出栏量, 头/羽
f_small (float): 第i类畜禽养殖户畜禽养殖第j项污染物产污系数(排污系数), kg/(头/羽)
Returns:
Q_livestock (float): 第i类畜禽养殖第j项污染物产生量(排放量), 吨.
i 指代畜禽类型,如生猪,奶牛,肉牛,蛋鸡
j 指代污染物类型,如氨氮,总磷,总氮,化学需氧量
注:当使用排放系数计算式,得到的是畜禽养殖的污染物排放量。当使用产污系数计算式,得到的是畜禽养殖的污染物产生量。
"""
Q_livestock = (q_large * f_large + q_small * f_small) / 1000
return Q_livestock
def calc_urbanNP_waste_discharged(C_f, phi, A, P, EMC):
"""
城镇面源污染排放量核算
Args:
C_f (float): 不产生地表径流的降雨进行校正的因子
phi (float): 径流区平均径流系数
A (float): 径流集雨面积, 公顷
P (float): 计算时段(t)内的降雨量, mm
EMC (float): 多场降雨事件的平均污染物的径流量加权平均浓度,mg/L
Returns:
Lt (float): 计算时段(t)内径流排放污染负荷,kg;
"""
Lt = C_f * phi * A * P * EMC * 100
return Lt
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment