🧄Python时间序列模型使用河流-地下水-泵水数据计算地下水位下降
Last updated
Was this helpful?
Last updated
Was this helpful?
Was this helpful?
Python | 时间序列 | 地下水 | 水文 | 泵水 | 河流水位 | 传递函数噪声建模 | 物理脉冲响应函数 | 降水过剩 | 河流水位变化 | 自回归噪音模型 | Makkink参考蒸发 | Hantush井函数 | 贝塞尔函数 | 赤池信息准则 | 畜水层 | 稳态回撤
import json
import logging
import sys
sys.path.insert(1, "..")
from util import get_pastastore
logging.basicConfig(level=logging.INFO)
with open("../tsa_settings.json", "r") as f:
settings = json.load(f)
name = f"frontiers_{settings['freq']}"
pstore = get_pastastore(name, settings["connector_type"])
conn = pstore.conn
prec = hpd.ObsCollection.from_knmi(xmid=oseries.x, ymid=oseries.y,
meteo_vars=["RH"], start=["1980-01-01"],
ObsClass=[hpd.PrecipitationObs])
evap = hpd.ObsCollection.from_knmi(xmid=oseries.x, ymid=oseries.y,
meteo_vars=["EV24"], start=["1980-01-01"],
ObsClass=[hpd.EvaporationObs])
# %% write to store
for o in tqdm(prec.loc[:, "obs"].values, desc="Add precipitation series"):
meta = _get_metadata_from_obs(o)
# convert to mm/day
conn.add_stress(o["RH"] * 1e3, o.name, "prec",
metadata=meta, overwrite=True)
for o in tqdm(evap.loc[:, "obs"].values, desc="Add evaporation series"):
meta = _get_metadata_from_obs(o)
# convert to mm/day
conn.add_stress(o["EV24"] * 1e3, o.name, "evap",
metadata=meta, overwrite=True)