CN116911071B - 一种液体氮肥施放方式的确定方法、装置、设备及介质 - Google Patents

一种液体氮肥施放方式的确定方法、装置、设备及介质 Download PDF

Info

Publication number
CN116911071B
CN116911071B CN202311147826.6A CN202311147826A CN116911071B CN 116911071 B CN116911071 B CN 116911071B CN 202311147826 A CN202311147826 A CN 202311147826A CN 116911071 B CN116911071 B CN 116911071B
Authority
CN
China
Prior art keywords
nitrogen element
preset
migration
soil
mode
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202311147826.6A
Other languages
English (en)
Other versions
CN116911071A (zh
Inventor
刘肖廷
刘琨
戴会超
马哲
秦明
李婉
程建
张璐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Gezhouba Electric Power Rest House
China Three Gorges Corp
Original Assignee
Beijing Gezhouba Electric Power Rest House
China Three Gorges Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Gezhouba Electric Power Rest House, China Three Gorges Corp filed Critical Beijing Gezhouba Electric Power Rest House
Priority to CN202311147826.6A priority Critical patent/CN116911071B/zh
Publication of CN116911071A publication Critical patent/CN116911071A/zh
Application granted granted Critical
Publication of CN116911071B publication Critical patent/CN116911071B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Fertilizing (AREA)
  • Cultivation Of Plants (AREA)

Abstract

本发明涉及工程仿真与数值模拟技术领域,公开了一种液体氮肥施放方式的确定方法、装置、设备及介质。本实施例首先通过室内土柱试验得到试验数据,并将试验数据与氮元素迁移模拟模型结合得到预设施放方式对应的氮元素的迁移过程;然后,通过迁移过程判断预设施放方式是否为目标施放方式,如果不是,则不断调整施放方式,反复进行试验,直至得到目标施放方式;最后,将液体氮肥通过目标施放方式等比例施放至目标沙漠地区中。本实施例通过室内土柱试验的方式确定目标施放方式不会造成不可逆的损失。而且通过真实的试验数据和预构建的氮元素迁移模拟模型模拟氮元素在沙漠地区中的迁移过程,使得试验结果准确性和真实性更强,更具说服力。

Description

一种液体氮肥施放方式的确定方法、装置、设备及介质
技术领域
本发明涉及工程仿真与数值模拟技术领域,具体涉及一种液体氮肥施放方式的确定方法、装置、设备及介质。
背景技术
随着光伏治沙模式的成功,沙漠地区绿色生态农牧业得到较快发展。而在荒漠中,通过使用土基材料,可以使得沙漠表层土壤化,从而增强沙漠表层对水分、养分和空气的存储能力,使作物得以生长。
但是由于沙漠地区水资源匮乏,土壤肥力差,因此,通常会采用液体氮肥以提升土壤中的氮元素。但是,如果施肥过度,又会导致土壤板结硬化。
因此,如何确定液体氮肥的施放量,能够在不造成土壤板结硬化的前提下,还能够保障土壤中的氮元素的含量,成为目前需要解决的问题。
发明内容
有鉴于此,本发明提供了一种液体氮肥施放方式的确定方法、装置、设备及介质,以解决由于施放方式不准确而导致土壤硬化板结硬化的问题。
第一方面,本发明提供了一种液体氮肥施放方式的确定方法,方法包括:
获取目标沙漠地区的垂向分布信息、环境信息、土基材料信息和沙漠沙信息;基于垂向分布信息、环境信息和土基材料信息和沙漠沙信息形成室内土柱,并将液体氮肥按照预设施放方式施放至室内土柱中,开展室内土柱试验,得到各预设位置分别对应的氮元素浓度;将各预设位置分别对应的氮元素浓度输入至预构建的氮元素迁移模拟模型中,得到与预设施放方式对应的氮元素的迁移过程;基于迁移过程判断预设施放方式是否为目标施放方式;当预设施放方式不是目标施放方式时,调整液体氮肥的施放方式,重新进行试验,直至基于迁移过程确定调整后的施放方式为目标施放方式时,结束试验,并将液体氮肥按照目标施放方式等比例施放于目标沙漠地区中。
本实施例提供的液体氮肥施放方式的确定方法,首先通过室内土柱试验得到试验数据,并将试验数据与预构建的氮元素迁移模拟模型结合,得到预设施放方式对应的氮元素的迁移过程;然后,通过迁移过程判断预设施放方式是否为目标施放方式,如果不是,通过不断调整施放方式,反复进行试验,得到目标施放方式;最后,将液体氮肥通过目标施放方式等比例施放至目标沙漠地区中。相比直接在目标沙漠地区进行试验,本实施例通过室内土柱试验的方式确定目标施放方式,不具破坏性,且不会造成不可逆的损失。而且通过真实的试验数据和预构建的氮元素迁移模拟模型模拟氮元素在沙漠地区中的迁移过程,使得试验结果准确性和真实性更强,更具说服力。
在一种可选的实施方式中,预构建的氮元素迁移模拟模型由双阶段的时间分数阶对流扩散反应方程和Hausdroff对流扩散反应方程耦合得到。
在一种可选的实施方式中,氮元素迁移模拟模型通过如下公式表示:
其中,为地表位置,/>为土基材料与沙漠沙接触位置,/>为沙漠沙的下边界,/>为垂向空间位置,/>为时间分数阶导数,/>中C为时间分数阶导数所属的类型,/>为时间分数阶阶数,/>为豪斯多夫导数,/>为豪斯多夫导数阶数,/>为液体氮肥中氮元素的浓度,/>为氮元素背景浓度值值函数,/>为土基材料的流动参数,/>为沙漠沙的流动参数,/>为土基材料的扩散系数,/>为沙漠沙的扩散系数,/>为硝化速率、/>为反硝化速率。
本实施例提供的液体氮肥施放方式的确定方法,结合土基材料和沙子中水氮迁移转化规律,采用由双阶段的时间分数阶对流扩散反应方程和Hausdroff对流扩散反应方程耦合形成的氮元素迁移模拟模型,精确高效的模拟了氮元素的迁移转化过程,从而准确评估出目标施放方式,防止了土壤板结硬化的发生,为沙漠地区绿色生态农牧业的发展提供技术支持。
在一种可选的实施方式中,将各预设位置分别对应的氮元素浓度输入至预构建的氮元素迁移模拟模型中,得到与预设施放方式对应的氮元素的迁移过程,包括:
将预设位置分别对应的氮元素浓度输入至预构建的氮元素迁移模拟模型中,采用最小二乘法进行反演,得到与预设施放方式对应的迁移系数;将迁移系数输入至氮元素迁移模拟模型中,得到与预设施放方式对应的氮元素的迁移过程。
本发明实施例得到的迁移过程,结合了预设施放方式对应的真实试验数据,因此,该迁移过程对氮元素的迁移转化过程刻画的贴切和真实,使得目标施放方式的确定也更加准确。
在一种可选的实施方式中,迁移过程为垂直方向上各位置处氮元素的浓度分布。
在一种可选的实施方式中,基于迁移过程判断预设施放方式是否为目标施放方式,包括:
基于迁移过程,确定土柱中土基材料部分对应的氮元素积蓄量;基于迁移过程,确定土柱中沙漠沙下边界的氮元素输出量;将氮元素积蓄量与第一预设浓度进行比较,得到第一比较结果;将氮元素输出量与第二预设浓度进行比较,得到第二比较结果;基于第一比较结果和/或第二比较结果,判断预设施放方式是否为目标施放方式。
在一种可选的实施方式中,基于第一比较结果和第二比较结果,判断预设施放方式是否为目标施放方式,包括:
当第一比较结果为氮元素积蓄量小于第一预设浓度,且,当第二比较结果为氮元素输出量小于第二预设浓度时,将预设施放方式确定为目标施放方式。
本实施例将获取的氮元素积蓄量与第一预设浓度进行比较,将氮元素输出量与第二预设浓度进行比较,通过两次比较的结果确定出一种既不会造成土基材料板结,又不会造成地下水污染的目标施放方式,从而为沙漠地区绿色生态农牧业的发展做出了贡献。
第二方面,本发明提供了一种液体氮肥施放方式的确定装置,装置包括:
获取模块,用于获取目标沙漠地区的垂向分布信息、环境信息、土基材料信息和沙漠沙信息;第一确定模块,用于基于垂向分布信息、环境信息和土基材料信息和沙漠沙信息形成室内土柱,并将液体氮肥按照预设施放方式施放至室内土柱中,开展室内土柱试验,得到各预设位置分别对应的氮元素浓度;第二确定模块,用于将各预设位置分别对应的氮元素浓度输入至预构建的氮元素迁移模拟模型中,得到与预设施放方式对应的氮元素的迁移过程;判断模块,用于基于迁移过程判断预设施放方式是否为目标施放方式;调整模块,用于当预设施放方式不是目标施放方式时,调整液体氮肥的施放方式,重新进行试验,直至基于迁移过程确定调整后的施放方式为目标施放方式时,结束试验,并将液体氮肥按照目标施放方式等比例施放于目标沙漠地区中。
在一种可选的实施方式中,第二确定模块中的预构建的氮元素迁移模拟模型由双阶段的时间分数阶对流扩散反应方程和Hausdroff对流扩散反应方程耦合得到。
在一种可选的实施方式中,第二确定模块中的氮元素迁移模拟模型通过如下公式表示:
其中,为地表位置,/>为土基材料与沙漠沙接触位置,/>为沙漠沙的下边界,/>为垂向空间位置,/>为时间分数阶导数,/>中C为时间分数阶导数所属的类型,/>为时间分数阶阶数,/>为豪斯多夫导数,/>为豪斯多夫导数阶数,/>为液体氮肥中氮元素的浓度,/>为氮元素背景浓度值值函数,/>为土基材料的流动参数,/>为沙漠沙的流动参数,/>为土基材料的扩散系数,/>为沙漠沙的扩散系数,/>为硝化速率、/>为反硝化速率。
在一种可选的实施方式中,第二确定模块,包括:
第一确定子模块,用于将预设位置分别对应的氮元素浓度输入至预构建的氮元素迁移模拟模型中,采用最小二乘法进行反演,得到与预设施放方式对应的迁移系数;第二确定子模块,用于将迁移系数输入至氮元素迁移模拟模型中,得到与预设施放方式对应的氮元素的迁移过程。
在一种可选的实施方式中,第二确定子模块中的迁移过程为垂直方向上各位置处氮元素的浓度分布。
在一种可选的实施方式中,判断模块,包括:
第三确定子模块,用于基于迁移过程,确定土柱中土基材料部分对应的氮元素积蓄量;第四确定子模块,用于基于迁移过程,确定土柱中沙漠沙下边界的氮元素输出量;第一比较子模块,用于将氮元素积蓄量与第一预设浓度进行比较,得到第一比较结果;第二比较子模块,用于将氮元素输出量与第二预设浓度进行比较,得到第二比较结果;判断子模块,用于基于第一比较结果和/或第二比较结果,判断预设施放方式是否为目标施放方式。
第三方面,本发明提供了一种计算机设备,包括:存储器和处理器,存储器和处理器之间互相通信连接,存储器中存储有计算机指令,处理器通过执行计算机指令,从而执行上述第一方面或其对应的任一实施方式的液体氮肥施放方式的确定方法。
第四方面,本发明提供了一种计算机可读存储介质,该计算机可读存储介质上存储有计算机指令,计算机指令用于使计算机执行上述第一方面或其对应的任一实施方式的液体氮肥施放方式的确定方法。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是根据本发明实施例的液体氮肥施放方式的确定方法的流程示意图;
图2是根据本发明实施例的沙漠地区垂向分布示意图;
图3是根据本发明实施例的另一液体氮肥施放方式的确定方法的流程示意图;
图4是根据本发明实施例的又一液体氮肥施放方式的确定方法的流程示意图;
图5是根据本发明实施例的液体氮肥施放方式的确定装置的结构框图;
图6是本发明实施例的计算机设备的硬件结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
随着光伏治沙模式的成功,沙漠地区的农牧业得到了长足的发展。通过在荒漠中使用土基材料,使得沙漠表层土壤化,从而促进作物的生长。但是,由于沙漠地区缺乏水资源,而且土壤肥力差,因此,需要通过施加液体氮肥提升土壤中的氮元素。但是,如果液体氮肥施放过量,就会导致土壤硬化板结。因此,本发明实施例提供了一种液体氮肥施放方式的确定方法,通过室内土柱试验采集的试验数据与预构建的氮元素迁移模拟模型确定出当前液体氮肥的施放方式对应的氮元素迁移过程,从而基于该迁移过程判断当前施放方式是否会造成土壤硬化板结,进而确定出一种既能保证土壤中保有充足氮元素,还不会造成土壤硬化板结的液体氮肥的施放方式。
根据本发明实施例,提供了一种液体氮肥施放方式的确定方法实施例,需要说明的是,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
在本实施例中提供了一种液体氮肥施放方式的确定方法,可用于计算机设备,图1是根据本发明实施例的液体氮肥施放方式的确定方法的流程图,如图1所示,该流程包括如下步骤:
步骤S101,获取目标沙漠地区的垂向分布信息、环境信息、土基材料信息和沙漠沙信息。
具体的,垂向分布信息即目标沙漠地区垂直方向上的层级分布,如图2所示,自上而下分别为地表、土基材料以及沙漠沙。
具体的,环境信息即目标沙漠地区所处环境的基本信息,环境信息包括环境温度、环境湿度、风速等。
具体的,土基材料信息包括土基材料的类别、配比、厚度、导水率、渗透系数、土基材料中污染物背景浓度值、土基材料中的初始氮元素浓度等。
具体的,沙漠沙信息包括沙漠沙导水率、渗透系数、沙漠沙中污染物背景浓度值、沙漠沙中的初始氮元素浓度。
示例性地,一般而言,氮元素的背景浓度值为连续的数据,因此,通常通过背景浓度值值函数来表示土基材料中的初始氮元素浓度和沙漠沙中初始氮元素浓度。需要说明的是,背景浓度值值函数是氮元素浓度关于垂向位置的函数。
示例性地,沙漠沙的导水率和渗透系数是通过分析沙漠沙样本中的粒径组成后,经过测量得到的。测量方式此处不作具体限定,本领域技术人员可以实现结果为依据自行选择。
示例性地,土基材料中污染物背景浓度值和沙漠沙中污染物背景浓度值,是通过在目标沙漠地区布置测站进行监测得到的。
步骤S102,基于垂向分布信息、环境信息和土基材料信息和沙漠沙信息形成室内土柱,并将液体氮肥按照预设施放方式施放至室内土柱中,开展室内土柱试验,得到各预设位置分别对应的氮元素浓度。
具体的,开展室内土柱试验,室内土柱的垂向分布与目标沙漠地区的垂向分布信息保持一致。且土柱中土基材料的选取和分布与土基材料信息保持一致,沙漠沙的选取和分布与沙漠沙信息保持一致。在形成稳定土层之后,控制室内土柱上边界温度、湿度和风速与沙漠地区的环境信息保持一致。
且在室内土柱的侧方设置有均匀分布的取样孔,此处的取样孔即预设位置。在各预设位置处布设水质分析仪,用于监测示踪剂浓度,从而测量该位置处氮元素的浓度。
具体的,预设施放方式的可以根据本领域技术人员的经验进行确定,此处不作具体限定。施放方式包括施放液体氮肥的浓度和液体氮肥的施放量。
具体的,在室内土柱试验的前期准备均完成之后,将液体氮肥按照预设施放方式施放至室内土柱中,等待一段时间后,采集各预设位置处氮元素的浓度。需要说明的是,预设位置一部分位于土基材料中,一部分位于沙漠沙材料中。
步骤S103,将各预设位置分别对应的氮元素浓度输入至预构建的氮元素迁移模拟模型中,得到与预设施放方式对应的氮元素的迁移过程。
具体的,氮元素的迁移过程与液体氮肥的施放方式对应,施放方式不同,迁移过程也不完全相同。
步骤S104,基于迁移过程判断预设施放方式是否为目标施放方式。
具体的,目标施放方式为既能保证土壤中氮元素保有量最多,又不会造成土壤硬化板结的施放方式。
步骤S105,当预设施放方式不是目标施放方式时,调整液体氮肥的施放方式,重新进行试验,直至基于迁移过程确定调整后的施放方式为目标施放方式时,结束试验,并将液体氮肥按照目标施放方式等比例施放于目标沙漠地区中。
具体的,当预设施放方式不是目标施放方式时,可以通过调整液体氮肥的施放量或者调整施放液体氮肥的浓度调整施放方式,并重新进行步骤S102至步骤S104的操作,直至调整后的施放方式为目标施放方式时,结束试验。
本实施例提供的液体氮肥施放方式的确定方法,首先通过室内土柱试验得到试验数据,并将试验数据与预构建的氮元素迁移模拟模型结合,得到预设施放方式对应的氮元素的迁移过程;然后,通过迁移过程判断预设施放方式是否为目标施放方式,如果不是,通过不断调整施放方式,反复进行试验,得到目标施放方式;最后,将液体氮肥通过目标施放方式等比例施放至目标沙漠地区中。相比直接在目标沙漠地区进行试验,本实施例通过室内土柱试验的方式确定目标施放方式,不具破坏性,且不会造成不可逆的损失。而且通过真实的试验数据和预构建的氮元素迁移模拟模型模拟氮元素在沙漠地区中的迁移过程,使得试验结果准确性和真实性更强,更具说服力。
在本实施例中提供了一种液体氮肥施放方式的确定方法,可用于计算机设备,图3是根据本发明实施例的液体氮肥施放方式的确定方法的流程图,如图3所示,该流程包括如下步骤:
步骤S301,获取目标沙漠地区的垂向分布信息、环境信息、土基材料信息和沙漠沙信息。详细请参见图1所示实施例的步骤S101,在此不再赘述。
步骤S302,基于垂向分布信息、环境信息和土基材料信息和沙漠沙信息形成室内土柱,并将液体氮肥按照预设施放方式施放至室内土柱中,开展室内土柱试验,得到各预设位置分别对应的氮元素浓度。详细请参见图1所示实施例的步骤S102,在此不再赘述。
步骤S303,将各预设位置分别对应的氮元素浓度输入至预构建的氮元素迁移模拟模型中,得到与预设施放方式对应的氮元素的迁移过程。
具体地,上述步骤S303包括:
步骤S3031,将预设位置分别对应的氮元素浓度输入至预构建的氮元素迁移模拟模型中,采用最小二乘法进行反演,得到与预设施放方式对应的迁移系数。
具体的,采用最小二乘法进行反演的过程已属于成熟的现有技术,此处不再赘述。
具体的,迁移系数包括,时间分数阶阶数、Hausdroff(豪斯多夫)导数阶数、土基材料的流动参数、沙漠沙的流动参数、土基材料的扩散系数、沙漠沙的扩散系数、硝化速率和反硝化速率。
具体的,预构建的氮元素迁移模拟模型由双阶段的时间分数阶对流扩散反应方程和Hausdroff对流扩散反应方程耦合得到。
示例性地,预构建的氮元素迁移模拟模型通过如下控制方程表示:
其中,为地表位置,/>为土基材料与沙漠沙接触位置,/>为沙漠沙的下边界,/>为垂向空间位置,/>为Caputo类型的时间分数阶导数,/>中C表示该类型导数为Caputo类分数阶导数,0表示初始时刻,/>表示时间分数阶阶数,代表土基材料整体对氮元素的滞留程度,/>,/>取值越低,滞留效果越强,/>为Hausdroff导数,定义为/>为Hausdroff导数阶数,由于沙漠沙无滞水效果,此处/>1,/>为液体氮肥中氮元素的浓度,/>为氮元素背景浓度值值函数,/>为土基材料的流动参数,/>为沙漠沙的流动参数,/>为土基材料的扩散系数,/>为沙漠沙的扩散系数,/>为硝化速率、/>为反硝化速率。
其中,时间分数阶导数定义为:
其中,时间Hausdroff导数的离散形式为:
此处,,/>表示第/>个时刻,每个时间短布点都为/>,总点数为/>
步骤S3032,将迁移系数输入至氮元素迁移模拟模型中,得到与预设施放方式对应的氮元素的迁移过程。
具体的,迁移系数与施放方式对应。
具体的,氮元素的迁移过程即将迁移系数输入至氮元素迁移模拟模型后,得到的由连续数据形成的控制方程。
具体的,迁移过程中包括氮元素在土基材料中的迁移过程,还包括氮元素在沙漠沙中的迁移过程。
步骤S304,基于迁移过程判断预设施放方式是否为目标施放方式。详细请参见图1所示实施例的步骤S104,在此不再赘述。
步骤S305,当预设施放方式不是目标施放方式时,调整液体氮肥的施放方式,重新进行试验,直至基于迁移过程确定调整后的施放方式为目标施放方式时,结束试验,并将液体氮肥按照目标施放方式等比例施放于目标沙漠地区中。详细请参见图1所示实施例的步骤S105,在此不再赘述。
本实施例提供的液体氮肥施放方式的确定方法,结合土基材料和沙子中水氮迁移转化规律,采用由双阶段的时间分数阶对流扩散反应方程和Hausdroff对流扩散反应方程耦合形成的氮元素迁移模拟模型,精确高效的模拟了氮元素的迁移转化过程,从而准确评估出目标施放方式,防止了土壤板结硬化的发生,为沙漠地区绿色生态农牧业的发展提供技术支持。
在本实施例中提供了一种液体氮肥施放方式的确定方法,可用于计算机设备,图4是根据本发明实施例的液体氮肥施放方式的确定方法的流程图,如图4所示,该流程包括如下步骤:
步骤S401,获取目标沙漠地区的垂向分布信息、环境信息、土基材料信息和沙漠沙信息。详细请参见图1所示实施例的步骤S101,在此不再赘述。
步骤S402,基于垂向分布信息、环境信息和土基材料信息和沙漠沙信息形成室内土柱,并将液体氮肥按照预设施放方式施放至室内土柱中,开展室内土柱试验,得到各预设位置分别对应的氮元素浓度。详细请参见图1所示实施例的步骤S102,在此不再赘述。
步骤S403,将各预设位置分别对应的氮元素浓度输入至预构建的氮元素迁移模拟模型中,得到与预设施放方式对应的氮元素的迁移过程。详细请参见图3所示实施例的步骤S303,在此不再赘述。
步骤S404,基于迁移过程判断预设施放方式是否为目标施放方式。
具体的,上述步骤S404包括:
步骤S4041,基于迁移过程,确定土柱中土基材料部分对应的氮元素积蓄量。
具体的,由于迁移过程为时间维度上,在垂直方向不同位置处氮元素的分布。因此,可以从迁移过程(即控制方程)中提取出土基材料所在高度不同位置处分别对应的氮元素浓度。将所有氮元素浓度求和,得到土基材料部分对应的氮元素积蓄量。
示例性地,以步骤S3031记载的实施例为例,土基材料所处高度在地表位置和土基材料与沙漠沙接触位置/>之间,因此,仅需要从第一个控制方程中获取不同位置处氮元素的浓度,通过相加求和得到土基材料对应的氮元素积蓄量。
步骤S4042,基于迁移过程,确定土柱中沙漠沙下边界的氮元素输出量。
具体的,基于氮元素在沙漠沙中的迁移过程,获取沙漠沙下边界位置处氮元素的浓度,并基于氮元素浓度和地下水的重量确定沙漠沙下边界的氮元素输出量。
具体的,沙漠沙下边界的氮元素输出量相当于地下水中氮元素输入量。氮元素输入量的确定方式还可以是,在室内土柱的下方放置细密纱网,并使其表层浸泡在盛水的容器中,容器下放置计重器,实验结束后通过计重器测量容器中氮元素浓度,那么,地下水中氮元素输入量即为/>,其中,G为容器中水体的重量。
步骤S4043,将氮元素积蓄量与第一预设浓度进行比较,得到第一比较结果。
具体的,第一预设浓度为使土基材料硬化的氮元素浓度的最小值。
具体的,当第一比较结果为氮元素积蓄量小于第一预设浓度时,表示土基材料中氮元素积蓄量未达到土基材料板结硬化的临界值,那么,如果按照当前施放方式给目标沙漠地区施放液体氮肥,不会造成土壤硬化板结。当第一比较结果为氮元素积蓄量大于或等于第一预设浓度时,表示土基材料中氮元素积蓄量已经超过土基材料板结硬化的临界值,如果按照当前施放方式给目标沙漠地区施放液体氮肥,那么会有很大可能造成土壤硬化板结。
步骤S4044,将氮元素输出量与第二预设浓度进行比较,得到第二比较结果。
具体的,第二预设浓度为造成地下水污染的氮元素浓度的最小值。
具体的,当第二比较结果为氮元素输出量小于第二预设浓度时,表示沙漠沙下边界氮元素输出量未达到污染地下水的程度,因此,如果按照预设施放方式给目标沙漠地区施放液体氮肥,则不会造成地下水污染。相反,当第二比较结果为氮元素输出量大于或等于第二预设浓度时,表示沙漠沙下边界氮元素输出量超过了污染地下水的氮元素浓度,那么,如果按照预设施放方式给目标沙漠地区施放液体氮肥,将会造成地下水污染。
步骤S4045,基于第一比较结果和/或第二比较结果,判断预设施放方式是否为目标施放方式。
如果,在液体氮肥的施放过程中需要及考虑地下水污染问题,还要考虑土基材料硬化板结问题时,目标释放方式的确定需要同时考虑第一比较结果和第二比较结果。即当第一比较结果为氮元素积蓄量小于第一预设浓度,且,当第二比较结果为氮元素输出量小于第二预设浓度时,将预设施放方式确定为目标施放方式。
如果,在液体氮肥的施放过程中仅需要考虑土基材料硬化板结的问题,那么,只需要基于第一比较结果判断预设施放方式是否为目标施放方式。即,当氮元素积蓄量小于第一预设浓度时,将预设施放方式确定为目标施放方式。
同理,如果在液体氮肥的施放过程中仅需要考虑是否会造成地下水污染时,那么只需要根据第二比较结果判断预设施放方式是否为目标施放方式。即,当第二比较结果为氮元素输出量小于第二预设浓度时,将预设施放方式确定为目标施放方式。
步骤S405,当预设施放方式不是目标施放方式时,调整液体氮肥的施放方式,重新进行试验,直至基于迁移过程确定调整后的施放方式为目标施放方式时,结束试验,并将液体氮肥按照目标施放方式等比例施放于目标沙漠地区中。详细请参见图1所示实施例的步骤S105,在此不再赘述。
作为本发明实施例中的一个最优实施例,下文将结合实际应用场景对本方案进行详细描述。
获取目标沙漠地区的垂向分布信息、环境信息、土基材料信息和沙漠沙信息。并按照沙漠沙地区的垂向分布信息、环境信息、土基材料信息和沙漠沙信息构建室内土柱,使得室内土柱所属的环境以及室内土柱土层的垂向分布等信息均与目标沙漠地区保持一致。
当室内土柱试验准备就绪后,按照预设施放方式将液体氮肥施放至室内土柱中,等待一段时间后,获取各预设位置分别对应的氮元素浓度。
将各预设位置分别对应的氮元素浓度输入至由双阶段的时间分数阶对流扩散反应方程和Hausdroff对流扩散反应方程耦合而成的氮元素迁移模拟模型中,通过最小二乘法反演出预设施放方式对应的迁移系数,并将迁移系数输入至氮元素迁移模拟模型中,得到与预设施放方式对应的氮元素的迁移过程。
根据氮元素在土基材料中迁移过程,确定土基材料部分对应的氮元素积蓄量,并将氮元素积蓄量与使土基材料板结硬化的氮元素浓度临界值(即第一预设浓度)进行比较,得到第一比较结果。
根据氮元素在沙漠沙中的迁移过程,确定沙漠沙下边界氮元素的输出量。将氮元素输出量与造成地下水污染的氮元素输出量临界值(即第二预设浓度)进行比较,得到第二比较结果。
当第一比较结果为氮元素积蓄量小于第一预设浓度,且,第二比较结果为氮元素输出量小于第二预设浓度时,表示预设施放方式既不会造成土基材料硬化板结,又不会造成地下水污染,此时,可将预设施放方式确定为目标施放方式。当第一比较结果和第二比较结果为除上述情况外的其它情况时,预设施放方式不能作为目标施放方式。此时,还需要重新调整施放方式,直至调整后的施放方式可以作为目标施放方式时,方可将液体氮肥以目标施放方式等比例施放于目标沙漠地区中。
本实施例通过将试验数据和预构建的氮元素迁移模拟模型结合,得到氮元素分别在土基材料和沙漠沙中对应的迁移过程,通过该迁移过程确定出一种既不会造成土基材料板结,又不会造成地下水污染的目标施放方式,从而为沙漠地区绿色生态农牧业的发展做出了贡献。
在本实施例中还提供了一种液体氮肥施放方式的确定装置,该装置用于实现上述实施例及优选实施方式,已经进行过说明的不再赘述。如以下所使用的,术语“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
本实施例提供一种液体氮肥施放方式的确定装置,如图5所示,包括:
获取模块501,用于获取目标沙漠地区的垂向分布信息、环境信息、土基材料信息和沙漠沙信息。
第一确定模块502,用于基于垂向分布信息、环境信息和土基材料信息和沙漠沙信息形成室内土柱,并将液体氮肥按照预设施放方式施放至室内土柱中,开展室内土柱试验,得到各预设位置分别对应的氮元素浓度。
第二确定模块503,用于将各预设位置分别对应的氮元素浓度输入至预构建的氮元素迁移模拟模型中,得到与预设施放方式对应的氮元素的迁移过程。
判断模块504,用于基于迁移过程判断预设施放方式是否为目标施放方式。
调整模块505,用于当预设施放方式不是目标施放方式时,调整液体氮肥的施放方式,重新进行试验,直至基于迁移过程确定调整后的施放方式为目标施放方式时,结束试验,并将液体氮肥按照目标施放方式等比例施放于目标沙漠地区中。
在一些可选的实施方式中,第二确定模块503中的预构建的氮元素迁移模拟模型由双阶段的时间分数阶对流扩散反应方程和Hausdroff对流扩散反应方程耦合得到。
在一种可选的实施方式中,第二确定模块503中的氮元素迁移模拟模型通过如下公式表示:
其中,为地表位置,/>为土基材料与沙漠沙接触位置,/>为沙漠沙的下边界,/>为垂向空间位置,/>为Caputo类型的时间分数阶导数,/>中C表示该类型导数为Caputo类分数阶导数,0表示初始时刻,/>表示时间分数阶阶数,代表土基材料整体对氮元素的滞留程度,/>,/>取值越低,滞留效果越强,/>为Hausdroff导数,定义为/>为Hausdroff导数阶数,由于沙漠沙无滞水效果,此处/>1,/>为液体氮肥中氮元素的浓度,/>为氮元素背景浓度值值函数,/>为土基材料的流动参数,/>为沙漠沙的流动参数,/>为土基材料的扩散系数,/>为沙漠沙的扩散系数,/>为硝化速率、/>为反硝化速率。
在一种可选的实施方式中,第二确定模块503,包括:
第一确定子模块,用于将预设位置分别对应的氮元素浓度输入至预构建的氮元素迁移模拟模型中,采用最小二乘法进行反演,得到与预设施放方式对应的迁移系数。
第二确定子模块,用于将迁移系数输入至氮元素迁移模拟模型中,得到与预设施放方式对应的氮元素的迁移过程。
在一种可选的实施方式中,第二确定子模块中的迁移过程为垂直方向上各位置处氮元素的浓度分布。
在一种可选的实施方式中,判断模块504,包括:
第三确定子模块,用于基于迁移过程,确定土柱中土基材料部分对应的氮元素积蓄量。
第四确定子模块,用于基于迁移过程,确定土柱中沙漠沙下边界的氮元素输出量。
第一比较子模块,用于将氮元素积蓄量与第一预设浓度进行比较,得到第一比较结果。
第二比较子模块,用于将氮元素输出量与第二预设浓度进行比较,得到第二比较结果。
判断子模块,用于基于第一比较结果和/或第二比较结果,判断预设施放方式是否为目标施放方式。
在一些可选的实施方式中,判断子模块,包括:
确定单元,用于当第一比较结果为氮元素积蓄量小于第一预设浓度,且,当第二比较结果为氮元素输出量小于第二预设浓度时,将预设施放方式确定为目标施放方式。
上述各个模块和单元的更进一步的功能描述与上述对应实施例相同,在此不再赘述。
本实施例中的液体氮肥施放方式的确定装置是以功能单元的形式来呈现,这里的单元是指ASIC(Application Specific Integrated Circuit,专用集成电路)电路,执行一个或多个软件或固定程序的处理器和存储器,和/或其他可以提供上述功能的器件。
本发明实施例还提供一种计算机设备,具有上述图5所示的液体氮肥施放方式的确定装置。
请参阅图6,图6是本发明可选实施例提供的一种计算机设备的结构示意图,如图6所示,该计算机设备包括:一个或多个处理器10、存储器20,以及用于连接各部件的接口,包括高速接口和低速接口。各个部件利用不同的总线互相通信连接,并且可以被安装在公共主板上或者根据需要以其它方式安装。处理器可以对在计算机设备内执行的指令进行处理,包括存储在存储器中或者存储器上以在外部输入/输出装置(诸如,耦合至接口的显示设备)上显示GUI的图形信息的指令。在一些可选的实施方式中,若需要,可以将多个处理器和/或多条总线与多个存储器和多个存储器一起使用。同样,可以连接多个计算机设备,各个设备提供部分必要的操作(例如,作为服务器阵列、一组刀片式服务器、或者多处理器系统)。图6中以一个处理器10为例。
处理器10可以是中央处理器,网络处理器或其组合。其中,处理器10还可以进一步包括硬件芯片。上述硬件芯片可以是专用集成电路,可编程逻辑器件或其组合。上述可编程逻辑器件可以是复杂可编程逻辑器件,现场可编程逻辑门阵列,通用阵列逻辑或其任意组合。
其中,存储器20存储有可由至少一个处理器10执行的指令,以使至少一个处理器10执行实现上述实施例示出的方法。
存储器20可以包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需要的应用程序;存储数据区可存储根据计算机设备的使用所创建的数据等。此外,存储器20可以包括高速随机存取存储器,还可以包括非瞬时存储器,例如至少一个磁盘存储器件、闪存器件、或其他非瞬时固态存储器件。在一些可选的实施方式中,存储器20可选包括相对于处理器10远程设置的存储器,这些远程存储器可以通过网络连接至该计算机设备。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
存储器20可以包括易失性存储器,例如,随机存取存储器;存储器也可以包括非易失性存储器,例如,快闪存储器,硬盘或固态硬盘;存储器20还可以包括上述种类的存储器的组合。
该计算机设备还包括通信接口30,用于该计算机设备与其他设备或通信网络通信。
本发明实施例还提供了一种计算机可读存储介质,上述根据本发明实施例的方法可在硬件、固件中实现,或者被实现为可记录在存储介质,或者被实现通过网络下载的原始存储在远程存储介质或非暂时机器可读存储介质中并将被存储在本地存储介质中的计算机代码,从而在此描述的方法可被存储在使用通用计算机、专用处理器或者可编程或专用硬件的存储介质上的这样的软件处理。其中,存储介质可为磁碟、光盘、只读存储记忆体、随机存储记忆体、快闪存储器、硬盘或固态硬盘等;进一步地,存储介质还可以包括上述种类的存储器的组合。可以理解,计算机、处理器、微处理器控制器或可编程硬件包括可存储或接收软件或计算机代码的存储组件,当软件或计算机代码被计算机、处理器或硬件访问且执行时,实现上述实施例示出的方法。
虽然结合附图描述了本发明的实施例,但是本领域技术人员可以在不脱离本发明的精神和范围的情况下做出各种修改和变型,这样的修改和变型均落入由所附权利要求所限定的范围之内。

Claims (11)

1.一种液体氮肥施放方式的确定方法,其特征在于,所述方法包括:
获取目标沙漠地区的垂向分布信息、环境信息、土基材料信息和沙漠沙信息;
基于所述垂向分布信息、所述环境信息和所述土基材料信息和所述沙漠沙信息形成室内土柱,并将液体氮肥按照预设施放方式施放至所述室内土柱中,开展室内土柱试验,得到各预设位置分别对应的氮元素浓度;
将各所述预设位置分别对应的氮元素浓度输入至预构建的氮元素迁移模拟模型中,得到与所述预设施放方式对应的氮元素的迁移过程;
基于所述迁移过程判断所述预设施放方式是否为目标施放方式;
当所述预设施放方式不是目标施放方式时,调整所述液体氮肥的施放方式,重新进行试验,直至基于所述迁移过程确定调整后的施放方式为所述目标施放方式时,结束试验,并将所述液体氮肥按照所述目标施放方式等比例施放于所述目标沙漠地区中;
所述预构建的氮元素迁移模拟模型由双阶段的时间分数阶对流扩散反应方程和Hausdroff对流扩散反应方程耦合得到;
所述氮元素迁移模拟模型通过如下公式表示:
其中,为地表位置,/>为土基材料与沙漠沙接触位置,/>为沙漠沙的下边界,/>为垂向空间位置,/>为时间分数阶导数,/>M为时间分数阶导数所属的类型,/>为时间分数阶阶数,/>为豪斯多夫导数,/>为豪斯多夫导数阶数,/>为液体氮肥中氮元素的浓度,/>为氮元素背景浓度值值函数,/>为土基材料的流动参数,/>为沙漠沙的流动参数,/>为土基材料的扩散系数,/>为沙漠沙的扩散系数,/>为硝化速率、/>为反硝化速率。
2.根据权利要求1所述的方法,其特征在于,所述将各所述预设位置分别对应的氮元素浓度输入至预构建的氮元素迁移模拟模型中,得到与所述预设施放方式对应的氮元素的迁移过程,包括:
将所述预设位置分别对应的氮元素浓度输入至预构建的氮元素迁移模拟模型中,采用最小二乘法进行反演,得到与所述预设施放方式对应的迁移系数;
将所述迁移系数输入至所述氮元素迁移模拟模型中,得到与所述预设施放方式对应的氮元素的迁移过程。
3.根据权利要求2所述的方法,其特征在于,所述迁移过程为垂直方向上各位置处氮元素的浓度分布。
4.根据权利要求2或3所述的方法,其特征在于,所述基于所述迁移过程判断所述预设施放方式是否为目标施放方式,包括:
基于所述迁移过程,确定所述室内土柱中土基材料部分对应的氮元素积蓄量;
基于所述迁移过程,确定所述室内土柱中沙漠沙下边界的氮元素输出量;
将所述氮元素积蓄量与第一预设浓度进行比较,得到第一比较结果;
将所述氮元素输出量与第二预设浓度进行比较,得到第二比较结果;
基于所述第一比较结果和/或所述第二比较结果,判断所述预设施放方式是否为所述目标施放方式。
5.根据权利要求4所述的方法,其特征在于,所述基于所述第一比较结果和所述第二比较结果,判断所述预设施放方式是否为所述目标施放方式,包括:
当所述第一比较结果为所述氮元素积蓄量小于所述第一预设浓度,且,当所述第二比较结果为所述氮元素输出量小于所述第二预设浓度时,将所述预设施放方式确定为所述目标施放方式。
6.一种液体氮肥施放方式的确定装置,其特征在于,所述装置包括:
获取模块,用于获取目标沙漠地区的垂向分布信息、环境信息、土基材料信息和沙漠沙信息;
第一确定模块,用于基于所述垂向分布信息、所述环境信息和所述土基材料信息和所述沙漠沙信息形成室内土柱,并将液体氮肥按照预设施放方式施放至所述室内土柱中,开展室内土柱试验,得到各预设位置分别对应的氮元素浓度;
第二确定模块,用于将各所述预设位置分别对应的氮元素浓度输入至预构建的氮元素迁移模拟模型中,得到与所述预设施放方式对应的氮元素的迁移过程;
判断模块,用于基于所述迁移过程判断所述预设施放方式是否为目标施放方式;
调整模块,用于当所述预设施放方式不是目标施放方式时,调整所述液体氮肥的施放方式,重新进行试验,直至基于所述迁移过程确定调整后的施放方式为所述目标施放方式时,结束试验,并将所述液体氮肥按照所述目标施放方式等比例施放于所述目标沙漠地区中;
所述第二确定模块中的所述预构建的氮元素迁移模拟模型由双阶段的时间分数阶对流扩散反应方程和Hausdroff对流扩散反应方程耦合得到;
所述第二确定模块中的所述氮元素迁移模拟模型通过如下公式表示:
其中,为地表位置,/>为土基材料与沙漠沙接触位置,/>为沙漠沙的下边界,/>为垂向空间位置,/>为时间分数阶导数,/>M为时间分数阶导数所属的类型,/>为时间分数阶阶数,/>为豪斯多夫导数,/>为豪斯多夫导数阶数,/>为液体氮肥中氮元素的浓度,/>为氮元素背景浓度值值函数,/>为土基材料的流动参数,/>为沙漠沙的流动参数,/>为土基材料的扩散系数,/>为沙漠沙的扩散系数,/>为硝化速率、/>为反硝化速率。
7.根据权利要求6所述的装置,其特征在于,所述第二确定模块,包括:
第一确定子模块,用于将所述预设位置分别对应的氮元素浓度输入至预构建的氮元素迁移模拟模型中,采用最小二乘法进行反演,得到与所述预设施放方式对应的迁移系数;
第二确定子模块,用于将所述迁移系数输入至所述氮元素迁移模拟模型中,得到与所述预设施放方式对应的氮元素的迁移过程。
8.根据权利要求7所述的装置,其特征在于,所述第二确定子模块中的所述迁移过程为垂直方向上各位置处氮元素的浓度分布。
9.根据权利要求7或8所述的装置,其特征在于,所述判断模块,包括:
第三确定子模块,用于基于所述迁移过程,确定所述室内土柱中土基材料部分对应的氮元素积蓄量;
第四确定子模块,用于基于所述迁移过程,确定所述室内土柱中沙漠沙下边界的氮元素输出量;
第一比较子模块,用于将所述氮元素积蓄量与第一预设浓度进行比较,得到第一比较结果;
第二比较子模块,用于将所述氮元素输出量与第二预设浓度进行比较,得到第二比较结果;
判断子模块,用于基于所述第一比较结果和/或所述第二比较结果,判断所述预设施放方式是否为所述目标施放方式。
10.一种计算机设备,其特征在于,包括:
存储器和处理器,所述存储器和所述处理器之间互相通信连接,所述存储器中存储有计算机指令,所述处理器通过执行所述计算机指令,从而执行权利要求1至5中任一项所述的液体氮肥施放方式的确定方法。
11.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机指令,所述计算机指令用于使计算机执行权利要求1至5中任一项所述的液体氮肥施放方式的确定方法。
CN202311147826.6A 2023-09-07 2023-09-07 一种液体氮肥施放方式的确定方法、装置、设备及介质 Active CN116911071B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311147826.6A CN116911071B (zh) 2023-09-07 2023-09-07 一种液体氮肥施放方式的确定方法、装置、设备及介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311147826.6A CN116911071B (zh) 2023-09-07 2023-09-07 一种液体氮肥施放方式的确定方法、装置、设备及介质

Publications (2)

Publication Number Publication Date
CN116911071A CN116911071A (zh) 2023-10-20
CN116911071B true CN116911071B (zh) 2024-01-26

Family

ID=88358690

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311147826.6A Active CN116911071B (zh) 2023-09-07 2023-09-07 一种液体氮肥施放方式的确定方法、装置、设备及介质

Country Status (1)

Country Link
CN (1) CN116911071B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2020101615A4 (en) * 2020-08-02 2020-09-10 Beijing Normal University A Method for Source Apportionment of PAHs in Roadway Sediments Coupled with Transport and Transformation Process
CN116539482A (zh) * 2023-05-12 2023-08-04 中国长江三峡集团有限公司 河流中磷元素迁移检测方法、装置、电子设备及存储介质
CN116628959A (zh) * 2023-05-12 2023-08-22 中国长江三峡集团有限公司 一种携沙河流污染物迁移检测方法、装置及电子设备

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2020101615A4 (en) * 2020-08-02 2020-09-10 Beijing Normal University A Method for Source Apportionment of PAHs in Roadway Sediments Coupled with Transport and Transformation Process
CN116539482A (zh) * 2023-05-12 2023-08-04 中国长江三峡集团有限公司 河流中磷元素迁移检测方法、装置、电子设备及存储介质
CN116628959A (zh) * 2023-05-12 2023-08-22 中国长江三峡集团有限公司 一种携沙河流污染物迁移检测方法、装置及电子设备

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
关中盆地包气带氮迁移转化数值模拟及预测;姜桂华;王文科;;西北大学学报(自然科学版)(05);全文 *
时间分数阶扩散方程的新型格式;张佳丽;中国优秀硕士学位论文全文数据库基础科学辑(第09期);全文 *

Also Published As

Publication number Publication date
CN116911071A (zh) 2023-10-20

Similar Documents

Publication Publication Date Title
Ritter et al. Using inverse methods for estimating soil hydraulic properties from field data as an alternative to direct methods
CN109451532B (zh) 一种基站位置的核查方法及装置
Ahuja et al. Parameterization of agricultural system models: Current approaches and future needs
CN109448124B (zh) 用于河道的水质模拟方法和装置
Hollaway et al. The challenges of modelling phosphorus in a headwater catchment: Applying a ‘limits of acceptability’uncertainty framework to a water quality model
CN105550784A (zh) 一种空气质量监测站优化布点方法
Cichota et al. Describing N leaching from urine patches deposited at different times of the year with a transfer function
Nilawar et al. Use of SWAT to determine the effects of climate and land use changes on streamflow and sediment concentration in the Purna River basin, India
CN109559098B (zh) 一种海绵城市试点区低影响开发设施模拟方法
Lilburne et al. Computer‐based evaluation of methods to sample nitrate leached from grazed pasture
CN114239278A (zh) 土壤重金属累积过程时空模拟模型构建方法
Miao et al. Simulating pesticide leaching and runoff in rice paddies with the RICEWQ–VADOFT Model
CN116911071B (zh) 一种液体氮肥施放方式的确定方法、装置、设备及介质
Turek et al. Criteria for the estimation of field capacity and their implications for the bucket type model
CN108733952A (zh) 一种基于序贯模拟的土壤含水量空间变异性三维表征方法
Vogeler et al. Effect of variability in soil properties plus model complexity on predicting topsoil water content and nitrous oxide emissions
Grossel et al. Simulating the spatial variability of nitrous oxide emission from cropped soils at the within-field scale using the NOE model
Sogbedji et al. N rate and transport under variable cropping history and fertilizer rate on loamy sand and clay loam soils: II. Performance of LEACHMN using different calibration scenarios
CN109472002A (zh) 污染物浓度分布图的生成方法及装置
CN115164943B (zh) 定位设备的测试方法、装置、电子设备及存储介质
CN116911072B (zh) 透镜体分布占比的确定方法、装置、计算机设备及介质
CN109507367A (zh) 确定大气污染物精细化分布的方法及装置
Xu et al. Nitrogen leaching in a typical agricultural extensively cropped catchment, China: experiments and modelling
CN114048692A (zh) 一种基于监测数据与模型融合的水体月尺度水质预判方法
Bouhachlaf et al. Monitoring soil elements for irrigation management using Internet of Things (IoT) sensors

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant