Python数据分析实例，利用Pandas建立流域三层蒸发和蓄满产流模型

最新推荐文章于 2024-06-03 19:22:50 发布
最新推荐文章于 2024-06-03 19:22:50 发布
阅读量4.3k 收藏 82
点赞数 12
分类专栏： 水文 文章标签： python 数据分析
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Joyonlyonly/article/details/104843095
版权

流域三层蒸发和蓄满产流模型

根据水量平衡方程推求土壤含水量

if i ==1:
	df["WU"][i] = WUM
	df["WL"][i] = WLM
	df["WD"][i] = WDM
              
else:
	dw=df["PE"][i-1]-df['R'][i-1]
	df["WU"][i] = df["WU"][i-1]+dw
	df["WL"][i] = df["WL"][i-1]
	df["WD"][i] = df["WD"][i-1]
	if df["WU"][i] < 0:
		df["WL"][i] = df["WL"][i-1]+df["WU"][i]
		df["WU"][i] = 0
		df["WD"][i] = df["WD"][i-1]
		if df["WL"][i] < 0:
			df["WD"][i] = df["WD"][i-1]+df["WL"][i]
			df["WL"][i] = 0
			df["WU"][i] = 0
			if df["WD"][i]<0:
				df["WL"][i] = 0
				df["WU"][i] = 0
				df["WD"][i] = 0
                                
	if df["WU"][i] > WUM:
		df["WL"][i] = df["WU"][i]-WUM+df["WL"][i-1]
		df["WU"][i] = WUM
		df["WD"][i] = df["WD"][i-1]
		if df["WL"][i] > WLM:
			df["WD"][i] = df["WL"][i]-WLM+df["WD"][i-1]
			df["WL"][i] = WLM
			df["WU"][i] = WUM
			if df["WD"][i] > WDM:
				df["WL"][i] = WLM
				df["WU"][i] = WUM
				df["WD"][i] = WDM

土壤三层蒸发模型，EU,EL,ED的计算

if df["WU"][i]+df["P"][i] >= df["EP"][i]:
	df["EU"][i]=df["EP"][i]
	df["EL"][i]=0
	df["ED"][i]=0
if df["WU"][i]+df["P"][i] < df["EP"][i] and df["WL"][i] >= C*WLM:
	df["EU"][i]=df["WU"][i]+df["P"][i]
	df["EL"][i]=(df["EP"][i]-df["EU"][i])*(df["WL"][i]/WLM)
	df["ED"][i]=0
if df["WU"][i]+df["P"][i] < df["EP"][i] and C*(df["EP"][i]-df["EU"][i])<= df["WL"][i] and df["WL"][i]<C*WLM:
	df["EU"][i]=df["WU"][i]+df["P"][i]
	df["EL"][i]=C*(df["EP"][i]-df["EU"][i])
	df["ED"][i]=0
if df["WU"][i]+df["P"][i] < df["EP"][i] and df["WL"][i]<C*(df["EP"][i]-df["EU"][i]):
	df["EU"][i]=df["WU"][i]+df["P"][i]
	df["EL"][i]=df["WL"][i]
	df["ED"][i]=C*(df["EP"][i]-df["EU"][i])-df["EL"][i]

产流R计算

a = WMM*(1-math.pow((1-df["W"][i]/WM),(1/(1+b))))
if a + df["PE"][i] <= WMM:
	df["R"][i] = df["PE"][i]+df["W"][i]-WM+WM*math.pow(1-((df["PE"][i]+a)/WMM),b+1)

else:
	df["R"][i] = df["PE"][i]+df["W"][i]-WM
if df["R"][i] < 0:
	df["R"][i] = 0

分水源Rg,Rd计算

if df["R"][i]>0:
	df['R/PE'][i]=(df["R"][i])/(df["PE"][i])
	if df['R/PE'][i] < 1:
		if df["PE"][i] <= FC:
			df["RD"][i] = 0
			df["RG"][i] = df["PE"][i]*df['R/PE'][i]
		else:
			df["RG"][i] = FC*df['R/PE'][i]
			df["RD"][i] = (df["PE"][i]-FC)*df['R/PE'][i]
		if df['R/PE'][i] == 1:
			if df["PE"][i] > FC:
				df["RG"][i] = FC
				df["RD"][i] = df["PE"][i]-FC
			else:
				df["RG"][i] = df["PE"][i]
				df["RD"][i] = 0

完整代码

import pandas as pd 
import numpy as np
import math
file_path = '/Users/joythegreat/Desktop/三层蒸发/'

#设定常量
WUM,WLM,WDM,b,C,FC= 20,60,40,0.3,(1/6),2.5
WM = WUM+WLM+WDM
WMM = WM*(1+b)

#读入表格
df = pd.read_excel(file_path + '大作业1.xlsx',index_col='日期')
df = df.fillna(0)


