一种基于自然单元的2.5维自然电场数值模拟方法技术

本发明专利技术提供了一种基于自然单元的2.5维自然电场数值模拟方法，具体步骤如下：S1、建立二维自然电场地电模型，依据模型参数设置离散点集；S2、构建2.5维自然电场地电模型边值问题，推导其基于自然单元的变分形式及其积分形式；S3、构建基于Laplace插值函数的自然单元法，推导2.5维自然单元形函数及其导数的基本方程；S4、基于S2、S3推导的公式，计算各波数条件下的总刚度矩阵；S5、处理边界条件及电源信息；S6、求解大型稀疏方程组得到各波数条件下的波数域自然电位剖面，最后将各波数域自然电位剖面通过傅里叶反变换得到空间域下的自然电位剖面。

【技术实现步骤摘要】
一种基于自然单元的2.5维自然电场数值模拟方法
本专利技术涉及自然电场数值模拟方法
，尤其涉及一种基于自然单元的2.5维自然电场数值模拟方法。
技术介绍
随着我国经济的持续高速发展，在工业化、城镇化过程中的环境污染问题日益严重。即时检测即时控制污染源的扩散是保护土壤和地下水不受污染或少受污染的积极方法。与钻孔测试、定期采样分析等传统环境监测手段相比，自然电场法具有效率高、成本低、快速无损、无二次污染等优点，且对多种与污染扩散伴生的微生物活动信号、氧化还原反应以及渗流运动相当敏感，适用于地下污染物的检测以及长周期实时监测。有针对性地开展自然电场数值模拟工作，有助于反演解释，提高自然电场法在工程与环境污染监测、检测等方面的应用效果。自然电场数值模拟所采用的常规数值方法有积分方程法、有限元法、有限差分法以及有限体积法等。这些数值方法各有优势，但也同时存在共同的缺陷，即需要构建网格单元以构造插值形函数，该过程不仅增加了前期处理工作量，也使得复杂模型数值模拟变得困难。而自然单元法节点布置灵活，形函数构造过程无需网格单元信息，前处理过程简单方便，模拟精度高，能有效实现自然电场等稳定电流场复杂模型的数值模拟。
技术实现思路
本专利技术目的在于提供一种基于Laplace插值构造自然单元形函数，实现一种基于自然单元的2.5维自然电场数值模拟方法。为实现上述目的，本专利技术提供了一种基于自然单元的2.5维自然电场数值模拟方法，其具体步骤如下：S1：建立二维自然电场地电模型，依据模型参数设置离散点集；S2：构建2.5维自然电场地电模型边值问题，推导其基于自然单元的变分问题及其积分形式；S3：构建基于Laplace插值函数的自然单元法，推导2.5维自然单元形函数及其导数的基本方程；S4：基于S2、S3推导的公式，计算各波数条件下的总刚度矩阵；S5：处理边界条件及电源信息；S6：求解大型稀疏方程组得到各波数条件下的波数域自然电位剖面，然后将各波数域自然电位剖面通过傅里叶反变换得到空间域下的自然电位剖面。作为本专利技术的进一步方案：所述步骤S1中的所述二维自然电场地电模型的尺度设置为20m*30m(深度*横度)；所述离散点集设置有651个节点，地下水位位于2m深度处，上层电阻率为50欧姆米，下层电阻率为100欧姆米。作为本专利技术的进一步方案：所述步骤S2中的所述构建2.5维自然电场地电模型边值问题的方法为采用第三类边界条件进行构建。作为本专利技术的进一步方案：所述步骤S3中的自然单元法为一种基于Voronoi图和Delaunay三角化几何结构，以自然邻点插值为试函数的一种无网格数值方法；作为本专利技术的进一步方案：所述步骤S4中对波数及其权值的选择具体如下：波数k＝[0.004758,0.0407011,0.1408855,0.393225,1.088038]；波数权值g＝[0.0099472,0.031619,0.0980327,0.2511531,0.7260814]。所述总刚度矩阵的具体过程如下：(1)对全域的各背景积分网格进行循环；(2)对各背景积分网格中的高斯积分点进行循环；(3)搜索各高斯积分点的自然邻点；(4)计算各自然邻点相关于高斯积分点的形函数及导数；(5)计算各高斯积分点的子刚度矩阵；(6)将各子刚度矩阵加载到总刚度矩阵中，直至全域高斯积分点循环完毕。作为本专利技术的进一步方案：所述步骤S6中各波数条件下的波数域自然电位剖面的具体计算过程如下：(1)在各波数条件下，求解得到当前波数对应的总刚度矩阵；(2)结合电源项，求解方程组，得到各波数条件下的波数域总自然电位剖面；(3)通过傅里叶反变换将各波数域自然电位剖面转换为空间域中的总自然电位剖面。应用本专利技术的技术方案，具有以下有益效果：(1)本专利技术所提供的方法不依赖于网格单元，节点布置灵活，计算精度高，适用于自然电场复杂模型的高精度数值模拟。(2)本专利技术提供的方法，通过基于Laplace自然邻点插值的二维自然单元法，经傅里叶正反变换，实现波数域到空间域的2.5维自然电场数值模拟。(3)本专利技术提供的方法，在模型边界处满足精确的线性插值，可以准确施加边界条件。(4)本专利技术提供的方法，所有计算最终都可转化为简单的代数运算，计算效率高，编程简便。(5)本专利技术提供的方法，能推动自然电场等多源稳定电流场复杂模型数值模拟工作的开展，为自然电场法应用到污染监测、检测以及其他工程与环境地球物理问题中提供数值计算基础。(6)本专利技术中步骤S1所建立的二维自然电场地电模型，该模型可近似模拟微生物降解地下有机污染物时，在地下水位附近发生的氧化还原反应所引起的自然电位异常。(7)本专利技术中步骤S3中所采用的自然单元法具有无网格法和经典有限元方法的优点，又克服了两者的缺陷；其形函数满足插值性质，可以像有限元法一样直接施加本质边界条件，不存在基于移动最小二乘拟合的无网格方法不能直接施加本质边界条件的难题。(8)本专利技术中步骤S3中的所述Delaunay三角化具有最大最小角和空外接圆两个重要性质，且在二维情况下，Delaunay三角形网格还具有所形成的三角形互不重叠和所形成的三角形可以覆盖整个平面的特点。(9)本专利技术中，自然邻点插值函数选用Laplace插值。Laplace插值的优点主要体现在：A、直接利用节点的一阶Voronoi结构的边长和点到Voronoi边的距离计算插值基函数，避免了Sibson插值所采用的Waston面积叠加算法需要反复进行迭代的麻烦，计算量大大降低；B、除凸边界外，在凹边界上也满足精确的线性插值，可以准确施加边界条件。(10)本专利技术中通过采用将各波数域自然电位剖面通过傅里叶反变换得到空间域下的自然电位剖面，其通过数值模拟结果表明，地质微生物降解地下有机污染物时所发生的氧化还原反应能在地表产生负的自然电位异常；同时表明本专利技术所提的数值模拟方法是有效的可行的。除了上面所描述的目的、特征和优点之外，本专利技术还有其它的目的、特征和优点。下面将参照图，对本专利技术作进一步详细的说明。附图说明构成本申请的一部分的附图用来提供对本专利技术的进一步理解，本专利技术的示意性实施例及其说明用于解释本专利技术，并不构成对本专利技术的不当限定。在附图中：图1是本专利技术基于自然单元的2.5维自然电场实施例地电模型；图2是本专利技术中节点A的Voronoi单胞；图3是本专利技术中7节点Voronoi图；图4是本专利技术中Delaunay外接圆的示意图；图5是本专利技术中Delaunay三角化的示意图；图6是本专利技术中7节点Voronoi图中的x计算点的示意图；图7是本专利技术中x计算点的自然邻点及其Laplace插值元素的示意图；图8是本专利技术中基于Laplace插值的自然邻点插值函数示意图；图9是本专利技术中自然邻点插值函数对本文档来自技高网...

