CN112733386B - 一种基于深度神经网络的多组分反应性溶质运移参数反演方法 - Google Patents
一种基于深度神经网络的多组分反应性溶质运移参数反演方法 Download PDFInfo
- Publication number
- CN112733386B CN112733386B CN202110100068.7A CN202110100068A CN112733386B CN 112733386 B CN112733386 B CN 112733386B CN 202110100068 A CN202110100068 A CN 202110100068A CN 112733386 B CN112733386 B CN 112733386B
- Authority
- CN
- China
- Prior art keywords
- neural network
- model
- parameter
- component reactive
- solute transport
- 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
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Data Mining & Analysis (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Computational Linguistics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种基于深度神经网络的多组分反应性溶质运移参数反演方法,该方法的包括以下步骤:根据观测资料,建立多组分反应性溶质运移模型,确定待反演参数和参数取值范围;准备训练样本集;数据归一化处理;训练正演神经网络;训练反演神经网络;利用反演神经网络估计模型参数;计算参数估计点处的正演神经网络误差;通过自适应更新正演神经网络,提升其局部精度,逐步降低反演结果的误差。本发明可以解决多组分反应性溶质运移参数反演问题,为多组分反应性溶质运移模拟技术的应用提供技术保障。
Description
技术领域
本发明涉及一种基于深度神经网络的多组分反应性溶质运移参数反演方法,属于地下水数值模拟反问题研究。
背景技术
多组分反应性溶质运移模拟是分析多孔介质水流系统及化学组分迁移规律的重要工具,在地下水污染、二氧化碳地质储存与核废料埋藏安全评估等方面有着广泛的研究应用。
在用数值模型解决实际问题的研究过程中,准确设定模型参数是数值模型获得可靠模拟结果的前提。在多组分反应性溶质运移模拟问题中,许多参数难以通过目前的测量手段直接获取。针对这一问题,利用实际观测数据,结合数据同化方法的反演计算是目前主要的解决方案。
发明内容
本发明要解决的技术问题是提供一种基于深度神经网络的多组分反应性溶质运移参数反演方法,可以用于估计和设置多组分反应性溶质运移模型中的关键参数,从而解决多组分反应性溶质模拟过程中的参数识别问题。
一种基于深度神经网络的多组分反应性溶质运移参数反演方法,该方法包括以下步骤:
步骤1:通过资料收集或实地采集等方式,得到各化学组分的实际浓度观测数据;根据研究对象实际情况,构建实际观测数据所对应的多组分反应性溶质运移模型;同时确定模型待反演参数和各参数的取值范围;
步骤2:根据模型参数取值范围,随机采样若干组参数样本,逐一代入步骤1中所建立的多组分反应性溶质运移模型,得到每组样本对应的模型响应,形成由“模型参数—模型响应”样本对构成的训练样本集;
步骤3:对步骤1中的实际观测数据和步骤2中的训练样本集作归一化处理;
步骤4:构建一个由模型参数作为输入,模型响应作为输出的正演神经网络;用步骤3中归一化处理后的训练样本集对本步骤中的正演神经网络进行训练;
步骤5:构建一个由模型响应观测数据作为输入,模型参数作为输出的反演神经网络;将该反演神经网络的输出作为步骤4中的正演神经网络的输入,得到一个串联式连接的神经网络架构;在此架构下用步骤3中归一化后的实际观测数据对反演神经网络进行训练;
步骤6:将步骤3中的归一化后的实际观测数据输入步骤5中的反演神经网络,网络输出得到归一化后的模型参数估计结果;将其反归一化,即为步骤1中所建立的待研究对象的多组分反应性溶质运移模型参数估计值;
步骤7:计算步骤6得到参数估计值处的正演神经网络相对于多组分反应性溶质运移模型的误差;
步骤8:判断步骤7中误差是否小于容许误差;如果是,则认为步骤6中所得到的模型参数估计值即为所求的反演结果;否则,执行后续步骤;
步骤9:如果步骤7中正演神经网络在参数估计点处的误差大于容许误差,则在步骤6中所得到的模型参数估计值附近增加n个参数样本,并逐一代入步骤1中所建立的地下水多组分反应性溶质运移模型,得到对应的模型响应,得到局部样本集;计算步骤2中的训练样本集中的所有样本点的参数与步骤6得到的参数估计值之间的欧氏距离,将欧式距离最大的n个样本删除,并加入局部样本集得到更新后的训练样本集;同时对新的样本集做归一化处理;
步骤10:利用步骤9中得到的归一化处理后的新训练样本集,用迁移学习的方法对步骤4中的正演神经网络进行更新,得到更新后的正演神经网络;再执行步骤5~步骤8中的操作;
步骤11:如果步骤8中的判断条件依然不满足,则继续执行步骤9~步骤10,直到循环次数达到最大终止条件,输出所有循环中近似误差最小的那一次参数估计结果作为最终的反演结果。
所述步骤1中多组分反应性溶质运移模拟程序为TOUGHREACT。
所述步骤1中的模型参数一般包括渗透率,阳离子交换系数,各化学组分初始浓度,各化学组分边界浓度和阳离子交换容量等。
所述步骤2中随机采样通过拉丁超立方采样方法实现。
所述步骤2中通过编写程序,实现TOUGHREACT程序的输入文件中模型参数改写,TOUGHREACT程序可执行程序的调用和输出文件中结果的读取,从而得到由模型参数和模型响应构成的训练样本集。
所述步骤3中的数据归一化通过“0-1归一化”方法进行处理;其中,参数归一化的方法如下:
模型响应的归一化方法如下:
所述步骤4和步骤5中的正演神经网络和反演神经网络均通过基于二维卷积的神经网络进行构建;神经网络包括一个前处理模块,三个卷积特征提取运算模块和一个全连接模块依次连接而成;其中,前处理模块将向量形式的输入数据映射为一个32×32×1的矩阵数据,每个卷积模块由一个卷积层、一个批标准化层和一个池化层组成,全连接模块由一个展平层和一个前连接层组成,连接各层的激活函数为Swish函数:
所述步骤4和步骤5中的正演神经网络和反演神经网络的搭建及训练均在pytorch框架下实现;
所述步骤4中正演神经网络训练的损失函数为基于L1范数的目标函数:
公式(5)中θForward表示正演神经网络中各层之间的连接权重参数;FForward表示正演神经网络;
所述步骤5中反演神经网络训练的损失函数根据贝叶斯定理推导而来,具体如下:
公式(6)中,θInverse表示反演神经网络中各层之间的连接权重参数;FInverse表示反演神经网络;表示正演神经网络输出向量中的第i个元素;yobsNorm[i]表示yobsNorm向量中第[i]个元素;
所述步骤6中的参数反归一化处理通过以下公式实现:
所述步骤7中,计算步骤6得到参数估计值处的正演神经网络相对于多组分反应性溶质运移模型的误差为基于归一化数据的均方根误差(RMSE),具体计算方法如下:
则RMSE通过如下公式计算得到:
所述步骤8中的容许近似误差根据不同的研究对象和求解精度需求通过人为设定。
所述步骤10中利用步骤9中得到的更新后的训练样本集,仅对步骤4中的正演神经网络的输出层的神经网络连接权重参数进行更新训练,对其余连接权重参数固定,从而得到更新后的正演神经网络。
所述步骤11中的最大循环次数通过人为设定。
本发明的有益效果:
1、本发明降低了反演结果的误差。
2、本发明解决了多组分反应性溶质运移参数反演问题及模型校正问题,为多组分反应性溶质运移模拟技术的应用提供技术保障。
附图说明
图1为步骤4中的由正演神经网络和反演神经网络构成的串联式架构。
图2为构成正演神经网络和反演神经网络的预处理模块,卷积特征提取运算模块和全连接模块的示意图。
图3为实施案例中各离子的模拟计算结果(a),加入噪声扰动后的结果(b)。
图4参数反演结果对应的拟合图。
具体实施方式
实施例:
本实施例是一个假想的多孔介质中反应性溶质运移模拟柱实验模型,长为0.08m。初始状态溶液中的溶质主要为1mmol/LNaNO3和0.2mmol/LKNO3。该柱子一端以流速0.1m/h注入0.6mmol/L的CaCl2溶液。多孔介质的孔隙度为0.3,渗透率为7.0×10-12m2,密度为2650kg/m3。该系统中主要存在化学反应为阳离子交换过程,其表达式以如下两个为主:
Na++K-X=K++Na-X;
Na++0.5Ca-X2=0.5Ca2++Na-X;
K+和Ca2+以Na+为参照的阳离子交换系数分别为:KNa/K=0.1995和KNa/Ca=0.3981,阳离子交换容量为0.01779meq/100g。
各离子的初始浓度和边界浓度如表1:
表1柱子系统种各化学组分的初始浓度和边界浓度
基于上述信息,用TOUGHREACT进行模型构建,并计算出柱子末端Ca2+,Na+,K+和Cl-四种离子在25个时间点的浓度值,对这些浓度值加入服从分布ε~N(1,0.052)的噪声扰动,每种离子的浓度分布情况见图3;
在本实施例中:
步骤2中各参数先验区间见表2;
步骤2中的样本采样数为200;
步骤8中设置的容许误差为0.003;
步骤9中的局部采样参数个数n为20;
步骤11中的最大循环次数设置为20;
本实施例,经过18次循环,达到循环终止条件;得到各参数的反演结果及相对误差见表2;
表2待反演参数的真实值,先验区间,反演估计值和相对误差
根据表3的结果,参数反演值与其真实值都非常接近,相对误差均小于5%;
再将参数代入步骤1中建立的多组分反应性溶质运移模型,得到对应的模拟结果,绘制拟合曲线,如图4所示;
根据图4,在将参数估计值代入模型后,得到的模拟曲线与观测点数据的分布基本拟合。
由此可见,本发明能解决参数识别及模型校正问题。
Claims (13)
1.一种基于深度神经网络的多组分反应性溶质运移参数反演方法,其特征在于:包括以下步骤:
步骤1:通过资料收集或实地采集方式,得到各化学组分的实际浓度观测数据;根据研究对象实际情况,构建实际观测数据所对应的多组分反应性溶质运移模型;同时确定模型待反演参数和各参数的取值范围;
步骤2:根据模型参数取值范围,随机采样若干组参数样本,逐一代入步骤1中所建立的多组分反应性溶质运移模型,得到每组样本对应的模型响应,形成由“模型参数—模型响应”样本对构成的训练样本集;
步骤3:对步骤1中的实际观测数据和步骤2中的训练样本集作归一化处理;
步骤4:构建一个由模型参数作为输入,模型响应作为输出的正演神经网络;用步骤3中归一化处理后的训练样本集对本步骤中的正演神经网络进行训练;
步骤5:构建一个由模型响应观测数据作为输入,模型参数作为输出的反演神经网络;将该反演神经网络的输出作为步骤4中的正演神经网络的输入,得到一个串联式连接的神经网络架构;在此架构下用步骤3中归一化后的实际观测数据对反演神经网络进行训练;
步骤6:将步骤3中的归一化后的实际观测数据输入步骤5中的反演神经网络,网络输出得到归一化后的模型参数估计结果;将其反归一化,即为步骤1中所建立的待研究对象的多组分反应性溶质运移模型参数估计值;
步骤7:计算步骤6得到参数估计值处的正演神经网络相对于多组分反应性溶质运移模型的误差;
步骤8:判断步骤7中误差是否小于容许误差;如果是,则认为步骤6中所得到的模型参数估计值即为所求的反演结果;否则,执行后续步骤;
步骤9:如果步骤7中正演神经网络在参数估计点处的误差大于容许误差,则在步骤6中所得到的模型参数估计值附近增加n个参数样本,并逐一代入步骤1中所建立的地下水多组分反应性溶质运移模型,得到对应的模型响应,得到局部样本集;计算步骤2中的训练样本集中的所有样本点的参数与步骤6得到的参数估计值之间的欧氏距离,将欧式距离最大的n个样本删除,并加入局部样本集得到更新后的训练样本集;同时对新的样本集做归一化处理;
步骤10:利用步骤9中得到的归一化处理后的新训练样本集,用迁移学习的方法对步骤4中的正演神经网络进行更新,得到更新后的正演神经网络;再执行步骤5~步骤8中的操作;
步骤11:如果步骤8中的判断条件依然不满足,则继续执行步骤9~步骤10,直到循环次数达到最大终止条件,输出所有循环中近似误差最小的那一次参数估计结果作为最终的反演结果。
2.根据权利要求1所述的一种基于深度神经网络的多组分反应性溶质运移参数反演方法,其特征在于:所述步骤1中用TOUGHREACT程序构建地下水多组分反应性溶质运移模型。
3.根据权利要求1所述的一种基于深度神经网络的多组分反应性溶质运移参数反演方法,其特征在于:所述步骤1中的模型参数包括渗透率、阳离子交换系数、各化学组分初始浓度、各化学组分边界浓度和阳离子交换容量。
4.根据权利要求1所述的一种基于深度神经网络的多组分反应性溶质运移参数反演方法,其特征在于:所述步骤2中随机采样通过拉丁超立方采样方法实现。
5.根据权利要求1所述的一种基于深度神经网络的多组分反应性溶质运移参数反演方法,其特征在于:所述步骤2中通过编写程序,实现TOUGHREACT程序的输入文件中模型参数改写,TOUGHREACT程序可执行程序的调用和输出文件中结果的读取,从而得到由模型参数和模型响应构成的训练样本集。
6.根据权利要求1所述的一种基于深度神经网络的多组分反应性溶质运移参数反演方法,其特征在于:所述步骤3中的数据归一化通过“0-1归一化”方法进行处理;其中,参数归一化的方法如下:
模型响应的归一化方法如下:
8.根据权利要求1所述的一种基于深度神经网络的多组分反应性溶质运移参数反演方法,其特征在于:所述步骤4和步骤5中的正演神经网络和反演神经网络的搭建及训练均在pytorch框架下实现。
9.根据权利要求1所述的一种基于深度神经网络的多组分反应性溶质运移参数反演方法,其特征在于:所述步骤8中的容许近似误差根据不同的研究对象和求解精度需求通过人为设定。
12.根据权利要求1所述的一种基于深度神经网络的多组分反应性溶质运移参数反演方法,其特征在于:所述步骤10中利用步骤9中得到的更新后的训练样本集,仅对步骤4中的正演神经网络的输出层的神经网络连接权重参数进行更新训练,对其余连接权重参数固定,从而得到更新后的正演神经网络。
13.根据权利要求1所述的一种基于深度神经网络的多组分反应性溶质运移参数反演方法,其特征在于:所述步骤11中的最大循环次数通过人为设定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110100068.7A CN112733386B (zh) | 2021-01-25 | 2021-01-25 | 一种基于深度神经网络的多组分反应性溶质运移参数反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110100068.7A CN112733386B (zh) | 2021-01-25 | 2021-01-25 | 一种基于深度神经网络的多组分反应性溶质运移参数反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112733386A CN112733386A (zh) | 2021-04-30 |
CN112733386B true CN112733386B (zh) | 2021-08-17 |
Family
ID=75593970
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110100068.7A Active CN112733386B (zh) | 2021-01-25 | 2021-01-25 | 一种基于深度神经网络的多组分反应性溶质运移参数反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112733386B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113177330B (zh) * | 2021-05-27 | 2022-07-22 | 吉林大学 | 一种瞬变电磁快速统计学反演方法 |
CN113537354B (zh) * | 2021-07-19 | 2022-07-12 | 吉林大学 | 一种基于深度学习的含水层结构阶段式随机反演识别方法 |
CN113553773B (zh) * | 2021-08-16 | 2023-01-24 | 吉林大学 | 基于贝叶斯框架结合神经网络的地空电磁数据反演方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080071709A1 (en) * | 2006-08-22 | 2008-03-20 | Kjt Enterprises, Inc. | Fast 3D inversion of electromagnetic survey data using a trained neural network in the forward modeling branch |
CN111462230A (zh) * | 2020-02-18 | 2020-07-28 | 天津大学 | 一种基于深度强化学习的台风中心定位方法 |
CN111639746A (zh) * | 2020-05-14 | 2020-09-08 | 中国科学院国家空间科学中心 | 一种基于cnn神经网络的gnss-r海面风速反演方法及系统 |
CN111812732A (zh) * | 2020-06-29 | 2020-10-23 | 中铁二院工程集团有限责任公司 | 一种基于卷积神经网络的大地电磁非线性反演方法 |
CN112149353A (zh) * | 2020-09-24 | 2020-12-29 | 南京大学 | 基于卷积神经网络识别dnapl污染物在地下含水层分布的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111126591B (zh) * | 2019-10-11 | 2023-04-18 | 重庆大学 | 一种基于空间约束技术的大地电磁深度神经网络反演方法 |
-
2021
- 2021-01-25 CN CN202110100068.7A patent/CN112733386B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080071709A1 (en) * | 2006-08-22 | 2008-03-20 | Kjt Enterprises, Inc. | Fast 3D inversion of electromagnetic survey data using a trained neural network in the forward modeling branch |
CN111462230A (zh) * | 2020-02-18 | 2020-07-28 | 天津大学 | 一种基于深度强化学习的台风中心定位方法 |
CN111639746A (zh) * | 2020-05-14 | 2020-09-08 | 中国科学院国家空间科学中心 | 一种基于cnn神经网络的gnss-r海面风速反演方法及系统 |
CN111812732A (zh) * | 2020-06-29 | 2020-10-23 | 中铁二院工程集团有限责任公司 | 一种基于卷积神经网络的大地电磁非线性反演方法 |
CN112149353A (zh) * | 2020-09-24 | 2020-12-29 | 南京大学 | 基于卷积神经网络识别dnapl污染物在地下含水层分布的方法 |
Non-Patent Citations (3)
Title |
---|
Statistic inversion of multi-zone transition probability models for aquifer characterization in alluvial fans;Lin Zhu 等;《STOCHASTIC ENVIRONMENTAL RESEARCH AND RISK ASSESSMENT》;20150612;第30卷(第3期);第1005-1016页 * |
不同碎屑矿物CO2参与的水-岩作用效应数值模拟;李凤昱 等;《石油学报》;20160930;第37卷(第9期);第1116-1128页 * |
基于深度学习的地下水模拟高维不确定性分析和反演;莫绍星;《中国博士学位论文全文数据库 工程科技Ⅰ辑》;20210115;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112733386A (zh) | 2021-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112733386B (zh) | 一种基于深度神经网络的多组分反应性溶质运移参数反演方法 | |
CN110780164B (zh) | 基于yolo的绝缘子红外故障定位诊断方法及装置 | |
CN109828089B (zh) | 一种基于dbn-bp的水质参数亚硝酸氮的在线预测方法 | |
CN110807554B (zh) | 基于风电/光伏经典场景集的生成方法及系统 | |
CN110232471B (zh) | 一种降水传感网节点布局优化方法及装置 | |
CN106022954B (zh) | 基于灰色关联度的多重bp神经网络负荷预测方法 | |
CN104899448B (zh) | 一种集合卡尔曼滤波静态局地化方案的自适应补偿方法 | |
CN109872415B (zh) | 一种基于神经网络的车速估计方法及系统 | |
CN108710974A (zh) | 一种基于深度置信网络的水体氨氮预测方法及装置 | |
CN112132875B (zh) | 一种基于面特征的多平台点云匹配方法 | |
CN107622276B (zh) | 一种基于机器人仿真与物理采样结合的深度学习训练方法 | |
CN113344956A (zh) | 基于无人机航摄三维建模的地物轮廓提取分类方法 | |
CN112949089B (zh) | 一种基于离散卷积残差网络的含水层结构反演识别方法 | |
CN112949944B (zh) | 一种基于时空特征的地下水位智能预测方法及系统 | |
CN107909192A (zh) | 土壤重金属含量的估测方法及装置 | |
CN115797808A (zh) | 一种无人机巡检缺陷图像的识别方法、系统、装置及介质 | |
CN115952685A (zh) | 基于集成深度学习的污水处理过程软测量建模方法 | |
CN117933095B (zh) | 一种基于机器学习的地表发射率实时反演及同化方法 | |
CN114972339A (zh) | 用于推土机结构件生产异常检测的数据增强系统 | |
CN114034334B (zh) | 岩溶管道污染源和流量的识别方法 | |
CN117313564B (zh) | 反演海洋-大气光学参数的方法、装置及存储介质 | |
CN114638358A (zh) | 消光系数转颗粒物浓度方法、设备以及计算机存储介质 | |
CN108334977B (zh) | 基于深度学习的水质预测方法及系统 | |
CN115618988A (zh) | 基于深度神经网络的海水温度、盐度与流速的三维时空场联合预测方法 | |
CN117057221A (zh) | 一种基于机器学习的滨海非均质含水层刻画实现方法及装置 |
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 |