for i in range(1,154):
        
        #计算WU,WL,WD
        if i ==1:
                df["WU"][i] = WUM
                df["WL"][i] = WLM
                df["WD"][i] = WDM
                
               
        else:
                dw=df["PE"][i-1]-df['R'][i-1]
                df["WU"][i] = df["WU"][i-1]+dw
                df["WL"][i] = df["WL"][i-1]
                df["WD"][i] = df["WD"][i-1]
                if df["WU"][i] < 0:
                        df["WL"][i] = df["WL"][i-1]+df["WU"][i]
                        df["WU"][i] = 0
                        df["WD"][i] = df["WD"][i-1]
                        if df["WL"][i] < 0:
                                df["WD"][i] = df["WD"][i-1]+df["WL"][i]
                                df["WL"][i] = 0
                                df["WU"][i] = 0
                                if df["WD"][i]<0:
                                        df["WL"][i] = 0
                                        df["WU"][i] = 0
                                        df["WD"][i] = 0
                                
                if df["WU"][i] > WUM:
                        df["WL"][i] = df["WU"][i]-WUM+df["WL"][i-1]
                        df["WU"][i] = WUM
                        df["WD"][i] = df["WD"][i-1]
                        if df["WL"][i] > WLM:
                                df["WD"][i] = df["WL"][i]-WLM+df["WD"][i-1]
                                df["WL"][i] = WLM
                                df["WU"][i] = WUM
                                if df["WD"][i] > WDM:
                                        df["WL"][i] = WLM
                                        df["WU"][i] = WUM
                                        df["WD"][i] = WDM
                                
                
                
                            
    
        #计算E
        if df["WU"][i]+df["P"][i] >= df["EP"][i]:
                df["EU"][i]=df["EP"][i]
                df["EL"][i]=0
                df["ED"][i]=0
        if df["WU"][i]+df["P"][i] < df["EP"][i] and df["WL"][i] >= C*WLM:
                df["EU"][i]=df["WU"][i]+df["P"][i]
                df["EL"][i]=(df["EP"][i]-df["EU"][i])*(df["WL"][i]/WLM)
                df["ED"][i]=0
        if df["WU"][i]+df["P"][i] < df["EP"][i] and C*(df["EP"][i]-df["EU"][i])<= df["WL"][i] and df["WL"][i]<C*WLM:
                df["EU"][i]=df["WU"][i]+df["P"][i]
                df["EL"][i]=C*(df["EP"][i]-df["EU"][i])
                df["ED"][i]=0
        if df["WU"][i]+df["P"][i] < df["EP"][i] and df["WL"][i]<C*(df["EP"][i]-df["EU"][i]):
                df["EU"][i]=df["WU"][i]+df["P"][i]
                df["EL"][i]=df["WL"][i]
                df["ED"][i]=C*(df["EP"][i]-df["EU"][i])-df["EL"][i]
                
        #计算E,PE,W
        df["E"][i] = df["EU"][i]+df["EL"][i]+df["ED"][i]
        df["PE"][i] = df["P"][i] - df["E"][i]
        df["W"][i] = df["WU"][i]+df["WL"][i]+df["WD"][i]

        #计算产流
        a = WMM*(1-math.pow((1-df["W"][i]/WM),(1/(1+b))))
        if a + df["PE"][i] <= WMM:
                df["R"][i] = df["PE"][i]+df["W"][i]-WM+WM*math.pow(1-((df["PE"][i]+a)/WMM),b+1)

        else:
                df["R"][i] = df["PE"][i]+df["W"][i]-WM
        if df["R"][i] < 0:
                df["R"][i] = 0
                
        #分水源计算
        if df["R"][i]>0:
                df['R/PE'][i]=(df["R"][i])/(df["PE"][i])
                if df['R/PE'][i] < 1:
                        if df["PE"][i] <= FC:
                                df["RD"][i] = 0
                                df["RG"][i] = df["PE"][i]*df['R/PE'][i]
                        else:
                                df["RG"][i] = FC*df['R/PE'][i]
                                df["RD"][i] = (df["PE"][i]-FC)*df['R/PE'][i]
                if df['R/PE'][i] == 1:
                        if df["PE"][i] > FC:
                                df["RG"][i] = FC
                                df["RD"][i] = df["PE"][i]-FC
                        else:
                                df["RG"][i] = df["PE"][i]
                                df["RD"][i] = 0


                        
df.loc['和'] = df.apply(lambda x: x.sum())
pd.set_option('display.max_columns',None)
pd.set_option('display.width',1000)

#保留一位小数
for i in range(len(df)):
    for j in range(len(df.iloc[i])):
        df.iloc[i][j]=round(df.iloc[i][j],1)
#写入CSV文件
df.to_csv(file_path + '大作业.csv')

Python数据分析萌新，用pandas完成了一次专业相关的作业，有点兴奋，如有错误，还请指正~