【技术保护点】
1.一种基于自然单元的2.5维自然电场数值模拟方法，其特征在于，具体步骤如下：/nS1：建立二维自然电场地电模型，依据模型参数设置离散点集；/nS2：构建2.5维自然电场地电模型边值问题，推导其基于自然单元的变分问题及其积分形式；/nS3：构建基于Laplace插值函数的自然单元法，推导2.5维自然单元形函数及其导数的基本方程；/nS4：基于S2、S3推导的公式，计算各波数条件下的总刚度矩阵；/nS5：处理边界条件及电源信息；/nS6：求解大型稀疏方程组得到各波数条件下的波数域自然电位剖面，然后将各波数域自然电位剖面通过傅里叶反变换得到空间域下的自然电位剖面。/n

【技术特征摘要】
1.一种基于自然单元的2.5维自然电场数值模拟方法，其特征在于，具体步骤如下：
S1：建立二维自然电场地电模型，依据模型参数设置离散点集；
S2：构建2.5维自然电场地电模型边值问题，推导其基于自然单元的变分问题及其积分形式；
S3：构建基于Laplace插值函数的自然单元法，推导2.5维自然单元形函数及其导数的基本方程；
S4：基于S2、S3推导的公式，计算各波数条件下的总刚度矩阵；
S5：处理边界条件及电源信息；
S6：求解大型稀疏方程组得到各波数条件下的波数域自然电位剖面，然后将各波数域自然电位剖面通过傅里叶反变换得到空间域下的自然电位剖面。


2.根据权利要求1所述的基于自然单元的2.5维自然电场数值模拟方法，其特征在于：所述步骤S2中的所述构建2.5维自然电场地电模型边值问题的方法为采用第三类边界条件进行构建。


3.根据权利要求1所述的基于自然单元的2.5维自然电场数值模拟方法，其特征在于：所述步骤S3中的自然单元法为一种基于Voronoi图和Delaunay三角化几何结构，以自然邻点插值为试函数的一种无网格数值方法。


4.根据权利要求1所述的基于自然单元的2.5维自然电场数值模拟方法，其特征在于：所述Voronoi图和Delaunay三角化的概念来自于计算几何结构，所述的几何结构是由一组不规则点定义的最基本的几何结构。
...

【专利技术属性】
技术研发人员：崔益安，谢静，张丽娟，郭友军，柳建新，
申请(专利权)人：中南大学，
类型：发明
国别省市：湖南;43