CN112507596B - 地下水水位的分析方法、装置、电子设备和存储介质 - Google Patents
地下水水位的分析方法、装置、电子设备和存储介质 Download PDFInfo
- Publication number
- CN112507596B CN112507596B CN202011529864.4A CN202011529864A CN112507596B CN 112507596 B CN112507596 B CN 112507596B CN 202011529864 A CN202011529864 A CN 202011529864A CN 112507596 B CN112507596 B CN 112507596B
- Authority
- CN
- China
- Prior art keywords
- groundwater
- model
- movement model
- motion model
- flux
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
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)
- Geophysics And Detection Of Objects (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本申请实施例提供的地下水水位的分析方法、装置、电子设备和存储介质,通过根据多点通量逼近算法对第一地下水运动模型进行数值离散处理,得到第二地下水运动模型,基于第二地下水运动模型,确定待评估区域的地下水水位,实现了通过多点通量逼近算法对第一地下水运动模型进行离散以确定地下水水位,提高了对地下水水位的分析精度。
Description
技术领域
本申请实施例涉及地下水技术领域,尤其涉及一种地下水水位的分析方法、装置、电子设备和存储介质。
背景技术
随着地下水科学和计算机科学的发展,地下水数值模拟技术得到了快速发展和广泛应用,利用数学模型对地下水流等问题进行模拟成为地下水研究领域最有效、应用最广泛的方法之一。
由于绝大部分数学模型是无法用解析法求解的,通过数值化就是将数学模型转化为可解的数值模型是求解数学模型的常用手段。现有技术中,通过采用有限单元法(finiteelement method,FEM)或者有限差分法(finite difference method,FDM)对的地下水水流的数学模型进行离散,得到地下水水流的数值模型,从而得到数值解,达到分析地下水的目的。
然而,由于有限单元法和有限差分法的网格剖分主要以三角形、四边形或者两种的混合为主,而实际地下水的边界情况通常都比较复杂,因此,采用现有技术分析地下水时存在准确度不高的问题。
发明内容
本申请提供一种地下水水位的分析方法、装置、电子设备和存储介质,以解决现有技术中存在的准确度不高的问题。
第一方面,本申请实施例提供一种地下水水位的分析方法,包括:
根据多点通量逼近算法对第一地下水运动模型进行数值化处理,得到第二地下水运动模型;
基于所述第二地下水运动模型,确定待评估区域的地下水水位。
可选地,所述根据多点通量逼近算法对第一地下水运动模型进行数值化处理,得到第二地下水运动模型,包括:
对待评估区域进行网格化处理,得到待评估区域的二维平面网格;
根据多点通量逼近算法,确定所述二维平面网格中每个单元网格的通量;
将每个单元网格的通量代入所述第一地下水运动模型中,得到所述第二地下水运动模型。
可选地,所述基于所述第二地下水运动模型,确定待评估区域的地下水水位,包括:
获取所述待评估区域地下水流场的特征参数;
根据所述特征参数和所述第二地下水运动模型,确定所述待评估区域的地下水水位。
可选地,所述根据多点通量逼近算法对第一地下水运动模型进行数值化处理,得到第二地下水运动模型之前,所述方法还包括:
对地下水二维流运动模型进行预处理,得到所述第一地下水运动模型。
可选地,所述对地下水二维流运动模型进行预处理,得到所述第一地下水运动模型,包括:
根据高斯定理对所述地下水二维流运动模型进行变换处理,得到控制体格式的地下水二维流运动模型;
对所述控制体格式的地下水二维流运动模型中的时间项进行离散,得到所述第一地下水运动模型。
可选地,所述根据所述特征参数和所述第二地下水运动模型,确定所述待评估区域的地下水水位,包括:
对所述第二地下水运动模型进行矩阵化处理,得到矩阵格式的第二地下水运动模型;
根据所述特征参数求解所述矩阵格式的第二地下水运动模型,得到所述待评估区域的地下水水位。
第二方面,本申请实施例提供一种地下水水位的分析装置,包括:
第一处理模块,用于根据多点通量逼近算法对第一地下水运动模型进行数值化处理,得到第二地下水运动模型;
第二处理模块,用于基于所述第二地下水运动模型,确定待评估区域的地下水水位。
可选地,所述第一处理模块具体用于:
对待评估区域进行网格化处理,得到待评估区域的二维平面网格;
根据多点通量逼近算法,确定所述二维平面网格中每个单元网格的通量;
将每个单元网格的通量代入所述第一地下水运动模型中,得到所述第二地下水运动模型。
第三方面,本申请实施例提供一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如上述第一方面所述的地下水水位的分析方法。
第四方面,本申请实施例提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现如上述第一方面所述的地下水水位的分析方法。
本申请实施例提供的地下水水位的分析方法、装置、电子设备和存储介质,通过根据多点通量逼近算法对第一地下水运动模型进行数值离散处理,得到第二地下水运动模型,基于第二地下水运动模型,确定待评估区域的地下水水位,由于多点通量逼近算法在对第一地下水运动模型数据值化的过程中,根据待评估区域地下水流场的边界等实际情况,进行任意多边形的剖分,因此,采用本申请技术方案的方法分析得到的地下水水位更能反映待评估区域地下水的实际情况,从而提高了对地下水水位的分析精度。
附图说明
图1为本申请实施例一提供的地下水水位的分析方法的流程示意图;
图2为本申请实施例一提供二维平面网格的结构编号示意图;
图3为本申请实施例一提供的利用多点通量逼近算法求解节点O0的通量的结构编号示意图;
图4为本申请实施例二提供的地下水水位的分析装置的结构示意图;
图5为本申请实施例三提供的一种电子设备的结构示意图。
具体实施方式
下面结合附图和实施例对本申请作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释本申请,而非对本申请的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本申请相关的部分而非全部结构。
现有技术中在求解地下水数学模型的过程中,采用有限单元法或者有限差分法进行离散,如常用的GMS、VisualMODFLOW、VisualGroundwater等地下水模拟软件均是采用这种方法,而这些方法在进行网格剖分时主要以三角形、四边形或者两种的混合为主,而实际地下水的边界情况通常都比较复杂,因此,采用这些方法进行地下水的分析和评估时存在准确度不高的问题。如何实现任意多边形的网格剖分成为我国地下水数值模拟水平提升过程中有待突破的瓶颈,本申请发明人在长期进行地下水研究的过程中发现,多点通量逼近(multi-point flux approximations,MPFA)算法是一种只关注过节点的边的数量,而对网格单元的边数不做具体要求的离散方法,因此,适用于不规则扩展网格单元,通过将多点通量逼近算法应用到地下水数学模型的求解过程中,可以解决现有技术中存在的技术问题。
本申请技术方案的主要思路:基于现有技术对地下水水流分析时存在准确度不高的技术问题,本申请提供一种地下水水位的分析方法,通过多点通量逼近算法对地下水水流的数学模型进行离散,并对离散后地下水水流的数值模型进行求解,得到待评估区域的地下水水位,由于本申请的技术方案在离散过程中可以根据待评估区域地下水的实际流场和边界等情况,进行任意多边形的剖分,因此,得到的地下水水位更能反映待评估区域地下水的实际情况,提高了了对地下水水位的分析精度,提高了我国地下水的模拟水平。
本申请实施例提供一种分析地下水水位的技术方案,适用于对潜水含水层及承压含水层中的地下水水位进行分析。下面将以潜水含水层中地下水水位的分析为例,对本申请的技术方案加以描述。
实施例一
图1为本申请实施例一提供的地下水水位的分析方法的流程示意图,该方法可以由具有地下水水位分析功能的电子设备如终端设备、服务器等来执行,如图1所示,示例性地,本实施例的地下水水位的分析方法,包括:
S110、对地下水二维流运动模型进行预处理,得到第一地下水运动模型。
本步骤中,为了将多点通量逼近算法应用到潜水含水层的地下水二维流运动模型的求解过程中,需要先对地下水二维流运动模型进行适当的变形处理。具体地:
t时刻,潜水面上任意位置(x,y)处的渗流规律,可以用地下水二维流运动模型表示如下:
其中,Sy为给水度;K的水力传导系数;WM为源汇项;h为压力水头;H为总水头,H=h+z(z为位置水头)。
对公式(1)进行积分,并利用高斯定理得到反映一个控制单元的地下水渗流规律的控制体格式的地下水二维流运动模型为:
其中,Sa为控制单元的面积;Sk为控制单元的边,k为角标,表示边的序号;▽H为水力梯度,n为边的外法线向量。右边第一项表示一个控制单元的所有边的通量代数和,用fk表示第k条边的通量,则公式(2)变为:
进一步地,对公式(3)中的时间项进行离散,得第一地下水运动模型,表示如下:
公式(4)中m为时间步。
需要说明的是,第一地下水运动模型,即公式(4),是通过对公式(3)采用全隐式的离散方法得到的。
本步骤中,通过对地下水二维流运动模型进行预处理,引入了通量的概念,由公式(4)可知,对于某一个控制单元来说,通过计算其所有边的通量代数和就可以确定该位置的水位。
S120、根据多点通量逼近算法对第一地下水运动模型进行数值化处理,得到第二地下水运动模型。
为了实现对第一地下水运动模型的求解,本步骤中,通过多点通量逼近算法对第一地下水运动模型进行离散,即数值化处理,从而得到第二地下水运动模型。具体地,本步骤中,通过对待评估区域进行网格化处理,得到待评估区域的二维平面网格,根据多点通量逼近算法,确定二维平面网格中每个单元网格的通量,将每个单元网格的通量代入第一地下水运动模型中,得到第二地下水运动模型。
其中,网格化处理是指对待评估区域进行网格剖分,由于本申请中采用多点通量逼近算法对第一地下水运动模型进行离散,因此,在剖分过程中可以根据水流场的实际分布和边界情况等,采用任意边数的多边形,不必局限于现有技术中的三角形或四边形。
二维平面网格由多个网格单元组成,一个网格单元即为一个控制单元,每个网格单元由多个节点和多条边构成。下面将引入网格单元的通量和节点的通量的两个概念。
网格单元的通量是指过该网格单元的所有边的通量的代数和。
节点的通量是指过该节点的所有半边的通量的代数和。
由于网格单元、节点和边之间的关系,申请实施例中,通过确定每个网格单元包含的节点的通量,确定每个网格单元的通量。
本实施例中,在进行待评估区域的地下水水位之前可以由测试人员根据待评估区的边界和形状等编写相关的文本文件,并将其输入到进行地下水水位分析的电子设备中,当需要进行地下水水位分析时,电子设备在获取相应的文本文件进行网格剖分,从而得到能够反应地下水实际流场情况的二维平面网格。示例性地,表1示出了本实施例提供的一种进行网格剖分的文本文件,其中可以包括网格单元的单元编号、边起点坐标、边终点坐标、边界条件等。
表1
示例性,图2为本申请实施例一提供二维平面网格的结构编号示意图,以节点O0为例,如图2所示,过节点O0的边一共有五条,分别O0P1、O0P2、O0P3、O0P4和O0P5,O0为所有边的初端编号,P1-P5为边的末端编号,T1、T2、T3、T4、T5分别为对应边的中点,O1、O2、O3、O4、O5分别为对应的网格单元的单元编号。在图2中,网格单元O1由四条边(O0P1、O0P2、P1M1、P2M1)和四个节点(O0、P1、P2、M1)构成,网格单元O3由五条边(O0P4、O0P3、P3M2、P4M3和M2M3)和五个节点(O0、P3、P4、M2、M3)构成。
可知,在图2中节点O0的通量即为过半边O0T1、半边O0T2、半边O0T3、半边O0T4和半边O0T5的通量之和,因此,为得到节点O0的通量就需要先分别确定半边O0T1、O0T2、O0T3、O0T4和O0T5的通量。
下面将详细介绍如何通过多点通量逼近算法,确定二维平面网格中每个节点的通量。
示例性地,图3为本申请实施例一提供的利用多点通量逼近算法求解节点O0的通量的结构编号示意图,在图3中,同样地,O0为边O0Pk、O0Pk+1的初端编号,Pk、Pk+1分别为边O0Pk、O0Pk+1的末端编号,Tk为O0Pk边的中点,Tk+1为O0Pk+1边的中点,Ok为网格单元编号,为Tk处的水头,/>为Tk+1处的水头,Hk为网格单元k的中心水头,/>为O0Tk边的外法线向量,/>为O0Tk+1边的外法线向量。图3中的阴影部分为第k个单元网格的一部分。
假设网格单元k中总水头H都是线性分布的,则根据地下水二维流运动模型可知,网格单元k的中心水头的梯度可以表示为:
其中,Ak为网格单元k中阴影的面积;和/>分别为点Tk,Tk+1和O0的横坐标,R是旋转矩阵,/>
半边O0Tk,O0Tk+1的通量可分别用如下公式表示:
公式(6)中hk为网格单元k的压力水头。
将公式(5)中的▽Hk代入公式(6)中,并将其改写成矩阵的形式,得:
其中,为/>和/>的系数,用通用公式表示为:
由半边法向流的局部守恒性,可知:
对于公式(9),将Hk、Hk-1移到等式右边,并将其写成矩阵的形式得:
其中,为边的总水头矩阵,表示为/>H为中心的总水头矩阵,表示为/>M为边的总水头的系数矩阵,表示为/>N为中心的总水头的系数矩阵,表示为/>M和N可以由下式给出:
由公式(10)可知,网格单元k边上的总水头可以中心水头表示为:
将公式(12)代入公式(7)中,可得只用中心的总水头表示的第k条边的半边通量:
公式(13)表示过节点O0的一条边的半边通量,假设过节点O0的边有n0条,则过节点O0的所有半边通量用矩阵表示为:
f'=CH (14)
其中,f'为节点O0的通量矩阵,表示为C为局部矩阵,表示为
通过剖分得到的二维平面网格中有多个节点,根据式(14)可计算得到过每个节点的半边通量,进而根据节点、边、半边和网格单元之间的数学关系,通过一定的运算,就可以可得到每个网格单元的所有边的通量。
假设剖分得到的二维平面网格中有q个网格单元,分别将将每个网格单元的通量代入第一地下水运动模型(公式(4))中,就可以得到第一地下水运动模型的离散格式,也就是第二地下水运动模型。
第二地下水运动模型用公式表示为:
其中,表示第q个网格单元的给水度,/>表示第q个网格单元的面积,Hq,m表示第q个网格m时刻的总水头,Hq,m-1表示第q个网格单元m-1时刻的总水头,/>表示第q个网格单元中第k条边m时刻的通量,/>表示第q个网格单元的源汇项。
本步骤中,通过多通量逼近算法实现了对第一地下水运动模型的离散,为后续确定待评估区域的地下水水位奠定了基础。
S130、基于第二地下水运动模型,确定待评估区域的地下水水位。
本步骤中,根据S120中得到的第二地下水运动模型和待评估区域的特性就可以确定待评估区域的地下水水位。具体地,本步骤中可以通过获取待评估区域地下水流场的特征参数,根据特征参数和第二地下水运动模型,确定待评估区域的地下水水位。
其中,特征参数用于反映待评估区域的地下水流场特性,本实施例中,获取的特征参数可以包括如下参数至少一项:含水层厚度、水力传导系数、给水度、初始水头和源汇项。
一种可能的实现方式中,本实施例中,在进行待评估区域的地下水水位之前可以由测试人员采集或收集反映待评估区域水文特征的相关特征参数,生成文本文件,并将其输入到进行地下水水位分析的电子设备中,当需要进行地下水水位分析时,电子设备获取该文本文件,并根据文本文件中的具体参数求解第二地下水运动模型,从而确定待评估区域的地下水水位。示例性地,表2示出了本实施例提供的一种记录特征参数的文本文件。
表2
可选地,为便于进行第二地下水运动模型的求解,可以先对第二地下水运动模型进行矩阵化处理(将公式(15)表示成矩阵的形式),得到矩阵格式的第二地下水运动模型,再根据特征参数求解矩阵格式的第二地下水运动模型,得到待评估区域的地下水水位。
矩阵格式的第二地下水运动模型用公式表示为:
DH=B (16)
其中,D为总体矩阵,H为待求的水头(即水位),B为右向量。
为提高第二地下水运动模型的求解精确度,本实施例中可以通过迭代算法求解矩阵格式的第二地下水运动模型,可选地,本实施中可以采用迭代算法有:皮卡逐次逼近法(Picard)或牛顿-拉夫森(Newtonz-Raphson)迭代法。
本实施例中,通过根据多点通量逼近算法对第一地下水运动模型进行数值化处理,得到第二地下水运动模型,基于第二地下水运动模型,确定待评估区域的地下水水位,由于多点通量逼近算法在对第一地下水运动模型数据值化的过程中,根据待评估区域地下水流场的边界等实际情况,进行任意多边形的剖分,因此,采用本申请技术方案的方法分析得到的地下水水位更能反映待评估区域地下水的实际情况,从而提高了对地下水水位的分析精度。
可以理解的是,对于承压含水层来说,只要将相应的地下水水流模型替换为承压含水层中的地下水水流模型,采用相似的方法分析和求解,就可以实现对承压含水层中的地下水水位的评估。
实施例二
图4为本申请实施例二提供的地下水水位的分析装置的结构示意图,如图4所示,本实施例中地下水水位的分析装置10包括:
第一处理模块11和第二处理模块12。
第一处理模块11,用于根据多点通量逼近算法对第一地下水运动模型进行数值化处理,得到第二地下水运动模型;
第二处理模块12,用于基于所述第二地下水运动模型,确定待评估区域的地下水水位。
可选地,第一处理模块11具体用于:
对待评估区域进行网格化处理,得到待评估区域的二维平面网格;
根据多点通量逼近算法,确定所述二维平面网格中每个网格单元的通量;
将每个网格单元的通量代入所述第一地下水运动模型中,得到所述第二地下水运动模型。
可选地,第二处理模块12具体用于:
获取所述待评估区域地下水流场的特征参数;
根据所述特征参数和所述第二地下水运动模型,确定所述待评估区域的地下水水位。
可选地,第一处理模块11还用于:
对地下水二维流运动模型进行预处理,得到所述第一地下水运动模型。
可选地,第一处理模块11具体用于:
根据高斯定理对所述地下水二维流运动模型进行变换处理,得到控制体格式的地下水二维流运动模型;
对所述控制体格式的地下水二维流运动模型中的时间项进行离散,得到所述第一地下水运动模型。
可选地,第二处理模块12具体用于:
对所述第二地下水运动模型进行矩阵化处理,得到矩阵格式的第二地下水运动模型;
根据所述特征参数求解所述矩阵格式的第二地下水运动模型,得到所述待评估区域的地下水水位。
本申请实施例所提供的地下水水位的分析装置前述方法实施例所提供的地下水水位的分析方法,具备执行方法相应的功能模块和有益效果。本实施例的实现原理和技术效果与上述方法实施例类似,此处不再一一赘述。
实施例三
图5为本申请实施例三提供的一种电子设备的结构示意图,如图5所示,该电子设备20包括存储器21、处理器22及存储在存储器上并可在处理器上运行的计算机程序;电子设备20处理器22的数量可以是一个或多个,图5中以一个处理器22为例;电子设备20中的的处理器22、存储器21可以通过总线或其他方式连接,图5中以通过总线连接为例。
存储器21作为一种计算机可读存储介质,可用于存储软件程序、计算机可执行程序以及模块,如本申请实施例中的第一处理模块11和第二处理模块12对应的程序指令/模块。处理器22通过运行存储在存储器21中的软件程序、指令以及模块,从而执行设备/终端/服务器的各种功能应用以及数据处理,即实现上述的地下水水位的分析方法。
存储器21可主要包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需的应用程序;存储数据区可存储根据终端的使用所创建的数据等。此外,存储器21可以包括高速随机存取存储器,还可以包括非易失性存储器,例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中,存储器21可进一步包括相对于处理器22远程设置的存储器,这些远程存储器可以通过网格连接至设备/终端/服务器。上述网格的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
实施例四
本申请实施例四还提供一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序在由计算机处理器执行时用于执行一种地下水水位的分析方法,该方法包括:
根据多点通量逼近算法对第一地下水运动模型进行数值化处理,得到第二地下水运动模型;
基于所述第二地下水运动模型,确定待评估区域的地下水水位。
当然,本申请实施例所提供的一种包计算机可读存储介质,其计算机程序不限于如上所述的方法操作,还可以执行本申请任意实施例所提供的地下水水位的分析方法中的相关操作。
通过以上关于实施方式的描述,所属领域的技术人员可以清楚地了解到,本申请可借助软件及必需的通用硬件来实现,当然也可以通过硬件实现,但很多情况下前者是更佳的实施方式。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如计算机的软盘、只读存储器(Read-Only Memory,ROM)、随机存取存储器(RandomAccess Memory,RAM)、闪存(FLASH)、硬盘或光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网格设备等)执行本申请各个实施例所述的方法。
值得注意的是,上述地下水水位的分析装置的实施例中,所包括的各个单元和模块只是按照功能逻辑进行划分的,但并不局限于上述的划分,只要能够实现相应的功能即可;另外,各功能单元的具体名称也只是为了便于相互区分,并不用于限制本申请的保护范围。
注意,上述仅为本申请的较佳实施例及所运用技术原理。本领域技术人员会理解,本申请不限于这里所述的特定实施例,对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本申请的保护范围。因此,虽然通过以上实施例对本申请进行了较为详细的说明,但是本申请不仅仅限于以上实施例,在不脱离本申请构思的情况下,还可以包括更多其他等效实施例,而本申请的范围由所附的权利要求范围决定。
Claims (8)
1.一种地下水水位的分析方法,其特征在于,包括:
根据多点通量逼近算法对第一地下水运动模型进行数值化处理,得到第二地下水运动模型;
基于所述第二地下水运动模型,确定待评估区域的地下水水位;
所述根据多点通量逼近算法对第一地下水运动模型进行数值化处理,得到第二地下水运动模型之前,所述方法还包括:
对地下水二维流运动模型进行预处理,得到所述第一地下水运动模型;
所述第一地下水运动模型,表示如下:
其中,/>为控制单元的面积,Sy为给水度,H为总水头,m为时间步;k为角标,表示边的序号;fk表示第k条边的通量,WM为源汇项;
所述对地下水二维流运动模型进行预处理,得到所述第一地下水运动模型,包括:
根据高斯定理对所述地下水二维流运动模型进行变换处理,得到反映一个控制单元的地下水渗流规律的控制体格式的地下水二维流运动模型;
对所述反映一个控制单元的地下水渗流规律的控制体格式的地下水二维流运动模型中的时间项进行离散,得到所述第一地下水运动模型。
2.根据权利要求1所述的方法,其特征在于,所述根据多点通量逼近算法对第一地下水运动模型进行数值化处理,得到第二地下水运动模型,包括:
对待评估区域进行网格化处理,得到待评估区域的二维平面网格;
根据多点通量逼近算法,确定所述二维平面网格中每个单元网格的通量;
将单元网格的通量代入所述第一地下水运动模型中,得到所述第二地下水运动模型。
3.根据权利要求1所述的方法,其特征在于,所述基于所述第二地下水运动模型,确定待评估区域的地下水水位,包括:
获取所述待评估区域地下水流场的特征参数;
根据所述特征参数和所述第二地下水运动模型,确定所述待评估区域的地下水水位。
4.根据权利要求3所述的方法,其特征在于,所述根据所述特征参数和所述第二地下水运动模型,确定所述待评估区域的地下水水位,包括:
对所述第二地下水运动模型进行矩阵化处理,得到矩阵格式的第二地下水运动模型;
根据所述特征参数求解所述矩阵格式的第二地下水运动模型,得到所述待评估区域的地下水水位。
5.一种地下水水位的分析装置,其特征在于,包括:
第一处理模块,用于根据多点通量逼近算法对第一地下水运动模型进行数值化处理,得到第二地下水运动模型;
第二处理模块,用于基于所述第二地下水运动模型,确定待评估区域的地下水水位;
所述第一处理模块还用于对地下水二维流运动模型进行预处理,得到所述第一地下水运动模型;
所述第一地下水运动模型,表示如下:
其中,/>为控制单元的面积,Sy为给水度,H为总水头,m为时间步;k为角标,表示边的序号;fk表示第k条边的通量,WM为源汇项;
所述第一处理模块还具体用于:
根据高斯定理对所述地下水二维流运动模型进行变换处理,得到反映一个控制单元的地下水渗流规律的控制体格式的地下水二维流运动模型;
对所述反映一个控制单元的地下水渗流规律的控制体格式的地下水二维流运动模型中的时间项进行离散,得到所述第一地下水运动模型。
6.根据权利要求5所述的装置,其特征在于,所述第一处理模块具体用于:
对待评估区域进行网格化处理,得到待评估区域的二维平面网格;
根据多点通量逼近算法,确定所述二维平面网格中每个单元网格的通量;
将每个单元网格的通量代入所述第一地下水运动模型中,得到所述第二地下水运动模型。
7.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现如权利要求1-4中任一所述的地下水水位的分析方法。
8.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现如权利要求1-4中任一所述的地下水水位的分析方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011529864.4A CN112507596B (zh) | 2020-12-22 | 2020-12-22 | 地下水水位的分析方法、装置、电子设备和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011529864.4A CN112507596B (zh) | 2020-12-22 | 2020-12-22 | 地下水水位的分析方法、装置、电子设备和存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112507596A CN112507596A (zh) | 2021-03-16 |
CN112507596B true CN112507596B (zh) | 2023-09-29 |
Family
ID=74922915
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011529864.4A Active CN112507596B (zh) | 2020-12-22 | 2020-12-22 | 地下水水位的分析方法、装置、电子设备和存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112507596B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114757093A (zh) * | 2022-03-28 | 2022-07-15 | 中国水利水电科学研究院 | 一种地下水基准水位确定方法和装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109299428A (zh) * | 2018-10-10 | 2019-02-01 | 深圳市沃泰凯尔科技开发有限公司 | 应用运动波的水流计算方法及系统 |
CN111723505A (zh) * | 2020-06-19 | 2020-09-29 | 中国水利水电科学研究院 | 一种流域水质水量监测系统 |
-
2020
- 2020-12-22 CN CN202011529864.4A patent/CN112507596B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109299428A (zh) * | 2018-10-10 | 2019-02-01 | 深圳市沃泰凯尔科技开发有限公司 | 应用运动波的水流计算方法及系统 |
CN111723505A (zh) * | 2020-06-19 | 2020-09-29 | 中国水利水电科学研究院 | 一种流域水质水量监测系统 |
Non-Patent Citations (2)
Title |
---|
Anis Younes等.Monotonicity of the cell-centred triangular MPFA method for saturated and unsaturated flow in heterogeneous porous media.Journal of Hydrology.2013,第504卷第132-141页. * |
陈智伟等.利用新版GRACE时变模型反演珠江流域陆地水储量变化.大地测量与地球动力学.2020,第40卷(第3期),第305-310页. * |
Also Published As
Publication number | Publication date |
---|---|
CN112507596A (zh) | 2021-03-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Marco et al. | Exact 3D boundary representation in finite element analysis based on Cartesian grids independent of the geometry | |
CN108959328B (zh) | 知识图谱的处理方法、装置及电子设备 | |
Kim et al. | Isogeometric analysis with trimming technique for problems of arbitrary complex topology | |
Piggott et al. | Anisotropic mesh adaptivity for multi-scale ocean modelling | |
Lacasta et al. | An optimized GPU implementation of a 2D free surface simulation model on unstructured meshes | |
Guarguaglini et al. | Fast reaction limit and large time behavior of solutions to a nonlinear model of sulphation phenomena | |
Bazyar et al. | Transient seepage analysis in zoned anisotropic soils based on the scaled boundary finite‐element method | |
Boddula et al. | A moving least squares based meshless local petrov-galerkin method for the simulation of contaminant transport in porous media | |
CN106934093A (zh) | 模拟三维地下水流运动的三重网格多尺度有限元方法 | |
Karahan et al. | Transient groundwater modeling using spreadsheets | |
CN110083853B (zh) | 模拟地下水流运动的有限体积Yeh多尺度有限元法 | |
Ackerman et al. | Computational modeling for climate change: Simulating and visualizing a resilient landscape architecture design approach | |
CN112507596B (zh) | 地下水水位的分析方法、装置、电子设备和存储介质 | |
Veneziano et al. | Scaling and fractals in hydrology | |
Lu et al. | A reduced-order fast reproducing kernel collocation method for nonlocal models with inhomogeneous volume constraints | |
Behrens et al. | Efficiency considerations in triangular adaptive mesh refinement | |
Of et al. | On the numerical solution of minimal energy problems | |
Karahan et al. | Time‐dependent groundwater modeling using spreadsheet | |
Zhang et al. | A graphics processing unit-accelerated meshless method for two-dimensional compressible flows | |
Hwang | Development of a compact and accurate discretization for incompressible navier-stokes equations based on an equation-solving solution gradient, part I: kernel scheme fundamentals | |
Raghavan et al. | Implicit constraint handling for shape optimisation with pod-morphing | |
Shin et al. | C0-discontinuous Galerkin methods for a wind-driven ocean circulation model: Two-grid algorithm | |
Morales et al. | New additional conditions for the numerical uniqueness of the Boussinesq and Timpe solutions of elasticity problems | |
Ruffin et al. | A normal ray refinement technique for Cartesian-grid based Navier–Stokes solvers | |
Oner et al. | A parallel adaptive viscoelastic flow solver with template based dynamic mesh refinement |
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 |