Joyonlyonly
关注
专栏目录
[Python从零到壹] 十七.可视化分析之Matplotlib、Pandas、Echarts入门万字详解
杨秀璋的专栏
08-24 3万+
前一篇文章讲述了数据预处理、Jieba分词和文本聚类知识，这篇文章主要介绍Matplotlib和Pandas扩展包绘图的基础用法，同时引入Echarts技术，该技术主要应用于网站可视化展示中。本文内容以实例为主，给读者最直观的图形感受。两万字基础文章，希望对您有所帮助。
Python实战 | 使用 Python 的日志库（logging）和 pandas 库对日志数据进行分析
Web3 & Basketball
11-08 4458
在 Python 中，实现日志收集和分析的方法有很多，这里我为您介绍一个简单的示例，使用 Python 的日志库（logging）和 pandas 库对日志数据进行分析。以上代码将模拟的日志数据保存到 CSV 文件，并使用 pandas 对其进行简单的统计和分析。这些日志记录将分别对应不同的日志级别，从低到高依次为：DEBUG、INFO、WARNING、ERROR、CRITICAL。4. 运行上述代码后，您将在当前目录下看到一个名为 app.log 的日志文件，其中包含了您刚刚记录的日志。
水文预报新安江模型原理及Matlab代码
最新发布
JoeyChiu_的博客
06-03 2527
详细介绍水文预报中的新安江模型原理，并附带Matlab完整代码
三层蒸发蓄满产流模型（新安江模型）python计算程序.py
06-04
这是一个水文水资源学上的三层蒸发模式的蓄满产流模型（新安江模型）的python计算程序，本程序的运行需要结合本人上传的PEdata.xlsx文件，该excel文件包含了计算蓄满产流所需要的降水、蒸发能力数据。
三层蒸发模型
03-18
三层蒸发模型
matlab建立三层蒸发模型,水文模型三层蒸发模型的问题（初学者）
weixin_33501343的博客
03-17 2247
按土壤垂向分布的不均匀性将土层分为三层，用三层蒸散发模型计算蒸散发量。参数有流域平均张力水容量WM(mm)，上层张力水容量WUM(mm)，下层张力水容量WLM(mm)，深层张力水容量WDM(mm)，深层蒸散发扩散系数C。具体计算为若P+WU>=EP,则EU=EP,EL=0,ED=0; 若P+WU若WL>C*LM,则WL=(EP-EU)WL/LM,ED=0;若WL=C*(EP-EU),则...
python数据分析基础007 -利用pandas带你玩转excel表格（中上篇）
苏凉.py的博客
04-07 9832
利用pandas带你玩转excel，建议收藏！！
100天精通Python（数据分析篇）——第75天：Pandas数据预处理之数据标准化
热门推荐
努力让自己发光，对的人才能迎着光而来
02-13 1万+
专栏导读 1. 数据标准化是什么？ 2. 数据标准化的作用 3. 数据标准化的方法 4. 离差标准化 5. 标准差标准化 6. 小数定标标准化
三层蒸发_三层蒸发_
10-01
这是新安江模型的三层蒸发模型计算，计算蒸发
蓄满产流分层计算模型
11-10
蓄满产流分层计算模型蓄满产流分层计算模型蓄满产流分层计算模型
python数据分析基础010 -利用pandas带你玩转excel表格（终篇）
苏凉.py的博客
04-14 4192
带你玩转excel表格，建议收藏！！
《Pandas数据分析》（六）——模型的建立
qq_41430142的博客
08-25 1178
第三章 模型搭建和评估–建模 经过前面的两章的知识点的学习，我可以对数数据的本身进行处理，比如数据本身的增删查补，还可以做必要的清洗工作。那么下面我们就要开始使用我们前面处理好的数据了。这一章我们要做的就是使用数据，我们做数据分析的目的也就是，运用我们的数据以及结合我的业务来得到某些我们需要知道的结果。那么分析的第一步就是建模，搭建一个预测模型或者其他模型；我们从这个模型的到结果之后，我们要分析我的模型是不是足够的可靠，那我就需要评估这个模型。 我们拥有的泰坦尼克号的数据集，那么我们这次的目的就是，完成泰坦
Python 建立流域三层蒸发和蓄满产流模型(二水源划分)
m0_52964920的博客
03-11 4489
python pandas建立流域三层蒸发和蓄满产流模型(二水源划分)。
Python 水文预报三水源新安江模型构建
m0_52964920的博客
03-18 7206
Python构建新安江模型，包括三层蒸发模式、产流计算、分水源计算、汇流计算，既能学习原理又能引导写出自己的代码。
时段单位线法汇流计算python程序
qq_33015657的博客
06-09 3070
产汇流理论是水文学的重要理论之一。本文给出了汇流计算的经典方法——时段单位线法的python计算程序。需要的基础数据是时段单位线和地面净雨过程。 '''这是一个时段为6h,时段个数为15的单位线推流计算程序''' import numpy as np num=15 ur=10#单位地面净雨量，通常取10mm t=[i for i in range(1,num+1)]#时段序列 h=[24,23,3.2]#地面净雨过程 uh=[0,2,15,35,41,25,15,9,6,4,3,2,1,0]#时段单位线
水文预报三层蒸发模型与径流划分python程序
qq_54521224的博客
04-11 474
df = pd.read_csv("10月.csv", index_col='月份')df.to_csv('10月.csv')# 计算WU,WL,WD。# 计算E,PE,W。
整个建模流程
qq_41721660的博客
02-16 300
逻辑回归建模流程代码
评论 7
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值